Holographic thermalization with a chemical potential in Gauss-Bonnet gravity

Holographic thermalization is studied in the framework of Einstein-Maxwell-Gauss-Bonnet gravity. We use the two-point correlation function and expectation value of Wilson loop, which are dual to the renormalized geodesic length and minimal area surface in the bulk, to probe the thermalization. The numeric result shows that larger the Gauss-Bonnet coefficient is, shorter the thermalization time is, and larger the charge is, longer the thermalization time is, which implies that the Gauss-Bonnet coefficient can accelerate the thermalization while the charge has an opposite effect. In addition, we obtain the functions with respect to the thermalization time for both the thermalization probes at a fixed charge and Gauss-Bonnet coefficient, and on the basis of these functions, we obtain the thermalization velocity, which shows that the thermalization process is non-monotonic. At the middle and later periods of the thermalization process, we find that there is a phase transition point, which divides the thermalization into an acceleration phase and a deceleration phase. We also study the effect of the charge and Gauss-Bonnet coefficient on the phase transition point.


Introduction
AdS/CFT correspondence [1] has been proved to be an effective tool to deal with the strong coupled systems. Especially in recent years, this correspondence has also been used extensively in non-equilibrium strong coupled systems [2,3,4,5,6,7]. One of the biggest motivation to use this gauge gravity duality to study the non-equilibrium strongly coupled systems is that the thermalization time of quark gluon plasma produced in RHIC and LHC experiments predicted by the perturbation theory is longer than the experiment result. The reason for this difference is that the perturbation theory [8] treats the thermalization process as a nearequilibrium process, where the static black hole in the bulk is dual to the boundary system in equilibrium with finite temperature. However, the sudden injection of energy in the thermalization process is a far-from-equilibrium behavior of strongly coupled systems, the holographic bulk thus should be a highly dynamical spacetime, which can be described as black hole formation or black hole merger. Based on this paradigm, there have been many models to study the far-from-equilibrium thermalization behaviors [9,10,11,12,13,14,15,16,17] recently. Especially in [18,19], Balasubramanian et al. find that one can use the two-point correlation function, Wilson loop, and entanglement entropy, which can further be evaluated in the saddle point approximation in terms of geodesic, minimal surface, and minimal volume individually, to detect the thermalization, where the initial state in the conformal field theory is dual to the AdS boundary in a higher dimensional space time, the sudden injection of energy is dual to the collapse of a thin shell of dust, and the final equilibrium state is dual to a static black brane. It is found that the holographic thermalization always proceeds in a top-down pattern, namely the UV modes thermalize firstly, followed by the IR modes 1 . They also find that there is a slight delay in the onset of thermalization and the entanglement entropy thermalizes slowest, which sets a timescale for equilibration. Later, such an investigation is generalized to the bulk geometry with electrostatic potential [20,21,22] and high curvature corrections [23,24,25] to see how the chemical potential and correction parameter affect the thermalization time in the boundary field theory, other extensions on this topic please see [26,27,28,29]. Made by Balasubramanian et al., there are further some elegant extensions on holographic thermalization very recently. Firstly the time-dependent spectral functions in conformal field is found to be tractable [30]. With it, many quantities of interest in the thermalization process, e.g., the time dependence of the occupation number of field modes in the boundary gauge theory, can be calculated. Secondly, the holographic thermalization is extended to the inhomogeneous case by considering an AdS 4 weak field inhomogeneous collapse [31,32]. It is found that the AdS description of the early time evolution is well-matched by free streaming, and the stress tensor approaches that of second order hydrodynamics near the end of the early time interval.
The purpose of this paper is to investigate holographic thermalization in the bulk with curvature corrections and a gauge potential. By holography, curvature correction corresponds to 1 N or 1 λ correction [33,34] to the boundary field theory 2 , and the gauge potential corresponds to a chemical potential in the boundary field theory. In Einstein gravity, there have been some works to study the effect of the chemical potential on the thermalization [20,21]. And there are also some works to study the effect of correction parameter of the high order curvature on the thermalization [23,24,25]. Here we want to explore whether the chemical potential and the correction parameter have the same effect on the thermalization time. If not, our model will provide theoretically a wider range of the thermalization time. To observe the thermalization process in the dual conformal field theory, we take the two-point correlation function and expectation value of Wilson loop as thermalization probes to study the thermalization behavior 3 . According to the AdS/CFT correspondence, this process equals to probing the evolution of a shell of charged dust that interpolates between a pure AdS and a charged Gauss-Bonnet AdS black brane by making use of the renormalized geodesic length and minimal area surface. Concretely we first study the motion profile of the geodesic and 1 On the weakly coupled side classical calculations have shown that the thermalization process is of the bottom-up type, i.e. low energetic modes reach thermal equilibrium first [8]. 2 According to the AdS/CFT correspondence, it is well known that IIB string theory on AdS5 ×S 5 background is dual to D = 4 N = 4 SU (Nc) super Yang-Mills theory. In the limits of large Nc and large 't Hooft coupling, the SYM theory is dual to IIb supergravity which is low energy effective theory of superstring theory. So the higher curvature terms in the bulk theory can arise as next to leading order corrections in the 1/N (large N) expansion of the boundary CFTs in the strong 't Hooft coupling limit. 3 Usually, one also can use the entanglement entropy to detect the thermalization process. Recently there have been many works to study holographic entanglement entropy in higher derivative gravitational theories [35,36,37,38,39,40,41]. Especially in [25], the author found that the entanglement entropy has similar behavior as the other observables during the thermalization process. For simplicity, we will not study entanglement entropy in this paper. minimal area, and then the renormalized geodesic length and minimal area surface in the charged Gauss-Bonnet Vaidya AdS black brane. When we study the effect of the gauge potential on the thermalization process, the Gauss-Bonnet coefficient is fixed and when we are interested in the effect of Gauss-Bonnet coefficient, the gauge potential is fixed. Our result shows that larger the Gauss-Bonnet coefficient is, easier the dual boundary system thermalizes, while larger the charge is, harder the dual boundary system thermalizes. That is, the Gauss-Bonnet coefficient have an opposite effect on the thermalization time compared with the gauge potential. In addition, we also obtain the analytical functions of the renormalized geodesic length with respect to the thermalization time as well as the renormalized minimal area surface with respect to the thermalization time. Based on the functions, we get the thermalization velocity. Our result shows that the velocity is negative at the initial time and positive at the middle and later periods, which indicates that the thermalization process is non-monotonic. We find that there is a phase transition point for the thermalization velocity at a fixed Q and α, which divides the thermalization process into an acceleration phase and a deceleration phase. Effect of Q and α on the phase transition points is also studied. The remainder of this paper is organized as follows. In the next section, we shall provide a brief review of the charged Vaidya AdS black brane in Gauss-Bonnet gravity. Then the holographic setup for non-local observables will be explicitly constructed in Section 3. Especially, we show that the relation between the renormalized length and the two-point correlation function is still valid in the Gauss-Bonnet gravity. Resorting to numerical calculation, we perform a systematic analysis of how the Gauss-Bonnet coefficient and chemical potential affect the thermalization time and thermalization velocity in Section 4. We end up with some discussions in the last section.

