Interplay between the holographic QCD phase diagram and mutual & n-partite information

In an earlier work, we studied holographic entanglement entropy in QCD phases using a dynamical Einstein-Maxwell-dilaton gravity model whose dual boundary theory mimics essential features of QCD above and below deconfinement. The model although displays subtle differences compared to the standard QCD phases, however, it introduces a notion of temperature in the phase below the deconfinement critical temperature and captures quite well the entanglement and thermodynamic properties of QCD phases. Here we extend our analysis to study the mutual and n-partite information by considering n strips with equal lengths and equal separations, and investigate how these quantities leave their imprints in holographic QCD phases. We discover a rich phase diagram with n ≥ 2 strips and the corresponding mutual and n-partite information shows rich structure, consistent with the thermodynamical transitions, while again revealing some subtleties. Below the deconfinement critical temperature, we find no dependence of the mutual and n-partite information on temperature and chemical potential.


Introduction
Recent developments in string theory suggest that the idea of gauge/gravity duality [1][2][3] can shed new light on the intriguing connection between quantum information notions, quantum field theories and spacetime geometries [4][5][6][7][8][9][10][11]. At the heart of these advancements is the seminal work of Ryu-Takayanagi [4,5], which gave a holographic framework for calculating entanglement entropy. The Ryu-Takayanagi (RT) proposal relates the entanglement entropy of the boundary theory to the area of minimal surfaces, which are homologous to the boundary of the subsystem and extend into the bulk. The Ryu-Takayanagi entanglement entropy proposal is one of the most significant and useful suggestions that has emerged from the gauge/gravity duality, providing not only a deep connection between quantum information and geometry but also opens a new way to calculate and understand other information theoretic quantities such as the mutual or n-partite information [10][11][12].

