Singularity development and supersymmetry in holography

We study the effects of supersymmetry on singularity development scenario in holography presented in \cite{Bosch:2017ccw} (BBL). We argue that the singularity persists in a supersymmetric extension of the BBL model. The challenge remains to find a string theory embedding of the singularity mechanism.


Introduction
Horizons are ubiquitous in holographic gauge theory/string theory correspondence [2,3]. Static horizons are dual to thermal states of the boundary gauge theory [4], while their long-wavelength near-equilibrium dynamics encode the effective boundary hydrodynamics of the theory [5]. Typically, dissipative effects in the hydrodynamics (due to shear and bulk viscosities) lead to an equilibration of a gauge theory state -a slightly perturbed horizon in a dual gravitational description settles to an equilibrium configuration. If the initial perturbation away from thermal equilibrium is sufficiently strong, non-hydrodynamic modes participate in the equilibration process 1 1 Even when not excited, non-hydrodynamic modes are important as they determine the convergence properties of the all-order effective hydrodynamic description [6,7].
-perturbed bulk horizon relaxes via quasinormal modes [8][9][10][11][12]. We call horizons that relax via hydrodynamic modes that attenuate in space-time domain or via positively gapped quasinormal modes 2 stable horizons. Bulk holographic dual with (dynamically) stable horizons describe boundary gauge theory states which are stable with respect to sufficiently small fluctuations.
Translationary invariant horizon can suffer an instability when a hydrodynamic mode in a system becomes unstable. An example of such an instability is discussed in [13]. In this case, the expected end-point of the evolution is a new inhomogeneous phase of the system 3 . Alternatively, a horizon can be de-stabilized when a positively gapped non-hydrodynamic mode becomes unstable when one lowers energy (or temperature) [15,16]. This instability realizes a holographic dual of a spontaneous symmetry breaking in a mean-field approximation -the system undergoes a second-order phase transition towards a new stable phase with a finite density condensate of the originally unstable mode.
Above classification of (un)stable horizons in a holographic framework matches both the gravitational and the field theory intuition. However, there is more to the story.
In [17,18] a new instability of the translationary invariant holographic horizons was identified: while there is a linearized instability below some critical energy density (or temperature), triggered by a non-hydrodynamic mode spontaneously breaking a discrete symmetry, there is no candidate equilibrium state with a condensate of this unstable mode.
The dynamics of the model [17,18] was extensively studied in [1] (BBL): both the presence of the linearized instability and the absence of the suitable equilibrium state for the evolution below the criticality was confirmed dynamically. Moreover, it was argued that the gravitational system evolves to a region of arbitrary large curvatures in the vicinity of the horizon, asymptotically turning such region singular in a finite time with respect to the boundary theory. The area density of the apparent horizon, associated with the non-equilibrium entropy density of the boundary gauge theory [19,20] also diverges within a finite boundary time. This latter observation is significant as it precludes evolution of instabilities towards any finite entropy density spatially inhomogeneous equilibrium states as well -in other words, despite the fact that BBL dynamics occurs in a constrained phase space (spatial homogeneity and isotropy), the conclusion that the system evolves to a singularity, violating the weak cosmic censorship conjecture, is robust.
The BBL model is not a top-down holographic construction. Thus, one might wonder whether the phenomenon discovered in [1,17,18] occurs in real string theory holographic examples. A particular aspect of the model is the unboundedness of the bulk scalar field potential. It was argued in [1] that there is no weak cosmic censorship conjecture violation once the potential is bounded. On the other hand, unbounded scalar potentials do occur in supersymmetric top-down holographic models (as e.g., in [21]). Additionally, the exotic critical phenomena of [17] was identified in supergravity model [21] (DG) in grand canonical ensemble. In this paper we partly address above questions.
Unfortunately, construction of a top-down holographic model realizing the dynamics of BBL model remains open: while the DG model is "exotic" in grand canonical ensemble, we study in section 2 the equilibrium properties and dynamics of DG model in microcanonical ensemble and show that it realizes a standard spontaneous symmetry breaking instability as in [15,16]. In section 3 we present a supersymmetric extension of the BBL model, which exhibits the exotic phenomenon of [17,18]. This sBBL model does not have an equilibrium state below the criticality, and, as in BBL model, its homogeneous and isotropic states evolve towards asymptotically divergent expectation value of the symmetry breaking operator. The similarities and differences of the BBL and sBBL model are further highlighted in section 4. Some technical details are delegated to appendix A.