Charged Vaidya AdS black branes in Gauss-Bonnet gravity
The D(D ≥ 5) dimensional Einstein-Maxwell theory with a negative cosmological constant and a Gauss-Bonnet term 4 representing a quadratic curvature correction is given by where G D is the D-dimensional gravitational constant, R is the Ricci scalar, Λ is the negative cosmological constant, α is the Gauss-Bonnet coefficient, and From the action in Eq.(2.1), many exact solutions have been found [44,45,46,47,48,49,50,51], here we are interested in the D-dimensional charged black brane solution 3) 4 The action in Eq.(2.1) is not a consistent truncation of an effective higher derivative gravitational action, it is one special case where the higher order corrections of the Maxwell field have vanishing coefficients. One also can take into account the corrections to the Maxwell field such as in [42,43]. where where V n is the volume of the unit radius sphere S D−2 . The U (1) gauge potential reads in which r h is the event horizon radius that is characterized by F (r h ) = 0 and u is the electrostatic potential, which can be identified with the chemical potential of the dual field theory. The Hawking temperature of the charged Gauss-Bonnet AdS black brane reads which can be viewed as the temperature of the dual conformal field theory on the AdS boundary. On the other hand, as r approaches to infinity, one can see the above black brane metric changes into (2.9) Thus this black brane solution is asymptotically AdS with AdS radius ℓ ef f . To get a Vaidya type evolving black brane, we would like first to make the coordinate transformation z = ℓ 2 r , in this case the black brane metric in Eq.(2.3) can be cast into where Then by introducing the Eddington-Finkelstein coordinate system, namely Eq.(2.10) changes into Now the charged Gauss-Bonnet Vaidya AdS black brane can be obtained by freeing the mass parameter and charge parameter as arbitrary functions of v [52,53,54]. This black brane is a solution of the equations of motion of the total action I total = I + I matter , (2.14) where I matter is the action for external matter fields. As one can show, such a metric is sourced by the null dust with the energy momentum tensor flux and gauge flux [52,53,54] 8πG D T matter where the dot stands for derivative with respect to coordinate v, M (v) and Q(v) are the mass and charge of a collapsing black brane 5 . It is obvious that for the charged Vaidya AdS black branes in Gauss-Bonnet gravity, the energy-momentum tensor depends on not only M (v) but also Q(v). In this case, Eq. (2.13) describes the collapse of a thin-shell of charged dust from the boundary toward the bulk interior of asymptotically anti-de Sitter spaces as M and Q are substituted by M (v) and Q(v).

