Critical point in a holographic defect field theory

We study a holographic gauge theory dual to the D3/D5 intersection. We consider a pure gauge B-field flux through the internal two-sphere wrapped by the probe D5--brane, which corresponds to a non-commutative configuration of adjoint scalars. There is a domain wall separating the theory into regions with different ranks of the adjoint group. At zero temperature the theory is supersymmetric and at finite temperature there is a critical point of a second order phase transition. We study the corresponding critical exponents and find that the second derivatives of the free energy, with respect to the bare mass and the magnetic field, diverge with a critical exponent of -2/3.


Introduction
Phase transitions are ubiquitous in nature and key to the understanding of a vast amount of phenomena, ranging from the description of boiling water to the comprehension of exotic states of matter such as quark-gluon plasma or quantum spin liquids. The modern classification distinguishes two major types of phase transitions: First order phase transitions are characterised by the release or absorption of a latent heat. Typically these transitions are abrupt and powerful (large latent heat per volume). Across such a phase transition a first derivative of the thermodynamic potential has a discrete jump proportional to the latent heat. An important feature of the first order phase transitions is that their characteristics depend on the microscopic details of the system.
The second order phase transitions (also called continuous phase transitions) are characterised by a diverging susceptibility (a second derivative of the free energy), scale invariance and infinite correlation length. The scale invariance washes out the microscopic details of the system allowing a classification of the continuous phase transitions into universality classes characterised entirely by the critical exponents of the divergent thermodynamic quantities. These remarkable properties make second order phase transition the subject of a considerable interest.
In this paper we will uncover the existence of a second order critical point in a flavoured defect field theory holographically dual to the D3/D5-brane intersection. It is well known that at finite temperature this system exhibits a first order meson melting phase transition [1] (see also refs. [2] and [3]). This meson melting transition is a close analogue of the second ordered confinement/deconfinement phase transition in Quantum Chromodynamics. It is believed that the first order of the meson melting transition in many holographic systems (including the D3/D5 system) is an artefact of the large N limit. In the dual geometry this is represented by the existence of a topology changing merger transition of the probe D5-brane. In more details it is a transition between embeddings closing above the horizon of the geometry (referred to as Minkowski embeddings) and embeddings that reach all the way to the horizon (Black hole embeddings).
Having mentioned the appealing properties of second order phase transitions it is natural to attempt to amend the holgraphic set-up so that the meson melting transition is continuous. One way to achieve this is to place the theory at a finite baryon density [4], corresponding to the introduction of a 10D electric field. Charge conservation requires a number of fundamental strings to be attached to the Mikowski probe D-branes to carry the charge into the AdS black hole, which is dynamically unstable. Remarkably, [4] the probe D-brane develops throats/spikes mimicking fundamental strings and Mikowski embeddings are approximately realised within the black hole embeddings. The phase transition no longer corresponds to a topology change transition and a critical point of a second order phase transition naturally arises for a given parameter set.
A similar mechanism, which is a magnetic analogue of the baryon density case, was considered in ref. [5], where a D5-brane probe in global AdS 5 ×S 5 space-time was considered. The dual field theory has fundamental flavours confined to a maximal two-sphere within the three-sphere where the field theory lives. Introducing external magnetic field to the dual field theory corresponds to the introduction of a magnetic monopole [6] in the world volume of the probe D5-brane, which prohibits embeddings wrapping a shrinking two-sphere within the AdS 5 subspace, unless D3-branes sourcing the monopole are attached. Remarkably these D3-branes are dynamically realised as throats by Mikowski D5-brane embeddings. Once again a topology change transition is avoided and a (quantum) critical point of a second order phase transition is realised.
In this paper we consider yet another approach to avoid the topology change transition and this is to consider a pure gauge B-field flux through the internal S 2 wrapped by the D5-brane probe within the compact S 5 subspace. This setup was first considered in refs. [8,9], where the meson spectrum was obtained. Various aspects of this holographic set-up were studied in ref. [10] where it was argued that the magnetic flux corresponds to a non-commutative configuration of adjoint scalars. Furthermore, the fundamental fields are confined to a defect that serves as a domain wall separating the dual gauge theory into regions with different ranks of the gauge group. Extensive study of the phase diagram with a very general ansatz was considered in ref. [12]. Finally, quasinormal modes analysis at vanishing bare mass was performed in ref. [13].
We focus on the holographic renormalization of the set-up, the existence of a second order phase transition and the properties of the theory near criticality. Interestingly we uncover the same critical exponents as in ref. [5]. We also revisit the supersymmetric regime of the theory generalizing it to the case of a non-vanishing bare mass and employing kappa symmetry to confirm the non-broken supersymmetry. The paper is organised as follows: • In Section 2 we describe general properties of the holographic set-up with emphasis on Ramond-Ramond charge conservation and the existence of a domain wall separating regions with different ranks of the gauge group.
• In Section 3 we study the general supersymmetric embeddings. We use kappa symmetry to argue that all of the original supersymmetry of the D3/D5 system is preserved by the flux on the internal S 2 . We construct supersymmetric D5-brane embeddings interpolating between a stack of half D3-brane embeddings parallel to the D3-branes sourcing the geometry and regular D5-brane embeddings. We discuss the field theory interpretation of the flux and how the domain wall is realised.
• Section 4 studies the thermodynamics and holographic renormalization of the theory. We derive expressions for the fundamental condensates associated to the bare mass and the position of the defect. We discuss the phase structure of the theory with a focus on the existence of a critical point of a second order phase transition. We obtain numerically the critical exponents of the theory.
• Finally, we conclude with a brief discussion in Section 5.
2 Holographic Set-Up

