Thermal chaos of a charged dilaton-AdS black hole in the extended phase space

We have studied thermal chaotic behavior in the extended phase space for a charged dilaton-AdS black hole by Melnikov method and present the effect of dilaton parameter on the thermal chaos. Our result show that for the temporal perturbation the thermal chaos in the charged dilaton-AdS black hole occurs only if the perturbation amplitude is larger than certain a critical value, but for the spatially perturbation, the chaos always exists irrespective of perturbation amplitude. These behaviors are similar to those in other AdS black hole, which can be regarded as the common features of the static AdS black holes. Moreover, we also find that the critical temporal perturbation amplitude leading to chaos increases with the dilaton parameter and decreases with the charge. This means that under the temporal perturbation the presence of dilaton parameter makes the onset of chaos more difficult, which differs from that of the charge parameter.


I. INTRODUCTION
Chaos is a kind of very complicated and irregular motion, which appears only in the non-linear and nonintegrable dynamical systems. It is highly sensitive to the initial conditions and presents the intrinsic random in systems [1][2][3]. The chaotic phenomenon can be detected by many methods including the Poincaré surfaces of section, the Lyapunov characteristic exponents, the fractal basin boundaries, the bifurcation diagram, the Melnikov method and so on. Although the chaotic motion is very complex, it is very common in nature.
In black hole physics, the chaotic orbits of particles have been found to exist in multi-black hole spacetime [4,5], in the Manko-Novikov black hole spacetime [6][7][8][9][10], in the perturbed Schwarzschild spacetime [11][12][13][14] , in the spacetime of a black hole immersed an external magnetic field [15] and in the accelerating and rotating black holes spacetime [16]. Moreover, it is of interest to find that the chaotic behaviors appear in ring string dynamics in the Schwarzschild black hole spacetime [17] and in AdS black hole spacetimes [18,19].
The chaotic behavior in van-der Waals fluid [20] has been found by applying the Melnikov method [21].
The fluid system is initially supposed to be in the unstable spinodal region, which allows the existence of the homoclinic orbit in phase space. It is shown that for the temporal perturbation the thermal chaos occurs only if the perturbation amplitude is larger than certain a critical value, but for the spatially perturbation, the chaos always exists irrespective of perturbation amplitude. Recent investigations indicate that thermodynamic behavior of AdS black holes in the extended phase space bears high resemblance to that of van-der Waals fluid system. By treating the cosmological constant as a thermodynamic pressure and its conjugate quantity as a thermodynamic volume, Kubiznak et al. [22] found that in the system of a charged AdS black hole the smalllarge black hole phase transition possesses the same critical behavior of liquid-gas phase transitions in the van der Waals fluid. The similar P − V critical behaviors are disclosed in the other AdS black hole spacetimes [23][24][25][26]. These similar critical behaviors mean that there must be some deep connection between AdS black holes and van-der Waals fluid system. Thus, it is very natural to probe the thermal chaos in AdS black holes with phase structures similar to that of van-der Waals fluid. For a Reissner-Nordström-AdS black hole, it is found the similar chaotic behaviors existed in the system of van-der Waals fluid. The critical temporal perturbation amplitude depend on the black hole charge Q and the presence of charge makes the onset of chaos easier [27]. For a Gauss-Bonnet AdS black hole [28], the presence of charge is necessary for chaos under temporal perturbations. However, under spatial perturbations, chaotic behavior always exist irrespective of whether the black hole carries charge or not. Moreover, the chaos under temporal perturbation is found to occur easier in the case with larger Gauss-Bonnet coupling constant. The same qualitative properties of thermal chaos are also found in the AdS black hole spacetime with Born-Infeld electrodynamics [29], which shows that the effect of Born-Infeld parameter on thermal chaos is similar to that of black hole charge.
In this paper, we will study thermal chaos in the extended phase space for a charged dilaton-AdS black hole, which belongs to a family of solutions in Einstein-Maxwell-dilaton gravity theory [30][31][32]. The presence of the dilaton field changes the causal structure of the spacetime and modifies the thermodynamic properties of the black holes. In an extended phase space, one can find that the thermodynamic quantities depend on the dilaton parameter [33]. For example, the thermodynamic pressure P depends on both the cosmological constant and the coupling strengthen between the dilaton scalar and electromagnetic fields. Correspondingly, the thermodynamic volume V is also a function of the coupling constant, which is different from those in the usual AdS black hole spacetimes in which P is determined entirely by the cosmological constant and V is only a function of the event horizon radius. Although the thermodynamic quantities depend on the dilaton parameter, it is found that the P − V critical behaviors present the similar feature as in usual AdS one and the critical exponents are independent of the details of the dilaton system [33]. And then, it is very natural to ask whether dilaton parameter affects thermal chaos of AdS black holes in Einstein-Maxwell-dilaton gravity theory. The main purpose of this paper is to study thermal chaos of a charged dilaton-AdS black hole [30] and to probe the effects of dilaton parameter on thermal chaotic behaviors. This paper is organized as follows. In Sec.II, we will review briefly thermodynamics of a charged dilaton-AdS black hole [30] in the extended phase space and dependence of the critical temperature T c on black hole parameters. In Sec. III, we will study the thermal chaos of the charged dilaton-AdS black hole flow under thermal perturbations and probe the effects of dilaton parameter on such kind of thermal chaos. Finally, we present a brief summary.

