Unstable horizons and singularity development in holography

In holographic applications one can encounter scenarios where a long-wavelength instability can arise. In such situations, it is often the case that the dynamical end point of the instability is a new equilibrium phase with a nonlinear scalar hair condensate outside the black hole horizon. We here review holographic setups where symmetric horizons suffer from long-wavelength instabilities where a suitable equilibrium condensate phase does not exist. We study the dynamics of the simplest model in this exotic class, and show that it uncovers arbitrarily large curvatures in the vicinity of the horizon which asymptotically turn such region singular, at finite time with respect to the boundary theory.


Introduction
The string theory/gauge theory correspondence [1,2] is by now a mature framework exploited to address interesting questions in strongly coupled gauge theories that are often inaccessible with other theoretical tools. In a nutshell, this duality establishes a holographic correspondence (a dictionary) between two objects: a non-abelian gauge theory and a higher-dimensional gravitational theory/string theory in asymptotically anti de-Sitter space-time. One particularly appealing consequence of this correspondence is the fact that questions about the gauge theory in strongly coupled regimes are mapped onto questions in classical gravity. Likewise, a dual gauge theory perspective allows for different, and often intuitive, understanding of instabilities in black hole/black brane spacetimes.
Indeed, lets recall the physics of holographic superconductors [3,4]. Consider the four-dimensional effective gravitational action 1 in asymptotically AdS 4 (dual to a threedimensional conformal field theory CF T 3 ), (1.1) The four dimensional gravitational constant κ is related to a central charge c of the CF T 3 as

2)
F µν is a field strength of a global U (1) symmetry of the CFT, and φ is a (neutral) gravitational bulk scalar with which is dual to a dimension ∆ φ = 2 operator O φ of a boundary theory 2 . Note that there is Z 2 symmetry in the model, associated with this scalar, φ ↔ −φ. As it is wellknown, there are two phases of equilibrium states of this CF T 3 at a finite temperature T and a U (1) global symmetry chemical potential µ, distinguished whether O φ = 0 or O φ = 0. The O φ = 0 phase exists for arbitrary temperature T ≥ 0 and it is gravitationally described by Reissner-Nordstrom AdS 4 black brane with unbroken Z 2 symmetry, correspondingly φ ≡ 0. For sufficiently small T /µ this Z 2 symmetric phase becomes unstable [4]: on the gauge theory side of the correspondence the instability is a generic instability of the order parameter in the mean-field theory of thermal secondorder phase transitions; on the gravity side, this is a Gregory-Laflamme (GL) type instability [6] (in the sense of being unstable to long-wavelength perturbations) due to scalarization of the Reissner-Nordstrom AdS 4 black brane horizon. To understand the gravitational origin of the instability the authors of [4] noted that even though the scalar φ is above the AdS 4 Breitenlohner-Freedman (BF) bound as the Reissner-Nordstrom AdS 4 black brane becomes extremal (T /µ → 0), it develops AdS 2 × R 2 near horizon geometry with the curvature radius L 2 2 = L 2 6 . In this limit (1.5) 1 We set the radius L of an asymptotic AdS 4 geometry to unity. 2 φ has two alternative quantizations in AdS 4 [5]; our results do not depend on this choice. and the bulk scalar φ becomes unstable (the quasinormal frequency of its linearized fluctuations has Im[ω] > 0). The condensation of the gravitational scalar φ at low temperatures is dynamically saturated by nonlinear effects, spontaneously breaking Z 2 symmetry and leading to a new equilibrium phase of the CFT with O φ = 0.
There exist many studies and generalizations of the described phenomena in holography 3 . In this paper we focus on a less-known, exotic property of certain black brane/black hole horizons. As in the example of the holographic superconductor above, imagine a holographic 4 horizon with a discrete (or continuous) symmetry. Suppose that there is a critical energy 5 or energy density (for gauge theory states with translational invariance) below which the horizon becomes unstable with respect to symmetry breaking (GL) fluctuations. There is an equilibrium phase with spontaneously broken symmetry, branching off the GL onset of the instability, yet, this phase does not exist below the critical energy; moreover, it has lower entropy above the criticality than the symmetric phase. Thus, the horizon representing the symmetric thermal state is unstable, but it is unknown what the end point of its instability is.
To our knowledge, the first realization of the above exotic scenario appeared in [8] which was later found in a top-down holographic model in [9]. Finally, the same exotic physics is behind the leading instability of small black holes in AdS 5 × S 5 (dual to SO(6)-symmetric states of strongly coupled N = 4 SYM plasma) [10][11][12][13]. Here, we study the endpoint of this exotic horizon instability.
In the next section we briefly review the bottom-up model of [8]. We discuss the equilibrium states of the system, and the linearized instability of symmetric phase states at low energy densities. We construct the symmetry-broken phase of the system and demonstrate that it is never preferred dynamically. In section 3 we employ a characteristic formulation of the gravitational dynamics [14] in our exotic model. We confirm the onset of the GL instability dynamically, and compare the linear growth (below the criticality) and decay (above the criticality) rates of the symmetry breaking fluctuations with the corresponding quasinormal mode (QNM) computations of section 2. Next, we present results for the full-nonlinear evolutions of unstable horizons.
Details of the numerical implementation as well as the convergence and the validation of the code are delegated to Appendix A. We conclude and discuss open questions in 3 See [7] and references therein. 4 We point out this feature occurs in top down holographic models, and thus is of importance to issues of equilibration and thermalization in strongly coupled gauge theories. 5 As we study dynamical phenomena, we work in a microcanonical ensemble.