Probe configuration
In this section we outline the holographic set-up. 1 We consider an AdS 5 × S 5 black hole supergravity background with metric: Here u 0 is the radius of the black hole horizon related to the Hawking temperature of the background via u 0 = πR 2 T and R is the radius of the internal S 5 given by R 4 = 4πg s N c α 2 . It is convenient to use the following parametrisation of the unit 5-sphere: Next we consider a D5-brane embedding extended along the time direction t, two of the directions of the dual field theory x 1 and x 2 and localised along the third coordinate x 3 2 . The D5-brane also wraps the unit S 2 inside the S 5 part of the geometry parametrised by α and β. We let the D5-brane embedding develop a non-trivial profile θ(u) and x 3 (u). The dual defect field theory was first analysed in ref. [7] in the special case when the D5-branes have vanishing separation at infinity (corresponding to a zero bare mass) and the theory is superconformal. In this paper we are interested in a sort of deformation of that theory caused by turning on a non-trivial profile for the U (1) gauge field on the S 2 ⊂ S 5 . We consider the following ansatz: One can show that the ansatz (2.4) is consistent with the equations of motion of the D5-brane. The full (DBI + WZ) action of the D5-brane embedding is given by: where F (2) = −P[B (2) ] + 2πα F (2) , F (2) is the field strength of the U (1) gauge field living on the D5-brane and B (2) is the Kalb-Ramond B-field. The main contribution 3 to the Wess-Zumino term in (2.5) is given by: We will show that in addition to the leading order contribution from equation (2.2) the background C (4) has an additional contribution from the D3-brane charge density supporting the magnetic monopole implied by the ansatz for the U (1) gauge field (2.4).

Charge Conservation.
Let us consider the gauge transformation: where dΛ (3) has legs along the D5-brane world volume. For the variation of the WZ action (2.6) we obtain: where M 3 is a three dimensional slice of the D5-brane world volume at the point u min where the internal S 2 (parametrised by α and β) shrinks to zero radius 4 . As one can see from equation (2.8) there is a three dimensional boundary in the world volume of the D5-brane. In fact, this is the magnetic monopole (in six dimensions) supporting the ansatz (2.4). It is a standard result in D-brane dynamics that magnetic monopole in the world volume of a Dp-brane is sourced by a D(p-2)-brane ending on the Dp-brane. In our case we need to attach certain number n of D3-branes to the boundary M 3 (with the appropriate orientation) to cancel the variation of the WZ action in (2.8). Indeed, under the gauge transformation (2.7) the WZ action of this additional D3-branes would transform as: It is straightforward to verify that the condition δS W Z + δS D3 W Z = 0 is equivalent to: that is the flux of the U (1) field through the unit S 2 is quantised according to the Aharonov-Bohm effect. Even though we have ensured charge conservation at the boundary M 3 , we now have the problem of attaching the other end of the stack of n D3-branes. We have several options: we could let the new stack of D3-branes carry the charge to infinity, we could attach the new stack to the original stack of D3-branes (most likely unstable) or if at finite temperature we could let the new stack of D3-branes carry the charge to the horizon of the geometry. Note that these are only kinematic considerations and our goal is not to engineer any of the above brane configurations. Instead, in the next sections, we analyse the dynamics of the D5-brane probe and show that the additional stack of D3-branes is dynamically realised as a throat in the D5-brane worldvolume.
At zero temperature there are supersymmetric solutions, generalizing the ones obtained in ref. [10] and interpolating between D5-brane embeddings at large radial distance and "natural" (parallel to the stack sourcing the geometry) D3-brane probes. Note that in addition to solving the charge conservation problem, this introduces a domain wall separating regions with different ranks of the gauge group [10].
At finite temperature the only complete solutions that we uncover are black hole D5-brane embeddings. The Minkowski phase is realised within the black hole one with the D5-brane embeddings developing D3-brane throats carrying the charge to the horizon.
3 The SUSY case. A Domain Wall.
In this subsection we focus on the zero temperature case (u 0 = 0 in equation (2.1)). The massless case was first analysed in ref. [10]. In the following we will consider the general case (allowing a non-zero bare mass of the defect fields) and will supplement these studies with a kappa symmetry derivation of the supersymmetric embeddings. We begin by considering the effective action of the probe brane.