DG model in microcanonical ensemble
In [21] the authors studied equilibrium states of d = 3 N = 8 superconformal gauge theory dual to AdS 4 × S 7 at a finite temperature and a finite chemical potential with respect to a diagonal U(1) R ⊂ SO(8) global symmetry within different consistent truncations of D = 11 supergravity on S 7 , allowing for different patterns of spontaneous global symmetry breaking. We consider the following two consistent truncations, with the effective actions where • "DG-A" model Lagrangian (2 + 2 equal charged scalars model in [21]) is : • "DG-B" model Lagrangian (4 equal charged scalars model in [21]) is: For the maximally supersymmetric quantization the dimensions ∆ of the CFT opera- and to the bulk scalar φ (DG-B model) is In both models A 1 is the bulk gauge field dual to U(1) R global symmetry; the bulk scalars γ i and φ have the U(1) R charge 1, while σ is R-symmetry neutral.
Note that DG-A model is invariant under the Z 2 symmetry: Its truncation to a Z 2 -even sector produces DG-B model with the identification

Grand canonical ensemble (a review)
We focus on DG-A model; the construction for the DG-B model is a consistent truncation as explained in (2.9). Following [21], we set the chemical potential µ 1 = 1 for U(1) R symmetry, and for the other global U(1) (holographic dual to the bulk gauge fieldĀ 0 )μ 0 = 0, i.e., spatially homogeneous and isotropic thermal equilibrium states of the CFT are represented by a bulk geometry with the asymptotics: where r is a standard asymptotic-AdS radial coordinate (see (2.13) below) and {q 1 ,q 0 } determine the U(1) R × U(1) charge densities, see (2.16). Imposing the supersymmetric quantization (2.6), we identify the expectation values of the dual operators The background metric ansatz takes form with the following asymptotic expansions at infinity: where parameter m is related to the energy density of the state, see (2.16). Assuming that a regular Schwarzschild horizon is located at r = r h , we have asymptotic expansions in (r − r h ) > 0 as (2.15) Parameters in asymptotic expansions (2.10), (2.11), (2.14) and (2.15) determine the thermodynamic properties of an equilibrium CFT state, the energy density E, the pressure P = 1 2 E, the grand potential density Ω = −P , the temperature T , the entropy density s, and the U R (1) × U(1) charge density {ρ 1 ,ρ 0 }: (2.16) We use numerical shooting method developed in [22] with a radial coordinate to relate the boundary and the near-horizon asymptotics of the dual gravitational background. We reproduce the results reported in [21]: At any temperature T > 0 there is a phase of the CFT with zero condensates 4 (2.12) and a grand potential density Ω RN : There is a critical temperature T c in the system separating phases with nonzero condensates (2.12) i.e., model DG-B, exists at T ≥ T c . As its grand potential density Ω DG−B is DG-B model realizes the exotic thermodynamics discovered in [17] in grand canonical ensemble.
A phase with  i.e., model DG-A, exists at T ≤ T c . As its grand potential density Ω DG−A is DG-A model realized a standard holographic dual of a spontaneous symmetry breaking in grand canonical ensemble [15,16]. O γ 1 condensate in various phases of the CFT and the corresponding grand potential densities close to criticality, see (2.19), are presented in fig. 1.

Microcanonical ensemble
In previous section we reproduced some of the results of [21] to confirm that model DG-B indeed exhibits exotic thermodynamics in the spirit of [17] in grand canonical ensemble. The singularity mechanism identified in [1] hinges upon the persistency of exotic phase structure in microcanonical ensemble, i.e., in dynamical evolution with fixed energy density and charges. Unfortunately, the instability mechanism in both DG-A and DG-B models in microcanonical ensemble lead to standard picture of the spontaneous symmetry breaking 5 [15,16].
To study DG-A/B models in microcanonical ensemble, the only modification needed 5 We verified that this remains true for supersymmetry breaking assignment of dimensions of bulk operators dual to {γ i , σ}. is the change of the boundary conditions on gauge fields in (2.10); now we require where the "chemical potential" parameters {µ 1 ,μ 0 } are adjusted to keep fixed U(1) R charge density (parameter q 1 ), as well as the vanishing charge for the remaining U(1) (in model DG-A). For numerical analysis we set 6 ρ 1 = −q 1 = 1 2 -such a choice will put the phase transition numerically close to the one in grand canonical ensemble (2.19) with µ 1 = 1.
Phase diagram in microcanonical ensemble is presented in fig. 2. The dashed blue curve represents symmetry unbroken phase -the RN black brane. The latter curve ends at the extremal RN solution: in the microcanonical description of the CFT below the critical energy density.

