Holographic QCD3 and Chern-Simons theory from anisotropic supergravity
Abstract
Based on the gauge-gravity duality, we study the three-dimensional QCD (QCD3) and Chern-Simons theory by constructing the anisotropic black D3-brane solution in IIB supergravity. The deformed bulk geometry is obtained by performing a double Wick rotation and dimension reduction which becomes an anisotropic bubble configuration exhibiting confinement in the dual theory. And its anisotropy also reduces to a Chern-Simons term due to the presence of the dissolved D7-branes or the axion field in bulk. Using the bubble geometry, we investigate the the ground-state energy density, quark potential, entanglement entropy and the baryon vertex according to the standard methods in the AdS/CFT dictionary. Our calculation shows that the ground-state energy illustrates degenerate to the Chern-Simons coupling coefficient which is in agreement with the properties of the gauge Chern-Simons theory. The behavior of the quark tension, entanglement entropy and the embedding of the baryon vertex further implies strong anisotropy may destroy the confinement. Afterwards, we additionally introduce various D7-branes as flavor and Chern-Simons branes to include the fundamental matter and effective Chern-Simons level in the dual theory. By counting their orientation, we finally obtain the associated topological phase in the dual theory and the critical mass for the phase transition. Interestingly the formula of the critical mass reveals the flavor symmetry, which may relate to the chiral symmetry, would be restored if the anisotropy increases greatly. As all of the analysis is consistent with characteristics of quark-gluon plasma, we therefore believe our framework provides a remarkable way to understand the features of Chern-Simons theory, the strong coupled nuclear matter and its deconfinement condition with anisotropy.
Si-wen Li111Email: [email protected], Sen-kai Luo222Email: [email protected], Ya-qian Hu333Email: [email protected],
Department of Physics, School of Science,
Dalian Maritime University,
Dalian 116026, China
1 Introduction
While quantum chromodynamics (QCD) is the underlying theory to describe the strong interaction, it is usually very difficult to solve at low energy due to its asymptotic freedom, especially in the dense matter with finite temperature. Therefore it provides motivation to study the dynamics of strongly coupled non-Abelian gauge field theory via the gauge-gravity duality as an alternative option [1, 2]. On the other hand, the heavy-ion collision (HIC) experiments show that quark-gluon plasma (QGP) created in the collision is strongly coupled [3, 4] and anisotropic [5, 6, 7, 8], hence constructing the type IIB supergravity in order to investigate the anisotropic and strongly coupled QGP or Yang-Mills theory through gauge-gravity duality is naturally significant [9, 10, 11] since, as it has been well-known, the most famous example in the gauge-gravity duality is the corresponding between four-dimensional super Yang-Mills theory on D3-branes and type IIB super string theory on .
A remarkable work in the top-down holographic approach to study the anisotropy in gauge theory is [12] in which the black D3-brane solution in type IIB supergravity is anisotropic due to the presence of the axion field or dissolved D7-branes in the bulk. Following the AdS/CFT dictionary, the thermodynamics and transport properties in an anisotropic plasma are explored holographically by the gravity solution in [12] which attracts many interests [13, 14]. In particular, another concern in [12] is that the presented axion field leads to a theta term in the dual theory and the parameter is spatially dependent. Since -dependence involves the topological property in gauge theory, it also gets many attentions in theoretical and phenomenological researches [15]. Although the experimental value of the parameter is very small, it may influence many observable effects in gauge theories e.g. the deconfinement phase transition [16, 17], the glueball spectrum [18], the CP violation in hot QCD [19, 20], the chiral magnet effect [21, 22], the large N limit [23] and its holographic correspondence [24, 25, 26, 27]. Accordingly, the holographic duality proposed in [12] becomes a topical issue at one stage.
Keeping these in hand, in this work, we would like to study the holographic duality between the three-dimensional QCD (QCD3) and Chern-Simons theory based on [12]. As the parameter in the framework of [12] linearly depends on one of the three spatial coordinates, integrating by parts, one can get a three-dimensional Chern-Simons term as which accordingly is the part of the motivation for this work. Furthermore, QCD3 or the Chern-Simons theory involving fundamental matters with flavors and their large N ’t Hooft limit are also interesting topics especially in three-dimensional case [28, 29, 30, 31, 32, 33, 34], thus including flavors would also be our concern in this project. And the presented anisotropy might be more closed to the realistic physical situation in some materials.
However one of the key points here is to find a scheme to combine the gravity system in [12] with the three-dimensional theory in holography. Fortunately the answer could be found in the famous [35, 36] which provides the compactification in the D3-branes system in order to obtain a three-dimensional non-supersymmetric and non-conformal gauge theory, as it is successfully performed in the D4/D8 approach [37]. So by imposing the compactification method in [35, 36] to the supergravity system in [12], in this work we first obtain the bubble configuration of the bulk geometry which could be remarkably analytical if we take the compactification limit (i.e. the size of the compactification direction vanishes). Since the bubble configuration does not have a horizon, the dual theory is at zero temperature limit. And we examine the dual theory by introducing a probe D3-brane at the holographic boundary which exactly exhibits a Yang-Mills plus Chern-Simons theory as it is expected. A notable feature in our holographic setup is that the Chern-Simons level is naturally identified to the number of the D7-branes dissolved in the bulk which is automatically quantized. And we believe this provides a holographic proof to the quantization of the Chern-Simons level.
Afterwards, some of the observables are investigated by using the standard method according to the AdS/CFT dictionary in the bubble configuration of the bulk, specifically they are the ground-state energy density, quark potential, entanglement entropy and the baryon vertex. To simplify the calculation, we consider that the size of the compacted direction trends to be vanished in the bulk geometry throughout this work, so that the dual theory would become exactly three-dimensional. Then our results show that the ground-state energy is degenerate to the Chern-Simons coupling coefficient which is in agreement with the properties of the gauge transformation in the gauge Chern-Simons theory [38]. Besides, the behaviors of the quark potential and entanglement entropy depending on the position of the fundamental string or “slab” resultantly reveal that the confinement may be destroyed in hadron if the anisotropy becomes strong enough, because the entanglement entropy may also be a characteristic tool to detect the confinement [39, 40, 41, 42]. Moreover, we introduce a wrapped D5-brane on as the baryon vertex [43] in this geometry and study its embedding configuration as [44]. The numerical calculation confirms the wrapped configuration of the baryon vertex in this system and the D-brane force illustrates the bottom of the bulk is the stable position of a baryon vertex as it is expected to minimize its energy. Interestingly, the numerical calculation also displays the wrapped baryon vertex trends to become unwrapped by the increasing of the anisotropy which means the baryon vertex may not stably exist if the anisotropy becomes very large. And it is seemingly consistent with the analysis of the quark potential and entanglement entropy with respect to the confinement in this system i.e. strong anisotropy may destroy the confinement.
Last but not least, to explore the Chern-Simons topological feature involving the flavors in the dual theory, following [45, 46, 47], various D7-branes as flavor and Chern-Simons branes are introduced into the bulk bubble configuration as probes. In a transverse plane, the vacuum configuration of the D7-branes, which means the embedding function minimizes the energy of the D7-branes, is numerically evaluated and the calculation shows the vacuum structure is shifted by the presence of the axion field in the bulk. To further consider the spontaneous breaking of the flavor symmetry , we separate coincident of flavor branes living into the upper part of the transverse plane and the other coincident flavor branes living into the lower part of the plane while they extend to a same position at the holographic boundary. Therefore the interpretation of such configuration could be that the flavor symmetry spontaneously breaks down into at low energy in dual theory. Taking into account the contribution of the orientation of the flavor and Chern-Simons branes, we can get an effective flavor-dependent Chern-Simons level. Then evaluating the total energy including both flavor and Chern-Simons branes by counting the orientation in the effective Chern-Simons level, the associated topological phases in the dual theory can be obtained. A noteworthy conclusion here is that the total energy including flavor and Chern-Simons branes illustrates the topological phase transition may occur at a critical flavor mass which decreases due to the presence of the anisotropy or the axion field in the bulk geometry. By analyzing the phase diagram, it seemingly means the broken flavor symmetry , which may relates to the chiral symmetry, would become restored to if the anisotropy becomes sufficiently strong. And this behavior is also predicted by the numerical evaluation of the embedding of the flavor branes since the two branches of the flavor branes trend to become coincident when the anisotropy becomes large. Altogether, this framework may provide a holographic way to study the behavior of metastable vacua in large N QCD3 with a Chern-Simons term [48] and its deconfined condition with anisotropy, although the anisotropy is expected to be small for the numerical calculations in this project.
The outline of this manuscript is as follows. In Section 2, we briefly review the anisotropic black brane solution in the type IIB supergravity, then give our holographic setup to this work. In Section 3, we calculate several observables with respect to the constructed bulk geometry. In Section 4, we discuss the embedding of the flavor and Chern-Simons branes. In Section 5, we analyze the corresponding topological phase and its associated phase transition in holography. Summary and discussion are given in the final section. In addition, we list the relevant parts of the functions presented in the bulk geometry in the appendix which would be very useful to this work.
2 Holographic setup
2.1 Review of the anisotropic solution in type IIB supergravity
In this subsection, we review and collect the relevant content of the anisotropic solution in ten-dimensional type IIB supergravity in [12]. The remarkable anisotropic solution describes the bulk dynamics of D3-branes with D7-branes dissolved in the spacetime in the large limit and the D-brane configuration is given in Table 1.
Black brane background | ||||||
---|---|---|---|---|---|---|
D3-branes | - | - | - | - | ||
D7-branes | - | - | - | - |
As our concern would be the holographic duality, let us start with the type IIB supergravity action in string frame,
(2.1) |
where the index runs over 0 to 9, is the ten-dimensional gravitational coupling constant . To obtain an anisotropic solution, the associated equation of motion to (2.1) can be solved by the following anisotropic ansatz in string frame,
(2.2) |
where refers to the axion, dilaton and the unit volume form of a five-sphere . The parameters in the solution are given as follows,
(2.3) |
where represents the radius of the bulk, the string coupling and the ’t Hooft coupling constant respectively. The solution (2.2) describes the black branes with a horizon at and the anisotropy in direction. There would not be new field in the boundary which is located at because the D7-branes do not extend along the holographic direction . Since the dynamic of the axion , which magnetically couples to D7-branes, is taken into account, it is clear that the supergravity solution (2.2) includes the backreaction of D7-branes to the D3-brane bulk geometry. We note that the D7-branes are distributed along direction with the constant distribution density according to the solution for . So once the backreaction of D7-branes to the background geometry is included, it implies that is fixed in the large limit.
The regular functions in (2.2) depend on the holographic coordinate which must be determined by their equations of motion. However they are non-analytical in general. In order to avoid the conical singularities in the bulk, the Euclidean version of the bulk metric near the horizon,
(2.4) |
must impose the period to be . Hence it reduces to the formula of the Hawking temperature as,
(2.5) |
Suppose the temperature is sufficiently large (or equivalently ), the functions can be analytically written as the series of which are given in the Appendix. We note this high-temperature analysis would be remarkably useful in the following sections of this work.
2.2 Construction for the 2+1 dimensional theory
As it is known that the type IIB supergravity theory holographically corresponds to the super Yang-Mills theory on D3-brane, it would be very straightforward to construct the D3-brane configuration or the super Yang-Mills theory in order to obtain a non-supersymmetric and non-conformal dual theory by following the steps in the well-known [35, 36]. Specifically the first step is to take the spatial dimensions of the D3-brane to be compactified on a circle with a period . Therefore the dual theory is effectively three-dimensional below the Kaluza-Klein energy scale defined as . The second step is going to get rid of all massless fields other than the gauge fields, which is to impose respectively the periodic and anti-periodic boundary condition on bosonic and fermionic fields along . Afterwards the supersymmetric fermions and scalars acquire mass of order thus they are decoupled in the low-energy dynamics. So the dual theory below becomes three-dimensional pure gauge theory. By keeping these in mind, the next step is to identify the bulk geometry that corresponds to this gauge theory. The answer can be found by interchanging the roles of and i.e. performing a double Wick rotation to the metric presented in (2.2), which is
(2.6) |
where we have renamed after the double Wick rotation444Performing to the metric presented in (2.2) reduces to Then we rename by in order to obtain (2.6). . We note that in this notation the axion field becomes,
(2.7) |
while remains and now is periodic as,
(2.8) |
Using the formulas given in the appendix, we can obtain
(2.9) |
The solution (2.6) represents a bubble geometry of the bulk which is anisotropic on plane and defined only for . The D-brane configuration for the bubble solution (2.6) is given in Table 2.
Bubble background | ||||||
---|---|---|---|---|---|---|
D3-branes | - | - | - | - | ||
D7-branes | - | - | - | - |
Here we have renamed as since there is not a horizon in the bulk as it is illustrated in Figure 1.

