Holographic thermalization with a chemical potential from Born-Infeld electrodynamics

Giancarlo Camilo, Bertha Cuadros-Melgar and Elcio Abdalla

Instituto de Física, Universidade de São Paulo

C.P. 66318, CEP: 05315-970, São Paulo, Brasil

Escola de Engenharia de Lorena, Universidade de São Paulo,

Estrada Municipal do Campinho S/N, CEP: 12602-810, Lorena, Brasil

^{†}^{†}footnotetext: Email: gcamilo@usp.br, bertha@usp.br, eabdalla@usp.br

Abstract

The problem of holographic thermalization in the framework of Einstein gravity coupled to Born-Infeld nonlinear electrodynamics is investigated. We use equal time two-point correlation functions and expectation values of Wilson loop operators in the boundary quantum field theory as probes of thermalization, which have dual gravity descriptions in terms of geodesic lengths and minimal area surfaces in the bulk spacetime. The full range of values of the chemical potential per temperature ratio on the boundary is explored. The numerical results show that the effect of the charge on the thermalization time is similar to the one obtained with Maxwell electrodynamics, namely the larger the charge the later thermalization occurs. The Born-Infeld parameter, on the other hand, has the opposite effect: the more nonlinear the theory is, the sooner it thermalizes. We also study the thermalization velocity and how the parameters affect the phase transition point separating the thermalization process into an accelerating phase and a decelerating phase.

###### Contents

## 1 Introduction

Over the last years the AdS/CFT correspondence [1, 2, 3] has proved to be a powerful tool to describe gauge theories at strong coupling regime, where standard perturbative methods fail to work, by mapping them to a dual gravitational theory in higher dimensions whose treatment is tipically much easier. For that reason, it has found a wide range of applications in different areas of theoretical physics going from condensed matter systems to quantum chromodynamics. In particular, if the strongly coupled gauge theory is at zero temperature, it is well known that the dual gravity description involves a purely AdS space, while a finite temperature field theory requires the presence of an asymptotically AdS black hole. An interesting implication of this equivalence will be investigated in the present work, namely, the process of black hole formation in AdS space as the analog of nonequilibrium dynamics that leads a system towards thermal equilibrium after a sudden injection of energy.

One of the most interesting scenarios where these ideas have been applied is to describe properties of the quark-gluon plasma (QGP) formed in heavy ion colliders such as the RHIC and LHC. Recent results suggest that the QGP behaves as an ideal fluid with a very small shear viscosity over entropy density ratio () [4]. This implies that the QGP takes place at a strong coupling regime and therefore is amenable to a dual gravity treatment. Indeed, holographic calculations using the prototypical SU() supersymmetric Yang-Mills (SYM) theory at finite and its string theoretical AdS gravity dual show that there seems to be a small and universal lower limit for the ratio for all theories with gravity duals [5, 6]. This is one of the most prominent predictions of AdS/CFT at the moment. However, while the near-equilibrium dynamics (e.g., transport coefficients such as viscosity and electrical conductivity [7] and more general aspects of dissipative hydrodynamics [8, 9]) of the QGP is well known, the far from equilibrium process of formation of QGP after a heavy ion collision, often referred to as thermalization, is not well understood. The thermalization time scale observed at RHIC is considerably shorter than expected according to perturbative techniques [10], reinforcing the need of a strong coupling description of the thermalization process.

Some attempts were made in the literature to address this problem from a holographic point of view as a dual process of black hole formation via gravitational collapse in AdS space [11, 12, 13, 14, 15, 16, 17, 18, 19]. A slightly different and simpler model for holographic thermalization was introduced by Balasubramanian et al. in [20, 21], which despite its simplicity captures many important features of the thermalization process. It consists in the collapse of a thin shell of matter described by an AdS-Vaidya metric that interpolates between pure AdS space at early times and Schwarzschild-AdS black hole at late times. The authors used the dynamical background of the collapsing shell to study the time evolution of nonlocal thermalization probes of the boundary conformal field theory with well known dual gravity descriptions in terms of geometric quantities. They considered equal-time two-point correlation functions of local gauge invariant operators, expectation values of Wilson loop operators, and entanglement entropy which correspond in the gravity side to minimal lengths, areas, and volumes in AdS space, respectively. They found that the thermalization is a top-down process (i.e., UV modes thermalize first while IR modes thermalize later), in contrast to the predictions of bottom-up thermalization from perturbative approaches [22]. This has a clear and intuitive interpretation from the AdS/CFT perspective: UV modes correspond to small distance scales in the boundary of AdS space and, therefore, they do not capture much of the details of the collapse process happening deep into the bulk. IR modes, on the other hand, penetrate deeper into the bulk and for that reason are naturally more sensible to details of the bulk dynamics, therefore they should thermalize later. In addition, the authors found that the thermalization time scales typically as , where is the characteristic length of the probe.

A natural extension of this model was proposed in [23] (see also [24]) to include the effect of a non-vanishing chemical potential , which is usually the case in real heavy ion collision processes. The authors considered the collapse of a thin shell of charged matter in the bulk described by an AdS-Vaidya-like metric leading to a thermal equilibrium configuration given by a Reissner-Nordström-AdS black hole. They argue that varying the charge of the final state black hole from zero to the extremal value allows to explore the full range of chemical potential per temperature ratio in the dual strongly coupled QGP. Their main conclusion was that as the charge is increased, the thermalization time for renormalized geodesic lengths and minimal area surfaces becomes larger. Further investigations were made later to include Gauss-Bonnet higher curvature corrections [25, 26], angular momentum [27], noncommutative [28], Lifshitz and hyperscaling violating geometries [30, 29], de Sitter boundary field theories [31], as well as more elaborated dynamics for the collapsing shell [32, 33, 34] and a discussion on spectral functions of boundary two-point correlators [35]. Related work can also be found in [36, 37, 38, 39, 40, 41, 42].