Exotic hairy black holes at equilibrium
In this section we review the bottom-up holographic model of the exotic black holes presented in [8].
The effective four-dimensional gravitational bulk action, dual to a field-theoretic setup discussed in the introduction, takes the form where we split the action into (a holographic dual to) a CFT part S CF T ; its deformation by a relevant operator O r ; and a sector S i involving an irrelevant operator O i along with its mixing with O r under the renormalization-group dynamics. We take bulk quantization so that the scaling dimension of O r is ∆ r = 2; the scaling dimension of O i is ∆ i = 4 . In order to have asymptotically AdS 4 solutions, we assume that only the normalizable mode of O i is nonzero near the boundary.
The gravitational action (2.1) has Z 2 × Z 2 discrete symmetry that acts as a parity transformation on the scalar fields φ and χ. The discrete symmetry φ ↔ −φ is explicitly broken by the relevant deformation of the CFT, with Λ being the deformation mass scale, while the χ ↔ −χ symmetry is broken spontaneously. The mechanism for the long-wavelength instability at play in (2.1) was identified by Gubser [15] through the following observations: • consider the linearized dynamics of the χ-sector in the mass-deformed CFT dual to S CF T + S r in (2.1); • for the quartic coupling g < 0, the scalar χ has an effective mass • homogeneous and isotropic thermal equilibrium states of S CF T + S r at low temperature (energy densities) would result in large values of φ at the horizon of the dual gravitational description, thus driving m 2 χ below the effective BF bound.  given that the symmetric phase is perturbatively unstable with respect to linearized Z 2 -symmetry breaking fluctuations, the imaginary part of the frequency ω χ of these fluctuations at zero spatial momenta is positive (right panel).
A detailed analysis of the homogeneous and equilibrium states of the holographic model (2.1) in the canonical ensemble were presented in [8]. Here, we present results in the microcanonical ensemble. We omit all the technical details as the following discussion is a special case of the dynamical setup of section 3.
• There are two equilibrium phases of the holographic model (2.1), distinguished by the symmetry property under χ ↔ −χ: the symmetric phase with O i = 0, and the symmetry broken phase with O i = 0.
• The entropy density of the symmetric phase s sym as a function of the energy density E is presented in figure 1. While this phase is thermodynamically stable > 0, it is perturbatively unstable with respect to a linearized symmetry breaking fluctuations [16]: for E < E crit , with the quasinormal modes of the symmetry breaking linearized χ-fluctuations develop a positive imaginary part, Im(ω χ ) > 0. As emphasized in [16], this model is one of the explicit counterexamples of the Gubser-Mitra "correlated stability conjecture" [17,18].  on | k| (here E/E crit = 0.89780 (8)).
• Notice that there is a relation between the behavior of the unstable χ-mode and the 'GL' instability, in that the χ-instability requires long wavelength modes, i.e., the instability is cut-off at However, there is a clear distinction: in the GL instability [6] the unstable mode is hydrodynamic, while the χ-QNM behaves non-hydrodynamically away from the critical point, i.e., Im(ω χ ) = 0 as the spatial momentum (along the translationally invariant directions of the horizon) vanishes, | k| = 0.
• The expectation value of O i = 0 in the symmetry broken phase of the model as a function of the equilibrium energy density is presented in figure 3. This phase The equilibrium symmetry broken phase is never realized in a microcanonical ensemble as it has smaller entropy density compared to the symmetric phase for the same energy density. crit ,  figure 3, see also [8]. Notice that E vacuum sym → 0 in the conformal limit Λ → 0; to obtain the better estimate for E vacuum we extended the analysis of the symmetric phase to the lowentropy region, as shown in figure 5, and extrapolated the energy-entropy data to zero entropy density (indicated by the (red) dashed line). The right panel shows the dependence of the temperature T sym of the symmetric phase black holethe limit E → E vacuum appears to correspond to an extremal limit.
In this section we focused on the static phase diagram, along with the linearized (in)stabilities of these phases, of the holographic action (2.1), dual to non-conformal In what follows we will discuss the dynamical case. In section 4, we comment on properties of the model with QF T 3 residing in R×S 2 .
Additionally, we comment on the extension of the model (2.1) with the gravitational potential for the scalar χ bounded from below.