Since the wrap factor never goes to zero, the dual theory would exhibit confinement according to the behavior of the Wilson loop in this bulk geometry. Besides, we can notice that the gauge theory on D3-brane becomes purely three-dimensional theory if the compactified direction shrinks to zero i.e. or . And this limit exactly corresponds to the high temperature limit in the black brane solution (2.2) so that in the limit of , the functions in (2.6) are also analytical as they are given in the appendix. Accordingly we are going to consider the case in the limit throughout this work since our concern is the three-dimensional dual theory exactly555A safe statement is to further require in the construction of the bubble geometry since the limit of would lead to pure D7-brane background which may have many issues in a holographic approach. To avoid those issues, we may consider that with fixed in the large for the gravity background produced by multiple D-branes as in [59, 60] which means the anisotropy may not be very large in this setup..
2.3 The dual theory
The dual theory with respect to the bulk solution (2.6) is defined at the zero temperature limit since there is not a horizon i.e. . To examine the dual theory, let us introduce a single probe D3-brane located at the boundary . As we have discussed, in the bubble geometry (2.6), the low-energy modes in the dual field theory contain the gauge field only, so the effective action for such a probe D3-brane is given as,
(2.10) |
where the tension of the D3-brane is given as . And is the induced metric and the gauge field strength on the worldvolume of the D3-brane. Assuming does not have components along and does not depend on , then the quadratic expansion of the action (2.10) is
(2.11) |
where is the three-dimensional ’t Hooft coupling constant and
(2.12) |
is the Chern-Simons three-form. We use to denote the gauge potential and have imposed the boundary value to (2.11). Clearly the dual theory on the D3-brane is effectively three-dimensional Yang-Mills plus Chern-Simons theory below the energy scale . Notice that the dual theory is expected to be purely three-dimensional theory if we take the limit or equivalently .
It is remarkable to notice that in this holographic setup, the level number of the Chern-Simons term is integer automatically since the level number is exactly the number of D7-branes in the gravity side. This leads to a proof of the quantization of the Chern-Simons level via holography. On the other hand, when the backreaction of the D7-branes is included in the bulk, the dual theory is equivalently a topological massive theory. This can be confirmed once we derive the formula of the propagator with respect to action (2.11) which is [38],
(2.13) |
where and we use to denote the momentum in 2+1 dimensional spacetime. Therefore, the propagator defines the topological mass of the gauge field via the pole which is determined by the numbers of D7-branes. Thus the presented D7-branes involve the topological properties of the dual theory and, we will see, they contribute to the various vacuum configurations in the dual theory.
3 Observables
In this section, we will extract relevant information on the physics of the Yang-Mills theory with a Chern-Simons term dual to the anisotropic background given in (2.6). Using the standard holographic methods, we will focus on the ground-state energy, quark potential, entanglement entropy and baryon vertex in this system.
3.1 The ground-state energy
In the holographic dictionary, one of the basic entries is the relation between the renormalized on-shell supergravity action and partition function of the dual field theory [1, 2, 35, 36]. Therefore the ground-state energy density of the dual field theory can be obtained through the relation
(3.1) |
where refers to the infinite three-dimensional Euclidean spacetime volume. is the renormalized on-shell action for type IIB supergravity, which is given by
(3.2) |
where refers to the Gibbons-Hawking term and holographic counterterm for the type IIB supergravity. refers to the Euclidean version of (2.1) which in Einstein frame is given as,
(3.3) |
Since only has components on , the onshell action can be integrated out over to become an effective five-dimensional action as,
(3.4) |
where runs over 0 to 4. The cosmological constant is , is the five-dimensional scalar curvature and is the five-dimensional gravitational coupling constant. The action (3.4) is nothing but the five-dimensional axion-dilaton-gravity action. Hence the holographic counterterm can be chosen as [12],
(3.5) |
where refers to the boundary metric and refers to the conformal anomaly in the axion-dilaton-gravity system 666The counterterm in axion-dilaton-gravity system can also be found in [49, 50] . In this sense, the metric near the boundary is required to take the form as,
(3.6) |
in which the coordinate is the standard Fefferman-Graham (FG) coordinate. The relation between the coordinate presented in (2.6) and the FG coordinate is collected as,
(3.7) |
We also note that the standard Gibbons-Hawking term is given by the trace of the extrinsic curvature of the boundary as
(3.8) |
Plugging (3.4) (3.5) (3.8) into (3.2) and using the relation (3.1), the resultant free energy density of the dual theory is computed up to as,
(3.9) |
The free energy density depending on implies the ground-state is degenerate to the Chern-Simons coupling coefficient . Since is the distribution density of the D7-branes, the value of could be same for different . Therefore the (3.9) describes the vacuum energy in various branches characterized by its level number . To any given interval, we fix the length of for possible values of , then the ground-state free energy (3.9), which implies the vacuum with distinct , is shown in Figure 2.