In this paper we propose a further study of holographic thermalization with a chemical potential using Einstein gravity coupled to Born-Infeld (BI) nonlinear electrodynamics in the bulk. This is a natural generalization of the discussion initiated in [23], since it accomodates more elaborated dynamics for the gauge field including (all order) higher-derivatives of and, therefore, it may give rise to interesting effects on the chemical potential of the dual boundary theory that are not captured by the Maxwell description. Although BI electrodynamics has its origin a long time ago [43] as an attempt to obtain a finite self-energy of point-like charged particles, currently a renewed interest has been raised due to recent developments in superstring theory. In particular, it is well known that the low energy behavior of the vector modes of open strings is governed by the BI action [44, 45], while the low energy dynamics of D-branes is given by a similar non-Abelian version of the BI action [46] (see also [47, 48]). Consequently, BI electrodynamics provides a promising scenario to explore deviations from Maxwell electrodynamics, specially from the point of view of AdS/CFT calculations where string theory plays a prominent role (see [49, 50, 51, 52, 53, 54, 55, 56, 57, 58] for an incomplete list of previous works in this direction).

The paper is structured as follows. In section 2 we review the black hole solutions of Einstein-Born-Infeld theory in AdS space and construct their Vaidya-like extensions modelling the collapse of a thin charged shell, to be used in the sequence. Section 3 is devoted to the holographic setup for the non-local observables chosen as probes of thermalization, namely, the equal-time two-point correlators and the expectation value of Wilson loops. In Section 4 we give details of the numerical calculations and present all the results with the effects of the chemical potential and BI parameter on the thermalization curves and velocities. Finally, Section 5 contains our concluding remarks.

## 2 Vaidya AdS Black Hole solutions in Einstein-Born-Infeld theory

The starting point is the -dimensional Einstein gravity action with a negative cosmological constant (being the AdS curvature radius) minimally coupled to Born-Infeld electrodynamics

(1) |

where is given by

(2) |

The constant is the BI parameter with dimension of mass.^{1}^{1}1In the context of string theory, appears tipically in terms of the string parameter via . It is defined in such a way that the limit corresponds to the standard Maxwell Lagrangian. We choose units in which , being the Newton’s constant in dimensions.

The charged black hole solution to the equations of motion coming from the action (1), first obtained in [59] (see also [60]), reads

(3) |

where denotes the metric on the unit sphere and

(4) | |||||

In the above equation is the hypergeometric function and are integration constants related to the ADM mass and charge of the black hole via ^{2}^{2}2For simplicity, we will keep referring to and hereinafter simply as “mass”and “charge”parameters of the black hole without any risk of confusion.

being the volume of the . There is also a purely electric gauge field given by

(5) |

where is a constant corresponding to the electrostatic potential at , which will be related to the chemical potential in the dual gauge theory according to the AdS/CFT correspondence. It is defined such that the gauge field vanishes at the horizon, i.e.,

(6) |

The electric field associated to (5) is finite at the origin , which is a key feature of BI theories. The black hole function (4), on the other hand, is in general singular at the origin. Such a singularity is hidden behind an event horizon provided the free parameters are chosen so that the equation admits a real positive solution. We should also mention that taking the limit in (4) gives the well known Reissner-Nordström-AdS black hole studied in [61].

The solution (3) has the topology of at the AdS boundary . In the context of the AdS/CFT correspondence it is interesting to consider the limit where the AdS boundary is instead, since one is often interested in dual gauge theories living on flat space. This procedure is known in the literature as the “infinite volume limit”, and it arises only due to the presence of a negative cosmological constant [61]. The idea is to introduce a dimensionless parameter (which will be set to ) and rescale all dimensionful quantities as while at the same time blow up the as . This leaves the block of the metric almost invariant (except for the contribution of the constant term in (4)). Finally, taking yields

(7) |

where . Notice that now the horizon, defined by , is planar instead of spherical, so we should refer to (7) as a black brane instead of a black hole. In order to avoid the coordinate singularity at it will be interesting to express the metric in Eddington-Finkelstein coordinates by introducing a new time coordinate defined by , and also it will be convenient to work with an inverse radial coordinate such that the AdS boundary stays at while the singularity sits at infinity. The resulting metric is

(8) |

where we have defined

(9) |

Notice that near the AdS boundary .

The Hawking temperature of a black hole in the context of AdS/CFT can be viewed as the equilibrium temperature of the dual field theory living on the boundary. It is obtained as usual by continuing the black hole metric to its Euclidean version via and demanding the absence of conical singularities at the horizon. This results in a periodic Euclidean time whose period is identified with the inverse Hawking temperature. For the AdS Einstein-Born-Infeld black brane (7) this calculation gives

(10) |

This expression reduces to the Hawking temperature of the Reissner-Nordström-AdS black hole in the Maxwell limit [23]. When , the black brane is called extremal. If we think of all the parameters but the charge as fixed, then we can characterize the extremal black brane solution by a maximal value of charge given by

(11) |

According to the AdS/CFT dictionary, the asymptotic value of the time component of the gauge field at the AdS boundary (namely, the constant in equation (5)) corresponds to the chemical potential in the dual quantum field theory, . Actually, the precise relation should include some scale with length units since the chemical potential must have energy units (or [length]) while as defined by the action (1) is dimensionless. Hence, the chemical potential per temperature ratio of the boundary field theory is given by

(12) |

with and given by expressions (6) and (10), respectively. A remarkable feature is that if the horizon radius and the BI parameter are kept fixed, then by varying the charge from (vanishing ) to (vanishing ) it is possible to explore the whole range of values of the ratio in the dual field theory, i.e., from to .

A Vaidya-like extension of the BI AdS black brane metric (7) can be constructed by promoting the mass and charge to arbitrary functions and of the advanced time . The resulting dynamical metric has the same form as in (8) but now with instead of due to the time dependence introduced on the mass and charge. The same also holds for the gauge field (5), which now becomes . Such a spacetime describes the collapse of a thin-shell of charged dust from the boundary of the AdS space towards the bulk interior.

Of course such a metric is not a solution of the action (1) anymore: there must be some external matter action sourcing the time variation of and . If we take this external contribution into account, the Einstein-BI equations of motion become (we restore the factors of for a moment):

(13) | ||||

(14) |