JHEP04(2019)137
One of the main and original motivations of the gauge/gravity duality was to understand quantum chromodynamics (QCD) at strong coupling. Indeed, the lack of any non-perturbative theoretical tool as well as the large computational complexities and expenses involved in lattice simulations make the gauge/gravity duality the only reliable tool at our disposal to investigate QCD at strong coupling. The main idea here is to construct a gravity model whose dual boundary theory incorporates the essential features of QCDsuch as confinement/deconfinement transition, running of the coupling constant, temperature dependent Wilson and Polyakov loop expectation values, meson mass spectrum etcas accurately as possible. This research area, sometimes called AdS/QCD or holographic QCD, has attracted a lot of interest lately and by now many holographic models, both string theory inspired top-down as well as phenomenological bottom-up models, have been constructed which reproduce many QCD properties holographically .
Interactions in quantum field theories (QFT) via the entanglement in quantum states cause quantum information to be dispersed non-locally across space. It is therefore of great interest to examine how this structure of shared information changes with the size of the subsystem (the length scale), as it might provide important information about the confinement structure. With the exception of a few lattice related papers [44][45][46][47], the discussion of entanglement entropy in QCD like gauge theories has been rather limited. The conceptual as well as computational difficulties presented in the definition of entanglement entropy for interacting field theories make it extremely difficult to get any reliable nonperturbative estimate of the entanglement entropy relevant for QCD. On the other hand, the holographic RT proposal bypasses the technical difficulties presented in the computation of entanglement entropy of quantum field theories and therefore one can use this proposal to find the entanglement structure of QCD. This idea was first initiated in [48] in the top-down models of gauge/gravity, where a change in the entanglement entropy order was observed (from O(N 2 ) to O(N 0 ) or vice versa) as the size of the entangling region varied. In particular, a phase transition between connected and disconnected RT entangling surfaces in the confining background was obtained, causing a non-analyticity in the structure of entanglement entropy. This transition was suggested as an indication of (de)confinement in [48]. Importantly, the non-analytic behaviour of entanglement entropy in the confining background has received numerical confirmation from lattice papers as well [44][45][46]. The idea of [48] was then applied to many other town-down confining as well as soft wall models of holographic QCD [49][50][51][52][53][54][55][56][57]. Only recently the entanglement entropy computations for the phenomenological bottom-up models, which are somewhat more appropriate to model QCD holographically [58,59], were performed and the results were similar to those reported in [48].
Most holographic discussion concerning the entanglement structure of QCD has been restricted to the entanglement entropy only (with one subsystem). However, there are other information theoretic quantities such as the mutual and n-partite information that appear when two or more disjoint subsystems are considered [10,60,61]. These quantities do not suffer from the ambiguities associated with the entanglement entropy and can provide more information than the entanglement entropy alone. For example, they are finite and do not suffer from the usual UV divergences. These quantities, therefore, provide a UV JHEP04(2019)137 cutoff independent information as opposed to the entanglement entropy which explicitly contains the UV cutoff. Likewise, the tripartite information (n = 3), which quantifies the extensivity of the mutual information (n = 2), measures how much of the information that were presented in one part of the system can only be retrieved when having access to both parts of a bipartite system.
For two subsystems A 1 and A 2 , the mutual information is defined as the amount of information that A 1 and A 2 can share. In terms of the entanglement entropy it is written as where S(A 1 ), S(A 2 ) and S(A 1 ∪ A 2 ) are as usual the entanglement entropies of A 1 , A 2 and their union respectively. Form eq. (1.1) it is evident that I(A 1 , A 2 ) is zero for two uncorrelated subsystems. Moreover, the subadditivity property of the entanglement entropy also ensures that I(A 1 , A 2 ) is non-negative i.e. I(A 1 , A 2 ) provides an upper bound on the correlation functions between operators in A 1 and A 2 . For n disjoint subsystems, the above definition can be generalised to define n-partite information From the above definition, it is clear that n-partite information for n ≥ 2 is UV finite. Indeed, because of these desirable features, the mutual and n-partite information have been used, both from holography as well as from field theory point of view, to probe various interesting physics in a variety of systems, see for example [62][63][64][65][66][67][68][69][70][71][72][73][74][75][76][77][78]. The discussion of mutual and n-partite information with two or more disjoint intervals in holographic QCD is however relatively new. A partial discussion appeared in [79], where only the entanglement entropy phase diagram for n disjoint intervals in a top-down gauge/gravity confining model was discussed. However, as is well known, the top-down holographic QCD models usually face several limitations in mimicking real QCD. In particular, the boundary field theory of these top-down models generally contains additional Hilbert space sectors (arising from the KK modes of extra dimensions), the non-running coupling constant, undesirable conformal symmetries etc, whose analogue in real QCD do not exist [14][15][16][17][18][19]. On the other hand, the phenomenological bottom-up holographic QCD models, although often formulated in an ad-hoc manner to reproduce desirable features for the boundary QCD and lack solid gauge/gravity duality foundations, can overcome most of the difficulties presented in the top-down models . Since, the amount of correlation between two or more disjoint subsystems is actually characterised by the mutual or n-partite information, it is therefore of great interest to study them in a consistent bottomup holographic QCD model, in particular, to disclose further the additional information they can provide into the QCD vacuum structure and confinement mechanism.
In this work, our main aim is to fill the above mentioned gap by studying the mutual and n-partite information in a self-consistent bottom-up holographic QCD model by considering two or more subsystems. By now many phenomenological holographic QCD JHEP04(2019)137 models have been constructed, see for example [23][24][25][26][27][28][38][39][40][41][42][43], with each having own merits and demerits. Here, we consider a particular phenomenological Einstein-Maxwell-dilaton (EMD) holographic QCD model constructed in [80]. An important advantage of this model is that it can be solved exactly and the full gravity solution can be obtained analytically in terms of the gauge-dilaton coupling f (φ) and scale factor P (z) (see section 2 for more information). Moreover, by taking appropriate forms of these two functions, important real QCD properties like vector meson mass spectrum, confinement/deconfinement transition, Wilson loop area law etc can be realised holographically in this model as well. Importantly, by taking suitable forms of P (z) (see eq. (3.1) and (4.1)), this model not only predicts the standard confined and deconfined phases but also a novel specious-confined phase. This specious confined phase -although not exactly equivalent to the standard confined phase, however, shares many of its properties -is dual to a small black hole phase in the gravity side, and therefore has a notion of temperature. This, in turn, allow us to investigate temperature dependent profiles of various observables in the specious confined phase, which can be compared with lattice QCD confined phase predictions. Indeed, it was found that thermal behaviour of the quark-antiquark free energy, entropy and the speed of second sound in the specious-confined phase were qualitatively similar to lattice QCD predictions [80]. Because of these analytic and interesting features, the model of [80] has also been used in other holographic areas such as in holographic complexity [81,82], and here we consider it to investigate the mutual and n-partite information in holographic QCD.
In this work, following [80], we take two different forms P (z). Neither of these forms change the asymptotic structure of the boundary, however cause non-trivial modifications in the bulk spacetime. The first form P (z) = P 1 (z) gives thermal-AdS/black hole phase transition in the gravity side, which in the dual boundary theory corresponds to standardconfined/deconfined phase transition. We then examine the entanglement entropy, mutual and n-partite information in the obtained confined/deconfined phases by considering one or more strip geometries of length as the subsystems. In the confined phase, with one strip, the entanglement entropy again undergoes a connected to disconnected surface transition and exhibits non-analytic behaviour at the critical length c . These features are the same as suggested by [48]. However, with two or more strips, four minimal area surfaces {S A , S B , S C and S D } appear (see figure 6 for more details) which compete with each other and lead to an interesting phase diagram in the parameter space of and x (x being the separation length between the strips). In particular, two distinct tri-critical points appear where three minimal area surfaces coexist. Interestingly, depending on the surfaces involved, the order of the mutual and n-partite information may or may not change at the transition point. This is very different from the entanglement entropy behaviour, where the order always changes at the transition point. On the other hand, in the deconfined phase, only two minimal area surfaces (S A and S B ) appear with two or more subsystems. There is again a phase transition between these two surfaces, and this phase transition is always accompanied by a change in the order of mutual and n-partite information. We further find that these information theoretic results for the confined/deconfined phases qualitatively remain the same even when chemical potential is considered. We then discuss the holographic QCD phase diagram by examining the mutual and n-partite information JHEP04(2019)137 in the temperature-chemical potential plane. We find that these quantities capture the signature of thermal-AdS/black hole phase transition (or dual confined/deconfined phase transition), suggesting that these non-local observables, like the entanglement entropy, are also sensitive to the phase transition.
With the second form P (z) = P 2 (z), we instead find the small/large stable black hole phases, which on the dual boundary theory correspond to the specious-confined/ deconfined phases. The entanglement entropy computations of the specious-confined phase reveal a novel connected to connected surface transition (instead of a connected to disconnected transition), where the order of the entanglement entropy does not change. With two or more strips, unlike in the case of standard confined phase, the − x phase diagram in the specious-confined phase contains only two phases (S A and S B ) and the transition between S A and S B is again accompanied by a change in the order of mutual and n-partite information. Moreover, the mutual information also behaves desirably in the speciousconfined phase and satisfies non-negative property. Interestingly, the mutual and n-partite information of the deconfined phase here are similar to the deconfined phase mutual and n-partite information obtained using P 1 (z). Further, we investigate the dual speciousconfined/deconfined QCD phase diagram by studying the mutual and n-partite information in the temperature-chemical potential plane, and again find that these information theoretic quantities are sensitive to the phase transition.
The paper is organised as follow. In the next section, we briefly review our EMD gravity solution and then derive the necessary entanglement entropy formulae. In section 3, using the first form of P (z), we first examine the thermodynamics of the gravity solution and then discuss the entanglement entropy, mutual and n-partite information in the dual confined/deconfined phases. In section 4, we repeat the computations of section 3 with the second form of P (z). The last section is devoted to conclusions and an outlook to future research.

Holographic set up
In this section, we briefly describe the EMD gravity model as well as the holographic entanglement entropy and state only the useful expressions, which will be important for our investigation in later sections. The holographic EMD gravity model at finite and zero temperature as well various expressions for the entanglement entropy have been discussed in great detail in [58,80], and we refer the reader to [58,80] for more technical details.
The EMD action in five dimensions consists of Ricci scalar R, a field strength tensor F M N and a dilaton field φ, where the gauge kinetic function f (φ) represents the coupling between the gauge field A M and φ, G 5 is the Newton constant in five dimensions and V (φ) is the potential of the dilaton JHEP04(2019)137 field. Interestingly, using the following Ansätze, the equations of motion of the above EMD action can be explicitly solved analytically in terms of a scale function P (z) [39,40,58,80], The above (Einstein frame) gravity solution corresponds to a black hole with a horizon at z = z h . This solution is obtained by using the boundary condition that lim z→0 g(z) = 1 at the asymptotic boundary z = 0, and that g(z h ) = 0 at the horizon. Here µ is the chemical potential of the boundary theory, which is obtained from the asymptotic expansion of the gauge field. The form of coupling function f (z) in eq. (2.3) is also arbitrary and we chose f (z) = e cz 2 −P (z) so that the holographic meson mass spectrum of the boundary theory lies on a linear Regge trajectory, as governed by QCD phenomenology. Similarly, the magnitude of the parameter c = 1.16 GeV 2 is fixed by matching the holographic meson mass spectrum to that of lowest lying (heavy) meson states. Let us also note the expressions of black hole temperature and entropy, where V 3 is the volume of the three-dimensional plane. Another solution of EMD action can be obtained by taking the limit z h → ∞, which implies g(z) = 1. This solution corresponds to thermal-AdS (without horizon). The thermal-AdS solution has an asymptotically AdS structure at the boundary z = 0, however it can have, depending on P (z), a non-trivial structure in the bulk. The non-trivial structure of thermal-AdS in the IR region is in fact the same reason for having confinement behaviour in the dual boundary theory [20].