And the true vacuum should minimize the free energy (3.9).
Figure 2 also illustrates that the behavior of ground-state energy is similar to the same observable in QCD4 with finite theta angle, especially in some holographic approaches [24, 25, 26]. The reason is as follows. First, the black brane solution (2.2) and bubble solution (2.6) share same value of the onshell action in gravity side since the difference here between them is just the double Wick rotation. Then, on the other hand, the dual theory in the the black brane and bubble background is respectively four-dimensional theta-depended gauge theory and three-dimensional Yang-Mills-Chern-Simons given in (2.11). Therefore it is not surprised that their dual theories also share similar behavior of the ground-state energy according to the AdS/CFT dictionary (3.1) 777One may consider the dualities in string theory to study this similarity in another way, since the holographic investigation in [24, 25, 26] is based on the D4-brane approach in IIA string theory which is a T-duality version of IIB string theory.. In addition, we may find in (2.10) (2.11), the Chern-Simons three-form comes from the integral by part to the four-dimensional Wess-Zumino term in which the axion field plays exactly the role of the theta term as in QCD4. Thus the similarity in the behaviors of the ground-state energy with respect to QCD3 and QCD4 may also be understood by this integral relation.
Besides, the degeneracy to the Chern-Simons level number in the ground-state free energy via holography may also have an interpretation in terms of quantum field theory (QFT). As we have analyzed, the dual theory on D3-brane is the three-dimensional Yang-Mills-Chern-Simons whose action is given by (2.11). So under the local gauge transformation
(3.10) |
the Yang-Mills-Chern-Simons action presented in (2.11) transforms as,
(3.11) |
where is the winding number given by
(3.12) |
Therefore, under the gauge transformation, the partition function given by the Euclidean path integral is
(3.13) |
up to a phase . Comparing (3.13) with the AdS/CFT dictionary (3.1), it means the free energy must be degenerate to . This property of Chern-Simons theory also leads to an additive renormalization condition to the renormalized and bare Chern-Simons coupling coefficient
(3.14) |
up to the one-loop order calculation at least, according to [38]. It also implies the renormalized parameter (denoted by in our system) is degenerate to the bare parameter (denoted by ) up to a finite integer-valued shift.
3.2 Wilson loop and quark potential
In holography, the vacuum expectation value (VEV) of Wilson loop on a contour corresponds to the renormalized Nambu-Goto on-shell action of a fundamental open string whose endpoints span the contour [51] which is,
(3.15) |
The static quark-antiquark potential can be obtained by evaluating the Nambu-Goto action as,
(3.16) |
where refers to the bare mass of quark which is also the counterterm in the Nambu-Goto action. In the background (2.2), the bare mass can be evaluated by putting the fundamental open string extending along direction. Choosing , the induced metric on the worldsheet parametrized by is written as,
(3.17) |
Hence the bare mass is obtained by the Nambu-Goto action with respected to (3.17) as,
(3.18) |
Next in order to compute the quark-antiquark potential, we consider a rectangular contour with sides of length along and one spatial direction or . Notice the background metric is anisotropic in plane, the calculation of Wilson loop would be a little different when the string extends along and direction.
Parallel to the D7-branes
Let us first consider the situation that the open string extends parallel to the D7-branes as it is illustrated in Figure 3.

