Validity of Maxwell Equal Area Law for Black Holes Conformally Coupled to Scalar Fields in $\text{AdS}_5$ Spacetime

We investigate the $P-V$ criticality and the Maxwell equal area law for a five-dimensional spherically symmetric AdS black hole with a scalar hair in the absence of and in the presence of a Maxwell field, respectively. Especially in the charged case, we give the exact $P-V$ critical values. More importantly, we analyze the validity and invalidity of the Maxwell equal area law for the AdS hairy black hole in the scenarios without and with charges, respectively. Within the scope of validity of the Maxwell equal area law, we point out that there exists a representative van der Waals-type oscillation in the $P-V$ diagram. This oscillating part that indicates the phase transition from a small black hole to a large one can be replaced by an isobar. The small and large black holes share the same Gibbs free energy. We also give the distribution of the critical points in the parameter space both without and with charges, and obtain for the uncharged case the fitting formula of the co-existence curve. Meanwhile, the latent heat is calculated, which gives the energy released or absorbed between the small and large black hole phases in the isothermal-isobaric procedure.


Introduction
Since the seminal works by Hawking and Bekenstein on the radiation of black holes, the exploration of thermodynamic properties of black holes has received a wide range of attention [1][2][3] and also acquired great progress [4][5][6][7]. Of more particular interest is the thermodynamics of anti-de Sitter (AdS) black holes [8,9] where the AdS/CFT duality plays a pivotal role in recent developments of theoretical physics [10,11]. In the context of AdS/CFT correspondence [12,13], the Hawking-Page phase transition [14] of five-dimensional AdS black holes can be explained as the phenomenon of the confinement/deconfinement transition in the fourdimensional Yang-Mills gauge field theory [15]. Another archetypical example of the AdS/CFT correspondence is the holographic superconductor, which can be regarded as the scalar field condensation around a four-dimensional charged AdS black hole [16].
With the cosmological constant Λ being treated as a thermodynamic pressure variable [17][18][19][20][21][22] and its conjugate variable being considered as the thermodynamic volume, the thermodynamics in the extended phase space has been getting more and more attentions. In this paradigm, the mass of black holes is identified as the enthalpy rather than the internal energy. This idea has also been applied to other known parameters, such as the Born-Infeld parameter [23,24], the Gauss-Bonnet coupling constant [25], the noncommutative parameter [26], and the Horndeski non-minimal kinetic coupling strength [27], etc. All the parameters just mentioned can be regarded as a kind of thermodynamic pressure. Furthermore, there exists a similar situation in the exploration of charged AdS hairy black holes [28] of Einstein-Maxwell theory conformally coupled to a scalar field in five dimensions. The model's action has been given by [29][30][31][32] where κ = 16πG, G is the five-dimensional Newton constant, R the scalar curvature, F the electromagnetic field strength, g µν the metric with mostly plus signatures, and g = det(g µν ).
In fact, one can see that the above model is the most general scalar field/gravity coupling formulation whose field equations are of second order for both gravity and matter. Hence, we can say that this formulation is a generalization of the Horndeski theory [34] whose action contains a non-minimal kinetic coupling of a massless real scalar field and the Einstein tensor. Even more importantly, being a simple and tractable model, it provides a significant advantage for studying the phase transition of hairy black holes on the AdS spacetime where the back-reaction of the scalar field on the metric can be solved analytically in five dimensions. In the paradigm where there exists a complete physical analogy between the four-dimensional Reissner-Nordström AdS black hole and the real van der Waals fluid [19] in the phase transition, the recent research [28] shows that this charged AdS hairy black hole also exhibits the van der Waals-type thermodynamic behavior, and moreover, such a black hole undergoes reentrant phase transition which usually occurs in higher curvature gravity theory. We note that all the interesting results just mentioned are available by making the coupling parameter dynamical, i.e. treating the coupling parameter as a certain thermodynamic variable.
Based on the results pointed out above, we take advantage of the well-established Maxwell equal area law [35][36][37][38][39][40] to make a further investigation of the van der Waals-type phase transition, of the co-existence curve, and of the P − V critical phenomenon for this charged AdS hairy black hole. We give analytically the critical values of the charged hairy black hole and show in detail the behavior of the phase transition from a small black hole to a large one. We also give the distribution of the critical points in the parameter space of q (the coupling constant of the scalar field) and e (the electric charge of the black hole), and obtain for the uncharged case the fitting formula of the co-existence curve. More importantly, we analyze the validity and invalidity of the Maxwell equal area law for the five-dimensional charged AdS hairy black hole and determine the conditions of the law holding. Meanwhile, the latent heat is calculated, which gives the energy released or absorbed between the small and large black hole phases in the isothermal-isobaric procedure.
The paper is organized as follows. In section 2, we review the analytic solution of the charged AdS hairy black hole in D = 5 dimensions and some relevant thermodynamic quantities. In section 3, we calculate the P − V critical values and investigate the Maxwell equal area law for this five-dimensional charged AdS hairy black hole. This section contains two subsections which correspond to the scenarios without and with charges, respectively. Finally, we devote to drawing our conclusion in section 4.