JHEP04(2019)137
Let us now briefly discuss the holographic entanglement entropy and its relevant expressions in EMD gravity model. We concentrate only on one entangling surface, as the mutual and n-partite information can be obtained from it using eqs. (1.1) and (1.2). According to the RT prescription, the entanglement entropy of the subsystem A is given by the area of the minimal surface γ A which extends from the AdS boundary into the bulk and shares the same boundary ∂A as the subsystem A, Here we consider a strip of length as the subsystem, i.e. the strip domain − /2 ≤ y 1 ≤ /2, 0 ≤ y 2 ≤ L y 2 and 0 ≤ y 3 ≤ L y 3 defines the entangling surface on the boundary. With the strip subsystem, there are two local minima surfaces of eq. (2.5): a (U-shaped) connected and a disconnected surface [58]. The entanglement entropy of the connected surface is given by the following expression, where z * is the turning point of the connected minimal area surface and is related to the strip length in the following way .

(2.7)
On the other hand, the entanglement entropy for the disconnected surface is given by where z d = ∞ for thermal-AdS background and z d = z h for AdS black hole background. The first term in eq. (2.8) comes from the entangling surface along the horizon and therefore does not contribute for the thermal-AdS background. Hence, S EE discon is actually independent of for the thermal-AdS background. As found in [58], this behaviour of disconnected entanglement entropy provided several interesting features in the entanglement entropy phase diagram in the dual confined phase. As we will show shortly this behaviour provides even richer phase structure in the mutual and n-partite information. For the AdS black hole background, however, the first term provides a finite contribution to the holographic entanglement entropy.
3 Case I: the confined/deconfined phases As in [58,80], let us first consider the following simple form of the scale function P (z), P (z) = P 1 (z) = −āz 2 .  It is easy to observe that P 1 (0) → 0, asserting that spacetime asymptotes to AdS. Also, where m 2 = ∆(∆ − 4) with ∆ = 3, satisfying the well known relation of the gauge/gravity duality. The parameterā = 0.145 is fixed by requiring the transition temperature T c of the thermal-AdS/black hole (or the dual confinement/deconfinement) phase transition to be around 270 MeV at µ = 0 in the pure glue sector, as is observed in large N lattice QCD [83].

Black hole thermodynamics
The thermodynamics of the gravity solution with eq. (3.1) has been discussed in [80] and here we just briefly highlight its main features. The thermodynamic results are shown in figures 1 and 2. For small values of µ, there appears a minimum temperature T min below which no black hole solution exists whereas above T min two black hole solutions -a large stable black hole (marked by 1 ) and a small unstable black hole (marked by 2 )exist at each temperature. The entropy increases with temperature in the large black hole phase indicating its stable nature whereas the entropy decreases with temperature in the small black hole phase indicating its unstable nature. The phase structure is shown in figure 2, where a Hawking/Page type first order phase transition between large black hole and thermal-AdS phases can be observed. The phase transition takes place at T crit > T min .
For higher values of µ, the critical temperature of the Hawking/Page thermal-AdS/black hole phase transition however decreases and the phase transition stops at a critical chemical potential µ c . In particular, at µ c the unstable small black hole phase disappears and we have a single black hole phase which remains stable at all temperatures JHEP04(2019)137 (indicated by a magenta line in figures 1 and 2). For this model, we get µ c = 0.673 GeV. Since finding an exact value of the QCD critical point in the QCD (T, µ) plane, if existing, is extremely hard [84][85][86], a reasonable estimate is a few hundred MeV, so we observe our estimate for the critical point lies in the same range.
In [80], by calculating the free energy of the probe quark-antiquark pair, it was further shown that the above Hawking/Page phase transition on the gravity side is dual to the confinement/deconfinement phase transition on the dual boundary side. In particular, the large black hole phase was shown to be dual to the deconfined phase whereas the thermal-AdS phase was shown to be dual to the confined phase. Since the backreaction of the dilaton field is included in a self-consistent manner from the beginning in this model and T crit of the obtained confinement/deconfinement phase transition decreases with µ as well, a result again in line with lattice QCD, this model therefore provides a more realistic holographic QCD model compared to the soft and hard wall models. It is therefore of great interest to investigate how the information theoretic quantities behave in this more physical holographic QCD model.

Entanglement in holographic QCD phases with multiple strips
The aim of this subsection is to investigate the entanglement entropy, mutual and n-partite information in the above constructed confined/deconfined phases holographically. Unfortunately, analytic results are difficult to obtain and therefore we provide only numerical results here. We will first discuss the results in the thermal-AdS background and then discuss in the black hole background.