Taking the static gauge as , the Nambu-Goto action is given as,
(3.19) |
where the derivatives “ ′ ” are with respect to . Notice that the associated Hamiltonian to (3.19) is conserved i.e. a constant since the Lagrangian presented in (3.19) does not depend on explicitly. Accordingly we can reach
(3.20) |
Imposing the condition , the Hamiltonian reduces to
(3.21) |
which is equivalent to
(3.23) |
So according to (3.16), the quark-antiquark potential can be obtained by subtracting (3.18) from (3.23) as,
(3.24) |
where we have used to denote the quark potential with respect to the parallel case. The constants and depending on have to be calculated numerically and their behaviors with respect to are given in Figure 4 888One may find a tail in the behavior of and at in Figure 4 which seemingly implies they are non-monotonic functions of . The reason is that the formulas of the quark potential becomes divergent at , which is recognized as a IR divergence in the dual theory. As an effective description, we can introduce an IR cutoff by to remove the divergence (also in the numerical calculation), so that and could become monotonic functions of above the cutoff..

Perpendicular to the D7-branes
The second situation is that the open string extends vertically to the D7-branes as it is illustrated in Figure 3. In this case we take the static gauge as , then the Nambu-Goto action is given as,
(3.25) |
where the derivatives “ ′ ” are with respect to . And the associated Hamiltonian to (3.25) reduces to a constant which is,
(3.26) |
or equivalently
(3.27) |
As the analysis in the parallel case, the quark-antiquark potential can be obtained by plugging (3.27) into (3.25) then subtracting (3.18). The final result is given as,
(3.28) |
where the constant and as functions of are given in Figure 4.
As we have required that the size of the compactified direction is sufficiently small , so that the functions in (2.6) are analytical. Thus the quark tension can be evaluated analytically. In the large limit, to minimize its energy, the fundamental string trends to stretch as much as possible over . According to (3.25), its effective tension would be proportional to and the fundamental string would move approximately vertically up to UV cutoff around the extrema . Therefore in the large limit we can obtain
(3.29) |
which illustrates an area law of the Wilson loop. So the quark tension is evaluated as,
(3.30) |
As we can see the quark tension in the parallel case decreases while it increases in the perpendicular case. And this result is in agreement with our numerical calculation in the limit . While it is not strictly to discuss the case that the anisotropy becomes large, the quark tension could be vanished if increases, which may imply the deconfinement.
3.3 Entanglement entropy
Since the entanglement entropy is expected to be a probe to characterize the phase transition in the field theory, in this subsection let us compute the entanglement entropy in the anisotropic background (2.6).
We will take into account the region A and its complement, the region B, as two physically disjoint spatial regions in the dual theory. Based on the AdS/CFT correspondence [39, 52], in the dual theory the quantum entanglement entropy between region A and B is identified to be the surface stretched in the bulk whose boundary coincides with the boundary of A. In general the classical area of surface in the correspondence of AdSd+2/CFTd+1 is given as,
(3.31) |
where is the Newton constant in dimensional spacetime and is the induced metric on . In order to represent the entanglement entropy at a fixed time, surface must be space-like. A natural generalization of (3.31) to the ten-dimensional geometry in string theory is
(3.32) |
where the induced metric should be given in the string frame and the entanglement entropy can be obtained by minimizing the action (3.32). As the most simple case, we consider the “slab” geometry of A as , however it would be straightforward to realize that the resultant entanglement entropy given by (3.32) depends on the configuration of the slab as we have seen in the cases of studying the Wilson loop, since the background metric (2.6) is anisotropic in the plane. Therefore let us proceed the calculations to obtain the entanglement entropy in two cases: parallel and perpendicular case, which is similar as what we have analyzed with the setup of Wilson loop in the previous section999It also provides a parallel setup to compute the holographic entanglement entropy in anisotropic supergravity background in [53, 54]. .

Perpendicular case
Let us first deal with the case that the region A is perpendicular to the D7-branes as it is illustrated in Figure 5. The side of region reduces to a curved line in plane, hence becomes a function of in the induced metric. Impose (2.6) into (3.32), after simple calculations we obtain the action of as,
(3.33) |
where the derivatives “ ′ ” are with respect to , refers to the infinity volume of the two-dimensional worldvolume and
(3.34) |
The associated Hamiltonian to (3.33) must be a constant since the Lagrangian presented in (3.33) is independent on . Thus it leads to,
(3.35) |
i.e.
(3.36) |
where refers to the width of the region A. Obviously, the entanglement entropy given by (3.33) is divergent since it is proportional to the area of region A which therefore has to be renormalized. The “counterterm” can be obtained by evaluating the action of a surface extending along . Afterwards, the finite entanglement entropy follows the formulas as,
(3.37) |
By varying the value of , the relation of and , and is illustrated numerically in Figure 6 which displays the typical swallow-tail behavior with respect to and .




Notice that the critical value of , denoted by satisfying decreases when increases. Although our numerical calculation may be exactly valid for small anisotropy, it implies the deconfinement phase transition may trend to be vanished when increases if the entanglement entropy does characterize the confinement as in [39, 40, 41, 42].
Parallel case
Let us turn to the case that the region A is parallel to the D7-branes as it is illustrated in Figure 5. Using (3.32), the action of is given as,
(3.38) |
where the derivatives “ ′ ” are with respect to . Then we can obtain the relation
(3.39) |
where since the associated Hamiltonian is constant. Thus the the width is given as,
(3.40) |
By subtracting the divergence in (3.38), the finite part of the entanglement entropy is given as,
(3.41) |
And the relation of and , and is also illustrated in Figure 6. While the numerical calculation shows the typical swallow-tail behavior with respect to and , the associated phase transition trends to be vanished since there would not be a critical satisfying if anisotropy becomes sufficiently large. And again this conclusion is seemingly in agreement with the analysis of the Wilson loop if the entanglement entropy characterizes the confinement.
3.4 Baryon vertex
In the gauge-gravity duality, the baryon vertex is identified as a probe D-brane wrapped on the additional dimensions denoted by the spherical coordinates with open strings [43], and it can be treated as operator to create the baryon state101010In the gauge-gravity duality, the baryon state is created by quantizing the baryonic brane somehow. While it may be a little tricky in IIB string theory, one could review the quantization of the baryonic brane in IIA string theory by the approach of instanton [61] and matrix model [62].. Accordingly, in the type IIB supergravity on , the baryon vertex is a D5-brane wrapped on .
To search for the stable wrapped configuration of D5-brane in the background (2.6), it would be straightforward to investigate the condition of the force balance for the baryon vertex. To begin with, let us decompose the metric on by the coordinates of polar angle and , then define the radius coordinate as,
(3.42) |
Therefore the metric (2.6) becomes,
(3.43) |
where
(3.44) |
which implies and refers respectively to the radius and polar angle in plane. Since the baryon vertex D5-brane extends along the directions of , the induced metric on a probe D5-brane is,
(3.45) |
where the derivatives “ ′ ” are with respect to . To include the baryon potential, we turn on a single component of the gauge field on the D5-brane as , then the effective action for such a D5-brane is,
(3.46) |
where and
(3.48) |
the equation of motion for can be written as,




