Minimal Submanifolds asymptotic to AdS4xS2 in AdS5xS5

In this work, as an attempt in understanding the interplay between deformations of the external and internal spaces associated with the probe-brane embedding submanifold, we construct the zero-temperature phase diagram for the coupled phase between two two-dimensional defects stacked parallel in a four-dimensional ambient spacetime. Different UV parameters are turned on for different defects. We study the system in the quenched strong coupling limit, using holography with probe brane approximation, realized explicitly through the D3/D5 system. This coupled phase is holographic dual to the presence of spontaneous symmetry breaking of the individual ultraviolet flavor symmetries associated with the double heterostructure of the defect layers. We characterize this solution by its infrared geometric data, and present the numerical result showing a first-order phase transition between the asymmetrically coupled phase and the more mundane decoupled phase.

submanifold for both the external and internal embedding is possible. In fact, by stacking two such defects parallel, we show that, even when the two defect contents are in general different in their UV parameters, hence holographic dual to double heterostructure, there still exists an asymmetrically connected submanifold solution for the probe brane system. This configuration is dual to a coupled phase in the field theory side. We also find that such an asymmetric configuration is thermodynamically favorable upon some mild detuning of UV data, and there is a first-order transition from this asymmetrically coupled phases to the more mundane decoupled phase, with the increase of the difference of their UV parameters.
The organization of the paper is as follows: In Section 2 we first review the notion of the D7 and D5 embeddings in AdS 5 × S 5 , which also provides the mathematical setting. We also introduce the notion of single-sheeted, also later called decoupled phase, for the rest of the work. After that we present the calculation subtleties associated with the construction of the asymmetric joint-sheeted, also called nontrivial coupled phase. In Section 3 we present the numerical results, showing a first-order transition between these two aforementioned phases.

Methods
In this section, we introduce the notion used by Graham and Karch [15] for setting up our convention. They consider submanifolds asymptotic to AdS k+1 ×S l inside AdS n+1 ×S m . Specialized to the D3/D5 embedding, k = 3, l = 2, n = 4, m = 5, this probe brane system is holographic dual to two localized defects [5,7], with quark content of massive N = 2 N f -favor hypermultiplets in the fundamental representation, inserted into the background of N = 4 U (N c ) super-Yang-Mills gauge field theory, in the N c ≫ 1, large 't Hooft coupling limit. The two defects are separated by ∆x in the common transverse direction, and each sector is characterized by its own quark mass m L or m R in the UV Lagrangian. We will construct the phase diagram by tuning the dimensionless ratios, m L ∆x and m R ∆x.