With thermal-AdS background: one strip
Let us first discuss the results in the thermal-AdS background with one strip. This will set the notation and convention for the rest of the section. The results are shown in figures 3 and 4, where the variation of strip length with respect to the turning point z * of the connected entangling surface and the difference between connected and disconnected entanglement entropies ∆S EE = S EE con − S EE discon respectively are plotted. We find that there exist three RT minimal area surfaces: two connected and one disconnected. The two connected surfaces, shown by solid and dashed lines in figure 3, only exist below a maximum length max , and above this max only the disconnected entangling surface remains. Out of the two connected surfaces, the one that occurs at small z * (solid line) always has a lower the entanglement entropy than the one which occurs at large z * (dashed line). This suggests that the small z * connected surface solution is the true minima of eq. (2.6) for small .
Interestingly, a connected to disconnected entanglement entropy phase transition takes place as we increase the strip size . A pictorial illustration of this is shown in figure 5. In particular, ∆S EE = S EE con − S EE discon changes sign from a negative to positive value as increased. This implies that S EE discon has a lower entanglement entropy for large whereas S EE con has a lower entanglement entropy for small . The strip length at which this conncted/disconncted phase transition occur defines an crit (< max ). For c = 1.16, we find  Since the area of the disconnected surface in the thermal-AdS background is actually independent of , the corresponding entanglement entropy becomes independent of as well. These results can be summarized as, where N is the number of colors of the dual gauge group. The above kind of connected/disconnected phase transition between two entangling surfaces, where the order of the entanglement entropy changes at the critical point was first found in the top-down models of the gauge/gravity duality in [48,87]. The behavior that the entanglement entropy scales as N 0 for large and as N 2 for small , led the authors of [48] to interpret the subsystem length as the inverse temperature "T c ∝ 1 c ". Indeed, below the deconfinement critical temperature the color confined degrees of freedom count as order O(N 0 ) whereas above the critical temperature the deconfined gluon (colored) degrees of freedom count as order O(N 2 ), the same counting as estimated by the holographic JHEP04(2019)137 entanglement entropy. This further suggests that the entanglement entropy can act as a tool to diagnose confinement.
In [48], the above type of connected/disconnected entanglement entropy phase transition was suggested to be a characteristic feature of confining theories. Our analysis further validates this claim as we found identical results for the holographic entanglement entropy, however now in a genuine self-consistent bottom-up confining model where the running of the coupling constant is incorporated from the beginning. Moreover, a similar type of non-analyticity in the entanglement entropy has also been observed in SU(2) gauge theory using lattice simulations [44]. Therefore, such non-analyticity in the structure of entanglement entropy seems to be a universal feature of all confining theories. Interestingly, our gauge/gravity duality estimate for the length scale at which non-analyticity in the entanglement entropy appears ( crit ∼ 0.2 fm) is roughly of the same order as estimated by lattice simulations ( c rit ∼ 0.5 fm) [44,46], lending further support to the idea that the gauge/gravity duality can yield compelling predictions for QCD-like theories. As we will see shortly, this non-analytic behavior persists even when two or more subsystems are considered, albeit in those cases various other types of non-analyticity also emerge.

With thermal-AdS background: two strips
Having discussed the holographic entanglement entropy with one strip, we now move on to discuss it with two strips. For convenience, we concentrate only on equal size strips, where Interestingly, depending on the size and separation x between the two strips, there can now be four minimal area surfaces. These surfaces are shown in figure 6. We see that with two strips there can be connected (S A and S B ), disconnected (S D ) as well as a combination of connected and disconnected (S C ) surfaces.
The entanglement entropy expressions of these four minimal surfaces can be written down as, where as usual S EE con and S EE discon are the entanglement entropies of the connected and disconnected entangling surfaces with one strip.
With two strips, there can be various kinds of phase transitions between different entangling surfaces. The complete phase diagram in the phase space of ( , x) is shown in figure 7. We find that for small x, crit it is the S A phase which has the lowest entanglement entropy. The S B phase, however, becomes more favorable as increases. For a further increase in , keeping x( crit ) small, a phase transition from S B to S C takes place. For x = 0, this phase transition from S B to S C occurs at = crit /2, as can also be observed from eq. (3.4). In general, the phase transition line between S B and S C is given by 2 + x = crit . With a further increase in both x, crit the S D configuration ultimately becomes more favorable. In the near origin region, where both x, crit are small, the results are similar to the CFT expectation. This is interesting considering that QCD is expected to become conformal at extremely high temperature [36]. This further provides support to the analogy of strip length as the inverse temperature.  Illustration of the four different minimal area surfaces for two strips of equal length separated by a distance x in the thermal-AdS background. S A corresponds to minimal area surface for small and large x, S B corresponds to minimal area surface for small and small x, S c corresponds to minimal area surface for large and small x and S D corresponds to minimal area surface for large and large x. Interestingly two tri-critical points appear in the two strips phase space, which has no analogue in the one strip case. These two tri-critical points are shown by black dots in figure 7 and are denoted by c1 and c2 respectively. The first tri-critical point, where (S A , S B and S C ) phases coexist, occurs at ( = 0.35 GeV −1 , x = 0.27 GeV −1 ) whereas the second tri-critical point, where (S B , S C and S D ) phases coexist, occurs at ( = 0.96 GeV −1 , x = 0.96 GeV −1 ). These two tri-critical points again suggest non-analyticity in the structure of entanglement entropy. Importantly, the order of the entanglement entropy does not change as we go from one phase to another via the first tri-critical point, however, the order may or may not change as we go from one phase to another via the second tri-critical point. In particular, if one of the phases involved in the transition is S D then only the order changes (from N 2 to N 0 or visa versa), otherwise, it does not change.
It is also interesting to analyze the structure of mutual information in the above phases. For this, let us first note the expressions of mutual information  which also implies We see that depending on the critical point the order of the mutual information may or may not change as we go from one phase to another. For example, going from S B to S A phase (by increasing x) causes a change in the mutual information order (from O(N 2 ) to O(N 0 )) whereas going from S B to S C phase (by increasing ) does not cause such change in order. If we presume that the mutual information, like the entanglement entropy, also carries the information about the degrees of freedom of the system then our analysis suggests that the mutual information can also be used to extract useful information of QCD phases. Moreover, the mutual information also varies smoothly as we pass from one phase to another. For example, the mutual information connects smoothly between S B and S C phases as we pass through the S B − S C critical line. This result is shown in figure 8, where the variation of I and I [n=2] C as a function of is shown. Similarly, the mutual information also smoothly goes to zero as we approach S A (or S D ) phase from S B (or S C ) phase. This is shown in figure 9. These interesting results from holographic might have analogous realization in real QCD. Unfortunately, unlike the entanglement entropy, we do not have any lattice results for the QCD mutual information yet. In this regard, these results for the mutual information can be considered as a prediction from holography. x Figure 11. The n strip phase diagram of various entangling surfaces in the thermal-AdS background. The red, green, blue and brown curves correspond to n = 2, 3, 4 and n = 5 respectively. In units GeV.