(3.49) |
So the displacement is solved as,
(3.50) |
and the associated Hamiltonian to (3.47) is,
(3.51) |
Accordingly the force at the cusp of D5-brane is given as,
(3.52) |
where is the standard radius coordinate of the bulk and we have used the index to refer to the value of the variables at . For stable configuration located at , the D5-brane must satisfy the zero-force condition i.e. if there is no other probe brane. It implies the zero-force condition can be achieved only if since is always positive in the bulk. Therefore the stable position for the baryon vertex in this system is located at to consistently minimize its action. To confirm our analysis, we also solve numerically the equation of motion associated to the Lagrangian presented in (3.47) to obtain the embedding function of the probe D5-brane111111To obtain the equation of motion for the embedding function, we can set in the Lagrangian presented in (3.47).. For a stable solution, we impose the boundary condition to the equation of motion, then the numerical solution for the configuration of the baryon vertex in plane is given in Figure 7. As we can see, there always exists wrapped solution for the baryon vertex and the configuration of is nearly independent on the variation of . Thus the stable position for the D5-brane is indeed if the baryon vertex is the only probe brane.
To close this subsection, let us evaluate the baryon mass in this holographic system. Since the stable baryon vertex must be located at , its mass can be obtained by evaluate its onshell action (3.46) by setting . Thus the Euclidean action of the D5-brane and baryon mass are identified via holography as,
(3.53) |
Plugging the metric in (3.45) and the relation of in (2.9) into (3.53), the baryon mass is evaluated as,
(3.54) |
which is enhanced in the presence of the Chern-Simons term represented by . As we have specified in the Section 2, the dual theory is a topological massive theory due to the presence of the Chern-Simons term, so it would be easy to understand that baryon would also become topologically massive once the fundamental fermion is introduced. In addition, as the baryon mass is proportional to the worldvolume of the D5-brane (3.53), we find Figure 7 also illustrates the worldvolume is increased by the anisotropy which implies the increase of the baryon mass by the anisotropy.
4 Embedding of the D7-branes and the vacuum structure
In this section, we continue the holographic setup by introducing the various D7-branes to identify the dual field theory as QCD3 in large limit. We first address the D7-branes as the flavor degrees of freedom then take into account another D7-brane in which the low-energy effective theory on its worldvolume is expected to be pure Chern-Simons theory. And we will see, the topological property in the dual theory can be studied by analyzing the configuration and orientation of these distinct D7-branes.
4.1 The flavor brane
As most works about the gauge-gravity duality, the fundamental matter can be added to the D3-brane background by embedding flavor D7-branes as probe. In our setup, we add coincident copies of D7-branes as probes to the D3-brane background (2.6), transverse to the compactified direction, spanning the denoted by , the holographic direction denoted by and four of the five directions in as [45, 46]. The D-brane configuration including various D7-branes is illustrated in Table 3 where we have decomposed the directions of as and .
Bubble background | |||||||
---|---|---|---|---|---|---|---|
D3-branes | - | - | - | - | |||
D7-branes | - | - | - | - | - | ||
D7-branes | - | - | - | - | - | ||
CS D7-branes | - | - | - | - | - |
The leftover direction is transverse to both the color D3-branes and flavor D7-branes, so a bare mass for the flavors can be introduced by imposing a separation between color and flavor branes along direction at the UV boundary which breaks the parity in QCD3. In the D3-D7 approach, the direction corresponding to the scalar in the D7-brane worldvolume couples to the mass operator of the fermions. And according to the gauge-gravity duality, the profile along the transverse direction of the flavor branes corresponds to the meson operator in the dual theory.
Then let us investigate the embedding of the flavor branes with their effective action. To specify the embedding of the flavor branes with respect to the transverse direction , we first define a new radial coordinate as,
(4.1) |
Up to order of , the relation of and can be solved as,
(4.2) |
We note that the ambiguity in this relation can be omitted by choosing for . In this coordinate, the background metric (2.6) can be written as,
(4.3) |
where . Afterwards we impose the coordinate transformation i.e. where is one angular coordinate in , then the metric (4.3) takes the final form as,
(4.4) |
Since the flavor D7-branes extend along , the induced metric on the flavor branes is obtained as,
(4.5) |
where and the derivatives “ ′ ” are with respect to . So the action for a single flavor D7-brane is,
(4.6) |
The behavior of the embedding function of the flavor D7-branes can be obtained by solving the associated equation of motion to the action in (4.6), which is,
(4.7) |
Since a parity transformation acts as , for the massless case, we have to impose the boundary condition
(4.8) |
Keeping these in hand, we numerically solve the equation of motion (4.7) with various in the region which is illustrated in Figure 8.


The numerical calculation shows a very small shift with respect to the dependence on (we have enlarged the shift in the figure) and this behavior is opposite to the approach of the D3-D(-1) background [45]. We note here the two branches of the flavor branes trend to become coincident if increases. The full configuration of the embedding D7-branes is given in Figure 8 and we also calculate the massive case by setting the boundary condition in Figure 9.


4.2 The Chern-Simons brane
Due to the presence of the axion field in the bulk, the D3-D7 brane background (2.6) corresponds to the QCD3 with a Chern-Simons term in holography. Thus the present D7-branes should contribute to some vacuum properties of the dual theory. However, once we evaluate the embedding function of such a D7-brane as probe, its equation of motion implies the only solution for and its onshell action automatically vanishes. To include the set-up of a pure Chern-Simons theory, we follow the discussion in the D3-brane approach [46], to introduce coincident Chern-Simons (CS) as probe D7-branes where the configuration is given in Table 3. At very low energies, all other excitations on the Chern-Simons branes decouple and only a Wess-Zumino term is left as,
(4.9) |
Therefore we can see the gauge-gravity duality in this setup reduces to the well-known level/rank duality in quantum field theory (QFT) expectations precisely [45, 46, 56].
To obtain the embedding of the Chern-Simons brane, we need to evaluate the equation of motion by its effective action. Recall the induced metric on a D7-brane (4.5), the action for Chern-Simons brane takes the same formula as given in (4.6), thus its associated equation of motion is given in (4.7) while must be a constant for a Chern-Simons brane. Accordingly, we impose the following ansatz
(4.10) |
to the equation of motion of which reduces to a constraint for and as,
(4.11) |
where we have set . And this constraint can be solved numerically as it is illustrated in Figure 10.