II. THERMODYNAMICS OF A CHARGED DILATON-ADS BLACK HOLE IN EXTENDED PHASE SPACE
Let us now review briefly the thermodynamics of a charged dilaton-AdS black hole in extended phase space.
In a four-dimensional spacetime, the action of Einstein-Maxwell theory with a dilation field can be expressed as [30][31][32][33] A ν . The quantity b is a positive arbitrary constant and α is the coupling parameter between the dilaton and Maxwell fields. From the action (1), one can obtain a spherical symmetric black hole solution, whose metric has the form [30] with Here γ is related to the coupling parameter α by γ = α 2 α 2 +1 . The parameter m is associated with the ADM mass M of black hole (2) by m = 2(α 2 + 1)b −2γ M and q is the black hole charge. In the absence of the dilaton field (α = 0), the above metric describes the geometry of the usual Reissner-Nordström AdS black hole.
For the black hole (2), with the event horizon radius r + (i.e., the largest root of equation f (r) = 0), one can obtain Hawking temperature and entropy respectively. Introducing the thermodynamic pressure P , its conjugate quantity as volume V and the electric potential U one can find that the first law of thermodynamics and the corresponding Smarr formula can be expressed as respectively. In the limit α → 0, the above thermodynamic pressure P and volume V become which are consistent with that of Reissner-Nordström AdS black hole. Similarly, defining the specific volume one can find that thermodynamic pressure P can be rewritten as The rich phase structures for the black hole (2) have been studied in the extended phase space [33]. It is shown in Fig. (1) that there exists Large-Small black hole phase transition, which is qualitatively similar to the gas-fluid phase transition in the van-der Waals system. This is a second order phase transition occurred at the critical temperature Obviously, the critical temperature T c decreases with the parameters b and q. The change of T c with α depends on the value of q, which is shown in Fig.(2). The critical temperature T c for different b increases monotonically with α in the case with the larger q, but it first decreases and then increases in the case with the smaller q. The diagrams in Fig. (1) indicate that as the temperature of black hole T is less than the critical temperature T c the P − v curve can be divided into two stable regions and a unstable one. Two stable correspond to the small black hole region and the large black hole one, respectively. The two points v α and v β are determined by ∂P (v,T0)