Analytic solution in D = 5 dimensions
The model described by eqs. (1.1) and (1.2) admits [28,31] an exactly electrically charged solution in five dimensions, where the function f takes the form, is the metric of the three-dimensional surface with a constant curvature, the curvature is positive for k = 1, zero for k = 0, and negative for k = −1, and m and e are two integration constants corresponding to the mass and the electric charge of the black hole, respectively. Here the parameter l represents the curvature radius of the AdS spacetime, which is associated with the cosmological constant Λ whose role is analogous to the thermodynamic pressure, Moreover, the parameter q is characterized as the coupling constant of the scalar field, where ε = −1, 0, 1 and there exists an additional constraint: 10b 0 b 2 = 9b 2 1 , to ensure the existence of this black hole solution. These conditions imply that q only takes values 0, ±|q|. Meanwhile, the scalar field configuration takes the form, 5) and the Maxwell gauge potential reads where the field strength still takes the standard form: The location of this hairy black hole horizon is denoted by r h , which is taken to be the largest real positive root of f (r) = 0. That is, the horizon radius r h satisfies the following polynomial equation [31], For the horizon thermodynamic properties of this hairy black hole, some thermodynamic quantities have been calculated in refs. [28,31] and are listed below for our later use. The thermodynamic enthalpy M, temperature T , entropy S, and charge Q take the following forms, where ω 3(k) denotes the area of the compact three-dimensional manifold with the metric dΩ 2 3(k) . In order to develop the first law and the Smarr relation, the coupling parameter is dealt with [28] as a dynamical variable, which means that q is extended to be a continuous and real parameter. Thus, q should appear in the Smarr relation and its variation should be included in the first law of thermodynamics to make the first law of black hole thermodynamics be consistent with the Smarr relation.
Overall, the extended first law of thermodynamics can be written in terms of the thermodynamic quantities mentioned above as follows, where the thermodynamic volume V , the electric potential Φ, and the extensive variable K conjugate to the coupling parameter q have the forms, (2.14) Correspondingly, the extended Smarr relation can be deduced [28], Meanwhile, the Gibbs free energy and the equation of state for this hairy black hole can be written in terms of eqs. (2.8), (2.9), and (2.10) as follows, Next, we shall discuss the Maxwell equal area law for this hairy black hole in order to investigate its critical behavior and the co-existence curve of two phases resorting to eqs. (2.13) and (2.18) in the (P, V ) plane. Note that in the planar case, i.e. k = 0, according to eq. (2.4), we can see q = 0. It implies that there are no hairs. In addition, for the hyperbolic case, i.e. k = −1, there are no physically critical values as pointed out by ref. [28]. Hence, we shall focus only on the spherical case, i.e. k = 1, in the following context.