The numerical calculation also illustrates a very small shift with respect to the dependence on for the embedding of the Chern-Simons brane. Then the ground energy of the Chern-Simons brane can be evaluated by its action as
(4.12) |
The action is minimized at , so we can obtain the ground energy density of the Chern-Simons brane is obtained by
(4.13) |
5 The phase diagram involving the massive flavors
In the previous sections, we have evaluated the embedding of the Chern-Simons and flavor branes with massless boundary condition. Here we are going to study the configurations having both Chern-Simons and flavor branes with massive boundary condition since the vacuum of QCD3 with Chern-Simons term would include both of them in general. Then the phase diagram would be obtained by evaluating the energies of these D7-branes.
5.1 The energy of massive embedding flavor brane
In order to obtain the energy of the flavor branes with massive boundary condition, let us take a look at the asymptotic behavior of the embedding functions which satisfies the equation of motion (4.7), although the embedding configuration of the flavor branes has been illustrated in Figure 9. Since the flavor mass corresponds to the boundary value of , we can find the asymptotics of (4.7) at () as,
(5.1) |
where the relation (4.2) and boundary behavior of have been imposed. For massless case, the general form of the asymptotics at large for would be,
(5.2) |
where is a constant energy scale. So a simple way to obtain the asymptotics at large with massive boundary condition is to consider a very small variation of ,
(5.3) |
due to
(5.4) |
Keeping these in hand, then let us investigate the associated variation in the on-shell action of the flavor brane which is given by recalling (4.6) (4.7)
(5.5) |
Notice the relation of the flavor brane energy and on-shell action is , using (5.5), the contribution of the massive part to the flavor brane energy (density) is evaluated as,
(5.6) |
where we have introduced a energy scale and “” corresponds to the negative/positive mass of the flavor as it is illustrated in Figure 9. Afterwards, the total energy of flavor brane with massive boundary condition can be obtained by its massless part of energy plus the massive contribution as,
(5.7) |
Thus the flavor condensate in this system can be obtained as,
(5.8) |
for positive/negative mass.
To close this subsection, let us evaluated the total energy of the D-brane configuration that of flavor branes extend in the upper plane while the other flavor branes extend in the lower plane, as it is illustrated in Figure 11. Since the energy of each flavor brane should be equivalent, for the massless case, the total energy is given by
(5.9) |
Then for the massive case, suppose the flavor branes have a common mass , the degeneracy between the flavor branes extending in upper and down lower plane is lifted for , therefore we could get the total energy as,
(5.10) |
We note that this formula of the total energy in this D-brane configuration would be useful to study the various phase in the dual theory.
5.2 The topological phase
In order to identify the vacuum structure of the D7-branes with the various phases in the dual theory, we need to give a well-defined Chern-Simons level in which the Chern-Simons level depends on the number of the flavor branes. The details to obtain this goal have been discussed in [45, 46], so let us briefly outline the main idea and investigate the phase transition in our holographic setup.
First the effective Chern-Simons level according to the holographic duality is given as,
(5.11) |
where is a circle whose location is among all other coordinates except . Thus it is a fixed point in the plane. Then we define the “ sector”, that is a D-brane configuration with flavor branes extending in the upper plane and flavor branes extending in the lower plane. Afterwards, we take into account the contribution to by counting the number of the orientation in the plane of the D7-branes. For instance, the contribution to the effective Chern-Simons number of D7-branes with counterclockwise/clockwise orientation in the plane is positive/negative respectively, as it is illustrated in Figure 11 (while we show the massless case, it would be same for the massive case.).