Effective Action
It is convenient to rewrite the metric of AdS 5 × S 5 as: (3.1) where the connection to the coordinates in equation (2.1) is given by: l = u sin θ and ρ = u cos θ. Now considering an ansatz: l(ρ), x 3 (ρ) and A µ = 0 we obtain: Note that x 3 (ρ) is a cyclic variable and can be eliminated by a Legender transformation: where we have used that ∂L/∂x 3 (ρ) = const = 0 5 and Remarkably, the Legender transformed Lagrangian (3.3) is exactly the same as in the supersymmetric H = 0 case and the regular solution to the equation of motion for l is given by l(ρ) = m = const. In fact, in the next section we will use kappa symmetry to show that the solution: is supersymmetric.

Kappa Symmetry
The kappa symmetry matrix is given by: where Γ (0) is given by: Note that the matrices γ a are the pull back for the "flat" gamma matrices given by: where Eμ µ are the tetrads of the metric (3.1) and Γμ are the flat gamma matrices associated to the tetrads. Note also that the matrices γ a are obtained by lifting the index of γ a with the inverse of the induced metric. Next we proceed by considering the ansatz l(ρ) and x 3 (ρ), the goal is to obtain the solution (3.5) from the requirement to preserve kappa symmetry. With this ansatz we obtain: Now using the ansatz for the gauge field (2.4) we obtain: Substituting Γ (0) from equation (3.9) we obtain: Now using that the Killing spinor ε satisfies: 6 we obtain the kappa symmetry condition: Note that if we set l (ρ) = 0, x 3 (ρ) = 0 and H = 0, equation (3.13) reduces to: which reduces the supersymmetry of the background by half. Now let us restore a non trivial x 3 (ρ) and a non vanishing H while keeping the projection (3.14). We obtain: where we used that Γᾱβ ⊗ σ 3 ε = −Γρx 3 ε. It is easy to check that only for x 3 (ρ) given in equation (3.5) we have: and kappa symmetry is preserved without additional projection. Therefore, the solution (3.5) preserves N = 2 supersymmetry and the corresponding dual field theory is given by a generalization of the Lagrangian proposed in ref. [7] to a nonzero bare mass.