Dynamics of DG-B model
In the previous section we argued that DG-B model describes the dominant phase of the CFT below the critical energy density in microcanonical ensemble. Here we describe dynamics of spatially homogeneous and isotropic states of the model. We confirm that the criticality in the model is a standard holographic realization of the spontaneous symmetry breaking.
We begin with undoing the implicit bulk gauge symmetry fixing in the effective action of DG-B model (2.4): We further introduce a pair of scalar fields f i in lieu of {φ, θ}: Assuming translational invariance along the spatial directions, we take leading for the following equations of motion together with the constraint equations: The constraint equations are preserved by the evolution equations provided they are satisfied at a given timelike surface -which in our case is the AdS boundary.
The general asymptotic boundary (r → ∞) solution of the equations of motion, given by 7 is characterized by seven (generically time-dependent) parameters The last two constraints in (2.31) imply that only five of them are independent: provided f i,2 relate to f i,1 as follows: Additionally, µ 1 in (2.24) is identified with Finally, λ(t) is the residual radial coordinate diffeomorphisms parameter which can adjusted to keep the apparent horizon at a fixed location, which in our case will be r = 1: To initialize evolution at t = 0, we provide the bulk scalar profiles, Further details of the numerical implementation can be found in Appendix A.1.

QNMs and linearized dynamics of DG-B model
We discuss here the spectrum of QNMs associated with the spontaneous U(1) R symmetry breaking and the linearized dynamics in DG-B model.
The spectrum of the symmetry breaking QNM is presented in fig. 3. Note that ω I > 0, signalling the instability, once E < E crit , see (2.26).
As in [1], we verify our dynamical code for DG-B by "turning off" the scalar fields present the evolution of f 1,1 (t) and  fig. 6 presents, we find that this is indeed the case. Notice that approach to equilibrium is rather slow (t equilibration T 20 in this case ) -as emphasized in [12] this is expected in the sector of the gauge theory responsible for the critical behaviour close to transition.
We confirm the conclusion reached in section 2.2: even though DG-B model exhibits an exotic thermodynamics discovered in [17] in grand canonical ensemble, it represents the well-known holographic realization of the spontaneous symmetry breaking [15,16] in microcanonical ensemble.

Supersymmetric extension of BBL model
We would like to generalize BBL construction [1,17,18] to a "supersymmetric" model.
As we show, this is rather easy to achieve.
In four-dimensional gauged supergravity [24,25] the Lagrangian of the scalar fields coupled to gravity is restricted to where the kinetic term is and the potential is determined from the (real) superpotential W Assuming the metric ansatz, supersymmetric RG flow equations are obtained from the supersymmetry variation of a the fermions: • the first order RG flow equations are consistent with the second order EOMs derived from the Lagrangian; • because flow equations for scalars are of the first order, supersymmetry imposes specific quantization for the scaling dimensions for the dual operators (e.g., see (2.6) for DG-A model).
Recall, for the BBL model we have: Note that the potential is unbounded from below for the nonlinear coupling g < 0.
The latter is necessary for the exotic phase structure of [17,18]. For a supersymmetric generalization of the model, we introduce the superpotential W, leading to P sBBL = P BBL + O(φ n χ 6−n ) , n ≥ 0 , (3.8) i.e., the sBBL scalar potential captures the leading nonlinearity of the BBL model, but differs for higher order nonlinear interactions. As in BBL model, there is a single nonlinear coupling constant g. Explicitly, (3.9) Notice that additional terms in P sBBL (compare to P BBL ) are trivially unbounded from below when − 21 + 27 √ 77 104 < g < 0 . Notice that the quantization (3.11) is different from the one used in [1]; as a result, numerical implementation of sBBL model (see appendix A.2) is closer to that of DG-B model rather than the original BBL model.