Therefore the effective Chern-Simons level with respected to the orientation reads,
(5.12) |
where is the number of the Chern-Simons branes which is given as since both branches of flavor branes at the intersection point count one-half i.e. .
At low energy, the interpretation of such a D-brane configuration in holography is that the flavor symmetry spontaneously breaks down to . So Goldstone bosons are created and the associated target space is Grassmann,
(5.13) |
On the other hand, since the presence of the Chern-Simons leads to a level/rank duality , the dynamics of a sector at low-energy takes the symmetry,
(5.14) |
in which the vacuum of the dual theory is described by the sectors via holography. Afterwards, we can investigate the phase diagram by evaluating the total energy including both flavor and Chern-Simons branes. Recall (4.13) and (5.10), the total energy including flavor and Chern-Simons branes is collected as,
(5.15) |
The phase diagram can be obtained by minimizing the energy in (5.15). To find the dependence on the anisotropy (denoted by ), we require that the value of is fixed when we minimize . Then the results are collected as, for ,
(5.16) |
where refers to the symmetry group of the corresponding topological phase. And for , the associated topological phase and free energy are collected as,
(5.17) |
where refers to the critical value of the mass when the phase transition occurs, given as,
(5.18) |
Here a noteworthy feature is that the critical mass may become vanished if the axion field or the anisotropy in the bulk becomes very non-negligible i.e becomes sufficiently large. In this sense the Grassmann phase would not exist which seemingly means the broken flavor symmetry is restored. Although our setup may be exactly valid only for small anisotropy, this result is instructively suggestive to study the anisotropic behavior of the metastable vacua in QCD3 via holography.
6 Summary and discussion
In this work, we construct the anisotropic black D3-brane solution in IIB supergravity [12] then obtain the anisotropic bubble configuration for QCD3 with a Chern-Simons term due to the presence of the axion field. The analytical formulas for the the background geometry is available since the dual theory is exactly three-dimensional theory in the compactification limit, as it is expected. With this analytical bulk geometry, we investigate the ground-state energy density, quark potential, entanglement entropy and baryon vertex in the dual theory according to the AdS/CFT dictionary. Technically, we consider small anisotropy to avoid the difficulty in our numerical calculation. Then all the results show the dependence of the axion field or the anisotropy in bulk as it is expected. Afterwards we introduce various probe D7-branes as flavor and Chern-Simons branes to include flavor matters and topological numbers in the dual theory. By examining the embedding functions and counting the orientation of these D7-branes, we obtain the vacuum energy associated to the corresponding effective Chern-Simons level, hence the phase transition can be achieved by comparing the various vacuum energies.
To close this work, let us give some comments about this project. Due to the presence of the axion in bulk, the quark potential and entanglement entropy are shifted as some holographic studies in four-dimensional QCD with an axion e.g. [23, 26]. However our work additionally implies the quark tension and the potential phase transition illustrated in the behavior of the entanglement entropy could be destroyed in the presence of strong anisotropy. These can be found in Section 3.2 and Section 3.3: when the anisotropy increases, we can see the quark tension trends to become vanished and there would not be a critical value of satisfying the entanglement entropy . As the entanglement entropy could be a tool to characterize the confinement [39, 40, 41, 42], this behavior implies there would be no phase transition for i.e. no confinement for strong anisotropy. Besides, the baryon vertex also reveals the unwrapped trend when the anisotropy becomes large and the “unwrapped baryon vertex” also means deconfinement [44]. In a word, this holographic approach shows us the confinement can not maintain in an extremely anisotropic situation. Interestingly, this conclusion is in agreement with the fact that the QGP is anisotropic and deconfined, so it may provide a holographic way to understand the features of the strongly coupled matter with anisotropy.
On the other hand, the dependence on the axion or the anisotropy of the critical mass in the topological phase transition would also be a parallel computation to the D3-(D-1) approach [45] and the extension of [46] by including an dynamical axion field. Moreover, as the critical mass trends to be vanished when the anisotropy increases, it means the Grassmann target space would not exist. Namely the broken flavor symmetry would be restored if the anisotropy becomes sufficiently large. We note that this behavior is also illustrated in Figure 8. As we can see in Figure 8 the flavor branes in the upper and lower branches trend to be coincident if increases greatly, i.e. the flavor symmetry would be restored to if . Accordingly, this behavior implies the flavor symmetry, which is related to the chiral symmetry, would be restored when the anisotropy is extremely strong.
Combine the above together, this holographic system reveals a potential conclusion that is the confinement will not maintain and the flavor symmetry (or probably chiral symmetry) would be restored in an extremely anisotropic situation. This conclusion adheres to intuition, because if the anisotropy becomes very large for a fixed as , the dual theory depending on would include modes above the scale thus the dual theory is decompactified and non-confining according to [38, 39]. Remarkably, all the analyses are exactly coincident with the characteristic properties of QGP, particularly it is usually to be treated as the fundamental assumption to study the deconfined matter in holography [57, 58, 59, 60]. Therefore, while we can only work out numerically the case that the anisotropy is small in this project, our framework would be very instructive to study QCD and Chern-Simons theory with anisotropy.
Acknowledgements
We would like to thank Niko Jokela for helpful discussion. This work is supported by the National Natural Science Foundation of China (NSFC) under Grant No. 12005033, the research startup foundation of Dalian Maritime University in 2019 under Grant No. 02502608 and the Fundamental Research Funds for the Central Universities under Grant No. 3132022198.
Appendix: The analytical formulas for the anisotropic background
In the high temperature limit , . The functions presented in the anisotropic black brane background (2.2) can be written as a series of up to as [12],
(A-1) |
where
(A-2) |
We note that (A-1) and (A-2) is in fact a series of (or equivalently ), so the high temperature limit exactly refers to the case in [12] which corresponds to or in (A-1) (A-2). In the black brane background (2.2), refers to the horizon. In the bubble background (2.6), the double wick rotation reduces to the replacement in the black brane solution. Note that is replaced by in the bubble background (2.6) since refers to the bottom of the bulk instead of a horizon as it is illustrated in Figure 1. And the formulas of remain while we replace by . Clearly the high temperature limit in the black brane solution corresponds to the limit of dimension reduction in the bubble solution i.e. the limit for that the compactified direction shrinks to zero in (2.6) i.e. , or equivalently .
On the other hand, the dual theory on the D3-branes in the bubble background (2.6) is effectively three-dimensional below the energy scale . So if we take , the dual theory would become exactly three-dimensional theory. Therefore, the analytical formulas (A-1) (A-2) for functions can always be employed by replacing in the bubble background (2.6) since the dual theory is always expected to be exactly three-dimensional theory, which means is always expected even if becomes large but fixed. In a word, using (A-1) (A-2) in the bubble background means the dual theory is exactly three-dimensional theory. In this sense, we believe our analysis in this work with (A-1) (A-2) is also valid for large under .
References
- [1] O. Aharony, S. S. Gubser, J. M. Maldacena, H. Ooguri and Y. Oz, “Large N field theories, string theory and gravity”, Phys. Rept. 323 (2000) 183, arXiv:hep-th/9905111.
- [2] E. Witten, “Anti-de Sitter space and holography”, Adv.Theor.Math.Phys. 2 (1998) 253-291, arXiv:hep-th/9802150.
- [3] E. Shuryak, “Why does the quark gluon plasma at RHIC behave as a nearly ideal fluid?”, Prog.Part.Nucl.Phys. 53 (2004) 273-303, arXiv:hep-ph/0312227.
- [4] E. Shuryak, “What RHIC experiments and theory tell us about properties of quark-gluon plasma? ”, Nucl.Phys.A 750 (2005) 64-83, arXiv:hep-ph/0405066.
- [5] W. Florkowski, R. Ryblewski, “Highly-anisotropic and strongly-dissipative hydrodynamics for early stages of relativistic heavy-ion collisions”, Phys.Rev.C 83 (2011) 034907, arXiv:1007.0130.
- [6] R. Ryblewski, W. Florkowski, “Non-boost-invariant motion of dissipative and highly anisotropic fluid”, J.Phys.G 38 (2011) 015104, arXiv:1007.4662.
- [7] M. Martinez, M. Strickland, “Dissipative Dynamics of Highly Anisotropic Systems”, Nucl.Phys.A 848 (2010) 183-197, arXiv:1007.0889.
- [8] M. Martinez, M. Strickland, “Non-boost-invariant anisotropic dynamics”, Nucl.Phys.A 856 (2011) 68-87, arXiv:1011.3056.
- [9] D. Giataganas, U. Gürsoy, J. Pedraza, “Strongly-coupled anisotropic gauge theories and holography”, Phys.Rev.Lett. 121 (2018) 12, 121601, arXiv:1708.05691.
- [10] E. Banks, “Phase transitions of an anisotropic N=4 super Yang-Mills plasma via holography”, JHEP 07 (2016) 085, arXiv:1604.03552.
- [11] D. Ávila, D. Fernández, L. Patiño, D. Trancanelli, “Thermodynamics of anisotropic branes”, JHEP 11 (2016) 132, arXiv:1609.02167.
- [12] D. Mateos, D. Trancanelli, “Thermodynamics and Instabilities of a Strongly Coupled Anisotropic Plasma”, JHEP 07 (2011) 054, arXiv:1106.1637.
- [13] A. Rebhan, D. Steineder, “Violation of the Holographic Viscosity Bound in a Strongly Coupled Anisotropic Plasma”, Phys.Rev.Lett. 108 (2012) 021601, arXiv:1110.6825.
- [14] D. Mateos, D. Trancanelli, “The anisotropic N=4 super Yang-Mills plasma and its instabilities”, Phys.Rev.Lett. 107 (2011) 101601, arXiv:1105.3472.
- [15] E. Vicari, H. Panagopoulos, “Theta dependence of SU(N) gauge theories in the presence of a topological term”, Phys. Rept. 470, 93 (2009), arXiv:0803.1593.
- [16] M. D’Elia, F. Negro, “Theta dependence of the deconfinement temperature in Yang-Mills theories”, PhysRevLett.109.07200, arXiv:1205.0538.
- [17] M. D’Elia, F. Negro, “Phase diagram of Yang-Mills theories in the presence of a theta term”, Phys. Rev. D 88, 034503 (2013), arXiv:1306.2919.
- [18] L.Debbio, G. Manca, H. Panagopoulos, A. Skouroupathis, E. Vicari, “Theta-dependence of the spectrum of SU(N) gauge theories”, JHEP 0606:005,2006, [arXiv:hep-th/0603041].
- [19] D. Kharzeev, R.D. Pisarski, M. H. G. Tytgat, “Possibility of spontaneous parity violation in hot QCD”, Phys.Rev.Lett. 81 (1998) 512-515, arXiv:hep-ph/9804221.
- [20] K. Buckley, T. Fugleberg, A. Zhitnitsky, “Can theta vacua be created in heavy ion collisions?”, Phys.Rev.Lett. 84 (2000) 4814-4817, arXiv:hep-ph/9910229.
- [21] K. Fukushima, D. Kharzeev, H. Warringa, “The Chiral Magnetic Effect”, PhysRevD.78.074033, arXiv:0808.3382.
- [22] D. Kharzeev, “The Chiral Magnetic Effect and Anomaly-Induced Transport”, j.ppnp.2014.01.002, arXiv:1312.3348.
- [23] E. Witten, “Theta Dependence In The Large N Limit Of Four-Dimensional Gauge Theories”, Phys.Rev.Lett.81:2862-2865,1998, arXiv:hep-th/9807109.
- [24] C. Wu, Z. Xiao, D. Zhou, “Sakai-Sugimoto model in D0-D4 background”, Phys.Rev.D 88 (2013) 2, 026016, arXiv:1304.2111.
- [25] L. Bartolini, F. Bigazzi, S. Bolognesi, A. Cotrone, A. Manenti, “Theta dependence in Holographic QCD”, JHEP 02 (2017) 029, arXiv:1611.00048.
- [26] F. Bigazzi, A. Cotrone, R. Sisca, “Notes on Theta Dependence in Holographic Yang-Mills”, JHEP 08 (2015) 090, arXiv:1506.03826.
- [27] S. Li, “A holographic description of theta-dependent Yang-Mills theory at finite temperature”, Chin.Phys.C 44 (2020) 1, 013103, arXiv:1907.10277.
- [28] S. Giombi, S. Minwalla, S. Prakash, S. Trivedi, S. Wadia, X. Yin, “Chern- Simons Theory with Vector Fermion Matter”, Eur. Phys. J. C 72, 2112 (2012), arXiv:1110.4386.
- [29] O. Aharony, G. Gur-Ari, R. Yacoby, “d=3 Bosonic Vector Models Coupled to Chern- Simons Gauge Theories”, JHEP 1203 (2012) 037, arXiv:1110.4382.
- [30] C. -M. Chang, S. Minwalla, T. Sharma, X. Yin, “ABJ Triality: from Higher Spin Fields to Strings”, J.Phys.A 46 (2013) 214009, arXiv:1207.4485.
- [31] S. Jain, S. P. Trivedi, S. R. Wadia and S. Yokoyama, “Supersymmetric Chern-Simons Theories with Vector Matter”, JHEP 10 (2012) 194, arXiv:1207.4750.
- [32] O. Aharony, G. Gur-Ari and R. Yacoby, “Correlation Functions of Large N Chern- Simons-Matter Theories and Bosonization in Three Dimensions”, JHEP 02 (2013), arXiv:1207.4593.
- [33] G. Gur-Ari, R. Yacoby, “Correlators of Large N Fermionic Chern-Simons Vector Models”, JHEP 02 (2013) 150, arXiv:1211.1866.
- [34] S. Jain, S. Minwalla, S. Yokoyama, “Chern Simons duality with a fundamental boson and fermion”, JHEP 1311 (2013) 037, arXiv:1305.7235.
- [35] E. Witten, “Anti-de Sitter space, thermal phase transition, and confinement in gauge theories”, Adv. Theor. Math. Phys. 2 (1998), 505-532, arXiv:hep-th/9803131.
- [36] K. Becker, M. Becker, J. Schwarz, “String Theory and M-Theory, A Modern Introduction”, Cambridge University Press, Cambridge, 2007.
- [37] T. Sakai, S. Sugimoto, “Low energy hadron physics in holographic QCD”, Prog.Theor.Phys. 113 (2005) 843-882, arXiv:hep-th/0412141.
- [38] G. V. Dunne, “Aspects of Chern-Simons theory”, arXiv:hep-th/9902115.
- [39] I. R. Klebanov, D. Kutasov, A. Murugan, “Entanglement as a probe of confinement”, Nucl.Phys.B 796 (2008) 274-293, arXiv:0709.2140.
- [40] J. Knaute, B. Kämpfer, “Holographic Entanglement Entropy in the QCD Phase Diagram with a Critical Point”, Phys.Rev.D 96 (2017) 10, 106003, arXiv:1706.02647.
- [41] N. Jokela, J. G. Subils, “Is entanglement a probe of confinement?”, JHEP 02 (2021) 147, arXiv:2010.09392.
- [42] M. A. Akbari, M Lezgi, “Holographic QCD, entanglement entropy, and critical temperature”, Phys.Rev.D 96 (2017) 8, 086014, arXiv:1706.04335.
- [43] E. Witten, “Baryons And Branes In Anti de Sitter Space”, JHEP 9807:006 (1998), arXiv:hep-th/9805112.
- [44] Y. Seo, S. Sin, “Baryon Mass in medium with Holographic QCD”, JHEP 04 (2008) 010, arXiv:0802.0568.
- [45] S. Li, S. Luo, M. Tan, “Three-dimensional Yang-Mills-Chern-Simons theory from a D3-brane background with D-instantons”, Phys.Rev.D 104 (2021) 6, 066008, arXiv:2106.04038.
- [46] R. Argurio, A. Armoni, M. Bertolini, F. Mignosa, P. Niro, “Vacuum structure of large N QCD3 from holography”, JHEP 07 (2020) 134, arXiv:2006.01755.
- [47] D. Hong, H. Yee, “Holographic aspects of three dimensional QCD from string theory”, JHEP 05 (2010) 036, JHEP 08 (2010) 120 (erratum), arXiv:1003.1306.
- [48] A. Armoni, T. Dumitrescu, G. Festuccia, Z. Komargodski, “Metastable vacua in large-N QCD3” , JHEP 01 (2020) 004, arXiv:1905.01797.
- [49] D. Mateos, R. C. Myers and R. M. Thomson, “Thermodynamics of the brane”, JHEP 0705, 067 (2007), arXiv: hep-th/0701132.
- [50] I. Papadimitriou, “Holographic Renormalization of general dilaton-axion gravity”, JHEP 1108, 119 (2011), arXiv:1106.4826.
- [51] J. Maldacena, “Wilson loops in large N field theories”, Phys.Rev.Lett. 80 (1998) 4859-4862, arXiv:hep-th/9803002.
- [52] S. Ryu, T. Takayanagi, “Holographic derivation of entanglement entropy from AdS/CFT”, Phys.Rev.Lett. 96 (2006) 181602, arXiv:hep-th/0603001.
- [53] M. Rahimi, M. Ali-Akbari, “Holographic Entanglement Entropy Decomposition in an Anisotropic Gauge Theory”, Phys.Rev.D 98 (2018) 2, 026004, arXiv:1803.01754.
- [54] I. Aref’eva, A. Patrushev, P. Slepov, “Holographic entanglement entropy in anisotropic background with confinement-deconfinement phase transition”, JHEP 07 (2020) 043, arXiv:2003.05847.
- [55] B. Gwak, M. Kim, B. Lee, Y. Seo, S. Sin, “Holographic D Instanton Liquid and chiral transition”, Phys.Rev.D 86 (2012) 026010, arXiv:1203.4883.
- [56] P. Hsin, N. Seiberg, “Level/rank Duality and Chern-Simons-Matter Theories”, JHEP 09 (2016) 095, arXiv:1607.07457.
- [57] O. Bergman, G. Lifschytz, M. Lippert, “Holographic Nuclear Physics”, JHEP 11 (2007) 056, arXiv:0708.0326.
- [58] S. Li, A. Schmitt, Q. Wang, “From holography towards real-world nuclear matter”, Phys.Rev.D 92 (2015) 2, 026006, arXiv:1505.04886.
- [59] S. Li, T. Jia, “Dynamically flavored description of holographic QCD in the presence of a magnetic field”, Phys.Rev.D 96 (2017) 6, 066032, arXiv:1604.07197.
- [60] F. Bigazzi, A. Cotrone, “Holographic QCD with Dynamical Flavors”, JHEP 01 (2015) 104, arXiv:1410.2443.
- [61] H.Hata, T. Sakai, S. Sugimoto, S. Yamato, “Baryons from instantons in holographic QCD”, Prog.Theor.Phys. 117 (2007) 1157, arXiv: hep-th/0701280.
- [62] K. Hashimoto, N. Iizuka, P. Yi, “A Matrix Model for Baryons and Nuclear Forces”, JHEP 10 (2010) 003, arXiv:1003.4988.