Maxwell equal area law
The thermodynamic behavior of black holes in the AdS background is analogous to that of the real van der Waals fluid. As was known, the critical behavior of the van der Waals fluid occurs at the critical isotherm T = T c when the P − V diagram has an inflection point, When T < T c , there is an oscillating part in the P − V diagram. We have to replace this oscillating part by an isobar in order to describe it in such a way that the areas above and below the isobar are equal to each other. This treatment is based on the Maxwell equal area law. Thus, by making an analogy between the black hole in the AdS background and the real van der Waals fluid, we find that there also exists an oscillating part below the critical temperature T c in the P − V diagram, which indicates that the first order phase transition occurs from a small black hole to a large one. This isobar that satisfies the Maxwell equal area law represents the co-existence curve of small and large black holes [40]. Normally, the Maxwell equal area law is constructed in the (P, V ) plane for a constant temperature, and it can also be made in the (T, S) or (Φ, Q) plane. These constructions are equivalent. Theoretically, the law can be established from the variation of the Gibbs free energy defined by eq. (2.17), With resorting to the first law of thermodynamics (2.12) and keeping in mind that the coexisting phases have the same Gibbs free energy, one thus arrives at the Maxwell equal area law in the (P, V ) plane by integrating eq. (3.2) at constant T , Q, and q, P (r 1 , T ) = P (r 2 , T ) = P * , where P * stands for an isobar, V 1 and V 2 denote the thermodynamic volume defined by eq. (2.13) for the small and large black holes with the horizon radii r 1 and r 2 , respectively. Thanks to the Maxwell equal area law eq. (3.3), we can also obtain the latent heat which represents the amount of energy released or absorbed from one phase to the other in the isothermal-isobaric condition, where S 1 and S 2 denote the entropy defined by eq. (2.10) for the small and large black holes with the horizon radii r 1 and r 2 , respectively. It is worth mentioning that the co-existence curve of the two phases, i.e. the small and large black holes that are described by an isobar in the Maxwell equal area law, is governed by the Clausius-Clapeyron equation, In the following we specialize in the Maxwell equal area law, the critical values, and the co-existing phases for this AdS hairy black hole.

Uncharged case: e = 0
For this hairy black hole with the spherical symmetry in the absence of a Maxwell field, according to eqs. (2.18) and (2.17), we can see that the equation of state reduces to be and that the Gibbs free energy becomes 1

Critical values
The critical values defined by eq. (3.1) have been obtained in ref. [28], which exist only for q < 0, and the relevant ratio reads We observe that this ratio is a constant that does not depend on the parameter q, and that it is different from that of the real van der Waals fluid or the five-dimensional Reissner-Nordström AdS black hole [19].

Maxwell equal area law
Inserting eqs. (3.6) and (2.13) into eq. (3.3), we give 2 the Maxwell equal area law, (3.10) When taking the critical limit r 1 = r 2 = r c , we can see that eq. (3.10) turns back to eq. (3.8). Moreover, with the help of eqs. (2.10) and (3.10) we obtain the latent heat L between the small and large black hole phases, In order to highlight the outstanding thermodynamic properties of the uncharged case, the numerical calculations about the Maxwell equal area law eq. (3.10) in the (P, V ) plane are displayed in Table 1, meanwhile the P −V critical behavior and the Gibbs free energy described by eqs. (3.6), (3.7), and (2.13) are portrayed in Figure 1 for the specific value of q = −2.
From Table 1, we can see clearly that the horizon radii of the small and large black holes, r 1 and r 2 , shrink into the critical horizon radius r c = 2.15443 when the isotherm condition T = T c = 0.11081 is taken. It implies that no phase transitions occur and correspondingly the latent heat L is of course equal to zero. At this moment, the isobar P * becomes the critical pressure P c = 0.01543. With gradually decreasing of the temperature T from the critical value T c = 0.11081, it is clear that r 1 decreases while r 2 increases, and the Maxwell equal area law is always valid. Furthermore, the small black hole with the radius r 1 and the large one with the radius r 2 share the same Gibbs free energy, where the values of the Gibbs free energy are listed in the fifth column of Table 1. Quantitatively, the latent heat L between the two phases takes a sharp increase with the temperature decreasing. Qualitatively, our analysis is given below. As we have known, r 1 decreases while r 2 increases when T decreases. Combining eq. (3.10) with eq. (3.11), we can get L ∼ 3π 4 r 2 2 when T is decreasing. Due to r 2 increasing, the latent heat L increases in the quadratic form of r 2 . As a result, the qualitative analysis coincides with the quantitative one shown in Table 1.  .10), the Gibbs free energy eq. (3.7), and the latent heat eq. (3.11) for the hairy black hole with the spherical symmetry in the absence of a Maxwell field.
Next, we take a close look at Figure 1 which portrays the critical behavior of the uncharged case. The left diagram demonstrates the representative P − V critical curves. When the temperature T exceeds the critical temperature T c , see the black curve, there is no any criticality. The middle diagram illustrates the typical equal area law, where we indeed observe a van der Waals-type oscillation shown in the red curve. By analogy with the real van der Waals fluid, this oscillating part must be replaced by an isobar in the P − V diagram. Adopting the Maxwell equal area law, i.e., the areas of the regions surrounded by an oscillation (red dashed curve) and an isobar (black solid line) are equal to each other, we can find out this isobar P * and correspondingly the thermodynamic volumes V 1 and V 2 . In other words, the phase transition from the small black hole to the large one occurs in this situation, and these two types of black holes share the same Gibbs free energy displayed by the crossing point A in the right diagram.
With the aid of the Clausius-Clapeyron equation eq. (3.5) and the Maxwell equal area law eq. (3.3), we can obtain the co-existence curve of two phases. Considering that the co-existence curve for the real van der Waals fluid has a positive slope everywhere and terminates at the critical point, we try to fit it using a polynomial. Thanks to the simple forms of the critical point eq. (3.8), such a treatment can be realized. By introducing the reduced parameters, we eventually obtain the parametrization form of the co-existence curve, +5.05509t 7 − 2.53495t 8 + 0.241581t 9 + 0.178083t 10 , where t ∈ (0, 1). Then we plot the numerical values governed by the Clausius-Clapeyron equation eq. (3.5) and the Maxwell equal area law eq. (3.3), and also plot the fitting formula eq. (3.13) of the co-existence curve of the small and large black holes in Figure 2. We can see clearly that the numerical values and the fitting formula match well with each other. Both of them terminate at the critical point, i.e. the point C (1, 1). With the increasing of q, the critical values, see eq. (3.8), take an increasing trend shown in the right diagram of Figure 2. At last, we make a brief summary of this subsection. For the AdS hairy black hole with the spherical symmetry in the absence of a Maxwell field, the P − V critical behavior is available only under the condition of q < 0, and the Maxwell equal area law always holds under this condition. Meanwhile, we take the specific value of q = −2 as an example to highlight the salient features of this black hole in Table 1 and Figure 1. Furthermore, because the critical point eq. (3.8) has a simple form, we obtain the fitting formula (3.13) of the co-existence The co-existence curve for the five-dimensional AdS hairy black hole without charges, the blue points for the numerical values, and the red solid curve for the fitting formula. Right: The co-existence curve in the parameter q space. Each blue dotted curve indicates a co-existence curve at a fixed q and the red curve at the boundary corresponds to the critical points.
curve, see Figure 2 for the illustration. As to the charged case, it has a little complexity that is exhibited in the following subsection.