Domain Wall
In this section we show that the solutions constructed above describe a domain wall separating the Minkowski space into regions with different gauge group. Furthermore, for massless embeddings the theory remains conformal. Let us integrate the equation for x 3 in (3.5). We obtain: One can see that for large ρ x 3 approaches the constant value x 3 (∞) and the solution asymptotes to the "usual" D5-brane embedding (the one with H = 0) on the other hand in the limit ρ → 0 we have x 3 → −∞, while the internal S 2 shrinks to zero radius. Therefore, at ρ → 0 the solution looks like a D3-brane spanned along the t, x 1 , x 2 and x 3 directions. We can estimate the number of these D3-branes. Let us evaluate the contribution to the DBI and WZ actions of the probe D5-brane of the half Minkowsi space M spanning the t, x 1 , x 2 directions and the ray −∞ < x 3 < −Λ with Λ 1. To leading order we obtain: Comparing equations (3.18) and (3.19) to the contribution to the DBI and WZ actions of a probe D3-brane (positioned at ρ = 0, l = m) of the same half-space: we arrive at equation (2.10) relating the number of effective D3-branes to the flux of the U (1) field. Therefore, the D5-brane embedding interpolates between the solutions analysed in ref. [7] and a stack of k D3-branes parallel to the stack of D3-branes sourcing the geometry.
In figure 1 we have presented plots of the solutions in the x 3 , u plane. The left plot represents a solution with finite mass m. The solution start as a regular D5-brane embedding at large u and bends towards the negative x 3 axes as u → m. As equations (3.18) and (3.19) suggest for u → m and x 3 → −∞ the D5-brane embeddings approximate a stack of k D3-branes. As a results in the dual field theory there is a domain wall (sitting at x 3 (∞)) splitting the Minkowski space into regions with different gauge groups. For x 3 > x 3 (∞) the gauge group is the original SU (N c ), while for x 3 < x 3 (∞) the gauge group is enhanced to SO(N c ) × SO(k) due to the additional D3-brane charge. The right plot in figure 1 represents massless (m = 0) D5-brane embedding. In this case the effective stack of k D3-branes overlaps with the original D3-branes. This enhances the gauge group on the left side (x 3 < x 3 (∞)) of the domain wall to the full SU (N c + k) group. Furthermore, one can check [10] that when m = 0 the induced metric of the D5-brane can be written as: whereρ = √ 1 + H 2 ρ. This is the metric of an AdS 4 × S 2 space time (with different radii of the AdS 4 and S 2 manifolds). This suggests that for m = 0 the classical conformal symmetry of the dual theory is not broken by the domain wall.
We can use the N = 2 supersymmetry to generalize the field theory interpretation of ref. [10] to the non-zero mass case. Namely, that the flux (2.10) carried by the D5-branes corresponds to producing a noncommutative configuration of adjoint scalars in a SU (k) subgroup of the SU (N c ) × SU (k) gauge group on the relevant side of the domain wall. Now using supersymmetry one can match the profile of the D5-branes to the profile of the adjoint scalars X A , one has : where we have used equations (2.10) and (3.17) and for simplicity we have set the position of the domain wall x 3 (∞) to zero. We have also used that m q = m/(2πα ). Clearly for m q = 0 equation (3.23) reduces to the one given in ref. [10] and the gauge group on one side of the domain wall is promoted to SU (N c + k).

Thermodynamics and Holographic Renormalization
Various properties of the termodynamics of the D3/D5 system have been studied in refs. [11], [12] and [14][15][16][17][18][19][20][21][22][23]. In ref. [12] a very general anstaz has been considered and the author has uncovered a number of interesting properties of the system. In this work we uncover the existence of a critical point of a second order phase transition drawing a parallel with the studies of ref. [5].

Holographic Renormalization
We begin by writing the Lagrangian corresponding to the effective action (2.5) for the background and ansatz given in equations (2.1)-(2.4). In fact, we will work with the Euclidean version of the Lagrangian, which in our case differ by an overall minus sign. Furthermore, we define:ũ = u/u 0 ;R 2 = R 2 /u 0 ; (4.1) and consider the dimensionless Lagrangian: (4.2) Note that x 3 is a cyclic variable and hence ∂L E ∂x 3 (u) = const. One can show that the following choice for the constant of integration: is the only one leading to finite Lagrangian density and embeddings consistent with charge conservation. With this choice it is easy to obtain: where to be consistent with the zero temperature limit we have chosen the positive sign for x 3 (u). Using (4.4) it is straightforward to to write down the Legender transform of L E : Solving perturbatively at large u the equation of motion for θ(u) derived from (4.5), one obtains: Using (4.6) it is easy to show that there are no new divergences introduced by the U (1) gauge field. Therefore, we can regularize the Euclidean actions corresponding to L E andL E using the approach of [1]. For D5-brane embeddings reaching the horizon (Black Hole embedding) we can write: The corresponding thermodynamic potentials (or rather their two dimensional densities) are given by: We identify F as the free energy density satisfying: where: S is the entropy density, T is the temperature, C m is the fundamental condensate, m q = u 0m /(2πα ) is the fundamental mass, C ∆x 3 is the condensate of the operator obtained by varying the defect field theory action with respect to its position along x 3 and ∆x 3 ≡ x 3 (∞) − x 3 (u 0 ). Note that since x 3 is defined up to an additive constant we can always define x 3 (u) to satisfy x 3 (u 0 ) = 0 and have ∆x 3 = x 3 (∞).
To calculate C ∆x 3 we vary the free energy density F at fixed T and m, which is equivalent of keeping fixed u 0 and θ(u). It is easy to obtain: Therefore, (4.14) To relate the fundamental condensate to the supergravity parameters it is convenient to consider the "canonical" free energyF given in (4.11). It is easy to show that the following relations hold:F Now fixing T and C ∆x 3 and varying m, that is fixing u 0 and H and varying θ(u) with large u expansion δθ(u) = δm/u + O(1/u 2 ), we get: where we used thatm = 2πα m q /u 0 . Therefore,