III. CHAOS IN THE CHARGED DILATON-ADS BLACK HOLE FLOW UNDER THERMAL PERTURBATIONS
As in refs. [20,[27][28][29], for a sake of simplicity, the charged dilaton-AdS black hole flow is assume to move along x-axis in a tube of unit cross section of fixed volume. And then, the position of a fluid particle can be described by the Eulerian coordinate x. The mass M of a column of fluid of unit cross section between the reference fluid particle with Eulerian coordinate x * and a fluid particle with the coordinate x can be expressed where ρ(x, t) is the fluid density at spatial position x and time t. The relation (14) means that the position x of any particle can be expressed as a function of the mass M and the time t, i.e., x = x(M, t). From Eq. (14), it one can also define the velocity as u(M, t) ≡ x t (M, t). With these quantities, one can find that the charged dilaton-AdS black hole flow can be described by a dynamical system with the balance equation of mass and momentum [20,[27][28][29] ∂v ∂t = ∂u ∂M , where τ is the Piola stress. As in refs. [20,[27][28][29], we assume that the charged dilaton-AdS black hole fluid is thermoelastic, isotropic and slightly viscous, which means that the Piola stress can be expressed as [34] Here A is a positive constant, and µ is a small positive constant viscosity. With the above form of Piola stress τ , the balance equation (16) can be rewritten as One can suppose that a finite black hole fluid tube of unit cross section contains the total mass 2π/s in a volume 2πv 0 /s, where s is a positive parameter. Introducing the change of variablesM = sM,t = st,x = sx and µ = ǫµ 0 , one can find that the range of the massM becomesM ∈ [0, 2π], and then the equation (18) can be written as where ǫ is a perturbation parameter. Here the overbars are omitted for later convenience.
Let us now to study the temporal chaos in the spinodal region in which the small and large black hole phase coexist. The system is assumed to lie initially in the unstable equilibrium state associated with the specific where δ is amplitude of the perturbation relative to the small viscosity. Expanding P (v, T ) around the equilibrium point (v 0 , T 0 ) as in refs. [20,21,[27][28][29], we have with The coefficients P T T (v 0 , T 0 ), P vT T (v 0 , T 0 ) and P T T T (v 0 , T 0 ) vanish because the pressure (12) is a linear function of the temperature T . The absence of P vv (v 0 , T 0 ) is attributed to a fact that the thermodynamics system (12) satisfies ∂ 2 P (v,T0) i.e., [20,21,[27][28][29] v( Here x i (t) and u i (t) can be regarded as hydrodynamical modes which describe the deviation from the initial equilibrium state with v = v 0 . Although a full analysis should consider the effects from the higher modes, it is very difficult to find the homoclinic orbit in the case containing the higher modes, even in the three-mode one [20]. It is argued that the existence of homoclinic orbits for higher mode approximations to the original partial differential equation or for the infinite-dimensional problem itself remains an open problem [20]. The two-mode case is investigated for the van der Waals fluid [20] and for charged AdS black hole [27] and Gauss-Bonnet one [28]. The recent investigations [29] show that the effects of the second mode has no basic difference from that of the first mode. Therefore, we here consider only the first mode (x 1 (t), u 1 (t) ) as in ref. [29] and omit the subscript 1 in the later formulae. In doing so, the dynamical equation (19) can be simplified further asẋ = u, With the denotation z ≡ [x, u] T , the above equations can be rewritten aṡ with and where a 2 ≡ (P v − As 2 ). In the case without thermal perturbation (i.e., ǫ = 0), one can find an analytical solution for the equation (25) [29,35] It is a homoclinic orbit which joins a saddle equilibrium point to itself. The solution (28) owns the two branches, which correspond to two wings of the butterfly-like orbit shown in Fig.(3), respectively. Under the time-periodic thermal perturbation (20) (i.e.,ǫ = 0), the above homoclinic orbit may break so that the possible chaos could appear in this system. The existence of chaos is determined by the Melnikov function, which has a form [36,37] with Combining Eq. (29) with Eqs. (26) and (27), one can obtain the Melnikov function for the charged dilaton-AdS black hole flow With the residue theorem, the Melnikov function can be further expressed as The coefficients K and L are similar to those obtained in ref. [29], but they depend on the dilaton parameter α in this case. After a simple analysis, one can find that M (t 0 ) has a simple zero if the condition is satisfied [38]. This means that the chaos appears if the amplitude of the perturbation δ is larger than δ c = sµ0L ωK . In Fig.(4), we present numerically the time evolution of equations of motion (24) in the phase plane x(t) − u(t) for the charged dilaton-AdS black hole. It is shown that as δ < δ c , the perturbation decays with time and the system will finally approach to its original unstable equilibrium state with v = v 0 . As δ > δ c , one can find that the trajectories in the phase plane become irregular and complex, and the dynamical evolution of the system exhibits chaotic feature. The critical value δ c depends on various the parameters b, q and α. From Fig.(5), we find that δ c decreases with the charge q and the parameter b, but increases with the dilaton parameter α. This means that the larger b or q makes the onset of chaos easier, but the larger α makes chaos more difficult under the temporal perturbation.
We are now in the position to study the thermal chaos a charged dilaton-AdS black hole due to a small spatially periodic perturbation. Firstly, we assume that the black hole is in the equilibrium state with a subcritical temperature T 0 . From van-der Waals-Korteweg theory, the stress tensor without flow (17) becomes where the notation ′ denotes the derivative with respect to x. P is thermodynamic pressure (12) of the charged dilaton-AdS black hole and A is a positive constant. For a static equilibrium state with no body forces, the balance of linear momentum is τ ′ = 0, which means that τ = B = constant. This constant quantity B also represents the ambient pressure at the end of the tube. Thus, the equation (35) for a static equilibrium state can be expressed further as For any fixed temperature T 0 < T c , the non-linear systems in Eq.
one can find that the dynamical equation (36) becomes This second-order differential equation can be rewritten as a pair of first-order differential equations v ′ = u, For the dynamical equation (38), the general solutions describing the homoclinic or heteroclinic orbit can be expressed as and the corresponding functions f (z(x − x 0 )) and g(z(x − x 0 ), x) in Melnikov function becomes Thus, the Melnikov function in this cases can be simplified as Therefore, M (x 0 ) always possesses simple zeros for the arbitrary values of L and N . This means that there always exist chaos for thermodynamic system of the charged dilaton-AdS black hole suffered from a spatially periodic thermal perturbation. This is the same as that obtained in other AdS black holes [20,[27][28][29], which may be a common feature of such kind of thermodynamic systems. In Fig. 7, we chose the homoclinic orbit or heteroclinic orbit as the initial configures and plot the solutions of the perturbed equation (38) in the v − v ′ plane. It is shown that there is spatial chaos under perturbations irrespective of black hole parameters and the perturbation strength.

IV. DISCUSSION AND SUMMARY
We have studied thermal chaotic behavior in the extended phase space of a charged dilaton-AdS black hole by Melnikov method and present the effect of dilaton parameter on the thermal chaos. Under the temporal perturbation in the spinodal region, we find that the chaos occurs only if the perturbation amplitude δ is larger than a critical value δ c , which is similar to those in the RN-AdS, Born-Infeld-AdS and Gauss-Bonnet AdS black holes. The dependence of δ c on black hole parameters show that δ c decreases with the charge q and the parameter b, but increases with the dilaton parameter α. This means that the presence of dilaton parameter makes the onset of chaos more difficult, which differs from those arising from the parameters q and b. For the spatially periodic thermal perturbation, we find that there always exists chaos for thermodynamic system of the charged dilaton-AdS black hole. Comparing it with those of the RN-AdS, Born-Infeld-AdS and Gauss-Bonnet AdS black holes, it is easy to obtain that for the temporal perturbation the thermal chaos occurs only if the perturbation amplitude is larger than certain a critical value. However, for the spatially perturbation, the chaos always exists irrespective of perturbation amplitude. These behavior can be regarded as a common feature of such kind of static AdS black holes.