Charged case: e = 0
For the spherically symmetric AdS hairy black hole with the Maxwell field, the equation of state takes the form, and the Gibbs free energy reads (3.15)

Critical values
Substituting eq. (3.14) into eq. (3.1), we can obtain the critical radius which satisfies the following equation, and derive the critical temperature, , (3.19) which turns back to eq. (3.9), i.e. the uncharged case, when e = 0. Moreover, if q = 0, this ratio equals 5/16, which coincides with the ratio of the five-dimensional Reissner-Nordström AdS black hole [19]. Nonetheless, this ratio depends in general on the values of parameters e and q.
Let us try to solve eq. (3.16). At first, we set w(r c ) ≡ r 4 c + 5qr c − 15e 2 and take a look at its asymptotic behavior. When r c → 0, we get a negative value, w(r c ) → −15e 2 . When r c → +∞, we see w(r c ) → +∞. In addition, we notice that for q > 0 the function w(r c ) is monotone increasing at the interval [0, +∞), and for q < 0 it is monotone increasing at the interval [(−5q/4) 1/3 , +∞) but monotone decreasing at the interval [0, (−5q/4) 1/3 ]. Hence, we can conclude that the equation w(r c ) = 0 must have one and only one positive root, 3 where the plus sign corresponds to the case of q < 0 and the minus sign to the case of q > 0, and the newly introduced parameter t c is defined as
I. The case of q < 0 In this situation, the entropy of the hairy black hole with the Maxwell field, eq. (2.10), is always positive. There are no additional restrictive conditions for investigating the Maxwell equal area law and the phase transition. The treatment is the same as that of the scenario without charges. We can directly write down the condition under which the van der Waalstype phase transition exists and the Maxwell equal area law holds, i.e., the temperature takes the values T < T c , where T c is given by eq. (3.17) together with the constraint q < 0. Table 2 displays the numerical results of the Maxwell equal area law (eq. (3.21)), and Figure 3 depicts the critical behavior of the equation of state (eqs. (3.14) and (2.13)) and the Gibbs free energy (eq. (3.15)) for an example of the charged hairy black hole at q = −2 and e = 2. We can see the similar critical behavior to that of the uncharged case. When T > T c , there is no criticality; when T < T c , there exists a representative van der Waals-type oscillation in the P − V diagram. By using the Maxwell equal area law and replacing the oscillating part by an isobar, we observe the phase transition from the small black hole to the large one in the middle diagram of Figure 3, and find that the two phases share the same Gibbs free energy displayed by the crossing point A in the right diagram. Moreover, when T is decreasing, the horizon radius of the small black hole, r 1 , decreases while that of the large black hole, r 2 , increases, and the latent heat L between the two phases presents a sharp increasing tendency. This observation means that the phase transition needs more energy at the lower temperature of the isothermal-isobaric procedure.  II. The case of q > 0 In the situation of a positive parameter q, the analysis becomes complicated. Due to the entropy described by eq. (2.10) being probably negative, we have to impose an extra constraint to the horizon radius, The purpose is to avoid the negative entropy that appears if r h < r s , where the negative entropy is regarded as an unphysical variable at present. The Maxwell equal area law holds under the condition T < T c together with the constraint eq. (3.22), where T c is given by eq. (3.17) under the condition q > 0. Now we consider an extreme situation, i.e., r c = r s , which leads to the critical value of the electric charge, This critical value (r s or e s ) gives the boundary of violating or maintaining the P −V criticality and the Maxwell equal area law. That is, if |e| < |e s |, we have r c < r s , resulting in the violation of the P − V criticality (eqs. (3.16)-(3.18)) and of the Maxwell equal area law (eq. (3.21)); if |e| > |e s |, we have r c > r s , maintaining the validity of the P − V critical values, but we have to add the condition r 1 > r s in order to establish the Maxwell equal area law. A sample of the latter situation is depicted in Figure 4 for the specific fixing of q = 2 and e = 1.5.
In this sample, we have e > e s = 1.30766 and figure out r s = 1.70998 and r c = 1.94471. Hence, the validity of the Maxwell equal area law depends on the temperature with the  Figure 4 correspond to T = 0.14000 < T c = 0.14207 that leads to r 1 > r s , and illustrates the equal area law holding and the characteristic swallowtail behavior existing. The red dashed curve corresponds to the negative entropy which is unphysical, but it does not affect the application for the equal area law and the characteristic swallowtail structure owing to r 1 > r s . On the contrary, if r 1 < r s , see the right two diagrams of Figure 4, the red dashed curve corresponding to the negative entropy indeed violates the equal area law and the characteristic swallowtail structure. Let us take a close look at this kind of violation. We can still determine the isobar P * because the initial point at r 1 (or V 1 ) and the terminal one at r 2 (or V 2 ) are independent of eq. (3.22). However, when T = 0.13800 that brings about r 1 < r s , i.e. the breaking of eq. (3.22), we find that the branch of negative entropy terminates at r s (or V s ) that is larger than r 1 , which evidently leads to the violation of the Maxwell equal area law and the characteristic swallowtail structure.
In addition, we take a look at the co-existence curve governed by the Clausius-Clapeyron equation eq. (3.5) and the Maxwell equal area law eq. are much more complicated than eq. (3.8), it is hard to give a perfect fitting formula of the co-existence curve of the small and large black holes. Nonetheless, we can know through the above analysis that the thermodynamic behavior for the case with charges is the same as that of the case without charges within the scope of validity of the Maxwell equal area law, and that the co-existence curve always has a positive slope and terminates at the critical point. Figure 5 shows the critical points for different values of e and q, from which we can see that the critical values T c and P c take an increasing trend with an increasing q but a fixed e, while they take a decreasing trend with an increasing e but a fixed q.
Before the end of this section, let us take a look at the physical meaning of the appearance of negative entropy in the charged case with positive q. We illustrate this situation by setting q = 2 and e = 1.5, where the behavior of the Gibbs free energy (eq. (3.15)) can be classified into three types that are shown in Figure 6 and the left-down and right-down diagrams of Figure 4. We know that the stable black hole is thermodynamically favorable to the lower Gibbs free energy.
• If P * < P s < P t , see Figure 6, one can see that the global minimum of the Gibbs free energy emerges a discontinuous characteristic due to the negative entropy (see the red dashed curve). More precisely, when P s < P < P t , the global minimum is located on the curve a → b. When P * < P < P s , the global minimum is located on the curve c → d. When P < P * , the global minimum is located on the curve d → e. At P = P s , the Gibbs free energy has a discontinuous global minimum, which means the occurrence of the zeroth order phase transition as reported in ref. [28]. At P = P * , the global minimum of the Gibbs free energy has an inflection point, which implies a standard first order small/large black hole phase transition, i.e. the van der Waals-type phase transition, and the Maxwell equal area law holds.
• If P s > P t , see the left-down diagram of Figure 4, the global minimum of the Gibbs free energy is continuous and also has an inflection point (corresponding to the pressure P * ), indicating the existence of the van der Waals-type phase transition.
• If P s < P * , see the right-down diagram of Figure 4, the global minimum of the Gibbs free energy is continuous but has no inflection points, which indicates nonexistence of phase transitions and results in invalidity of the Maxwell equal area law.
All in all, it is the positivity/negativity of the entropy that affects the behavior of the global minimum of the Gibbs free energy. The discontinuity of the Gibbs free energy causes a zeroth order phase transition at the pressure P s ∈ [P * , P t ] under a certain temperature. Moreover, if P s < P * , the disappearance of inflection points leads to the failure of the Maxwell equal area law, so that no real physical isobar exists for phase transitions, as shown in the right-top diagram of Figure 4.