Phase structure
In this section we solve numerically the EOM for θ(ũ) obtained by varying the Lagrangian (4.5) and solve for x 3 (ũ) using equation (4.4). The standard classification of the possible D5-brane embeddings [1] is based on the their topology. In particular one distinguishes Black Hole embeddings characterised by reaching the horizon and Minkowski embeddings which close above the horizon by having the internal S 2 wrapped by the brane shrink to zero size. However, as we showed in section 2.1 charge conservation demands that a number of D3-branes is attached to the D5-branes at the point where the S 2 shrinks to support the required magnetic monopole. This observation renders the class of Minkowski embeddings incomplete since no stable configurations with attached D3-branes carrying away the Ramond-Ramond charge were identified. Instead, the Minkowski embeddings develop a narrow throat connecting them to the horizon of the geometry and become a subset of the Black hole embeddings. In figure 2 we have presented a 3D plot of such an embedding, where the horizon is represented by a two-sphere parametrised by θ and β, while the internal S 2 is represented by an S 1 parametrised by β. One can see that the Black hole embedding mimics a Minkowski embedding connected with a D3-brane throat to the horizon. Our next step is to explore the behaviour of the fundamental condensate C m given in equation (4.18) as a function of the bare mass m q at a fixed condensate C ∆x 3 defined in equation (4.14).
In figure 3 we have presented our results for various values of the parameter H. As expected when H = 0 we recover the first order phase transition pattern reported in ref. [1]. Note that the use of different colours for the equation of state curve represents the fact that the phase transition is a topology change transition between Minkowski (blue) and black hole (red) embeddings (corresponding to a confinement/deconfinement phase transition in the dual gauge theory).
Next we consider a small amount of C ∆x 3 condensate (H = 0.001) as one can see the phase transition pattern is extremely close to the one at vanishing C ∆x 3 condensate (H = 0) but the single colour represents the fact that black hole embeddings are the only stable embeddings and the phase transition takes place entirely within the deconfined phase of the theory.
In the next plot in Figure 3 we present the equation of state curve for a larger C ∆x 3 condensate (H = 0.03). As one can see from the plot the first order pattern is much less pronounced, that is: the multivalued region of the C m versus m q curve (or −c versusm) is very narrow. Note that the dashed vertical line in the first three plots represent the critical parameterm cr at which the phase transition takes place.
Finally, when we go to even larger values of the C ∆x 3 condensate (H = 0.06), the first order phase transition disappears as evident in figure 3. In figure 4 we present the corresponding plots of the quantityĨ E defined in equation (4.8) and related to the free energyF through equation (4.16). One can see that at H = 0.06 the first order phase transition is replaced by a crossover.
The observed behaviour suggests that at some H cr the first order phase transition ends on a critical point of a second order phase transition. In the next section we study numerically this critical point with a focus on the critical exponents of the fundamental condensate.

Critical Point
As the analysis of the previous subsection suggests, for some H cr (in the range 0.03 < H cr < 0.06) the first order phase transition ends on a critical point of a first order phase transition. Beyond this critical point we have a crossover and the the condensate versus bare mass curve is single valued. It is straightforward to determine numerically the value of H cr at which the multivaluedness of the condensate curve disappears, we obtained H cr ≈ 0.0443835.
In figure 5 we have presented our results for the condensate C m and the free energyF as a function of the bare mass. The dashed vertical curve represents the critical parameterm cr at which the second order phase transition takes place. If we denote byc cr the critical valuec cr =c(m cr ), and zoom in near the phase transition,   where ∂ 2F ∂m 2 is proportional is the mass susceptibility of the condensate C m and ∆ is the corresponding critical exponent. We are going to show numerically that ∆ = −2/3 with a very high confidence. In Figure 6 (a) and (b) we have presented a zoom in of the condensate and susceptibility near the the critical region. The dashed vertical line represents the critical parameterm cr = 1.19591, while the continuous curves represent fits with ∆ = −2/3 in equations (4.19) and (4.20). As one can see the fits are excellent. Next rather than feeding in the value of ∆ we extract its value using by taking a logarithm on both sides in equation (4.19) and using a linear regression. In Figure 6 (c) and (d) we have plotted log |m −m cr | versus log |c −c cr | for both m <m cr andm >m cr . By taking the mean of the two slopes we can estimate ∆ and its standard deviation. Our results is: which is extremely close to ∆ = −2/3. Interestingly, this is the same critical exponent as the one reported in ref. [5] 7 , where the same holographic technique (introducing a non-zero U (1) flux through a two sphere) was considered to deform the first order phase transition into a second order one.