Holographic thermalization
In this section, we are going to investigate the thermalization process of a class of strongly coupled system. According to the AdS/CFT correspondence, the rapid injection of energy on the boundary corresponds to the collapse of a black brane in the AdS space. So to describe the thermalization process holographically, one should choose the mass M (v) and charge Q(v) properly so that it can describe the evolution of the charged dust. It was found that this properties can be achieved by setting the mass parameter and charge parameter as [20,21], where θ(v) is the step function. In this case, in the limit v → −∞, the background corresponds to a pure AdS space while in the limit v → ∞, it corresponds to a charged Gauss-Bonnet AdS black brane. For the convenience of numerical calculations, M (v) and Q(v) are usually chosen as the smooth functions where v 0 represents a finite shell thickness.
Having the construction of a model that describes the thermalization process on the dual conformal field theory, we have to choose a set of extended non-local observables 6 in the bulk which allow us to evaluate the evolution of the system. In this paper, we shall focus mainly on the two-point correlation function at equal time and expectation value of rectangular space-like Wilson loop, which in the bulk correspond to the renormalized geodesic length and minimal area surface respectively. For simplicity but without loss of generality, we shall set the unit such that ℓ = 1 and r h = 1 in the later discussions. In addition, from Eq.(2.11), we know that the mass and the charge have the relation M = 1 + Q 2 , in this case we only need to change the charge to adjust the mass parameter during the numerical process.