Phase diagram and QNMs
Since BBL and sBBL models differ by higher-order nonlinearities in bulk scalar potentials, holographic renormalization for the models is the same. In particular, we can borrow the expressions for the thermodynamic quantities in [17], appropriate for the supersymmetric quantization (3.11). We present the results only -for detailed discussion follow [1,17,18]. • The entropy density of the symmetric phase s sym as a function of the energy density E is presented in figure 7 (left panel). While this phase is thermodynamically  • The spectrum of Z 2 symmetry breaking quasinormal modes in sBBL model is presented in fig. 8 as a function of energy density for g = −2 (red curve) and 8 √ −h is the induced metric on the cut-off surface ∂M 4 in the holographic renormalization. 9 It is easy to verify that for the supersymmetric RG flows (3.5) in sBBL model the (vacuum) energy density vanishes, as required by the boundary CFT supersymmetry.
The unusual feature of these modes, first observed in [18], is the fact that they have vanishing real part, i.e.,

Re(ω)
red & green = 0 . (3.14) Because of (3.14), these modes must disappear from the spectrum in the limit E Λ 3 → ∞, i.e., when sBBL approaches its UV (conformal) fixed point 10 . The blue curves in fig. 8 represent the sBBL QNM mode, which connects at asymptotically large energies to ∆ = 1 QNM of the CF T 3 . Notice that as E increases sufficiently far over the corresponding E crit , there is a "level crossing" (left panel), and the red and green curve QNMs cease to dominate the relaxation of the symmetry breaking fluctuations in the system -it is governed by the blue curve QNM. represents the symmetric phase of the system; the green and red curves are the exotic phases at g = −8 and g = −2 correspondingly, which exist only for E > E crit .
• As in BBL model, in sBBL model new phases with spontaneous Z 2 symmetry breaking appear in the microcanonical ensemble once E > E crit , fig. 9. Once again, these new phases have lower entropy density than the symmetric phase (blue curve) and thus never dominate the dynamics of the system.

Dynamics of sBBL horizons
In section 3.1 we established that the equilibrium physics of BBL and sBBL models in microcanonical ensemble is identical: in both models the Z 2 symmetric equilibrium phase becomes unstable below some critical energy density E crit ; there is no end-point for the instability and the 'hairy' phases bifurcate from the onset of the instability towards E > E crit . These new phases are exotic -they have lower entropy density than the symmetry preserving phase and thus never dominate dynamically.
We now study the dynamics of the sBBL model. We highlight the main results and refer the reader to [1] for further implementation details. Because the quantization of the operator O r dual to the bulk scalar φ is modified, see (3.11), the general asymptotic boundary (r → ∞) solution of the equations of motion is now given by It is characterized by two constants {p 2 , µ}, and three dynamical variables {p 1 (t), q 4 (t), λ(t)}.
These parameters have the following interpretation: • p 2 and p 1 (t) are identified with the deformation mass scale Λ and the expectation value of the relevant operator O r of the dual QF T 3 , • q 4 (t) is the normalizable coefficient of the bulk scalar χ, identified with the expectation value of the Z 2 -symmetry breaking irrelevant operator O i of the dual • µ is related to the conserved energy density E of the boundary QF T 3 as follows • λ(t) is the residual radial coordinate diffeomorphisms parameter which can adjusted to keep the apparent horizon at a fixed location, which in our case will be r = 1: To initialize evolution at t = 0, we provide the bulk scalar profiles, Some details for the modification of the numerical code of [1] are collected in appendix A.2.

Dynamics of symmetric sBBL sector and its linearized symmetry breaking fluctuations
Typical evolution in the symmetric sector (the bulk scalar χ identically vanishes) is shown in fig. 10. An important check on the consistency of the evolution is the fact that the expectation value of O r evolves at asymptotically late times to its equilibrium value at the corresponding energy density (here we take g = −8 and E = 0.42978(1)E crit ), computed independently in section 3.1.
One of the advantages of the holographic formulation of the dynamics of strongly interactive gauge theories is a natural definition of the non-equilibrium entropy density s(t), associated with the area density of the apparent horizon, The evolution of the entropy density with time for the same set of parameters as in fig. 10 is shown in fig. 11. The red dashed line (left panel) indicates the equilibrium value of the entropy density at the corresponding energy density, computed in section 3.1. From the plot it is clear that the entropy production rate is non-negative, i.e.,  s sym ≥ 0. In fact, using the gravitational bulk equations of motion we finḋ which can be analytically proven to be non-negative following [27]. Because a derivation of (3.23) involves gravitational bulk constraint equations not directly used in the numerical evolution code, an important (dynamical) consistency check on the code is the vanishing of (3.24) The right panel in fig. 11 monitors the quantity (3.24). (3.25)