Discussion
In this paper we consider the holographic gauge theory dual to the D3/D5-brane intersection with a U (1) flux on the transverse two-sphere wrapped by the D5-brane probes. We consider both the finite temperature and zero temperature phases of the theory.
In the zero temperature case Ramond-Ramond charge conservation causes the D5-brane embeddings to bend along the Neumann-Dirichlet direction (x 3 ) carrying the charge to infinity. Remarkably, these solutions preserve the N = 2 supersymmetry of the D3/D5 intersection and interpolate between D5-brane embeddings at large radial distance and D3-brane probes (parallel to the stack of D3-branes sourcing the geometry). This not only solves the charge conservation problem but splits the 1+3 dimensional Super Yang-Mills theory into regions with different ranks of the gauge group. That is the 1+2 dimensional defect (introduced by the D5-branes) behaves as a domain wall. In the case of zero bare mass the theories on both sides of the domain wall are [10] SU (N c + k) and SU (N c ), while in the case of finite bare mass we showed that the groups are SU (N c ) × SU (k) and SU (N c ). Given that the field content of the N = 2 theory is completely fixed by supersymmetry we extended the field theory interpretation of ref.
[10] to a non-zero bare mass. At finite temperature, a D3-brane probe is no longer "natural" in the sense that D3-brane probes parallel to the field theory directions are no longer stable solutions. As a result a solution interpolating between a D5-brane and a (parallel) D3-brane embedding is no longer possible. To avoid violation of Ramond-Ramond charge conservation the Mikowski embeddings develop a D3-brane throat connecting them to the black hole horizon, while still developing a profile in the Neumann-Dirichlet direction x 3 . Note that even though Black hole embeddings do not extend all the way to x 3 = −∞ the continuity of the horizon (the horizon is extended along x 3 ) requires the Ramond-Ramond charge to reach infinity [10] and the domain wall interpretation is valid at finite temperature too. Considering equation (3.23) it is natural to speculate that at finite bare mass the large gauge group is SU (N c )×SU (k) and only at vanishing bare mass the group is promoted to SU (N c + K).
Even though Mikowski embeddings violate charge conservation, they can still be approximately realised within the class of black hole embeddings and the first order phase transition between the two classes of embeddings is no-longer a topology change transition. 8 The latter opens up the possibility to deform the first order phase transition to a second order one. Furthermore, the non-trivial profile along x 3 is reflected in a non-zero condensate C ∆x 3 thermodynamically conjugated to the position of the defect and proportional to the U (1) flux on the transverse two-sphere. Increasing the flux (the C ∆x 3 condensate) lowers the latent heat of the first order phase transition and for sufficiently strong flux it ends on a critical point of a second order phase transition. Our analysis of the critical point showed that the mass susceptibility diverges at criticality with a critical exponent of ∆ = −2/3. Interestingly, this is the same critical exponent as the one reported in ref. [5], where the holographic set-up also includes a constant U (1) flux through a shrinking S 2 cycle.
Given that the dual field theory is 2 + 1 dimensional the critical point that we studied might be relevant to a condensed matter system exhibiting critical behaviour with the same critical exponents. One might be able to use the holographic set-up as a toy model to study the near critical regime of such condensed matter system.
Our work can be extended in several possible directions. The most straightforward extension is to supplement our thermodynamical analysis with a study of the quantum fluctuations/meson spectra of the theory. Another extension is considering the holographic set-up in the presence of an external electric field or R-charge chemical potential, which will enable us to study quantum phase transitions. Finally, it would be interesting to understand better the sumersymmetric phase of the theory and the observed domain wall.