Submanifold Extension: mathematical settings
The mathematical setting of this article is as the following: Considering the background spacetime being a product manifold X × K, where X is a n + 1 dimensional asymptotically hyperbolic manifold with n dimensional boundary thereof, M ≡ ∂X, and K is a m dimensional compact manifold. The metric of X × K is parametrized using the standard Fefferman-Graham form [16]: with boundary M located at r min = 0, hereafter referred as the UV end. The minimal submanifold is denoted as Z ⊂ X × K, with boundary ∂Z = N × S ⊂ M × K. We can parametrize M as (x α , u α ′ ), K as (s A , t A ′ ), such that: 1) The boundary of the submanifold N × S is given by u α ′ → 0, t A ′ → 0, and 2) t and u variables are "orthogonal" to the boundary, as in the following sense: In this article, we focus on submanifolds asymptotic to AdS k+1 × S l , embedding into AdS n+1 × S m . We will choose to embed S l into S m as follows: with S l sitting at the equator parametrized by the θ embedding function. This ansatz corresponds to the simplest symmetry-breaking pattern associated with the dual defect field theory.
Given the above parametrization, let us first consider the single-sheeted (later also called decoupled) phase of the minimal submanifold. In this phase, the minimal submanifold is not only asymptotic to AdS k+1 × S l , but when it extends into the bulk, the only change is the internal sphere radius, controlled by the θ embedding function, being endowed with nontrivial radial dependence, θ = θ(r). All externalspace embeddings are constant, hence no "bending" of the probe external location occurs along the radial direction. In fact, notice that radial dependence is actually the most generally allowed dependence if we insist on preserving the defect-translational invariance, as well as S l and S m−l−1 isometries for the dual field theory. With the standard Poincaré patch for AdS n+1 and static gauge for AdS l+1 , the area of the submanifold is then given by: where α is a generalization introduced to account for the possible difference between the curvature radii of internal compact space S m and the external AdS k+1 space, α ≡ (R S m /R AdS k+1 ). α = 1 in AdS 5 × S 5 , the standard near-throat limit of supergravity background generated by from the D3 brane. Given that the D3/D5 system is our primary concern, hereafter we commit ourselves to the α = 1 case. Also we have set R AdS k+1 to 1. The usual variational method yields the equation of motion for θ(r) as: Let's first consider the special case of the D3/D7 system, where the D7 probe brane is asymptotic to AdS 5 ×S 3 (k = 4, l = 3): This submanifold is completely filling the external space of the background spacetime, and we can easily check that we then have an one-parameter family of embeddings, θ = arcsin(mr), with m being the free parameter. θ is the so-called slipping mode in the literature. The same solution also holds for the D3/D5 system, with the probe brane embedding submanifold asymptotic to AdS 4 × S 2 (k = 3, l = 2). Such a simple analytical result is connected with the supersymmetry of the probe brane system [17]. In more general systems, however, one can only observe that θ being zero is still a solution, which extends into the one-parameter family of solutions by turning on the slipping mode with the initial slope, θ(r) = mr + O(r 2 ).
But one will need to numerically integrate out the equation to complete the profile.
However, we are naturally more interested in configurations with the external embedding being nontrivial, specifically the joint-sheeted configuration, corresponding to two defect stacked parallel to each other, with the submanifolds bending towards each other and eventually smoothly joined in the radial direction. An example of this connected configuration is the Wilson line used by Maldacena [18] to compute the quark-antiquark force: Near the boundary the submanifold is asymptotic to AdS 2 times a point.
However, deep inside the bulk, the submanifold is connected into a U-shaped configuration. Therefore, in this article, we parametrize such a connected submanifold as being the union of two different branches, , with each branch chosen to be N i × S i = R k × S l . As mentioned before, to preserve the dual field theory symmetries associated with the symmetries of both the defect space and the transverse space, only radial dependence is admissible for the embedding functions. For simplicity, instead of the simple internal global symmetry breaking pattern chosen by Eq.(2.3), we will also restrict ourselves by only turning on one nontrivial external embedding function, hereafter named x(r). Still working in the standard Poincaré patch for AdS n+1 and static gauge for AdS l+1 , the area of the submanifold is modified as: with the equations of motion given by: However, unlike the previous single-sheeted phase, these equations admit no analytical solution at least to us, and hence we are forced to resort to numerical methods for computing the profile.  coupled configuration is expressed by two separate but identical branches joining smoothly at the turning point r max , normalized to 1 in above. In Fig.1b, the dotted line stands for the asymptotically decoupled configuration with the same UV parameter, the asymptotic mass term m asymptotics ≡ θ ′ (r) r→0 . Notice we use the translational invariance to set x(r max ) = 0.