Renormalized geodesic length
The relation between the renormalized geodesic length and two-point correlation function at equal time has been discussed extensively [19,56]. Here we will give a review of the relation to check whether it is still valid in the Gauss-Bonnet gravity.
Start with the action of a complex scalar field of mass m in the bulk background gravitational field g ab which gives rise to the bulk propagator can be thought of as the Hamiltonian of a fictitious quantum mechanical model with T the proper time 7 . With this in mind, one can reformulate this quantum mechanical model in terms of path integral over the worldlines of a massive particle. In particular, the above bulk propagator can be expressed as In the saddle point approximation, i.e., Because the local observables such as the energy-momentum tensor and its derivatives can not explore deviations from thermal equilibrium in detail though it provides valuable information about the applicability of viscous hydrodynamics. 7 For the derivation of the Hamiltonian, please see appendix.
the bulk propagator can be evaluated as with X µ (τ ) the classical trajectory satisfying the above equation of motion. Obviously Eq.(3.7) is consistent with the formulation in [19,56] where ∆ is the conformal dimension of scalar operator O, which is similar tom in Eq.(3.7), andL ren indicates the renormalized length of the bulk geodesic between the points (t 0 , x n ) and (t 0 , x ′ n ) on the AdS boundary. In other words, the Gauss-Bonnet coefficient has not effect on the dual relation between the renormalized geodesic length and the two-point correlation function. In what follows, we will make use of Eq.(3.8) to explore how the Gauss-Bonnet coefficient and gauge potential affect the thermalization time. The reason that we use the renormalized geodesic lengthL ren in Eq.(3.7) is that the geodesic length is divergent near the boundary, which is related to the divergent part as (3.9) whereL = √ −dS 2 is the geodesic length between the points (t 0 , x n ) and (t 0 , x ′ n ) on the AdS boundary and z 0 is the IR radial cut-off. Next, we are concentrating on studyingL. Taking into account the spacetime symmetry of our Vaidya type black brane, we can simply let (t 0 , x n ) and (t 0 , x ′ n ) have identical coordinates except x 1 = −ℓ ef f l 2 ≡ −l 2 and x ′ 1 = −ℓ ef f l 2 ≡ −l 2 with l the separation between these two points on the boundary. In order to make the notation as simple as possible, we would like to rename this exceptional coordinate x 1 as x and employ it to parameterize the trajectory such that the proper length in the charged Gauss-Bonnet Vaidya AdS black brane can be given bỹ where the prime denotes the derivative with respect to x and Note that the integrand in Eq.(3.10) can be thought of as the Lagrangian L of a fictitious system with x the proper time. Since the Lagrangian does not depend explicitly on x, there is an associated conserved quantity In addition, based on the Lagrangian L of the system, we obtain So the equation of motion for z(x) can be written as Similarly, the equation of motion for v(x) can be solved as Next, we turn to studying the equations of motion in Eq.(3.15) and Eq.(3.16). Considering the reflection symmetry of our geodesic, we will use the following initial conditions As z(x) and v(x) are solved, we can get the IR radial cut-off and thermalization time by the boundary conditions as follows In addition, with the help of Eq.(3.12) and Eq.(3.18), the proper length of geodesic in Eq.(3.10) can be simplified asL which will be more convenient for us to obtain the renormalized geodesic length.

Minimal area surface
In this section, we are going to study the minimal area surface, which in the dual conformal field theory corresponds to the Wilson loop operator. Wilson loop operator is defined as a path ordered integral of gauge field over a closed contour, and its expectation value is approximated geometrically by the AdS/CFT correspondence as [19,57] where C is the closed contour, Σ is the minimal bulk surface ending on C withÃ ren its renormalized minimal area surface, and α ′ is the Regge slope parameter. In the Gauss-Bonnet gravity, we will assume Eq.(3.21) is still valid as that between the renormalized geodesic length and the two-point correlation function.
Here we are focusing solely on the rectangular space-like Wilson loop. In this case, the enclosed rectangle can always be chosen to be centered at the coordinate origin and lying on the x 1 − x 2 plane with the assumption that the corresponding bulk surface is invariant along the x 2 direction. This implies that the minimal area surface can be expressed as where we have set the separation along x 2 direction to be one and the separation along x 1 to bel with x 2 renamed as y and x 1 renamed as x. As before, from Eq.(3.22) we also can get a Lagrangian L and with it we can find a conserved quantity, i.e., which can simplify our equations of motion as Similarly, with the initial conditions as in (3.18) and the regularization cut-off as in (3.19), the renormalized minimal area surface can be cast intõ Next, we will investigate the evolution of the renormalized minimal area surface with respect to the thermalization time to explore how the gauge potential and Gauss-Bonnet coefficient affect the thermalization process.