The Vaidya-BI-AdS metric above-mentioned is a solution to these equations provided the external sources satisfy

(15) | ||||

(16) |

where the dot denotes . We notice that there is no dependence on above, and indeed this is exactly the same current found in [23] in the Vaidya-Reissner-Nordström-AdS case. , on the other hand, differs from the corresponding one in the Reissner-Nordström case due to the hypergeometric term (but naturally reduces to it in the limit, since ).

## 3 Holographic thermalization

In this section we are going to study the thermalization process of a strongly coupled quantum field theory whose bulk gravity dual corresponds to the Einstein-Born-Infeld system presented in the previous section. According to the AdS/CFT correspondence, a zero temperature state on the dimensional boundary theory is dual to pure in the bulk, while a thermal state corresponds to the BI-AdS black brane (8). Therefore, any dynamical bulk spacetime which interpolates between these two situations would be a natural candidate to holographically model the nonequilibrium process leading to thermalization of the boundary theory after a rapid injection of energy.

Following [20], we choose our dynamical bulk metric to be the Vaidya-BI-AdS metric discussed in Section 2, namely (we set the AdS radius hereafter)

(17) |

with the mass and charge functions given by (see [62] for an interesting discussion on the corresponding bulk null energy condition)

(18) | ||||

(19) |

Clearly, for we have pure AdS () and for we have and which is the BI-AdS black brane (8). Indeed, equations (18)-(19) are just smooth versions (convenient for the numerical analysis) of the step functions and , which represent a shock wave (a zero thickness shell of charged matter suddenly forming at ). The constant represents a finite shell thickness and for we go back to the step function.

The next step is to choose a set of observables to use as probes of thermalization. Since local observables in the boundary such as expectation values of the energy-momentum tensor are not sensitive to the thermalization process, one needs to consider extended non-local observables.^{3}^{3}3Holography provides a geometric intuition for why local operators are insensitive to details of the progress towards thermalization: being local, they are only sensitive to phenomena happening in their vicinity near the AdS boundary. Thus they are not aware of the details of phenomena occurring near the thermal scale. We need observables dual to AdS quantities that probe deeper into the bulk in order to see signals of thermalization. In this work we shall focus on equal time two-point correlation functions and expectation values of rectangular Wilson loops, which have well known holographic descriptions in the bulk in terms of renormalized geodesic lengths and minimal area surfaces, respectively. A third observable that could be used is the entanglement entropy of boundary regions, which have a very similar description in terms of minimal volumes of codimension-two surfaces in the bulk. However, as the results of entanglement entropy lead essentially to the same conclusions, we will not show them here in order to avoid unnecessary repetitions.

### 3.1 Renormalized geodesic lengths and two-point functions

We start with the equal-time two-point correlation functions of local gauge invariant operators of conformal dimension . The AdS/CFT correspondence provides a simple geometrical way to compute it in the bulk gravity dual when the operator is “heavy”, i.e., when [63]. Namely,

(20) |

where stands for the length of the bulk geodesic between the points and located on the AdS boundary.^{4}^{4}4If there is more than one geodesic we should sum over them on the right-hand side. Actually, one should be careful when doing such an approximation because the geodesic length above is divergent due to the contribution of the AdS boundary (because the AdS metric itself diverges at the boundary ). In order to extract a meaningful quantity we need to introduce a cutoff . It turns out that the divergent part of is universal and equals to [20]. Then, we define a (finite) renormalized geodesic length , which will be related to the renormalized two-point function just as in equation (20).

We choose the coordinate axes for the boundary directions such that the spatial separation between the points lies entirely over the direction, i.e., we consider space-like geodesics between the boundary points and , where the ellipsis denotes the remaining coordinates which are the same for both points. Then, by symmetry the geodesics cannot depend on coordinates other than , and we can use as the geodesic parameter (we call it simply hereinafter). The solutions of the geodesic equations are then given by a pair of functions and . The boundary conditions at the AdS boundary are

(21) |

The length functional between the referred points follows immediately from the line element (17) as being

(22) |

It clearly depends on the path taken from one point to the other. The geodesic corresponds to the functions and that minimize the length and can be found by standard methods. The geodesic length (which we will call ) is just the value of the length functional evaluated at the geodesic solution.

The variational problem simplifies by noticing that the integrand in (22) does not depend explicitly on and, therefore, there is a conserved Hamiltonian given by

(23) |

Using the conditions on the turning point of the geodesic,

(24) |

arising from the fact that the geodesic must be symmetric with respect to , the conservation equation simplifies to

(25) |

The advantage of working with the boundary conditions (24) instead of the original ones (21) is that we can use the conservation equation above to write the geodesic length (22) as

(26) |

Hence, after solving the equations of motion and finding the geodesic functions and for a given pair of initial conditions it is a trivial task to use the relations (21) to read the corresponding values of boundary separation and time , as well as obtaining the corresponding renormalized geodesic length by means of expression (26).

The hard part is to find the geodesic, which means to solve the Euler-Lagrange equations for and . With the help of the conservation equation (25) they can be written, respectively, as

(27a) | ||||

(27b) |

This is a set of coupled, highly nonlinear differential equations and for that reason it is quite hard to handle with analitical methods. However, for a given pair it is possible to find a numerical solution subject to the boundary conditions (24). Indeed, by solving for sufficiently many pairs of initial conditions (carefully chosen in order to give the same boundary separation and different times), we can track time after time the whole evolution of the geodesics in the Vaidya-BI-AdS spacetime. In particular, if we calculate the renormalized geodesic length of each of these solutions using (26) we will be able to see the full time evolution of towards thermalization. This will be done in Section 4, where we provide a detailed explanation of the numerical procedure as well as the choice of parameters and show our results.

### 3.2 Minimal area surfaces and Wilson loops

We now study a second class of thermalization probes, namely the expectation values of Wilson loop operators in the boundary field theory. The Wilson loop is a non-local gauge-invariant observable defined as the path-ordered integral of the gauge field over a closed path :

(28) |