With thermal-AdS background: n strips
As in the case of n = 2 strips, there will again be four minimal area surfaces for n > 2 strips. In fact, these can be the only minimal area surfaces have been proved in [79] as well. These surfaces for n = 3 are shown in figure 10 and can be easily generalized to higher n as well. The entanglement entropy expressions of these four minimal surfaces are now given by, Our results for different n are shown in figure 11, where it can be seen that the phase diagram is quite close to n = 2 case. In particular, there are again two tri-critical points where three different entangling surfaces coexist. The main difference from the n = 2 case arises in the size of region S B . As can be seen, the region in the parameter space where phase S B is most stable decreases as the number of strips increases. Although the location of the second tri-critical point ( c2 ) does not change with different n, however the first tri-critical point ( c1 ) moves more and more towards the origin. Therefore, the region S B gets smaller and smaller and eventually will disappear for n → ∞. The behaviour that S B shrinks to zero as n → ∞ can be seen analytically as well. Notice from eq. (3.7) that the transition line between S B and S C satisfies n + (n − 1)x = crit , which implies x, → 0 as n → ∞. Therefore, the size of S B goes to zero as n → ∞.
We now calculate the corresponding n-partite information. It is not hard to see that the n-partite information again vanishes in S A and S D phases. However, its computation in S B and S C phases is now more non-trivial as we need to be careful in evaluating the contributions of (n − 1, n − 2, . . . , 1) intervals entanglement entropy to the n-partite information. Therefore, depending on the values of and x, the explicit expression of npartite information will vary. For example, we have the following expression for 3-partite  information at x = 0.1 GeV −1 , The 3-partite information as a function of and x is shown in figures 12 and 13. Interestingly, unlike the mutual information, the 3-partite information is always non-positive. This indicates the monogamy of mutual information. At this point, it is instructive to point out that in quantum theories without gravity duals the 3-partite information can be negative, positive or zero. However, in the context of field theories with gravitational dual the 3-partite information is always negative, which points to the monogamous nature of the mutual information in holographic theories [10]. Here, we find a similar result of the 3-partite information in holographic QCD theories. Moreover, since the entanglement entropy phase diagram and the transition lines for various entangling surfaces for n = 2 and n = 3 are different, they generate various non-analyticities in the structure of 3-partite information. This should be contrasted with the mutual information where no such non-analytic behaviour was seen. Furthermore, we find that the length at which non-analyticity in the 3-partite information appears, increases with x.
Similar results can be obtained for other values of n as well. We find that the 4partite information, on the other hand, is always non-negative. This is also in line with the holographic suggestion that the n-partite information is positive (negative) for even (odd) n [68]. The 4-partite information also exhibits non-analytic behaviour at various places. Interestingly, the number of points where non-analyticity appears, increases with n. This is an interesting new result, and it would be interesting to find an analogous realization in real QCD using lattice simulations.

With black hole background: 1 strip
Having discussed the holographic entanglement entropy in the dual confined phase, we now move on to discuss it in the dual deconfined phase. Let us consider n = 1 case first. The results are shown in figures 14 and 15.
In the deconfined phase, the entanglement entropy behavior is quite different. In particular, no maximum length ( max ) appears and the connected entangling surface solution of eq. (2.7) persists for all i.e now a one-to-one relation between and z * appears. This is shown in figure 14. This also indicates that the turning point (z * ) of the connected surface shifts more towards the horizon z h as the subsystem size increases. Moreover, no phase transition between disconnected/connected surfaces appears as well. In particular, ∆S EE is always greater than zero, indicating S EE con ≤ S EE discon . It is the first term of eq. (2.8) that makes S EE con ≤ S EE discon . The equality sign here is realized only when the size of the subsystem A approaches the full system size, i.e. when → ∞. Not surprisingly, in the limit → ∞, the entanglement entropy reduces to the Bekenstein-Hawking black hole entropy, which is expected from the general property of the entanglement entropy that it reduces to the thermal entropy at finite temperature when the size of the subsystem A approaches its full system size. Moreover, for the AdS black hole background we always have, Similar results for the entanglement entropy appear for finite values of µ as well. In particular, we again have S EE con ≤ S EE discon , indicating no phase transition between disconnected/connected surfaces as varies.

With black hole background: 2 strips
Since the entanglement entropy of the connected entangling surface is always favored in the black hole background, S A and S B are now the only phases which appear in the deconfinement phase with two strips. Correspondingly, the phase transition appears only between S A and S B connected phases as opposed to the confined phase where the phase transition to disconnected phases (S C and S D ) also occurred. The two strip phase diagram of S A and S B at µ = 0 for various temperatures is shown in figure 16. Again, as in the confined phase, S A is more favorable at larger x whereas S B is more favorable at smaller x. Moreover, we find that S A becomes relativity more favorable than S B as we increase the temperature. This can be observed by comparing red and blue lines of figure 16, where T = 1.2 T c and T = 2.0 T c respectively are used.
The above two strips phase diagram persists for a finite µ as well. This is shown in figure 17, where different values of µ at constant temperature T = 1.2 T c are considered. Interestingly, higher values of µ instead try to make S B more favorable.
Next we analyse the behavior of mutual information in S A and S B . By definition I gets higher for higher temperatures. This is shown in figure 18. Importantly, since I at the transition line. Although not shown here for brevity, the mutual information also varies smoothly as we pass from S B to S A by changing x.
The mutual information also behaves smoothly when finite chemical potential is considered, and most of the above mentioned results remain true with chemical potential as well. The main difference appears in the large asymptotic value of I