Dynamics of the exotic unstable horizons
In this section we discuss dynamical properties of the holographic model (2.1), with the boundary QF T 3 formulated in R 1,2 . We follow closely the holographic numerical framework in the characteristic formulation as described in, e.g., [14].

Dynamical setup
We assume translational invariance along the spatial directions of the boundary. The relevant fields are described by Einstein equations define the following evolution equations of motion: together with the constraint equations: where ≡ ∂ r and d + ≡ ∂ t + A ∂ r . The constraint equations are preserved by the evolution equations provided they are satisfied at a given timelike surface (e.g., [19][20][21]) -which in our case is the AdS boundary.
The general asymptotic boundary (r → ∞) solution of the equations of motion, given by is characterized by two constants {p 1 , µ}, and three dynamical variables These parameters have the following interpretation: • p 1 and p 2 (t) are correspondingly the non-normalizable and normalizable coefficients of the bulk scalar φ, 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, Details of the numerical implementation, specific choices of the initial conditions (3.11) used, and code convergence tests can be found in Appendix A.

Dynamics of the symmetric sector
To study dynamics in the symmetric sector, we adopt initial conditions as described in Results of a typical evolution are presented in figures 6 and 7. Here, the energy density is below the critical one (2.5), The left panel of figure 6 shows time evolution of the expectation value of O r . Within a time scale t ∼ Λ −1 the system equilibrates. The equilibrium expectation value, defined is represented by a (green) dashed line. We used the value of O e r obtained from the independent analysis of the static configurations, reported in section 2, evaluated at the energy density (3.13). Consistency of (3.14) is an important check of our evolution.
The right panel of figure 6 illustrates the system's approach to equilibrium, which displays a typical φ quasinormal mode ring-down of the exotic black hole horizon.
The entropy density is an intrinsically equilibrium concept in QFTs. One benefit of the holographic framework is that it provides a well-motivated notion of the non-equilibrium (even far from equilibrium) entropy. Following [22,23] we identify nonequilibrium entropy density s with the Bekenstein-Hawking entropy corresponding to the apparent horizon (see (3.10)) area density The right panel of figure 7 shows the evolution of thus defined dynamical entropy density. Notice that in line with the second law of thermodynamics,ṡ(t) ≥ 0 and approaches at late times the equilibrium value s sym , computed independently for the static configuration with the energy density (3.13).
While we can study the non-equilibrium dynamics in Z 2 -symmetric sector of the holographic model (2.1), completely suppressing the χ-scalar fluctuations as in (3.12), at E < E crit , and in particular in the discussed example (3.13), this is an unphysical approximation -in realistic settings the fluctuations of the χ scalar will always be present, and would destabilize this Z 2 -symmetric dynamics. In the right panel of , (3.16) relative to the AdS 4 Kretschmann scalar K AdS 4 (recall K AdS 4 = const = 24), to emphasize the fact that even if the symmetric sector is unstable, its bulk dynamics is weakly curved. Thus, higher derivative supergravity and string corrections are arguably irrelevant for the onset of the Z 2 symmetry breaking instability of the exotic black hole horizons.