where is the number of colors and is the non-abelian gauge field. Wilson loops contain useful information about the non-perturbative behavior of non-abelian gauge theories, such as whether they exhibit confinement or not. It is possible, in principle, to express all gauge-invariant functions of in terms of Wilson loops by appropriate choices of the path , but unfortunately they are in general hard to compute.

The AdS/CFT correspondence again provides an elegant way to compute the expectation value of Wilson loops of a strongly coupled gauge theory with a gravitational dual in terms of a geometrical quantity in the bulk [64]:

(29) |

where is the inverse string tension, denotes the minimal area bulk surface whose boundary is the original contour , and is the area of that surface. will be a solution of the bosonic part of the string action (the Nambu-Goto action), which is nothing but the area of the classical world-sheet with as its boundary. Indeed, equation (29) has its origin in a saddle-point approximation of the string theory partition function around the classical solution, which is only valid when the dual gauge theory has strong coupling.

Now we proceed to compute the minimal area surfaces in the Vaydia-BI-AdS spacetime as described before just as we did for the geodesic lengths. We will focus on a spacelike rectangular Wilson loop on the boundary. The rectangle can always be chosen to lie on the - plane, centered at the origin, with sides on the direction and on the direction. One also assumes translational invariance along , such that the shape of the bulk surface depends only on and again we can use to parametrize the functions and that characterize the surface. The boundary conditions at the AdS boundary are again given by equations (21).

Using the Vaidya-BI-AdS metric (17), the Nambu-Goto action (or area functional divided by ) becomes

(30) |

Notice that an obvious consequence of our assumption of translational invariance is that the length factorizes. Since we are interested just in the dependence, we can study instead of itself and forget about in what follows. As for the geodesics, the pair of functions that minimizes the Nambu-Goto action will be the minimal surface . The on-shell value of the Nambu-Goto action (i.e., ), which we call , will be our object of interest.

The subsequent calculation is closely analogous to the geodesics case. There is again a conserved Hamiltonian associated to (30) and we can introduce the alternative boundary conditions on the turning point of the minimal surface, which are the same as equations (24), to find an expression similar to (25) for the conservation equation. Replacing this back into (30), the on-shell Nambu-Goto action becomes

(31) |

which can be used to easily obtain the minimal area surface once we have solved the equations of motion and found the functions and for a given pair of initial conditions . The corresponding values of boundary separation and time can be read from the original conditions (21) as well. Here again we have to face the problem that the area diverges due to the contribution near the AdS boundary, but we can regularize the divergent part using again a cutoff and subtract it from to define the renormalized minimal area [20].

The functions and are found by solving the Euler-Lagrange equations coming from the action (30). The simplified equations of motion (after using the conservation equation) are, respectively

(32a) | ||||

(32b) |

As before, for a given pair we can solve this set of equations numerically subject to the boundary conditions (21). Thus, for some chosen length of the Wilson loop, by iterating for enough pairs we can track the whole time evolution of the minimal surfaces and, in particular, of their renormalized areas towards thermalization. This will be the aim of Section 4.

## 4 Numerical results

### 4.1 Renormalized geodesic lengths

In this section, we numerically solve the geodesic equations of motion (27) in order to find how the geodesic length evolves with time. Afterwards, we explore how the charge of the black hole and BI parameter affect the thermalization time. For the latter it will be convenient in the numerical calculations to use an inverse BI parameter instead of the original , such that the Maxwell limit is and increasing accounts for increasingly nonlinear electrodynamics.

First of all, we fix the free parameters. We will take the shell thickness and AdS space UV cut-off to be and , respectively. Since the effect of the number of spacetime dimensions and boundary separation on the thermalization probes has already been analyzed in previous works [20, 21, 23], we focus here in the case (namely, AdS space, which is dual to a dimensional gauge theory) and a fixed boundary separation . The mass of the final state black brane can be expressed in terms of the radius of its event horizon using the definition of , i.e., the largest solution of . Then, if we choose to fix the horizon at ,^{5}^{5}5This will be interesting to compare the thermalization for different values of and , since the black hole will always form at the same location for all . the mass is given by

(33) |

which of course only holds provided that and take values consistent with the existence of an event horizon. For a given this means that is allowed to take values from to the extremal value (11), which now reads

(34) |

As we pointed out in Section 2, for each , considering values of charge we can study all the range in the dual gauge theory. It is instructive to illustrate this fact here by looking at the form of expression (12) in the present case (we choose the scale for simplicity), namely,

(35) |

A plot of this as a function of the charge for distinct values of is shown in Figure 1, showing that indeed all the range of is covered. We also notice that for small values of charge (up to ) the BI parameter has no effect on the chemical potential since all the curves agree, so we will concern only about charges above this value in what follows.

Now we proceed to investigate the effects of and on the thermalization process. We choose the test values and for the numerical analysis. The procedure is the following: for a given pair we solve the geodesic equations (27) subject to the boundary conditions (24) characterized by a pair of values . We do this recursively for various pairs of initial conditions and collect from these just those yielding a boundary separation .^{6}^{6}6There is a small subtlety here. is determined numerically, via equation (21), after we have solved the equations using the modified boundary conditions (24), hence we must establish a criterion for what we mean by “”. We adopt the convention of tolerance, meaning that “”here corresponds to . Each of these collected solutions correspond to a different stage of the motion of the geodesics, as we may check by computing the corresponding boundary time via equation (21). We then calculate the renormalized geodesic length of each of these collected solutions and construct a list of points which contains all the information about the time evolution of the renormalized geodesic length. Actually, it will be convenient to divide all the lengths by in order to obtain a dimensionless, independent quantity . In addition, we subtract from this the final (thermal) value just to force all thermalization curves to end at zero, so that the quantity to be plotted is versus .

Before showing the thermalization curves, in Figure 2 we present an intuitive simple view of the effect of the BI parameter on the thermalization. It consists of a sequence of snapshots of the time evolution of geodesic profiles as well as the shell of charged dust described by the Vaidya-BI-AdS metric to form a black brane at at late times, for different values of the inverse BI parameter and fixed charge . Each column, top to bottom, corresponds to the time evolution for a given value of . It is found that at the early stages of the evolution, up to , the value of has little effect on the dynamics. After that, plays a crucial role in the evolution. We see that the bigger is (i.e., columns to the right), the sooner the black hole is formed. This is clear from the bottom line of the picture, corresponding to boundary time , where the black brane has just formed while the one is about to form and the one still needs some time to do so. This is a hint that the thermalization of the dual boundary field theory occurs sooner as increases, what indeed will be confirmed in Figure 4.