The Regularity Constraint and the Cascaded Integration Scheme
To solve the embedding equations Eq.(2.7), the boundary conditions, required to generate the joint-sheeted configurations, are given by requiring the connection of the two branches inside the bulk space being smooth.
This however renders x ′ (r) r→rmax divergent and further denies us control over the boundary specification at r max , the turning point of the joint-sheet submanifold. Nevertheless, this seeming difficulty can be resolved by a closer inspection of the regularity requirement: Notice that the Lagrangian density is cyclic in the  this coupled configuration is expressed by two separate and different branches joining smoothly at turning point r max , normalized to 1 in above. In Fig.2b, the dotted lines stand for the corresponding asymptotically decoupled configuration for each branch, with different mass terms on the different side of the lobe. We can generate such graph within the x-parametrization from the IR end, and cascade with the r-parametrization when approaching the UV end. Notice we use the translational invarience to set x(r max ) = 0.
x-parameter, hence conserving the conjugated momentum thereof, Therefore, we can express x ′ in term of this integral of motion Π x , (2.9) The regularity requirement therefore translates into the following two scenarios: First, we have the denominator equal to zero, providing a condition relating all the parameters at the IR end: Second, we have the numerator equal to infinity, signaling at the IR end sitting not only a divergent x ′ (r max ) but also a divergent θ ′ (r max ) as well: One may expect the first scenario to be more relevant, given the smoother behavior of θ(r) r→rmax . But further analysis, by expanding the submanifold in power series given the regular behavior of θ(r) around r max , reveals that the first scenario is actually too restrictive. It constrains the two branches of the submanifold to be exactly the same, ending with only symmetric joint-sheet configurations. Therefore, to include the most general, asymmetric joint-sheeted configuration, the second scenario and hence singular behavior of θ(r) around r max will need to be considered.
To numerically generate the solution given by the second scenario for the boundary conditions, we adopt the following cascaded integration scheme: First, around r ∼ r max , we choose to parametrize the solution in terms of the x-variable, with the equations of motion given by: Choosing the turning point to be at x 0 = 0, the solution is uniquely determined by the "IR" data: terms of x-variable at the boundary, to approach the boundary with more numerical stability, we will then switch over to the r-parametrization before carrying out the integration to the boundary. With such a cascaded scheme, we can numerically find the most general, asymmetric connected configuration.

Extraction Scheme
The great virtue of the probe brane approximation is that various physical quantities of interest can be obtained without re-solving the Einstein equation with the probe source, as long as those quantities can be calculated using the free energy and thermodynamical relations, as already explained in [19]. Therefore, to determine the phase diagram, we can follow this tenet and simply compute the free energy using the onshell action for the given configuration, and determine which phase, single-sheeted (also called decoupled) or joint-sheeted (also called coupled), is thermodynamically favorable. However, the boundary divergence of AdS space renders the first step rather laborious: the canonical approach is to deploy the holographic renormalization [20], which carefully reconstructs the diffeomorphism-invariant counterterms by examining the divergent structure associated with the tentative cut-off plane, dual of the UV regulator for field theories.
However, in this study, given that we eventually only focus on the difference of free energies between the decoupled and coupled configurations, we will adopt the following extraction scheme (also known as background subtraction): given any coupled configuration, we first construct the dual disconnected configuration with the same UV parameter m asymptotic . Then the difference of these two configurations is, by construction, vanishing near the boundary, since all coupled configurations are asymptotic to the decoupled solution in the UV region (See Fig.1-2).
The rationale behind our extraction scheme is that the to-be-constructed holographic counterterms can only depend on the UV behaviors but not IR physics, which is already present in the decoupled configuration once the only relevant information, m L or m R is extracted. This scheme is also more in tune with the pre-Wilsonian renormalization philosophy, "sweeping under the rug", that no divergence should be present if every physical prediction is expressed with physical quantities.
In practice, given we can only work with the numerically generated configuration, there is a potential caveat associated with such an extraction scheme: the UV parameter m L/R is located at the singular point of the equations of motion, and we only extract this information numerically up to a small cut-off distance due to the inevitable numerical instability. Such a seemingly simple numerical recipe can be subject to more elaborate modification: One can construct the analytic solution expanded around the singular point, and extract the parameters by fitting at a point away from boundary, with now more numerical control.
However, given the phase transition we find is of first order in nature, and in the precision we are working on, the difference thus introduced is found to be numerically negligible with no qualitative change of the conclusion. configuration; Black(thin) line, the difference between the previous two, with the signed area being the free energy difference. One can observe the divergence term is indeed subtracted out when approaching the boundary, and the difference is only due to the back reaction inside the bulk. Notice in this specific case the decoupled configuration can extend deeper into the bulk even after the coupled counterpart already terminates (r max = 1 by normalization). The reverse situation also is possible, where the decoupled surface can terminate before the connected one.