With black hole background: n strips
Let us now briefly discuss the phase diagram with n > 2 number of strips. The results for different n are shown in figure 20, where it can be observed that the phase diagram is quite similar to n = 2 case. In particular, again a phase transition between S A and S B occurs as the separation between the strips varies. Due to numerical limitations, it is difficult to exactly establish the phase diagram for large , however, the numerical trend suggests that the critical separation length x crit approaches a constant value i.e. independent of n, for large .
The behavior of tri-partite information (I again indicates the monogamous nature of mutual information, however now in the deconfined phase. We further find that these results hold for finite chemical potential as well. Moreover, I [n=3] also exhibits non-analytic behaviour. The separation length x at which non-analyticity appears decreases with temperature whereas it increases (only slightly) when higher values of chemical potential are considered. Similarly, the 4-partite information also behaves desirably in the deconfined phase. We find that it is always non-negative and exhibits non-analyticities at various places. We further investigate how the 4-partite information varies with temperature and chemical potential. Again, we find that, like the 3-partite information, the separation length at which non-analyticity in 4-partite information appears decreases with temperature whereas it enhances with chemical potential.
At this point, it is instructive to point out that the above results for n-strip phase diagram and n-partite information in the dual deconfined phase of our gravity model are qualitatively similar to what one gets in the dual deconfined phase of the AdS-Schwarzschild black hole. This suggests that, as far as the entanglement structure is concerned, the excited profile of the dilaton field does not lead to a significant effect in the deconfined phase. As we will show shortly, the above mentioned entanglement features of the dual deconfined phase will remain true even when other scale factors are considered. This suggests some type of universality in the entanglement structure of holographic deconfined phases.

Thermal-AdS/black hole phase transition and mutual information
In recent years the holographic entanglement entropy has been used to probe and investigate black hole phase transitions. The main idea here is that since the entangling surface propagates from asymptotic boundary into bulk it therefore might be able to probe a change in the spacetime geometry which occurs during the phase transition. This idea has become a fruitful arena of research lately and has been applied in many different contexts, let us just mention a few [88][89][90][91]. In [58], we performed a similar analysis for the EMD gravity model and found that the holographic entanglement entropy does indeed capture the essence of thermal-AdS/black hole phase transition (discussed earlier in this section). However, one might also wonder whether other information quantities like mutual or npartite information can similarly be used as a diagnostic tool to probe black hole phase is the mutual information of the thermal-AdS space, which is independent of temperature and chemical potential and is constant for a fixed and x. Here, we have used fixed = 0.2 GeV −1 and x = 0.1 GeV −1 so that S B is the most stable phase. We find that, just like the entanglement entropy, the structure of mutual information also displays a striking similarity with the Bekenstein-Hawking thermal entropy. In particular, for small µ, there are again two branches in ∆I B and these two branches exist only above T > T min . The negative slope branch in figure 22 corresponds to the unstable solution whereas the positive slope branch corresponds to the stable solution.
We see that, like the Bekenstein-Hawking entropy, ∆I B also displays double valuedness for µ < µ c -an indication of black hole phase transition -whereas this double valuedness disappears for µ > µ c . An analogous similarity between Bekenstein-Hawking entropy and entanglement entropy has been used in recent years to advocate that the entanglement entropy can be used as a diagnostic tool to probe black hole phase transition [88][89][90][91]. We find that the similarity with the Bekenstein-Hawking entropy goes beyond the entanglement entropy and even the mutual information exhibits similar features in the T − µ plane. In figure 22, we have used = 0.2 GeV −1 and x = 0.1 GeV −1 however we have checked that physical quantities like T min and µ c do not change even when other values of and x are considered. Moreover, although not shown here for brevity, we find that similar results hold for 3 and 4-partite information as well. Our analysis therefore not only confirms the suggestions of [88][89][90][91] in a more advanced holographic bottom-up model but also put further weight on the expectation that other information theoretic quantities can also be used to investigate phase transitions.

Case II: the specious-confined/deconfined phases
In [80], by taking the following second form for the scale factor, a novel specious-confined phase on the dual boundary side was revealed. This speciousconfined phase did not strictly correspond to the standard confined phase, but exhibited many properties which resembled quite well with the standard QCD confined phase. The novelty of the specious-confined phase lies in the fact that it is dual to a non-extremal small black hole phase in the gravity side. It therefore has the notion of temperature, which in turn allows to investigate temperature dependent properties of the dual speciousconfined phase. Indeed, it was shown in [80] that thermal behaviour of the quark-antiquark free energy and entropy, as well as the speed of second sound in the specious-confined phase, were qualitatively similar to lattice QCD results. It is important to mention that for both choices of scale factor P 1,2 (z), the potential is bounded from above by its UV boundary value i.e. V (z) ≤ V (0). Therefore, for both choices P 1,2 (z) the EMD model under consideration satisfies the Gubser criterion to have a well-defined dual boundary theory [92].
With the above scale factor, the metric again asymptotes to AdS at the boundary z → 0. However, it causes non-trivial modifications in the bulk IR region which in turn greatly modifies thermodynamic properties of the system. As in the case of first scale factor P 1 (z), the parameters a = c/9 and b = 5c/16 in P 2 (z) are again fixed by demanding the specious-confined/deconfined phase transition to be around 270 M eV at zero chemical potential, as is observed in lattice QCD for pure glue sector.

Black hole thermodynamics
As we can see from figures 23 and 24, the thermodynamic behaviour of EMD gravity model gets greatly modified with P 2 (z). In particular, on the top of a large stable black hole phase (marked by 1 ) and an unstable black hole phase (marked by 2 ), now a new stable phase appears at low temperatures. This new stable phase is marked by 3 in figure 23 and corresponds to a small black hole phase (large z h ). Importantly, with P 2 (z) at least one black hole phase always exists at all temperature. Apart from these three black hole phases, there also exists a thermal AdS phase. However, we find that the free energy of this thermal-AdS phase is always greater than the stable black hole phases, indicating that it is thermodynamically unfavourable at all temperature. The normalised free energy behaviour, plotted in figure 24, further suggests a first order phase transition between small and large hole phases as the Hawking temperature varied. For µ = 0, the small/large phase transition occurs at T = 0.276 GeV. Therefore, the large black hole phase is thermodynamically favoured at T > T c whereas the small black hole phase is favoured at T < T c .
The small/large black hole phase transition persists for finite chemical potential as well. The complete phase diagram and the dependence of T c on µ can be found in [80], where it was shown that T c decreases with µ for µ < µ c = 0.312 GeV, mimicking yet another important feature of lattice QCD. At µ c , the small/large black hole phase transition ceases to exist and we have a single stable black hole phase which exists at all temperatures (shown by cyan curve in figure 24). Overall, this gravity model shows a Van der Waals type black hole phase transition, however now with a planar horizon instead of a spherical horizon [93][94][95][96].
In [39,40], the above small/large black hole phases were suggested to be dual to confined/deconfined phases in the dual boundary theory. In particular, the small black hole phase was suggested to be dual to the confined phase whereas the large black hole phase was suggested to be dual to the deconfined phase. However, as was pointed out in [80], the small black hole phase does not strictly correspond to the confined phase, as the Polyakov and Wilson loop expectation value do not strictly exhibit the standard behaviour. Interestingly, the dual boundary theory does however exhibit properties, such as the quark-antiquark free energy and entropy etc, which are qualitatively similar to the standard QCD confined phase. For this reason, the dual boundary theory of the small black hole phase was named as specious-confined phase.
Since the specious-confined phase has the notion of temperature (thereby allowing us to study temperature dependent properties of many important observables) and shares many interesting lattice QCD properties, it becomes important to investigate mutual and n-partite information in this phase as well.  Here µ = 0 and red, green and blue curves correspond to T /T c = 0.9, 0.8 and 0.7 respectively. In units GeV.
confined phase. In particular, the connected surface now exists for all . Further, the vs z * behaviour is now divided into three regions (instead of two as in the case of standard confined phase): first increases with z * then decreases and finally increases again. These three regions are marked by 1 , 2 and 3 respectively in figure 25. Importantly, between min and max , these three solutions for the connected entangling surface coexist for a given . This makes the transition between different entangling surfaces more non-trivial in the specious-confined phase.
The difference between connected and disconnected entanglement entropy is shown in figure 26. It turns out that for all , the area of the connected surface is always smaller than the disconnected surface. It implies that connected surface is always more favourable, and hence no connected/disconnected phase transition takes place in the specious confined phase. Subsequently, the entanglement entropy is always of order O(N 2 ) at any subsystem size. This is one of the biggest differences between standard confined and specious confined phases. Interestingly, however, now a new type of connected/connected surface phase transition appears in the specious confined phase. This connected/connected surface phase transition is shown in figure 26, where one can clearly observe a transition between connected surfaces 1 and 3 . The critical strip size at which this phase transition occurs is indicated by c . It is important to emphasise that this connected/connected surface transition is very different from the connected/disconnected surface transition observed in the standard confined phase. In particular, the order of the entanglement entropy does not change at the connected/connected critical point, This important result further confirms that a non-trivial interpretation of small black hole phase as the gravity dual of standard confined phase is not entirely correct, as otherwise mentioned in [39,40]. In order to further highlight the subtle relation between standard confined and specious-confined phases, we also like to emphasize that although ∂S EE ∂ in the specious confined phase is not strictly zero, however it is very small. For example, for larger , the entanglement entropy depends only mildly on (as can be seen from figure 26) and is practically independent of it. When going from let's say to /2, the change in the magnitude of entanglement entropy occurs only at the fifth decimal place. This feature of entanglement entropy in the specious-confined phase again resemble approximatelyhowever, not exactly -to the standard confined phase for which ∂S EE ∂ = 0 for large , highlighting once more the non-trivial similarities as well as differences between standard confined and specious-confined phases as was first pointed out in [80].
Interestingly the entropic C-function, which quantifies the number of degrees of freedom in a system at length scale (or the energy scale), also behaves desirably in the speciousconfined phase. In particular, the C-function decreases monotonically from UV to IR in the specious-confined phase as well. This is shown in figure 27, where a sharp drop in its magnitude (shown by vertical solid lines) is explicitly evident. Interestingly, our holographic estimate for the length scale c at which C-function drop sharply is of the same order as was observed in lattice QCD. For example, at vanishing T and µ we find an estimate c = 0.931 GeV −1 whereas SU(3) lattice gauge setup suggested c ≈ 4.3 GeV −1 . Moreover, the temperature dependent behaviour of c also qualitatively matches with lattice prediction. The (T, c ) holographic phase diagram is shown in figure 28, and can be compared with the SU(2) lattice gauge theory conjecture of [44], figure 8.
We further like to emphasise that the above discussed richness in the structure of entanglement entropy in the specious-confined phase remains true with finite µ as well. In particular, there are again novel connected/connected surface transitions with a sharp decrease in the magnitude of C-function at the critical point c . With finite µ, the main difference appears in the magnitude of c , which attains a higher value for higher µ, and moreover approaches µ dependent constant value at low T .

With small black hole background: two strip
Since there is no connected to disconnected surface transition in the specious-confined phase, the corresponding phase diagram for two and higher strips is much simpler than the standard confined phase. In particular, there will not be any disconnected phases like S C and S D (see figure 6) and only the connected phases like S A and S B remain. The entanglement entropy expressions for S A and S B are again given by eq. (3.7).
The phase diagram with two strips at µ = 0 is shown in figure 29. We again find a phase transition between S A and S B phases. In particular, for a given , S A phase has the lowest entanglement entropy for large x whereas S B phase has the lowest entanglement entropy for small x. It is interesting to observe that this phase diagram is quite similar to the two strip standard confined phase diagram if we remove the disconnected phases from the latter. In particular, the transition line between S A /S B is again almost constant for large x, l > c . Moreover, the nature of the critical point c , where non-analyticity in the entanglement entropy appears, is also quite similar to the second tricritical point c2 of the standard confined phase. This once again emphasizes the closeness of specious-confined phase with the standard confined phase. However here, as opposed to the standard confined phase, entanglement entropy order always changes at the transition line. Our analysis further suggests only a mild dependence of the phase diagram on temperature. In particular, S A /S B transition line at different temperatures almost overlap with each other.
In figure 29, two strip phase diagram with finite chemical potential at T = 0.9 T c is shown. We again find a similar type of S A /S B transition as x and are varied. The main difference arises in the magnitude of c , which only gets enhanced with µ. The higher µ therefore increases the parameter space of S B . Although not presented here for brevity, similar results occur at other temperatures as well. It is also instructive to investigate the mutual information of the specious-confined phase and compared it with the standard confined phase. The results for various temperatures are shown in figure 31. I was instead found to approach a temperature dependent constant value (see figure 18). For completion, we have also included zero temperature I [n=2] B behaviour. We find that I profile for various temperatures overlap with each other, both in small as well as in large regions, thereby suggesting its non-thermal nature in the specious-confined phase. This interesting new result has not been discussed in lattice QCD community yet, and it would be interesting to perform a similar temperature dependent analysis of the mutual information using lattice simulations and compare the corresponding lattice results with the holographic prediction.
We further find that I [n=2] B smoothly goes to zero as the S A /S B transition line is approached. Moreover, the order of mutual information also changes as we go from S A to S B and visa versa. The phase transition in the specious-confined phase is therefore always accompanied by a change in the order of mutual information, as opposed to the standard confined phase where its order may or may not change depending on the nature of the transition line.
The effect of chemical potential on I

With small black hole background: n > 2 strips
The entanglement phase diagram of n > 2 strips is very similar to n = 2 strips. Again, only S A /S B type of entangling surface phase transition appears. At a fixed T and µ, larger n only increases the parameter space of S B in the small region whereas it remains almost constant in the large region.
The 3-partite information as a function of x and are shown in figures 33 and 34. We find that 3-partite information, like mutual information, shows almost no dependence on temperature and chemical potential. In particular, various temperature and chemical potential dependent profiles of the 3-partite information overlap with each other, and even the length scale where non-analyticity in the 3-partite information appears does not change with temperature and chemical potential. The temperature independent behaviour of 3partite information in the specious confined phase is therefore very different from the deconfined phase where 3-partite information was instead found to vary with temperature (see figure 21).
Similarly, we find that even the 4-partite information does not show any dependence on temperature and chemical potential. Although it is hard to explicitly establish this result for a general n, however the structure of n-strip phase diagram and the general trend do suggest that the corresponding n-partite information is independent of temperature and chemical potential as well. Again, this behaviour of 4-partite information should be contrasted from the deconfined phase 4-partite information where it does depend on temperature and chemical potential.

With large black hole background: n strips
Now, we will briefly mention the results for the dual deconfined phase which corresponds to the large black hole phase. The results for one strip (or the entanglement entropy) are shown in figures 35 and 36. There is again a one to one relation between and z * , with neither max nor crit exist. The connected entangling surface moves more and more towards the horizon as the subsystem size increases, thereby probing deeper spacetime structure, and have lower entanglement entropy than the disconnected surface at all . Therefore, no connected/disconnected phase transition and crit exist in the deconfined phase. These results are similar to the deconfined phase results obtained in the previous section using the scale factor P 1 (z). In fact, we have checked by taking other forms of the scale factor P (z) as well that these results for the entanglement entropy remain the same in the dual deconfined phase. Moreover, even in inconsistent models like soft and hard walls, similar results in the deconfined phases can be obtained. This suggests a universality in the entanglement structure of the dual deconfined phase.
Similarly, for n ≥ 2 strips, we do not find many differences from the deconfined phase results discussed in section 3.2. S A and S B phases and the corresponding phase diagram display the same features as were observed previously. In the phase diagram, higher temperature again tries to enhance the parameter space of S A whereas higher chemical potential tries to enhance the parameter space of S B . The phase transition is again accompanied by a change in the order of mutual information in the deconfined phase. The mutual information also behaves desirably and display the same features as were observed previously using P 1 (z). In particular, I again asymptotically approaches to a temperature and chemical potential dependent constant value (i.e. its large asymptotic value gets enhanced with both temperature and chemical potential). The closeness of P 2 (z) deconfined phase with the deconfined phase of P 1 (z) goes beyond the mutual information and we find that even the 3 and 4 partite information exhibit similar features in these deconfined phases.

Small/large hole phase transition and mutual information
We close this section by analysis the mutual information in T − µ plane. The objective here is to see whether the 3-partite information, like the entanglement entropy, captures the small/large black hole (or the dual specious-confined/deconfined) phase transitions just as it did for the thermal-AdS/black hole phase transition. Our results for the mutual information are shown in figure 37 for a fixed = 0.2 GeV −1 and x = 0.1 GeV −1 , although the main results of our investigation remain unchanged even for other values of and x. We again find that the structure of mutual information resembles remarkably well with the Bekenstein-Hawking entropy (see figure 23). In particular, the mutual information also exhibits three branches in T − µ plane. The branches with positive slope are stable whereas the branch with a negative slope is unstable. These three branches exist only when µ < µ c , and for µ > µ c the negative slope branch ceases to exist and we have only one branch. Importantly, a critical feature like µ c does not change when other values of and x are considered. The mutual information therefore again experiences the deviations in the spacetime geometry caused by the black hole phase transition. Moreover, we have checked that similar results also hold for 3-and 4-partite information. Our analysis therefore once again suggests that the other information theoretic quantities can also be used to probe black hole (or dual specious-confined/deconfined) phase transitions.

Conclusions
In this paper, using the holographic RT prescription of the entanglement entropy, we have investigated the mutual and n-partite information of a strongly coupled QCD theory whose dual gravitational theory is described by a consistent phenomenological bottomup Einstein-Maxwell-dilaton model. This is an extension of our previous work on the entanglement entropy [58], and the objective here was to further investigate how the mutual and n-partite information can shed new light on the confinement mechanism. The excited dilaton field in this model allowed us to modify the nature of the gravity solutions by JHEP04 (2019)137 choosing an appropriate form of the scale factor. With one form factor (P 1 (z)), we found a thermal-AdS/black hole phase transition which on the dual boundary side corresponds to the standard confined/deconfined phase transition. We then studied the entanglement entropy phase diagram by considering n disjoint intervals. In the confining background, with n ≥ 2, we found a rich phase diagram consisting of four distinct connected and disconnected surfaces in the parameter space of and x. With n ≥ 2, unlike n = 1, the order of the entanglement entropy may or may not change as the transition point is crossed. We then analyzed the mutual and n-partite information in the confining background and found that the mutual information is monogamous and that the n-partite information exhibits non-analyticity in its structure. In the deconfining background, on the other hand, only two connected phases appeared, leading to a much simpler phase diagram. We further found that higher temperature makes the parameter space of S B phase smaller whereas higher chemical potential makes the parameter space of S B larger, suggesting that one should look for the low temperature/high chemical potential region of the QCD phase diagram to find a non-trivial profile of the mutual information. Moreover, the separation length at which non-analyticity in 3-partite information appears is found to decrease with temperature whereas it enhances with chemical potential.
The second form factor (P 2 (z)) instead leads to the small/large black hole phase transition which in the dual boundary side corresponds to the specious-confined/deconfined phase transition. In this case, the entanglement structure of the deconfined phase is found to be similar to the deconfined phase entanglement structure obtained using P 1 (z). However, the entanglement structure of the specious-confined phase displayed many dissimilarities with the standard confined phase. In particular, a novel connected/connected (instead of a connected/disconnected) surface phase transition appeared with n = 1, which greatly modified its n ≥ 2 entanglement entropy phase diagram compared to the standard confined phase. The small black hole phase also allowed us to probe the effect of temperature and chemical potential on the entanglement phase diagram. We found that n ≥ 2 entanglement entropy phase diagram is almost independent of temperature and chemical potential in small < c region whereas the phase space of S B gets slightly enhanced with chemical potential in large > c region. Interestingly, the mutual information and n-partite information although behaving desirably in the specious confined phase, however, exhibit no dependence on temperature and chemical potential. This is a new and interesting prediction from holography, but unfortunately, unlike the entanglement entropy, we do not yet have any corresponding independent lattice result to compare our holographic result.
We further studied the Hawking/Page and small/large black hole phase transitions using the mutual and n-partite information and found that, just like the entanglement entropy, these quantities also capture the essence of black hole phase transitions. Since the imprints of phase transition were also seen on the mutual and n-partite information, our analysis therefore extended the number of boundary observables that can be used to probe these phase transitions.
There are several directions in which the analysis of this paper can be expanded. In present work, we have considered the oversimplified case of equal strip and separation lengths. This makes the phase space effectively two dimensional. However, for the case JHEP04(2019)137 of n unequal strips and separations lengths the phase space would be 2n − 1 dimensional. Analysis of this multi-dimensional phase space structure would although be bit tedious, however it might shed new light on the QCD phases, especially in the confining background. Another direction to extend our work is to include a background magnetic field and use the entanglement structure to investigate (inverse) magnetic catalysis. We hope to come back to these issues soon.