In Figures 3 and 4 we show the thermalization curves for the renormalized geodesic lengths with varying charge at fixed and varying at fixed charges, respectively. Instead of plotting point by point all the results obtained as described above, we find more instructive to fit those points using some polynomial function and plot the resulting curve. Details about the fits will be given below. We use dashed curves in all the plots hereinafter to highlight the Maxwell limiting case studied in [23]. The thermalized (final) state corresponding to the completely formed black brane is reached in each case when the curve touches the zero point of the vertical axis. The effect of the charge on the thermalization is clear from Figure 3. As grows, the thermalization time increases, meaning that the dual field theory thermalizes later. This had already been pointed out in [23] for the case of Maxwell electrodynamics (), and now we show that the same holds for BI nonlinear electrodynamics. Since the charge corresponds to the chemical potential, this means that the smaller the chemical potential is, the faster is the pair production and the screening effect takes over in an easier way. This is compatible with lower dimensional models, where screening effects are known to prevail over confinement [65, 66]. The second, more interesting result, is the effect of the inverse BI parameter shown in Figure 4. As one can see, increasing decreases the thermalization time, which means that the more nonlinear the bulk theory is, the sooner its dual field theory thermalizes. This confirms our intuition coming from the analysis of the geodesic profiles and shell motion in Figure 2. Such a behavior is similar to the effect of the Gauss-Bonnet parameter on the thermalization reported in [25]. As we discuss in the conclusions, this seems to be a general feature of introducing extra derivatives in the bulk theory. The numerical values obtained for the thermalization times are summarized in Table 1.

– | – | ||

– | – |

Having smooth fit functions for all the sets of numerical data we can use them to study the thermalization velocities aiming for more details of the nonequilibrium process. These are plotted in Figure 5 (only for the cases fixed and fixed, respectively, to avoid unnecessary repetitions). We notice from the velocity curves the existence of a phase transition point at the middle stage of the thermalization, which divides the process into an accelerating and a decelerating phase. Furthermore, we see that the phase transition point is shifted depending on the values of and . Figure (5) shows that increasing the value of causes a delay in the phase transition point, meaning that the accelerated phase lasts longer for the theory. This is to be contrasted with the fact that the theory is the first to thermalize, indicating that the dynamical process in this case consists of a slowly accelerating phase followed by a quick deceleration towards the equilibrium state. On the other hand, Figure (5) shows that the charge has the opposite effect, i.e., as increases the phase transition point arrives earlier. In other words, for large values of (or in the boundary field theory) the thermalization process consists of a quick accelerating phase followed by a slowly decelerating phase to the final state.

It should be stressed that in this work, in contrast to the authors of [26], we find no evidence for a negative thermalization velocity at initial times. They argue that the velocity should be negative in the very beginning of the evolution corresponding to a“quantum”stage of the nonequilibrium process, which soon becomes“classical”once the velocity becomes positive. In our case (Figure 5) all the thermalization velocities start from zero and increase monotonically until the phase transition point, indicating that nothing particularly odd seems to happen at the initial stages of the thermalization process. This is also clear from the comparison of the numerical results with the fitting functions presented in Figure 6. The zoomed region shows that the numerical points sit all over a horizontal line for initial times (up to ) and therefore there is no reason for a non-vanishing slope at such stage. For that reason, we use for our fit functions degree polynomials with the first powers of () set to zero in order to ensure the strictly constant behavior up to . This allows us to make an accurate fit of the whole set of numerical points, which after all are the ones carrying the physical information.

### 4.2 Renormalized minimal area surfaces

In this section, we numerically solve the equations of motion (32) for the minimal area surfaces in order to track their time evolution. The strategy will follow closely that of the renormalized geodesic lengths done in the previous subsection, so we will not repeat all the details on the numerical procedure as well as the fixing of free parameters since they are essentially identical.

Again we choose the test values and for the numerical analysis. The procedure is the same as before, i.e., for a given pair we solve the geodesic equations (32) subject to the boundary conditions (24) characterized by a pair of values . We repeat this for various pairs of and collect just those yielding a Wilson loop with side (keep in mind that the other side does not influence in our analysis).^{7}^{7}7Remember that is determined from the numerics, via equation (21), so we again use he same criterion of for what we mean by “”. Each of the collected solutions will correspond to a different stage of the time evolution determined by calculating the boundary time via equation (21). Then we integrate each of the collected solutions using equation (31), subtract the universal divergent part to obtain , and construct a list of points to be plotted against . Here is a dimensionless quantity independent of the dimensions of the boundary Wilson loop and is the corresponding thermal value.

A sequence of snapshots of the time evolution of the minimal area surfaces as well as the shell of charged dust described by the Vaidya-BI-AdS metric is shown in Figure 7 for different values of the inverse BI parameter and a fixed charge . Each column, from top to bottom, follows the time evolution for a given value of . Again we can see that up to the value of has little effect on the dynamics, while at the final stages of the evolution plays a decisive role. This is clear from the bottom row at , where we see that the black brane has already formed while the other two are about to form. That illustration suggests that increasing the value of makes the black hole form earlier, which indeed will be confirmed below.

The thermalization curves for the renormalized minimal area surfaces are shown in Figures 8 and 9 for varying at fixed values of and varying at fixed charges, respectively. All the curves are polynomial fits of the numerical points (see below) and the zero point of the vertical axis corresponds to the final state of the process, i.e., the static Einstein-BI black brane fully formed. We immediately notice that all the effects are less evident than those displayed before for the geodesic lengths (the reason why we show the insets in Figure 9) due to our choice of as the characteristic scale in the boundary, in contrast to the used for the geodesics. This just illustrates the argument made in the beginning that the holographic thermalization is a top-down process. Figure 8 shows that for a given the effect of the charge is to delay the thermalization process as it is increased, that is to say, as the chemical potential in the dual field theory grows, the time needed to reach the thermal state also raises. This reinforces the conclusions drawn from the analysis of the renormalized geodesic lengths in the previous subsection. The effect of the inverse BI parameter for fixed charges can be inferred from Figure 9, namely, the larger is, the shorter the thermalization time is. This implies that the boundary field theory is easier to thermalize in the nonlinear case, which is in perfect agreement with our results from the previous subsection. The numerical values obtained for the thermalization times are summarized in Table 2.