Symmetric Coupled Phase
We first restrict to the symmetric connected configurations only. Fig.4 shows the on-shell area difference, which is also the free energy difference, between the symmetric coupled and decoupled configurations as a function of the normalized mass term ξ normalized : For every value of ξ normalized we find two coupled solutions. At very low ξ normalized , one of the symmetric coupled configuration is the thermodynamically favored phase, compared with its decoupled counterpart(with negative free energy difference); as ξ normalized increases, the area difference begins to shrink, and finally a first order phase transition is reached when ξ normalized reach the value of 0.165, after which the decoupled dual becomes more favorable ones; however, as ξ normalized keeps increasing and eventually above 0.315, the coupled configuration disappears entirely, and only the decoupled phase exists as the only allowed solution for the minimal submanifold.

Asymmetric Coupled Phase
Using the cascade evolution scheme, we can also compute the phase diagram for the asymmetric connected configuration. Fig. 5 shows the phase diagram of the first order phase transition for the minimal submanifold between the asymmetric coupled and decoupled configurations as a function of two physical parameters: with where we use L and R to denote the different branch under study. Notice that due to the reflection symmetry of the parameter space along the axes of both ξ L = ξ R and ξ L = −ξ R , we only need to scan the portion of Fig. 5. In Fig.5b, the partial scanning of the free energy density difference is also shown, which explicitly demonstrates that the phase transition is also of the first-order.

Conclusion and Future Research
In this article we investigate the submanifolds asymptotic to AdS 4 ×S 2 embedded into AdS 5 ×S 5 . We restrict to the largest unbroken residual symmetry. We find the asymmetric joint-sheeted configuration, with both external embedding and internal embedding radially deformed. These are dual to the coupled phase between two two-dimensional defects stacked parallel, with different UV parameters. While similar phases have been studied previously in the probe brane systems, what is novel about our work is that the asymmetric probe brane system is holographic dual to the double heterostructure in condensed matter systems, with different UV data in the two layers. The asymmetrically coupled phase is holographic dual to the spontaneous symmetry breaking of the original decoupled ultraviolet U (N f ) L × U (N f ) R flavor symmetries, associated with flavor rotation for individual defect content, into the diagonal subgroup U (N f ) due to the non-vanishing condensate developed at infrared, ψ i a (e A ) a b ψ j,b . This order parameter has also been studied in [21] for the chiral condensate of Saki-Sugimoto model. Solutions are shown completely parametrized by the infrared geometric data, being the IR value and slope of the internal deformation, upon which the entire profile is numerically constructed by the cascading integration scheme we adopted in this paper. Aimed with the complete numerical solution, we map out the phase diagram between this asymmetrically coupled phase and the competing mundane decoupled phase, and we find that the coupled phase is dominating at zero temperature, even when the UV parameters are mildly detuned, and with large enough difference the system undergoes a first-order phase transition to the decoupled phase.
Notice that more general deformations of the submanifolds are possible, corresponding to finer symmetry breaking pattern in the dual field theory side. Such details present no challenge for profile construction using the cascading integration scheme. However the corresponding phase diagram is computationally difficult to enumerate: For example, in the cases of deforming two internal embedding functions, ie. θ = θ(r) and ψ = ψ(r), we find that this solution can be constructed with 4 degrees of freedom, being their IR values and slopes of the two internal coordinates. The phase diagram will be 4 dimensional, which present itself an expensive numerical barrier. Tracing out this space to locate the transition boundary and hence its transition nature will be left for future investigation.