Conclusion
In this paper, we investigate the P − V criticality, the Maxwell equal area law, and the coexistence curve for the spherically symmetric AdS black hole with a scalar hair [28,31] both in the absence of and in the presence of a Maxwell field, respectively. Especially in the charged case, we give the exactly analytical P − V critical values, see eqs. (3.16), (3.17), (3.18), and (3.20). Meanwhile, we provide the conditions of validity of the Maxwell equal area law for the hairy black hole without and with charges, respectively. Our results can be summarized as follows: • The case of q < 0 -Scenario without charges: the Maxwell equal area law holds under the conditions T < T c and eq. (3.8).
-Scenario with charges: the Maxwell equal area law holds under the conditions T < T c and eq. (3.17).
• The case of q > 0 -Scenario without charges: the Maxwell equal area law is violated.
-Scenario with charges: besides the conditions T < T c and eq. (3.17) together with the extra constraint eq. (3.22), whether the Maxwell equal area law holds or not depends on the charge and the relation between the temperature and the horizon radius, which can be classified into the following two situations, * when |e| < |e s |, the law does not hold. * when |e| > |e s |, the law holds if r 1 > r s , but fails if r 1 < r s , where r 1 depends on the temperature T shown in Figure 4.
Within the scope in which the Maxwell equal area law holds, we point out that there exists a representative van der Waals-type oscillation in the P − V diagram. This oscillating part that indicates the phase transition from a small black hole to a large one can be replaced by an isobar and the small and large black holes share the same Gibbs free energy. These salient features have been illustrated in Tables 1 and 2, and Figures 1 and 3, from which we conclude that when the temperature T is decreasing, the horizon radius of the small black hole r 1 decreases while that of the large black hole r 2 increases, and moreover, the latent heat L between the two phases presents a sharp increasing tendency. This observation means that the phase transition needs more energy at the lower temperature of the isothermal-isobaric procedure. Furthermore, for the uncharged case we obtain the fitting formula (3.13) of the co-existence curve depicted in Figure 2 due to the simple form of the critical point eq. (3.8). For the charged case, we give the distribution of the critical points described by eqs. (3.17) and (3.18) in the parameter space of q and e in Figure 5. Finally, we point out that the positivity/negativity of the entropy makes effect on the global minimum of the Gibbs free energy. The inflection point of the Gibbs free energy leads to the van der Waals-type phase transition, but the Maxwell equal area law does not hold and the real physical isobar does not exist for phase transitions if no inflection points exist as shown in Figure 4. The discontinuity of the Gibbs free energy causes a zeroth order phase transition as shown in Figure 6.

B. Root of Eq. (3.16)
For a special quartic equation, such as x 4 + bx − c 2 = 0, we set where α, β, and γ are parameters to be determined. At first, we write the four roots of the above quartic equation by solving the two quadratic equations separated from the quartic one, Next, expanding the right hand side of eq. (B1) and comparing the terms with the same power to x, we have Solving eqs. (B3) and (B4), we first determine two of the three parameters, α and β, In order to fix the parameter γ, we then insert eq. (B6) into eq. (B5) and introduce the new parameter y ≡ γ 2 , which can change the sixth-order equation with respect to γ into a cubic equation with respect to the new parameter y, For this kind of cubic equations, there is only one real root that is given by a hyperbolic form of Viète's solution, which actually determines the last parameter γ. As a result, eqs. (B2), (B6), and (B8) give the exactly analytical roots of the original quartic equation.