– | – | ||

– | – |

Just as we did before, we also use the fit functions to study the thermalization velocities for the minimal area surfaces, , which are plotted in Figure 10 for the cases fixed and fixed, respectively. The plots confirm the results obtained from the geodesics in what concerns the different stages of the dynamical process. Namely, there is always a phase transition instant at the middle stage of the evolution where the process changes from an accelerating phase to a decelerating phase. Such a phase transition point is reached later as we increase for a fixed charge (see Figure 10), or sooner as the charge is increased for a given (see Figure 10). In other words, the thermalization of expectation values of Wilson loops in the dual field theory consists in a slow (quick) accelerating phase followed by a quick (slow) decelerating phase towards the final state as the nonlinearity parameter (the charge , or chemical potential ) is increased.

Once more we point out that our velocity plots do not provide any evidence of a negative thermalization velocity at initial times, in contrast to [26]. This can also be seen from Figure 11, where we show a comparison of the numerical results with the fitting functions. The zoomed region shows that at the initial times, up to , all the numerical points lie over a horizontal line and thus the curve must have vanishing slope in that region. We should mention that the fitting functions used here were degree polynomials with the linear term () set to zero to ensure the strictly constant behavior up to as required by the numerical data.

## 5 Conclusions

The effect of the chemical potential and the inverse BI parameter on the thermalization time of the dual boundary field theory is studied using the Vaidya-like toy model of a collapsing shell of charged dust. The bulk spacetime is dynamical, constructed to interpolate between a pure AdS space at initial times and a charged Einstein-Born-Infeld AdS black brane at late times. We use as thermalization probes the equal-time two-point functions and expectation values of Wilson loops, which have well defined dual gravity descriptions in terms of renormalized geodesic lengths and minimal area surfaces in the bulk. Another class of nonlocal observables, the entanglement entropy of boundary regions, which also have well known holographic descriptions in terms of minimal volumes of codimension-2 surfaces in the bulk, can also be used as a third thermalization probe. However, as the results are similar, we focus only on the first two since they are enough to capture all the relevant effects.

We conclude that as the charge (or, equivalently, the chemical potential) grows, the thermalization time is also increased. The inverse BI parameter, on the other hand, has the opposite effect, i.e., the larger is the value of , the shorter the thermalization time is. At the initial stages of the thermalization we find that and have little effect on the time evolution, becoming important only from the middle stage on. We arrive at the same results independently for both the renormalized geodesic lengths and minimal area surfaces. In each case they can be seen just by looking to a sequence of snapshots of the motion profiles of the geodesics and minimal surfaces as the time goes by or, alternatively, from the thermalization curves obtained from the full numerical analysis. The effect of the charge happens to be the same found in [23] using Einstein-Maxwell theory and [26], where Gauss-Bonnet curvature corrections were included. The outcome of introducing a nonvanishing , on the other hand, is a new result. In particular, the effect is the same as that of the Gauss-Bonnet parameter reported in [26]. Since BI electrodynamics consists essentially of higher order derivatives of the gauge field, we suggest that this may be a general feature of introducing extra derivatives in the bulk, i.e., that boundary gauge theories whose gravity dual carry more than two derivatives tend to thermalize sooner than two-derivative theories after a sudden injection of energy. It would be interesting to explore this idea with other theories possessing this characteristic in order to have definite answers.

Moreover, by fitting the numerical data with smooth functions we were also able to study the thermalization velocities associated to each curve. Although this is not rigorous, it reveals some interesting features of the dynamical process of thermalization and how they are affected by and . We notice the existence of a phase transition point at the middle stage of the thermalization, which divides the process into an accelerating and a decelerating phase. The phase transition point is shifted depending on the values of and . Namely, as increases, the phase transition point arrives earlier. This indicates that for large values of chemical potential the thermalization process consists of a quick accelerating phase followed by a slowly decelerating phase to the final state. Increasing the value of , oppositely, causes a delay in the phase transition point. This is to be contrasted with the fact that nonlinear theories thermalize first, indicating that the dynamical process for non-vanishing consists of a slowly accelerating phase followed by a quick deceleration towards the equilibrium state. We also show from the velocity plots that the thermalization process is monotonic and the velocity is always positive at the initial stages, in contrast to the claim by the authors of [26] that there should be a negative thermalization velocity at the very beginning of the evolution.

Acknowledgments

We would like to thank Gastón Giribet for valuable discussions at the early stages of this work. G. Camilo and E. Abdalla thank CNPq for financial support. E. Abdalla also thanks the support from FAPESP.

## References