Unstable sBBL dynamics
Having reproduced the phase diagram and the linearized symmetry breaking dynamics of sBBL model in section 3.2.1, we now present fully nonlinear dynamical results. We focus on the unstable case only, as for simulations with E > E crit , after a brief non-linear regime, the system evolves to symmetric equilibrium configurations discussed above.
We consider first the evolution with g = −8 and E = 0.42978(1)E crit and E crit as in (3.13  i.e., with E < E crit . Different color coding represents different spatial resolutions of the numerical runs. relative to the AdS 4 Kretschmann scalar K AdS 4 (recall K AdS 4 = const = 24). Different color coding on the plots represents numerical runs from the same initial conditions, but with different spatial relation -the different number of collocation points:  [28] might be more suitable to capture high-gradients in the evolution.
Even more leverage could be achieved with adaptive mesh refinement as in [29].  i.e., with E < E crit . Different color coding represents different spatial resolutions of the numerical runs.
more dramatic growth of the entropy density s. Better implementation of the numerics could answer whether dynamics in sBBL model at g = −2 is physically different from that at g = −8, and whether this difference is attributed to the sign of E crit .

Conclusions
In [1] it was argued that the phenomenological holographic model introduced in [17] violates the weak cosmic censorship conjecture. The gravitational dual describes dynamics of certain QF T 3 with spontaneous symmetry breaking, but without an equi-librium ground state below the instability threshold 12 . An arbitrary weakly curved initial gravitational configuration was shown to evolve, in a finite boundary time, to a configuration with both the divergent area density of the apparent horizon and the Kretschmann scalar evaluated at the horizon. The divergence of the area density signals that the singularity mechanism is robust -while the evolution was restricted to spatially homogeneous and isotropic states of the QF T 3 , any finite entropy density state with broken (boundary) spatial translational invariance and/or rotations can not dominate late time dynamics. Following the gravity-fluid correspondence [5], the BBL holographic model points to a development of the singularities from regular initial conditions in corresponding relativistic fluid mechanics.
In this paper we attempted to address two questions: (1): can BBL scenario be realized in a top-down string construction? (2): what is the role of supersymmetry on the singularity development?
Concerning (1), it was pointed out in [1] that consistent truncations studied in [21] exhibit the exotic thermodynamics of [17] in grand canonical ensemble. We studied here the corresponding models in details and established that there is no run-away instability in DG models in microcanonical ensemble -rather, we found yet another realization of the mean-field spontaneous symmetry breaking mechanism of [15,16].
Concerning (2), we showed that it is straightforward to modify BBL model to mimic the structure of the bulk scalar coupling to gravity ubiquitous in gauged supergravity consistent truncations [24,25]. The scalar potential of the "supersymmetric" generalization of the BBL model (called sBBL here) remains unbounded from below. The equilibrium phase diagram and the linearized symmetry breaking dynamics in BBL and sBBL models are conceptually identical. Further studies (and a better numerical code) are needed to firmly establish whether sBBL model also evolves to a geometry with divergent area density of the apparent horizon and the curvature. As we pointed out, the latter might depend on details of the scalar superpotential, i.e., the nonlinear coupling g.
Finally, the challenge remains to find embedding of BBL mechanism in string theory. 12 Potentially related phenomenon was reported in [30].

A Numerical setup
We adapt numerical code developed in [1] to study dynamics of spatially homogeneous and isotropic states in DG and sBBL models.

A.1 DG-B model
We introduce a new radial coordinate maintaining ′ ≡ ∂ x and˙≡ ∂ t , and redefine the fields as follows Using (2.32) and (2.37), we find the asymptotic boundary expansion x → 0 + for all the fields: The rest of the code implementation is as in [1].