Numerical results
In this section, we will solve the equations of motion of geodesic length and minimal area surface numerically, and then explore how the chemical potential and Gauss-Bonnet coefficient affect the thermalization time. Since there have been many works to study the effect of the space time dimensions and boundary separation on the thermalization probes [18,19,20,21], to avoid redundancy, we mainly discuss the case D = 5 and a fixed boundary separation in this paper. During the numerics, we will take the shell thickness and UV cut-off as v 0 = 0.01, z 0 = 0.01 respectively.
According to the AdS/CFT correspondence, we know that the electromagnetic field in the bulk is dual to the chemical potential in the dual quantum field theory, so we will use the electromagnetic field defined in Eq.(2.6) to explore the effect of the chemical potential on the thermalization process in the AdS boundary. However, as stressed in [20,58,59], the chemical potential has energy units in the dual field theory ([u] = 1/[L]) while A µ as defined in Eq.(2.1) is dimensionless, thus one has to redefine the electromagnetic field asÃ µ = A µ /p, where p is a scale with length unit that depends on the particular compactification. In this case,Ã µ and u have the same unit and the chemical potential can be expressed as .
In addition, due to the conformal symmetry on the boundary, the quantity which is physically meaningful is the ratio of u/T in asymptotically charged AdS space time, namely the chemical potential measured with the temperature as the unit. Therefore, through this paper we will use the following ratio to check the effect of the chemical potential on the thermalization time. Obviously, for the case r h = 1, Eq.(4.2) shows that u/T changes from 0 → ∞ provided Q changes from 0 → D−1 D−3 . In other words, to adjust the change of the ratio u/T in all the range, we only need to change Q from 0 → D−1 D−3 . For the case D = 5 in this paper, we will choose Q = 0.00001, 0.5, 1 in our numerical result. On the other hand, we will also consider the effect of the Gauss-Bonnet coefficient on the thermalization time, as in [23] we will take α = −0.1, 0.0001, 0.08 as the constraint of causality of dual field theory on the boundary is imposed 8 .
As the boundary conditions in Eq.(3.18) is adopted, the equations of motion of the geodesic in Eq.(3.15) and Eq.(3.16) for different α and different Q can be solved numerically. When we are interested in the effect of Q on the motion profile of the geodesic, the Gauss-Bonnet coefficient α is fixed, and when we are interested in the effect of Gauss-Bonnet coefficient α, the charge Q is fixed. Since different initial time v ⋆ corresponds to different stage of the motion of the geodesics, we also discuss the effect of v ⋆ on the motion profile. Figure ( Figure (2), we know that as the initial time increases, the shell of the dust approaches to the horizon of the charged Gauss-Bonnet black brane, which means that for the larger initial time, the thermalization has been behaved longer. For different initial time, we also keep a watchful eye on how the Gauss-Bonnet coefficient and charge affect the thermalization time, which are listed in Table (1). From this table, we can observe that for both the two different initial times, as the Gauss-Bonnet coefficient grows, the thermalization time decreases, which means that the quark gluon plasma in the dual conformal theory is easier to be thermalized. But as the charge grows, we find for different initial times, the thermalization time has different variation trends. For the case v ⋆ = −0.856, it is found that for a fixed Gauss-Bonnet coefficient, as the charge grows the thermalization time has a little difference, while for v ⋆ = −0.456, the thermalization time increases as the charge grows. That is, the charge has little effect on the thermalization time at the initial stage of the thermalization. We also can observe the effect of the Gauss-Bonnet coefficient and charge on the motion profile of the geodesics. From Figure (2), we find for a fixed charge, e.g. Q = 0.5, as the Gauss-Bonnet coefficient grows from α = −0.1 to α = 0.08, the shell will approaches to the horizon of the black brane and drops into the horizon lastly. In other words, for α = −0.1, the quark gluon plasma in the conformal field theory is thermalizing while for α = 0.08, it is thermalized. For a fixed Gauss-Bonnet coefficient, e.g. α = 0.08, we find at Q = 0.00001, the shell lies above the horizon while at Q = 1, the shell lies below the horizon.
Having the numerical result of z(x), we can study the renormalized geodesic length with the help of Eq.(3.9) and Eq. (3.20). As done in [20], we compare δL at each time with the final values δL CGB , obtained in a static charged Gauss-Bonnet AdS black brane, i.e. M (µ) = M , Q(µ) = Q. In this case, the thermalized state is labeled by the zero point of the vertical coordinate in each picture. To get an observable quantity that isl independent, we will plot the quantity δL = δL/l. Figure (3) gives the relation between the renormalized geodesic length and thermalization time for different charge at a fixed Gauss-Bonnet coefficient. In each picture, the vertical axis indicates the renormalized geodesic length while the horizontal axis indicates the thermalization time t 0 . For a fixed Gauss-Bonnet coefficient, e.g. α = 0.08, the thermalization time increases as Q raises. This phenomenon has been also observed previously when we study the motion profile of the geodesic. In [20], the effect of charge on the thermalization time is investigated in Einstein gravity, it was shown that there is an enhancement of the thermalization time as the ratio of chemical potential over temperature increases. Obviously, in the Gauss-Bonnet gravity, this phenomenon is not changed. In addition, From the same color line, e.g. green line, in (a), (b) and (c) in Figure (3), we know that as the Gauss-Bonnet coefficient increases, the thermalization time decreases for a fixed charge. In Figure (4), we plot this graphics for different charge. Note that [23] has investigated this phenomenon for the case Q = 0. It was found that for a fixed boundary separation there is always a time range in which the renormalized geodesic length takes the same value nearly. That is, during that time range, the Gauss-Bonnet coefficient has little effect on the renormalized geodesic length. Obviously, (a) in Figure (4) is consistent with their result. For Q = 0.5 and Q = 1, we also can observe this phenomenon, which are plotted in (b) and (c) in Figure (4).
Interestingly, we find the thermalization curve for a fixed charge and Gauss-Bonnet coefficient in Figure ( 3) For small time, the function is determined by the lower power of t 0 , while for large time it is determined by the higher power of t 0 . With the function, we can get the thermalization velocity, which is plotted in Figure (6). From this figure, we can observe two interesting phenomena for the thermalization process. One is that the thermalization velocity is negative at the initial time while it is positive at the middle and later periods, which means that the thermalization process is non-monotonic. In fact, this non-monotonic behavior has also been observed in [21]. The reason for this behavior is that at the initial time the thermalization is "quantum" while at the later periods, it is "classical". The author in [21] further argued that at the "quantum" stage the slope of the thermalization curve, namely the thermalization velocity, is negative and at the "classical" stage the slope is positive. Obviously, our result plotted in Figure (6) confirms their argument. The other is that there is a phase transition point for the thermalization velocity, which divides the thermalization into an acceleration phase and a deceleration phase. (a) and (b) in Figure (6) represent respectively the effect of Q and α on the phase transition points. From both figures, we know that the phase transition points vary for different Q and α. Especially in (b), we can observe clearly that as α increases, the phase transition points shift left. At a fixed Q and α, we can get the value of the phase transition point. For the case Q = 0.5 and α = 0.0001, it is easy to find that in the time range 0 < t 0 < 1.04978, the thermalization is an acceleration process while for t 0 > 1.04978 it is a deceleration process before it approaches to the equilibrium state. Surely as the slope of this velocity curve is produced, we also can get the values of the acceleration and deceleration.
Adopting similar strategy, we also can study the motion profile of minimal area as well  as the change of the renormalized minimal area surface to thermalization time. Based on the motion equations in (3.24) and the boundary conditions in (3.18), the numerical solution of z(x) can be produced. In this case, we can get the motion profile of minimal area for different charge Q and Gauss-Bonnet coefficient α, which is shown in Figure (7). From this figure, we know that for a fixed charge, e.g. the first row, as the Gauss-Bonnet coefficient increases, the shell surface approaches to the horizon surface step by step, which means the thermalization is faster. For a fixed Gauss-Bonnet coefficient, e.g. the third column, as the charge increases, the shell surface is removed from the horizon surface step by step, which means the thermalization is slower. The thermalization time for different α and Q have been listed in Table ( 2). It is shown that for a fixed charge the thermalization time decreases as α becomes larger, while for a fixed α, the thermalization time increases as Q becomes larger. That is, the charge has an inverse effect compared with the Gauss-Bonnet coefficient on the thermalization time. This phenomenon is similar to that of the geodesics.  and δA CGB is the renormalized minimal area surface for a charged Gauss-Bonnet AdS black brane. The relation between the renormalized minimal area surface and thermalization time for different charge Q is given in Figure (8) for a fixed Gauss-Bonnet coefficient α, in which the vertical axis indicates the renormalized minimal area surface while the horizontal axis indicates the thermalization time t 0 . For a fixed α, we find larger the charge Q is, longer the thermalization time is. That is to say, as the chemical potential in the dual field theory increases, the thermalization time rises too, which is the same as that obtained by studying the motion profile of the minimal area. In addition, for a fixed charge, we also study the effect of the Gauss-Bonnet coefficient on the thermalization time, which is shown in Figure  (9). It is obvious that larger the Gauss-Bonnet coefficient is, shorter the thermalization time is, which means that the quark gluon plasma is easier to thermalize. This behavior is similar to that of the geodesic which is given in Figure (3). As the case of the renormalized geodesic length, we find there is also an overlapped region for the case Q = 0.00001, 0.5, 1 respectively in Figure (9). We also can get the functions of the renormalized minimal area surface with respect to the thermalization time. At α = 0.08, the functions of the thermalization curve   for Q = 0.00001, 0.5, 1 can be expressed respectively as the thermalization velocity is plotted in Figure ( 11). From this figure, we also can observe the non-monotonic behavior and phase transition point of the thermalization process, which is similar to that of the renormalized geodesic length. We also can obtain the acceleration phase and deceleration phase for this case. By the slope of the thermalization velocity, we know that the acceleration range is 0 < t 0 < 1.106 and the deceleration range is 1.106 < t 0 < 1.481. Adopting the same strategy, we also can get the phase transition points for other Q and α.

Conclusions
Effect of the chemical potential and correction parameter on the thermalization in the dual boundary field theory is investigated by considering the collapse of a shell of charged dust that interpolates between a pure AdS and a charged Gauss-Bonnet AdS black brane. The two-point function and expectation values of Wilson loop are chosen as the thermalization probes, which are dual to the renormalized geodesic length and minimal area surface in the bulk. We first study the motion profiles of the geodesic and minimal surface and find that larger the Gauss-Bonnet coefficient is, shorter the thermalization time is, and larger the charge is, longer the thermalization time is. At the initial stage of the thermalization, we find that the charge has little effect on the thermalization time. We reproduce this result by studying the relation between the renormalized geodesic length and time as well as the renormalized minimal surface area and time respectively. In addition, we also find the functions of the thermalization probes with respect to the thermalization time by fitting the numerical result. Though this is a naive test, we still can get some useful information. Firstly, we get the thermalization velocity for a fixed charge and Gauss-Bonnet coefficient. From the velocity curve, on one hand, we know that the thermalization process is non-monotonic, and on the other hand we find there is a phase transition point, which divides the thermalization into an acceleration phase and a deceleration phase. Secondly according to the slope of the velocity curve, we also get the acceleration and deceleration of the thermalization. Thus our investigation provides a more accurate description of the thermalization process. Note that in our naive test, we can not see how the charge and Gauss-Bonnet coefficient affect the thermalization quantificationally since the coefficients of the functions are not determined. In future, we expect to find an analytical formulism to study the thermalization so that we can get more useful information on the thermalization process.
where ǫ is a small quantity, to be taken to go to zero in the later calculation. Now for convenience, we would like to resort to the Rienmann normal coordinate at x, where the behavior of metric is simplified as g µν (x) = η µν , ∂ ρ g µν (x) = 0, ∂ ρ ∂ σ g µν (x) = − (1 − ǫ ∂ ∂T + · · ·)U (x, y; T ).