- [1] J.M. Maldacena, The large- limit of superconformal field theories and supergravity, Adv. Theor. Math. Phys. 2 (1998) 231 [Int. J. Theor. Phys. 38 (1999) 1113] [hep-th/9711200].
- [2] S.S. Gubser, I.R. Klebanov, A.M. Polyakov, Gauge theory correlators from non-critical string theory, Phys. Lett. B 428 (1998) 105-114 [hep-th/9802109].
- [3] E. Witten, Anti-de Sitter space and holography, Adv. Theor. Math. Phys. 2 (1998) 253 [hep-th/9802150].
- [4] E. Shuryak, Why does the quark gluon plasma at RHIC behave as a nearly ideal fluid?, Prog. Part. Nucl. Phys. 53 (2004) 274 [hep-ph/0312227].
- [5] G. Policastro, D.T. Son, A.O. Starinets, Shear viscosity of strongly coupled supersymmetric Yang-Mills plasma, Phys. Rev. Lett. 87 081601 (2001) [hep-th/0104066].
- [6] D. T. Son, A. O. Starinets, Viscosity, Black Holes, and Quantum Field Theory, Ann. Rev. Nucl. Part. Sci. 57 (2007) 95-118 [arxiv:0704.0240 [hep-th]].
- [7] G. Policastro, D.T. Son, A.O. Starinets, From AdS/CFT correspondence to hydrodynamics, JHEP 09 (2002) 043 [hep-th/0205052].
- [8] S. Bhattacharyya, V.E. Hubeny, S. Minwalla and M. Rangamani, Nonlinear fluid dynamics from gravity, J. High Energy Phys. 0802 (2008) 04 [arXiv:0712.2456 [hep-th]].
- [9] P. K. Kovtun and A. O. Starinets, Quasinormal modes and holography, Phys. Rev. D 72 (2005) 086009 [hep-th/0506184].
- [10] U.W. Heinz, The quark-gluon plasma at RHIC, Nucl. Phys. A 721 (2003) 30-39 [nucl-th/0212004].
- [11] U.H. Danielsson, E. Keski-Vakkuri and M. Kruczenski, Black hole formation in AdS and thermalization on the boundary, JHEP 02 (2000) 039 [hep-th/9912209].
- [12] R.A. Janik and R.B. Peschanski, Gauge/gravity duality and thermalization of a boost-invariant perfect fluid, Phys. Rev. D 74 (2006) 046007 [hep-th/0606149].
- [13] R.A. Janik, Viscous plasma evolution from gravity using AdS/CFT, Phys. Rev. Lett. 98 (2007) 022302 [hep-th/0610144].
- [14] P.M. Chesler and L.G. Yaffe, Boost invariant flow, black hole formation and far-from-equilibrium dynamics in N = 4 supersymmetric Yang-Mills theory, Phys. Rev. D 82 (2010) 026006 [arXiv:0906.4426 [hep-th]].
- [15] S. Bhattacharyya and S. Minwalla, Weak field black hole formation in asymptotically AdS spacetimes, JHEP 09 (2009) 034 [arXiv:0904.0464 [hep-th]].
- [16] E. Caceres, A. Kundu, J.F. Pedraza, D.-L. Yang, Weak Field Collapse in AdS: Introducing a Charge Density, [arXiv:1411.1744 [hep-th]].
- [17] D. Garfinkle, L.A.P. Zayas, Rapid Thermalization in Field Theory from Gravitational Collapse, Phys. Rev. D 84 (2011) 066006 [arXiv:1106.2339 [hep-th]].
- [18] D. Garfinkle, L.A.P. Zayas, D. Reichmann, On Field Theory Thermalization from Gravitational Collapse, JHEP 02 (2012) 119 [arXiv:1110.5823 [hep-th]].
- [19] S. Lin and E. Shuryak, Toward the AdS/CFT gravity dual for high energy collisions. 3. Gravitationally collapsing shell and quasiequilibrium, Phys. Rev. D 78 (2008) 125018 [arXiv:0808.0910 [hep-th]].
- [20] V. Balasubramanian et al., Thermalization of strongly coupled field theories, Phys. Rev. Lett. 106 (2011) 191601 [arXiv:1012.4753 [hep-th]].
- [21] V.Balasubramanian et al., Holographic thermalization, Phys. Rev. D 84 (2011) 026010 [arXiv:1103.2683 [hep-th]].
- [22] R. Baier, A.H. Mueller, D. Schiff, D.T. Son, Bottom-up thermalization in heavy ion collisions, Phys. Lett. B 502 (2001) 51-58 [hep-ph/0009237].
- [23] D. Galante and M. Schvellinger, Thermalization with a chemical potential from AdS spaces, JHEP 07 (2012) 096 [arXiv:1205.1548 [hep-th]].
- [24] E. Caceres and A. Kundu, Holographic Thermalization with Chemical Potential, JHEP 09 (2012) 055 [arXiv:1205.2354 [hep-th]].
- [25] X.-X. Zeng, W.-B. Liu, Holographic thermalization in Gauss-Bonnet gravity, Phys. Lett. B 726 (2013) 481 [arXiv:1305.4841 [hep-th]].
- [26] X.-X. Zeng, X.-M. Liu, W.-B. Liu, Holographic thermalization with a chemical potential in Gauss-Bonnet gravity, JHEP 03 (2014) 031 [arXiv:1311.0718 [hep-th]].
- [27] I. Aref’eva, A. Bagrov and A.S. Koshelev, Holographic thermalization from Kerr-AdS, JHEP 07 (2013) 170 [arXiv:1305.3267 [hep-th]].
- [28] X.-X. Zeng, X.-M. Liu, W.-B. Liu, Holographic thermalization in noncommutative geometry, [arXiv:1407.5262 [hep-th]].
- [29] P. Fonda et al., Holographic thermalization with Lifshitz scaling and hyperscaling violation, JHEP 08 (2014) 051 [arxiv:1401.6088 [hep-th]].
- [30] M. Alishahiha, A.F. Astaneh, M.R.M. Mozaffar, Thermalization in backgrounds with hyperscaling violating factor, [arXiv:1401.2807 [hep-th]].
- [31] S-J. Zhang, B. Wang, E. Abdalla, E. Papantonopoulos, Holographic Thermalization in Gauss-Bonnet Gravity with de Sitter Boundary, [arXiv:1412.7073 [hep-th]].
- [32] W. Baron, D. Galante and M. Schvellinger, Dynamics of holographic thermalization, JHEP 03 (2013) 070 [arXiv:1212.5234 [hep-th]].
- [33] W.H. Baron and M. Schvellinger, Quantum corrections to dynamical holographic thermalization: entanglement entropy and other non-local observables, JHEP 08 (2013) 035 [arXiv:1305.2237 [hep-th]].
- [34] V.Balasubramanian et al., Inhomogeneous holographic thermalization, JHEP 10 (2013) 082 [arXiv:1307.7086 [hep-th]].
- [35] V. Balasubramanian et al., Thermalization of the spectral function in strongly coupled two dimensional conformal field theories, JHEP 04 (2013) 069 [arXiv:1212.6066 [hep-th]].
- [36] J. Abajo-Arrastia, J. Aparicio and E. Lopez, Holographic Evolution of Entanglement Entropy, JHEP 11 (2010) 149 [arXiv:1006.4090 [hep-th]].
- [37] J. Aparicio and E. Lopez, Evolution of Two-Point Functions from Holography, JHEP 12 (2011) 082 [arXiv:1109.3571 [hep-th]].
- [38] T. Albash and C.V. Johnson, Evolution of Holographic Entanglement Entropy after Thermal and Electromagnetic Quenches, New J. Phys. 13 (2011) 045017 [arXiv:1008.3027 [hep-th]].
- [39] V. Keranen, E. Keski-Vakkuri and L. Thorlacius, Thermalization and entanglement following a non-relativistic holographic quench, Phys. Rev. D 85 (2012) 026005 [arXiv:1110.5035 [hep-th]].
- [40] V.E. Hubeny, M. Rangamani and E. Tonni, Thermalization of Causal Holographic Information, JHEP 05 (2013) 136 [arXiv:1302.0853 [hep-th]].
- [41] S.A. Stricker, Holographic thermalization in N = 4 Super Yang-Mills theory at finite coupling, Eur. Phys. J. C 74 2727 (2014) [arXiv:1307.2736 [hep-th]].
- [42] H. Ebrahim and M. Headrick, Instantaneous Thermalization in Holographic Plasmas, [arXiv:1010.5443 [hep-th]].
- [43] M. Born and L. Infeld, Foundations of the New Field Theory, Proc. R. Soc. A 144 (1934) 425.
- [44] E.S. Fradkin and A.A. Tseytlin, Non-linear Electrodynamics From Quantized Strings, Phys. Lett. B 163 (1985) 123.
- [45] A. Abouelsaood, C.G. Callan, C.R. Nappi and S.A. Yost, Open Strings in Background Gauge Fields, Nucl. Phys. B 280 (1987) 599.
- [46] R.G. Leigh, Dirac-Born-Infeld action from Dirichlet sigma model, Mod. Phys. Lett. A 4 (1989) 2767.
- [47] S.P. de Alwis, Supergravity, the DBI action and black hole physics, Phys. Lett. B 435 (1998) 31 [hep-th/9804019].
- [48] G.W. Gibbons, Aspects of Born-Infeld theory and string/M-theory, Rev. Mex. Fis. 49S1 (2003) 19 [AIP Conf. Proc. 589 (2001) 324] [hep-th/0106059].
- [49] S. Banerjee, A Note on Charged Black Holes in AdS space and the Dual Gauge Theories, Phys. Rev. D 82 (2010) 106008 [arXiv:1009.1780 [hep-th]].
- [50] F. Nogueira, J.B. Stang, Density versus chemical potential in holographic field theories, Phys. Rev. D 86 (2012) 026001 [arXiv:1111.2806 [hep-th]].
- [51] R.-G. Cai and Y.-W. Sun, Shear viscosity from AdS Born-Infeld black holes, JHEP 09 (2008) 115 [arXiv:0807.2377 [hep-th]].
- [52] H.S. Tan, First order Born-Infeld hydrodynamics via gauge/gravity duality, JHEP 04 (2009) 131 [arxiv:0903.3424 [hep-th]].
- [53] Z. Zhao, Q. Pan, S. Chen, J. Jing, Notes on holographic superconductor models with the nonlinear electrodynamics, Nucl. Phys. B 871 (2013) 98 [arXiv:1212.6693 [hep-th]].
- [54] W. Yao, J. Jing, Holographic entanglement entropy in insulator/superconductor transition with Born-Infeld electrodynamics, JHEP 05 (2014) 058 [arXiv:1401.6505 [hep-th]].
- [55] W. Yao, J. Jing, Holographic entanglement entropy in metal/superconductor phase transition with Born-Infeld electrodynamics, Nuclear Physics B 889 (2014) 109 [arXiv:1408.1171 [hep-th]].
- [56] W. Yao, J. Jing, Analytical study on holographic superconductors for Born-Infeld electrodynamics in Gauss-Bonnet gravity with backreactions, JHEP 05 (2013) 101 [arXiv:1306.0064 [hep-th]].
- [57] Q. Pan, J. Jing, B. Wang, Holographic superconductor models with the Maxwell field strength corrections, Phys. Rev. D 84 (2011) 126020 [arXiv:1111.0714 [hep-th]].
- [58] Y. Liu, B. Wang, Perturbations around the AdS Born-Infeld black holes, Phys. Rev. D 85 (2012) 046011 [arXiv:1111.6729 [hep-th]].
- [59] T.K. Dey, Born-Infeld black holes in the presence of a cosmological constant, Phys. Lett. B 595 (2004) 484 [hep-th/0406169].
- [60] R.-G. Cai, D.-W. Pang and A. Wang, Born-Infeld black holes in (A)dS spaces, Phys. Rev. D 70 (2004) 124034 [hep-th/0410158].
- [61] A. Chamblin, R. Emparan, C. Johnson and R. Myers,Charged AdS black holes and catastrophic holography, Phys. Rev. D 60 (1999) 064018 [hep-th/9902170].
- [62] E. Caceres, A. Kundu, J.F. Pedraza, W. Tangarife, Strong Subadditivity, Null Energy Condition and Charged Black Holes, JHEP 01 (2014) 084 [arXiv:1304.3398 [hep-th]].
- [63] V. Balasubramanian and S.F. Ross, Holographic particle detection, Phys. Rev. D 61 (2000) 044007 [hep-th/9906226].
- [64] J.M. Maldacena, Wilson loops in large-N field theories, Phys. Rev. Lett. 80 (1998) 4859 [hep-th/9803002].
- [65] E. Abdalla, R. Mohayaee, A. Zadra, Screening in two-dimensional QCD, Int.J.Mod.Phys. A12 (1997) 4539-4558 [hep-th/9604063].
- [66] E. Abdalla, R. Banerjee, Screening in three-dimensional QED, Phys.Rev.Lett. 80 (1998) 238-240 [hep-th/9704176].