Long-wavelength (GL-type) instability of the symmetric sector
In this section we study linearized fluctuations of the Z 2 symmetry breaking operator O i in the symmetric phase of the holographic model (2.1). We initialize the symmetric sector of the model as explained in section 3.2 for energy densities above/below the critical one. The bulk scalar χ, dual to an irrelevant operator O i , is initialized as (in field χ and which determines its dynamics -which is kept unchanged.   (7) .
This is in excellent agreement with the independent computation of the χ-scalar QNM frequencies reported in figure 1 at energy density (3.18):   Notice that the fluctuations of χ do not oscillate (both in the stable, ie. those that give rise to equilibrium, and unstable scenarios), i.e., We believe this is a reflection of the spontaneous character of the symmetry breaking due to these fluctuations at the horizon 6 together with the boundary conditions adopted.

Fully non-linear evolutions of stable and unstable black holes
We now turn our attention to the fully non-linear behavior. Thanks to the simulations' ability to account for the backreaction of the field χ a rich phenomenology is uncovered.
To aid in the interpretation of the results, we monitor several quantities: • The dynamical behavior of p 2 and q 4 .
• The area of the Apparent and Event horizons (see Appendix A.5).
• The behavior of the Kretschmann curvature scalar K = R abcd R abcd (normalized by the value of K for pure AdS).
As a first case of study, we confirm that for E > E crit the behavior observed is consistent with that captured by the linearized analysis described in section 3.3. For this case, the system asymptotically approaches a stationary hairy black hole which is evidenced by a non-zero asymptotic value of p 2 as illustrated in figure 10 as well as the behavior of the normalized curvature scalar K shown in figure 11. This figure, also shows that at late times the event and apparent horizon coincide and remain stationary.
On the other hand, the case where E < E crit -identified in the previous section as unstable-leads to a markedly different behavior. For concreteness, we concentrate on the particular case defined by the following configuration.
Energy density (with Λ = 1) All these studies confirm the observed behavior that we describe next.
As the field φ "rolls down" the unbounded potential, the energy gained impacts the dynamics of its normalizable coefficient (p 2 ) as well as the normalizable coefficient of χ which grows without bounds as shown in figure 12. This behavior is evidenced in the black hole, which grows fast and eventually reaches the AdS boundary in finite asymptotic time as illustrated in figure 13. The figure shows both the apparent horizon (AH), and the event horizon (EH) as well as the curvature scalar evaluated on them.
Clearly, as time progresses, the AH approaches the EH and both asymptote to infinite size in a finite amount of time. This asymptotic behavior can be fit by the expression This fit indicates a finite time divergence at t ≈ 4.30. We find a similar asymptotic behaviour for the Kretschmann scalar evaluated at the horizon with K Σ 3 EH , as seen in figure 14. Thus, at late times, (3.26) consequently, K EH diverges in finite time at the boundary of AdS. Additionally, the (normalized) scalar curvatures on the AH and EH diverge with K EH ≤ K AH . Naturally, the code is eventually unable to keep up with the radically rapid dynamics which requires ever smaller timesteps to capture the following behavior. Nevertheless, we have been able to extract convergent solutions up to a sufficiently late stage to understand the behavior and fate of the spacetime. The picture that arises is that the spacetime explores arbitrarily large curvatures in finite time, and outgoing null geodesics emanating from such regions reach the boundary of AdS in finite asymptotic time as indicated in figure 15. This behavior would violate the spirit of the weak cosmic censorship conjecture, in that far observers can be reached by signals emanating from arbitrarily curved spacetime regions, and is similar to that recently reported in [26].

Conclusions
In this work we studied an interesting instability of the black hole horizons, observed first in [8]: below some critical energy density the horizon is unstable with respect to fluctuations spontaneously breaking a discrete symmetry. However, there is no static end-point associated with the nonlinear build-up of the symmetry breaking condensate (scalar hair at the horizon), as opposed to typical constructions of holographic superconductors [3]. The instability is perturbative in nature (i.e., describing a second order phase transition), and is triggered by an arbitrary small amplitude of the symmetry breaking mode, provided the conserved energy density of the state E is below a critical energy density. As a result, the onset of the instability, and dynamics close to it, can not be affected by higher-order nonlinearities in the gravitational scalar potential as long as the amplitude of unstable modes remain small. Additionally, the instability initiates in the long wavelength regime, i.e., at small bulk curvature, and thus can not be removed by higher derivative corrections to the gravitational effective action.
We studied the future development of the instability, using a characteristic formulation of asymptotically anti-de Sitter gravitational dynamics [14], and argued that, While the analysis of this work is focused on a specific phenomenological model of gauge theory/gravity correspondence, represented by the effective gravitational action (2.1), the phenomena described, i.e., horizon instability without the static end point, is realized within bona fide holographic correspondence scenarios (e.g., [9,13]). We expect that curvature singularities also arise in those models as well 7 . The observation that curvature singularities might arise dynamically implies that consistent truncations of string theory and supergravity, while suitable to address static states in the theory, may fail in dynamical settings -when the evolution enters the regime of highly curved geometry; and in such cases stringy corrections will be important.
What are the holographic implications of a singularity developing evolution for the boundary gauge theory? A standard lore is that states of a closed non-integrable interacting system with large number of degrees of freedom should dynamically equilibrate [28]. In the context of holography, early indications supported this for generic far-from-equilibrium, arbitrary low-energy states of strongly coupled conformal gauge theories, even when the dynamical evolution of these states was artificially restricted to symmetric submanifolds of the full phase space of the theory 8 [29,30]. Shortly there- 7 It is a straightforward exercise to examine this in a holographic model [9]; dynamics of small black hole localization in AdS 5 × S 5 is much mode difficult [27]. 8 The states in question were spatially isotropic, and invariant under all global symmetries, i.e., the R-symmetry.

EoS
AdS Boundary AH EH High K region In the diagram, "EoS" refers to the "End of the Simulation" while the star refers to the blow up of the Kretschmann at the boundary in finite time.
after it was argued [31][32][33] that in fact symmetric phase space of holographic conformal field theories has islands of stability that never equilibrate. In this study we identified yet another possibility: initial states of holographic strongly coupled gauge theories, well described classically in the gravitational frame, evolve to a singularity in finite time. Singularity is a signature for a breakdown of an approximation, and we see two possible reasons. First, it is possible the singularity is an artifact of our restriction of the state evolution to symmetric submanifolds of the full phase space of the theory and that symmetry breaking modes would allow for a smooth evolution. As we discuss below, hydrodynamic modes in the system can be gapped, leading to the same qualitative behaviour. It is more difficult to argue for the absence of light modes spontaneously breaking internal symmetries -the singularity observed might be an indication that some of these symmetries must be dynamically broken during the evolution (similar ideas were proposed in [27,34]). Second, the state evolution in the gravitational frame of the holography might not be always semiclassical.
We find it important to discuss another possible limitation of the study carried out here and its conclusions. We restricted the dynamics in our model to preserve boundary homogeneity and isotropy. One might argue that the physical phenomena discussed here arise as a consequence of such a restriction, and that a sufficiently generic initial state would smoothly evolve to an end point where these symmetries are spontaneously broken. We do not have a full answer to this question -gravitational simulations in the presence of spatial inhomogeneities are beyond the scope of this paper. At the very least, in the holographic model studies in [35] there is an exotic instability discussed, without the spatially modulated endpoint as well 9 . Since the instability and the evolution towards the singularity in our model can be triggered by arbitrarily small amplitude fluctuations of the χ-mode, i.e., energetically arbitrary close to the critical point, the potentially physics-modifying hydrodynamics modes can be gapped, rendering them irrelevant to the question as to whether or not the singularity observed is physical. To demonstrate this, we modified our model with a boundary with topology R 2 → S 2 . All the main features described in the former case remain in the latter, in particular: there is a horizon instability, there is no static end point below some critical energy density associated with the onset of the perturbative instability, the low-energy SO(3)-invariant states evolve to a singular solution. It would be interesting to explore in details the role of additional massless fields at the threshold of instability, and their effect on the singularity development.
A feature of the bulk scalar potential of our holographic model (2.1) is that it is unbounded in the χ-direction (recall that the nonlinear coupling g < 0). We study in Appendix B a modification of the model which "bounds" the χ-potential with higher order, nonlinear in χ, interactions 10 . Of course, the linearized instability is unaffected; likewise, the unstable phase with O i = 0 for E > E crit is unchanged qualitatively (close to E crit the higher-order nonlinear terms in the gravitational potential are suppressed).
However, we find a new static black hole phase with O i = 0 , that exists both for E < E crit and E > E crit for the bounded potentials and, at least in the vicinity of perturbative instability, has higher entropy density than the symmetric phase. Unlike the exotic branch of the black holes, this new phase does not bifurcate from the onset of long-wavelength instability of the symmetric phase. For E < E crit this new static phase is always the end point of the evolution; for E > E crit the new symmetry broken phase can only be reached if the initial amplitude of the symmetry breaking fluctuations is sufficiently large -the symmetry broken phase is a potential barrier separated from the 9 We would like to thank Ben Withers for bringing the reference to our attention. 10 We would like to thank Jorge Santos for raising the issue of the unboundedness of the scalar potential in our model with its potential effect on the singular evolution that prompted this analysis. symmetric phase whenever E > E crit . Our model (2.1) is a phenomenological example of the holographic correspondence, thus one might worry whether curvature diverging scenarios described here is realized in genuine (top-down) holographic dualities. We believe the answer to the question is in the affirmative: First, the unbounded potentials are ubiquitous in holography -a typical example is a well-studied N = 2 * holography [36][37][38], where the bulk gravitational scalars {α, χ} effective action takes form, (4.1) The reason why the scalar potentials in supergravity constructions can be unbounded comes from the fact that they arise from the superpotential as (for the N = 2 * example (4.1))

(4.2)
It is the − 1 3 W 2 contribution to V in (4.2) that is responsible for the unboundedness of V .
Second, the scalar potential in the top-down embedding of the exotic black hole phenomena [8] constructed in [9] (see eq.(2.32) there) is unbounded from below: V (ϕ) = −2 (2 + cosh(2ϕ)) . In new variables (A.3), the equations of motion used to evolve the system take form: Numerical code is organized as follows. • The value a h is determined from the stationarity of the apparent horizon at x = 1 as explained in the following subsection.
• A.2 Apparent horizon and the boundary condition for a Our numerical implementation requires an independent computation of a h ≡ a(t, x = 1) (see (A.14)), given radial profiles {p, p , q, q , s, s , dp, ds, dq} and the diffeomorphisms parameter λ at time t. Following [14], this is done by enforcing the time-independent location of the horizon. Apparent horizon is located as x = x h such that Assuming x h = 1, dx h dt = 0, and using equations of motion (A.5)-(A.9) we compute a h from Denoting p h , dp h , q h , dq h , σ h ≡ { p , dp , q , dq , σ } (t,x=1) (A. 19) we find (A.20)
Recall that both φ and χ are left invariant under the reparametrization transformations: To maintain this invariance we specify initial conditions for {p, q} (in λ 0 -invariant way) in terms of two amplitudes {A p , A q }: We then proceed as follows 12 : • given {A p , A q } we set λ 0 = 0 and perform [Step 2] (A.11) and [Step 3] (A.12); • having enough data, we follow (A.3) to compute the profile d + Σ(t = 0, x); • we find numerically the root x = x 0 of the equation • we set the trial value of λ 0 as 12 For this procedure the integration range over the radial coordinate x might exceed unity. tΛ Radius x Figure 17: Illustration of representative null rays traced to find the event horizon. The inset shows a zoom-in at late times which aids to visualize how null rays starting at different locations converge as they are traced backwards in time.

A.5 Event Horizon Finder
To find the event horizon we trace null geodesics at late times back in time and determine the surface R(t) where they converge. To do so, we start from g ab n a n b = 0 , (A. 26) where n a denotes the null tangent vector to the geodesics. Using (3.1) and the field redefinition (A.1) and (A.2) this relation implies which we solve numerically using either a RK4 integrator or an second order implicit integrator. The results obtained with both methods converge and are in excellent agreement. As described briefly above, we consider a collection of starting points at different radii and bisect the resulting behavior to home-in on R(t). The left (green) vertical dashed line corresponds to the energy density (3.18), E green,lef t < E crit . The right panel in figure 18 shows the order parameter, O i , for the new symmetry breaking phase as a function of the energy density E.  As shown in figure 18, the new symmetry broken phase has an interesting feature in the narrow energy range: Here, the symmetric phase is perturbatively stable, but the new symmetry broken phase is nonetheless entropically favorable; thus one excepts that the broken phase can be reached dynamically only if the amplitude of the initial symmetry breaking fluctuations is sufficiently large. We find that this is indeed the case. For the results presented in figure 21, O i t=0 = 0.0252 Λ 4 and the system equilibrates to a (metastable) symmetric phase. Figure 22 represents results of the simulation for the initial condition with