A Perturbative RS I Cosmological Phase Transition

We identify a class of Randall-Sundrum type models with a successful first order cosmological phase transition during which a 5D dual of approximate conformal symmetry is spontaneously broken. Our focus is on soft-wall models that naturally realize a light radion/dilaton and suppressed dynamical contribution to the cosmological constant. We discuss phenomenology of the phase transition after developing a theoretical and numerical analysis of these models both at zero and finite temperature. We demonstrate a model with a TeV-Planck hierarchy and with a successful cosmological phase transition where the UV value of the curvature corresponds, via AdS/CFT, to an $N$ of $20$, where 5D gravity is expected to be firmly in the perturbative regime.


Introduction
Randall-Sundrum (RS) models [1] offer an attractive solution to the hierarchy problem, and put the cosmological constant problem [2] into a new perspective [3,4]. In RS models, the warping of AdS space geometrically generates large hierarchies. Interestingly, the effective 4D value of the cosmological constant is a sum of terms involving the bulk 5D cosmological constant, and two brane tensions associated with the UV and IR branes. The tiny observed value of the cosmological constant is obtained by separately tuning the UV brane tension against the bulk cosmological constant, and the IR brane tension against the same bulk cosmological constant.
This "double" fine-tuning in the 5D theory is necessary to force a flat direction for the location of the branes, for which the potential would otherwise cause either collapse of the geometry or a run-away. The Goldberger-Wise stabilization mechanism offers a solution to this tuning problem, with a 5D scalar field developing a spatially varying vacuum expectation value in the bulk of the extra dimension, and leading to a non-trivial potential for the location of the IR brane [5,6], stabilizing the "radion." However, this solution relied upon the mistune in the brane tensions being small to begin with, so that the bulk scalar field vev did not deform the geometry significantly from AdS and that the backreaction of the scalar field on the geometry remained small. Thus a degree of tuning remained, as naive dimensional analysis (NDA) from consideration of graviton loops suggests that the mistune be parametrically larger, with natural values for a quartic coupling for the radion being λ ∼ O [(4π) 2 ].
In addition to this naturalness issue, a more phenomenological and pressing problem plagues these models: The phase transition during which the vacuum expectation value for the radion develops is first order, and estimates of bubble nucleation rates in early universe cosmology strongly suggested that a RS phase transition would not proceed to completion due to Hubble expansion out-pacing true-vacuum bubble creation. In the region of parameter space where nucleation is fast enough, perturbativity of the 5D gravity theory is right on or past the threshold of being lost [7][8][9]. In this work, we address whether a recently studied class of geometries which contain a light radion mode but deviate far from AdS achieve a better transition rate while remaining perturbative.
There is strong motivation for considering such models. In terms of the AdS/CFT correspondence, the double tuning of RS in the absence of a stabilization mechanism has a natural interpretation [10][11][12][13][14]. The tuning of the UV brane tension against the bulk cosmological constant is viewed as a tuning of the bare cosmological constant in a nonsupersymmetric CFT very close to zero. This is required as a cosmological constant term would explicitly break conformal invariance, yet there is no supersymmetry to enforce this cancellation. The second tuning of the IR brane tension is interpreted as a tuning associated with the scale invariant quartic associated with an order parameter associated with spontaneous breaking of the CFT. The flat direction for the "radion" degree of freedom in RS appears as a tuning of this allowed parameter in the CFT to zero. If non-zero, such a quartic coupling would forbid the generation of a condensate that spontaneously breaks the CFT [15]. A solution to this problem appears if one allows a deformation of the CFT, i.e. by the introduction of a near-marginal operator. This small scale dependence effectively deforms the scale invariant quartic into a more generic potential that may have non-trivial minima away from the origin. The Goldberger-Wise stabilization mechanism is a dual to this prescription, but, as noted above, has tuning issues as well as cosmology problems.
It had long been thought that this fine tuning is unavoidable, as it reflects a coincidence problem in the 4D CFT dual -a flat direction in the theory that happens to coincide with a near-zero in the β functions for the theory [16][17][18]. However, it has been shown that if the scalar potential has only a soft dependence on φ, with the coefficients of the higher order interaction terms in the GW bulk potential being small, then the scalar field enters a significant back-reaction regime before the higher curvature terms come to dominate and perturbative control is lost. It has been shown that despite this large back-reaction, the dual theory is still conformal, and there is still a light dilaton that realizes scale invariance non-linearly. These "soft-wall" scenarios are the models that are of interest in this work. This ansatz for this type of bulk scalar potential is equivalent, via the AdS/CFT dictionary, to having a beta function in the CFT that remains small for a large range of the coupling. With this type of presumed dynamics, the coupling explores a large range of values during the running, and the scale invariant quartic could potentially find a zero, essentially finding a flat direction dynamically, and permitting a condensate that spontaneously breaks the approximate conformal invariance without fine tuning [19][20][21]. Other holographic studies of this scenario show that the dilaton mass in such models is suppressed relative to the breaking scale, and the cosmological constant is also parametrically suppressed [22][23][24][25].
In this work, we explore aspects of the cosmological phase transition in these softwall models. Efforts first are focused on a clear exposition of the theory of the dilaton effective potential at vanishing temperature. We then perform numerical calculations of the zero temperature potential for various ranges of the free parameters. Next, we study the theory of the model at finite temperature, and again perform detailed numerics of the finite temperature potential. Finally, we put the above results to work on the problem of the early universe conformal phase transition, finding an enhanced nucleation rate in soft-wall dilaton scenarios, and a phase transition that completes for much smaller (and thus perturbative) values of the curvature corresponding to a larger N dual CFT. Finally, we comment on the potential for the early universe conformal phase transition to be observed as a stochastic gravitational wave background signal due to the dynamics of bubble collisions.

Zero-Temperature Dilaton Effective Theory
We consider classical solutions to theories with a real 5D scalar field minimally coupled to gravity. The action we consider has a bulk contribution given by: where κ −2 ≡ 2M 3 * , with M * being the 5D planck scale. 1 We consider metric solutions with flat 4D slices: which can equivalently be expressed in coordinates y = A( y) that we find particularly convenient for this work: There are branes at orbifold fixed points taken to reside at y = y 0 , and y = y 1 . The scalar field has brane localized potentials at these points: Utilizing˙to represent derivatives with respect to y, the Einstein and scalar field equations can be written as: The Einstein equations can be used to eliminate G in the scalar field equation of motion: The total value of the classical action can be expressed as a pure boundary term. In particular, after substituting for the kinetic and potential terms for φ using the Einstein field equations, and taking into account contributions from singular terms in the scalar curvature at the orbifold fixed points, the resulting 4D effective potential is given by [20] Since the effective action is a pure boundary term, the 4D potential depends only on the asymptotic behavior of the geometry and the scalar field.
In the next two subsections we discuss the application of these equations first to the case of constant bulk potential V (φ) = − 6k 2 κ 2 , which is a review of previous results in the literature placed in the context of the motivation for this work. This is then extended to more general potentials that correspond to our weak φ dependence ansatz where we perform a numerical analysis over a broad range of parameter space. The constant potential case corresponds via AdS/CFT to an undeformed CFT. In the section on general potentials, we add a term to the 5D action that corresponds to sourcing a marginally relevant operator that stabilizes the pure scale-invariant dilaton quartic coupling typical for conformal field theories.

Example: Constant Bulk Potential
The case of constant potential can be solved analytically [26], and the result for φ is given by 12 κ 2 log e 4(y−yc) 1 + 1 + e 8(yc−y) . (2.10) The integration constant y c is chosen so as to correspond to the value of y for which the behavior of φ changes qualitatively from φ ≈ constant = φ 0 to a behavior that is linear in We can also evaluate the expression for G(y) exactly. Taking V = − 6k 2 κ 2 , and defining f ≡ ke −yc , and µ ≡ ke −y we have (2.11) With the above information we can extract the dilaton potential. We take the "stiff wall" limit where we presume that the boundary potential fixes φ at particular values on the branes: φ(y 0 ) = φ UV and φ(y 1 ) = φ IR . We also take the potentials in this limit to be pure tensions: V 0 (φ UV ) = Λ 0 and V 1 (φ IR ) = Λ 1 . Defining µ 0 = ke −y 0 , µ 1 = ke −y 1 we have: The scalar boundary conditions determine a combination of the free parameters µ 0 , µ 1 , f , and φ 0 . We fix µ 0 by matching the 4D observed Planck scale, and for the purposes of this discussion, we will hold f fixed. The UV boundary condition essentially sets φ 0 ≈ φ UV , up to terms of order (f /µ 0 ) 4 . The IR boundary condition fixes the ratio so that µ 1 might be replaced by a function of f in the expression for the effective potential.
Neglecting terms of order (f /µ 0 ) 8 induced by the explicit breaking of conformal invariance associated with sourcing 4D gravity at the scale µ 0 , the effective potential as a function of f can be written as (2.14) The first term in this expression is the contribution to the bare cosmological constant. This is expected to be either tuned to zero by choosing Λ 0 = 6k κ 2 , or made vanishing by the introduction of additional UV symmetries such as supersymmetry. The second is the contribution to the cosmological constant via the spontaneous breaking of conformal symmetry, or in other words, the dilaton quartic. As before, the quartic is a sum of two terms, one from the IR brane tension and the other from the bulk cosmological constant. The former is suppressed by the hierarchy in φ, such that for very large separation between φ UV and φ IR , a very large negative value of Λ 1 is required to cancel the positive contribution from the bulk geometry. 2 Note that we expect higher curvature operators that are expected to be induced by quantum corrections to give contributions to the dilaton potential, but the form is still be that of a scale invariant quartic plus derivative terms unless a non-trivial scalar potential is included.
The interpretation of this result is that even when the IR brane does not play a major role, there is the notion of a breaking scale of conformal symmetry given by f ≡ ke −yc . This scale corresponds in 5D to a position in the extra dimension at which the leading behavior of the curvature, or equivalently the scalar field evolution, makes a transition from one type of behavior to another. In Figure 1, we show schematically the behavior of the scalar Near AdS Region Condensate y=y c ⇠ Constant / y Figure 1: This cartoon shows the evolution of the field φ. It begins in a region where φ is nearly constant, and the geometry is nearly AdS. At a critical value of the coordinate, y = y c which is close in proper distance to a curvature singularity, φ then begins linear evolution, as shown and the curvature quickly grows large. The gap of the theory is set by the position of this "soft-wall." field evolution along the y coordinate. The behavior begins with slow evolution where the geometry is nearly AdS and then transitions to linear behavior where the curvature becomes large. Without an IR brane, a singularity at finite proper distance from the UV brane terminates the geometry. We term this region in y where φ is linear the "condensate" region. The effective potential for f is precisely what is expected for an approximately conformal theory with explicit breaking manifest in the form of a bare CC, and from the introduction of the Planck brane itself, making the position of this turnover of the 5D behavior of the scalar-gravity background a candidate for the dilaton.
Note that this "soft-wall dilaton" is unstable, just as in the original RS I model. Unless Λ 1 is large and negative quartic coupling is large and positive, driving f to zero in the absence of a stabilization mechanism. This means that the effective potential is minimized when the conformal symmetry is unbroken. Alternatively Λ 1 can be tuned to make the quartic vanish, and give the dilaton a flat direction. Further, the ansatz of flat 4D metric slices is only valid in the case that the total cosmological constant vanishes, or when all terms in the effective potential are arranged so as to exactly cancel each other.
The tuning of the bare CC, and the tuning of the dilaton quartic are precisely the two tunings that are required in two-brane RS models.
That there is a lack of stability of the constant potential case with soft wall breaking of conformal symmetry comes as little surprise. Typical conformal theories without supersymmetry do not support spontaneous conformal breaking due to the presence of the scale invariant quartic (in other words, the lack of scalar flat directions). 3 A deformation of the CFT, or in other words a departure from conformality, is required to stabilize a VEV against the scale invariant quartic. In the next two sub-sections, we demonstrate how deformations of the CFT (introduced in the AdS dual by considering a nontrivial bulk scalar potential) can stabilize the soft-wall dilaton.

Non-Constant Bulk Potentials
To stabilize the dilaton, we consider adding a deformation to the bulk scalar potential. For example, a mass term for φ could be added: Note that is defined to be dimensionless. The non-zero mass term for φ corresponds via AdS/CFT to a non-trivial quantum scaling dimension for the CFT operator that maps to φ via the duality [29]. If φ takes a non-trivial value on the boundary of AdS, then this operator is sourced in the dual (approximate) CFT, contributing as a small explicit violation of conformal invariance. In the constant potential limit, this operator is precisely marginal, and does not deform the CFT. When is negative, which is the ansatz we will take in our work, the scalar field is tachyonic, and tends to grow with increasing y. This is dual to sourcing a near-marginal relevant deformation of the CFT. This is somewhat similar to what occurs in QCD or technicolor-like theories, although in that case conformal invariance is badly broken in the infrared when the coupling becomes strong. We make the assumption that the φ-dependent terms in the potential remain small in the region of large φ, or in the approximate CFT dual, that the β-functions remain small even when the coupling becomes large. Condensation is triggered not by strong coupling, as happens in QCD, but rather by a coincidence of the coupling constant having a value associated with a near zero in the effective dilaton quartic. The slow running of the coupling over a large range of coupling values allows the theory to explore the landscape of quartics until a near zero is found and the theory condenses.
In the region of large back-reaction, the equations of motion do not admit analytic solutions to the equations of motion when is non-zero. The equations can be solved approximately using the method of boundary layer matching, however this approximation begins to break down when φ is large, but its behavior is not yet governed by the IR condensate asymptotics. Due to these difficulties, we resort to numerical solutions to the equations of motion to study the behavior of this system.
We presume throughout the rest of this work that the brane localized potentials enforce the stiff wall limit, and the branes are thus localized at positions where the scalar field φ is equal to φ 0 and φ 1 in the UV and IR, respectively. We do not expect any aspects of the analysis to change much if this condition is relaxed.
The UV brane tension is fine tuned to enforce the condition that the cosmological constant vanish in the limit when the conformal breaking scale goes to zero. This is a tuning of the bare cosmological constant -the remaining cosmological constant at nonzero values of f is due purely to the conformal symmetry breaking condensate. It is this remaining dynamical contribution which is suppressed by the small value of .
The value of f depends then depends on initial conditions for the derivative of φ, with the behavior of φ in the UV encoding the information regarding the scale of symmetry breaking. The brane potentials select the value of φ at the position of the branes, but this derivative of φ is not fixed. Varying this derivative is equivalent to varying over the value of the breaking scale, f . For a given IR brane tension, there is a given value for f at the minimum of the effective potential that we derived above -equivalently, minimization of the action fixes the derivative of φ at the position of the UV brane. The value for f in terms of the geometry is We have verified numerically that the masses of resonances obtaining their mass from conformal symmetry breaking (e.g. masses of gauge boson KK modes) track almost exactly linearly with the above definition for f , even in cases where the backreaction is very large. This definition agrees in the limit of vanishing backreaction with the usual definition, f ≈ ke −y 1 , but differs substantially from it in the regions of interest in this study.
We discuss results in terms of the dimensionless quantity N , where N is expressed in terms of the 5D curvature and Planck scale: For perturbativity of the 5D gravity model, N must be taken to be somewhat large. Note that the effective 4D Planck scale, M 2 Pl = N 2 16π 2 k 2 , should be held fixed, so that a particular value of N corresponds to a given value of k. It is also convenient to work in terms of a rescaled G: G = k 2 G. Finally, we also rescale both of the brane tensions: Λ 0,1 = 6 k κ 2 Λ 0,1 . With these rescaled parameters and functions, and with the expression for f in Eq. (2.16), we can express the effective dilaton potential in terms of dimensionless quantities as and we can write For numerical evolution of the scalar equation of motion in a manner that is independent of N , we define a dimensionless scalar field φ = κφ, in which case the scalar equation of motion is¨ For initial and final values of the scalar field φ 0 and φ 1 , we employ a rescaling that is common in the literature: where v 0,1 are the dimensionless values of the field on the branes, set in the stiff wall limit. We report results as a function of v 0,1 , N , , and Λ 1 .
In Figure 2 properties of the zero temperature dilaton potential are displayed for various values of the bulk mass term, , and for various values of the IR brane tension. The hierarchy between the 4D effective Planck scale and the dilaton vev, f , are shown. In addition, color shading indicates the value of the dilaton potential at the minimum (effectively the contribution to the cosmological constant from the dilaton vev). We show in two columns the dependence on the IR brane value of the bulk scalar field. On the left are plots that correspond to small back-reaction, which is the usual hard wall RS model where v 1 = 1. On the right, we display results for the soft-wall model where v 1 = 10, and backreaction in the IR region is large.
The values of Λ 1 that accomplish breaking of conformal symmetry are much larger than in the hard-wall model for small values of N . One might ask whether this is consistent with estimates arrived at using the tools of naive dimensional analysis. One must be careful in applying these tools in the soft-wall: Low mass graviton KK-modes do not penetrate the soft-wall to reach the IR brane, whereas higher mass KK modes have more significant overlap with the IR brane. This means that the local cutoff associated with violation of unitarity in KK-graviton scattering is much higher as the self coupling of the gravitons is very small in this region. It is not clear that the values for small N are completely reasonable, but we are most interested anyway in larger values of N , corresponding to perturbative 5D gravity models. As N increases, the values of Λ 1 rapidly approach much smaller values that are easily consistent with NDA expectations.
-V min

Finite Temperature
The phase structure of near-conformal 4D theories is of interest both as a theoretical question and one of phenomenology. If naturalness of the electroweak symmetry breaking sector is due to strong near-conformal dynamics, it is important to study the cosmology of such theories. Studies of the RS1 phase transition indicate that it is strongly first order, with a critical temperature suppressed in comparison with the value of the condensate. This is due to the presence of a near flat direction at the minimum of the dilaton potential. It is this which allows for the light dilaton, and also for a suppression in the contribution of condensates to the effective IR value of the cosmological constant. At finite temperature, such non-compact flat directions are lifted, sending the dilaton field value to the origin, thus evaporating the condensate.
In order to study the theory at finite temperature, the class of geometries we study is opened up to include the possibility of a horizon (or "black brane") at some point y = y h in the 5D coordinate [30][31][32]. In AdS space, the Hawking radiation from such a black hole allows the black hole to reach equilibrium with the thermal bath. The action associated with the classical solution corresponds to the thermodynamical free energy of the system. The geometry we study has metric function The presence of a horizon is associated with a zero in the horizon function h(y) at position y h . As we are considering a thermal partition function, we work in Euclidean metric signature, with the time coordinate compactified on a circle: t ∈ [0, 1/T ).
The equations of motion for the metric functions h and G, and for the scalar field φ are given byḧ The effective potential is still given by a pure boundary term, although the singular terms due to orbifold boundary conditions at a putative black hole horizon require special treatment, as we discuss later in this section.
The bulk contribution to the effective potential arises from using the equations of motion to express the bulk action as a total 5-derivative: 25) The curvature tensor has singularities at the orbifold fixed points that give additional contributions to the effective action. Integrating the action over these singularities at the UV and IR branes gives the following contribution to the effective potential: Note that the equation of motion for h enforces an exact cancellation between the twoḣ terms.
In summary, adding together the contributions to the potential when there is no black hole horizon, including the two brane localized potentials which each contribute g ind (y i )V i , are: Our goal is to replace the IR brane at y 1 with a black hole horizon at y h , such that h(y h ) = 0 [13], however due to the structure of the manifold near the horizon, one cannot simply take h(y h ) = 0 in the above equation. The reason for this is that the manifold near the horizon is typically singular, with a cone feature appearing in a given t − y slice of the geometry, as shown in Figure 3.
In order to study such a horizon for generic bulk scalar potential, we presume that the horizon function has a zero for some finite y = y h . We further presume that h(y) is analytic, such that it has a Taylor expansion in the vicinity of the horizon. In this case, we have˙h h ≈ 1 y−y h , with the sign determined by the fact that h is positive in the physical region y < y h , and that it is passing through zero.
This behavior of the horizon function determines a boundary condition for φ that arises from taking the near-horizon limit of the scalar field equation of motion: This boundary condition enforces regularity of the solution for φ at the horizon -without this condition, φ diverges in the approach to the horizon [33].
To compute the effective potential when the IR is screened by a black hole horizon, we need to pay closer attention to the treatment of singular terms at the orbifold fixed points at y = y 0 and y = y h . The scalar curvature is singular in both places. In the UV, the singular terms can be treated as before, yielding a contribution to the effective potential that is given by The IR contribution is calculated via a proper regularization of the 2D conical singularity. There is a conically singular geometry near the black hole horizon corresponding to a system that is out-of-equilibrium. Quantum effects will generally cause the singularity to emit radiation until it reaches equilibrium with the surrounding thermal bath, at the minimum of the free energy of the thermodynamical system.
If a theory admits solutions to the h function which vanish at some finite value of y h , then we can study such systems in the near-horizon limit. Considering the near-horizon limit of the metric, where h ≈ḣ(y h )(y − y h ), we have (displaying only the dt and dy components of the metric): .
We now go to "good" coordinates, where the metric is manifestly that of a cone: (3.31) The opening angle of the cone is given by with the overall minus sign ensuring positivity of the angle sinceḣ is negative at the horizon. By capping the cone with a sphere of radius r, which has constant curvature 2/r 2 , the contribution to the action is rendered finite and r independent, allowing a sensible r → 0 limit: Note that a factor of two has been included as the integral is over the entire S 1 space in the S 1 /Z 2 orbifold, leading to a double-copy of the spherical cap, one on each side of the orbifold fixed point. The singular IR contribution to the 4D effective potential energy is then given by The first term cancels exactly the corresponding UV term, and we can write the complete effective potential in the presence of the black hole horizon as The value of y h that minimizes the free energy as a function of the horizon location is obtained by inverting the following relation: The right hand side of this equation for arbitrary y h is interpreted as the temperature of the black hole. The value of the free energy at the minimum is Up to terms that violate conformal invariance due to the introduction of the Planck brane or the Goldberger-Wise potential, the equilibrium temperature that minimizes the effective potential as a function of y h is associated with the value of y h that removes the conical singularity. We can use the h equation of motion to express this equilibrium temperature in terms of the near AdS-Schwarzchild UV behavior of h and G:ḣ(y 0 ) ≈ −4e 4(y 0 − y h ) and G(y 0 ) ≈ k 2 . Note that y h is the position where the horizon would be if there were no deformation of the geometry due to the varying φ field. In the absence of scalar backreaction, y h = y h . From the equations of motion one finds that the presence of the back-reaction delays the onset of the horizon, establishing the inequality y h ≤ y h .
As the position of the horizon, y h , is greater than y h , the equilibrium temperature is larger than it would be in the absence of scalar backreaction. This is potentially problematic, as this would mean that the temperature is not necessarily a monatonic function of the position of the horizon. The temperature would in fact grow when the backreaction becomes sizable, causing a deviation between y h and y h . We see that the temperature grows with increasing y h when dy h d y h > 4/3. Note however, that the entropy S = 4π κ 2 e −3y h is monotonically decreasing with increasing y h . These high temperature solutions with low entropy are disfavored relative to those of equal temperature but small y h and thus larger entropy.

Constant Bulk Potential at Finite Temperature
In the case of V (φ) = − 6k 2 κ 2 , with no dependence on φ, the scalar field equation of motion has a significantly simplified relationship to h: This scalar field equation of motion is integrable, and we find that the solution is given by: where C l is an integration constant. We note that this equation immediately excludes the case of constant bulk potential as a candidate for a spontaneously broken CFT at finite temperature, or where h = 0 for some finite y in a non-trivial scalar field configuration. Clearly, if h is vanishing, but C l is finite then φ must be divergent at the position of the horizon, and the horizon boundary condition Eq. (3.28) cannot be satisfied.
The equations can be satisfied for one particular value: C l = 0, which corresponds to φ = constant. Solving the Einstein equations for this case yields h = 1 − e 4(y h −y) and G(y) = k 2 , corresponding to the AdS-Schwarzchild geometry. This configuration is dual to an unbroken exact CFT at finite temperature.

Generic Potential at Finite Temperature
We now calculate the results for the free energy when a non-trivial bulk potential is considered. Using the potential explored in the zero-temperature analysis, Eq. (2.15), we numerically solve the coupled scalar and Einstein equations for a range of temperature and the free parameters of the model. At high temperatures the theory is in a quasi-AdS-Schwarzchild geometry, with a free energy given by Eq. (3.35), while at low temperatures, the theory transitions to the zero temperature geometry studied in Section 2.
Geometries that minimize the free energy can be found for a large range of temperatures for each configuration of parameters. The analysis reveals that there are striking differences between the free energy as a function of the temperature when one includes or does not include the effects of backreaction on the metric. In Figure 4, we show the results of the numerical analysis in terms of the value of the free energy at the minimum, having extremized over the position of the black hole horizon. The curve corresponding to = 0 gives the free energy for the AdS-Schwarzchild solution, or equivalently, the free energy where backreaction is neglected. The remaining curves have non-trivial scalar field profiles due to non-vanishing values of the bulk scalar mass term, and affect the free energy. Backreaction effects generally increase the value of the free energy at the minimum for a given temperature, meaning that it will be easier to make the transition: the critical temperature when the value of the free energy is equal to the minimum of the zero temperature effective dilaton potential is higher.

Phase Transitions
The phase transition connecting the finite temperature approximate CFT to the zero temperature soft-wall description dual to spontaneously broken near-conformal symmetry is first order, and thus proceeds via bubble nucleation [8,34,35]. In the early universe, bubble nucleation competes with Hubble expansion, with the phase transition proceeding to completion only if bubble creation outpaces Hubble dilution [36,37].
This requirement can be phrased as the necessity of having one or more bubbles in one unit of Hubble 4-volume. If there is a rate of bubble nucleation per unit volume given by Γ/V , then model parameters that satisfy at some time period in the early universe will correspond to a successful phase transition.
In the hot conformal phase, the universe is radiation dominated, and we have The decay rate is proportional to Γ/V ∝ e −S E , where S E is the euclidean action associated with the fields evolving from their initial to their final values at nucleation. The coefficient of proportionality is difficult to calculate in general, involving an "obdurate" functional determinant, but dimensional analysis says that this factor should be of order f 4 . Similarly, the nucleation temperature is not expected to be drastically different from f , and later analysis in this section confirms this, with results shown in Figure 7. Putting this together, the criterion for bubble nucleation is roughly Many terms have been left out that are subdominant for large M Pl /f hierarchies. The large gap between f and M Pl along with the exponential sensitivity to S E justifies using simple dimensional analysis on the coefficient of the rate.
The Euclidean action is associated with the path in field space from the initial state far away from the bubble to the final state associated with its interior. At finite temperature, the geometry is compactified along the time direction, with compactification radius given by the inverse temperature: t ∈ (0, 1/T ]. Two types of bubbles are possible -those with O(3) symmetry whose radius is large in comparison with 1/T , and those with small radius, where bubbles exhibit invariance under the full O(4). The one with lower action is the one that will determine the decay rate.
At low temperatures, or at the interior of a bubble in the cooling universe, there exists a warped extra dimensional solution with time compactified and no black hole horizon. At high temperatures, or far away from the bubble, a near-AdS-Schwarzchild solution (AdS-S) with horizon is seen, as discussed in the previous section. To determine the rate at which the phase transition proceeds, one requires the action associated with moving from one phase to the other.
Comparison of the minimum total free energy for each geometry identifies the preferred vacuum at various temperatures, and the structure of the potential that interpolates between the AdS-S minimum and the conformal breaking minimum specifies the dynamics that interpolates between the two phases. Our analysis parallels earlier work on conformal phase transitions [7,8], but with emphasis on the particulars of the soft-wall light dilaton construction.

Nucleation of soft-wall bubbles
The dynamics of the phase transition interpolate between the black hole solution and the zero-temperature soft-wall dilaton geometry. We assume that Hubble expansion adiabatically cools the finite temperature solution, with the black hole horizon position tracking the minimum of the free energy. Eventually, when the temperature is such that the criteria for tunneling are satisfied, bubbles in the black hole horizon form, with the interior of the bubbles containing the brane that cuts off the zero temperature geometry in the infrared region.
The criteria for tunneling are two-fold. First, the transition must be energetically favorable, with the minimum of the black hole solution free energy being greater than the minimum of the zero-temperature effective dilaton potential. This defines the critical temperature, T c . However, the rate of nucleation may not yet be high enough to overcome Hubble dilution, which is the second criterion. It is only when the bubble action reaches the critical value in Eq. (4.43) that bubbles begin to nucleate. The temperature associated with the critical action crossing is the nucleation temperature, T n , but the action is minimum not for f at the bottom of the effective potential, but rather for smaller values of f . We denote the value of f inside a bubble as the nucleation scale, f n . After the phase transition is completed, the dilaton will oscillate and decay down to the true minimum.
The full action interpolating between the black hole and stabilized dilaton solutions is not accessible in this calculation without a UV completion, as the black hole solution at large y h and the zero temperature small f regions both involve trans-planckian excursions of the bulk curvature and scalar field, and also it is not clear how to properly normalize fluctuations in the position of the black hole horizon. However one can estimate the bubble action in several hypothesis scenarios that depend on the size of the bubble radius relative to the inverse temperature, and the maximum size of the gradient of the fields as one moves from the interior to the exterior of a bubble. We presume that the contribution to the bubble action from evolution on the black hole side of the transition is small, and the finite temperature calculation serves only to give the proper nucleation temperature. The size of the bubble determines whether the bubble has O(3) or O(4) symmetry in the Euclidean action, and the steepness of the bubble wall determines whether a thick [38] or thin [39] wall approximation is a better estimate for the minimum action.
Thick walls typically dominate when the latent heat associated with the phase transition is comparable to the barrier height separating the minima. While the barrier height cannot be calculated due to loss of control of the theory in the small T h , small f regions, it appears in numerical simulation that the trend is to maintain a shallow potential. We The Euclidean action associated with a bubble during a finite temperature phase transition is calculated on a geometry where the time coordinate is compactified on a circle of radius 1/T . Small bubbles with R < 1/T exhibit O(4) symmetry, and the action above reduces to a radial integral along which f varies from its nucleation value out to zero at the boundary of the bubble: For the larger O(3) symmetric bubbles which wrap the time direction, one obtains N is a normalization factor associated with canonically normalizing the fluctuations of the soft wall dilaton.
The bubble action can be approximated by [38] where ∆R is the region where f is changing significantly, andV is the volume averaged value of the potential inside the bubble. For the case of the dilaton potential, which is typically shallow for small values of , the thick wall action is smaller than the thin wall, and so the phase transition is driven by thick wall bubbles with f varying throughout. Thus, we can take ∆R to be the same as R, and minimize the bubble action over R. This yields: The normalization factor N we take to be N = 3N 2 /2π 2 . This is likely larger than the actual normalization, which is affected by backreaction, and thus we expect our values of the action to be conservatively large. For the average value of the potential, we use the difference between the finite temperature potential minimum and the value of the soft-wall potential at the nucleation value of the dilaton, f n : The properties of the finite temperature near-conformal phase transition can be calculated as a function of the input parameters in the 5D model. In Figure 5,we display in each panel the value of the minimum bubble action (both for O(3) and O(4) bubbles) as a function of the dimensionless IR brane tension Λ 1 . Each point on the Λ 1 axis has an associated value of f at the minimum of the zero temperature potential, and a value of the potential itself at the minimum. On additional horizontal axes on the top of each plot, we display the hierarchy between f and the 4D Planck scale as well as |V min |/f 4 for reference. In the figure, we show both "hard" and soft-wall values of φ 1 , with the hard wall calculation displaying excellent agreement with analytical calculations in the literature [8]. As N increases, the bubble action increases, making it more difficult for the phase transition to complete for a fixed ratio of f /M Planck .
We note that the phase transition completes over a much wider range of parameter space, which is a strong success of the soft-wall models. In the hard wall description, smaller values of N are typically necessary, and perturbativity of the 5D gravity theory is not guaranteed.
In Figure 6, we show the values for the bubble action as a function of the bulk scalar mass, . Values of Λ 1 are chosen so as to center the plots with = 0.1 corresponding to a ratio f /M Planck of order the TeV-Planck hierarchy.
In Figure 7, for the softwall scenario where v 1 = 10, we display the values of the temperature at which nucleation begins to occur for each scenario, as well as the value of f n that the system initially tunnels to. Both O(3) and O(4) bubbles are shown. Of note is the behavior of the nucleation temperature and nucleation value of the condensate. Relative to the value of f at the minimum, both of these quantities become small as the value of N is increased. Finally, in Figure 8, we display the bubble action and nucleation properties for the largest value of N found for which a hierarchy of TeV-Planck will complete for = −0.1. For this large value of N the values of T n /f and f n /f are quite small, 10 −4 and 10 −3 respectively.

Gravitational Waves
During a first order phase transition that proceeds by bubble nucleation, gravitational waves are sourced both by the collisions of the bubbles themselves, which break the spherical symmetry of the solutions discussed above, and also by turbulence in the finite temperature plasma as the bubbles move through it.
The stochastic spectrum is determined primarily by only coarse features of the phase transition. The latent heat difference determines a parameter α, which is the ratio of the latent heat compared with the energy density in the finite temperature false vacuum phase: Secondly, there is a parameter which describes the rate of variation of the bubble nucleation rate. This parameter, β, can be derived in terms of the variation of the bubble  action with respect to temperature at the bubble nucleation temperature: The two types of sources of gravitational waves generate spectra with different values of the frequency at the peak of the signal, and different power laws for the fall-off on the tails of the signal. The results of [40] are used for the purpose of calculating the characteristics of the signal. In Figure 10, we display the expected density spectrum of gravitational waves generated from the phase transition. The signal strength increases quickly with N , and as the nucleation temperature decreases with increasing N for fixed f , the higher values of N have a peak at lower frequency. Smaller values of N For N = 6 and N = 9 may be visible at LISA, and perhaps eLISA if our approximations are overly conservative. For N = 20, the frequency looks too small to be detectable at LISA, however proposed pulsar timing array experiments would probe this region with sufficient sensitivity.
Of course, the type of transition we discuss is not restricted to be at the TeV scale, although we would not be solving the electroweak hierarchy in this case. The scale f could instead be associated with other higher energy scales such as the GUT scale, or perhaps a Peccei-Quinn scale. In this case, the signal peak would be at higher frequency.  In this plot the signal data calculated in this work is overlaid on sensitivity curves derived from [41], and available at http://rhcole.com/apps/GWplotter/

Conclusions
This work explored solutions to 5D Einstein-Scalar theories with the goal of studying the cosmological phase transition of models with naturally large hierarchies of scale that are induced by geometric warping. In particular, we have examined scenarios where gravitational backreaction is large, yet the theory remains approximately conformally invariant. A numerical analysis of the finite and zero temperature potentials was performed in order to calculate the properties of the phase transition.
In early constructions of the Randall-Sundrum model, stabilized by a Goldberger-Wise mechanism with small gravitational backreaction, small values of N were required to accomplish a successful cosmological phase transition. Perturbativity of 5D gravity in these models with such small values of N is in doubt. In this work it was demonstrated that in the soft-wall construction where gravitational backreaction is taken into account, the bubble nucleation rate is faster for a given N , and that far larger values of N can be accommodated in the sense that there is a successful early universe first order phase transition. For the parameters we have considered, N = 20 is near the threshold beyond which bubbles will not nucleate.
The gravitational wave signals associated with a first order transition near the weak scale, generated by collisions and turbulance in the plasma, are potentially strong enough to be visible at next generation gravitational wave observatories such as LISA (and perhaps eLISA for larger N ) if the nucleation temperature corresponds to a peak in the spectrum not far from the targeted frequency range of these experiments. As the nucleation temperature can be far from the value of the condensate, f , the peak in the spectrum can vary significantly while holding the condensate value fixed, with this separation growing with increasing N . and if we then have then the superpotential equation is solved. Consistency of these two equations then gives an equation of motion for σ(φ): The solution for σ is a trivial constant if the bulk potential is constant, and for a small deformation, σ evolves slowly as a function of φ.
In terms of the superpotential, the dilaton effective potential has the following form where we can use the relation between φ and W to extract the hierarchy associated with the brane separation: The value of the condensate is then given by where y(φ) is given by the integral equation above.

B Appendix: Approximate Analytic Results
While a full numerical analysis was needed to draw accurate conclusions regarding the phase structure of the models under consideration, there are analytical expressions that give insight into the results. In this appendix, we present these approximate analytic solutions and visually compare them to full numerical results, and, for comparison, to solutions used in the literature for models that have only small backreaction.
To begin, we observe that for the case of a constant bulk potential, the equation of motion for φ reduces to¨ (B.10) It can be shown that starting in the UV and moving to the IR, this equation has asymptotic behavior whereφ ≈ 0 in the UV, andφ ≈ √ 12 in the IR.
For the case of a nontrivial bulk potential the equation of motion is For a bulk potential which depends only mildly on φ and is polynomial in φ, the potential term is always suppressed. For the case of a quadratic potential, V ( φ) = Λ(1 + 3 φ 2 ), when the backreaction and the value of φ are both small, the second derivative term goes like 2 , and can be ignored in a leading approximation. When φ is large, the term is again suppressed as the term goes like 1/ φ in this limit, and¨ φ is again especially small. The equation is then relatively simple in cases where the potential term is small: For the case of the quadratic potential, we have˙ φ − φ 1+ 3 φ 2 ≈ 0 which has an exact solution: φ(y) UV = 3 P L φ 2 UV 3 exp (6(y − y 0 ) + φ 2 UV ) (B.12) where we have imposed the boundary condition φ(y 0 ) UV ≡ φ U V and PL is the product log function. This UV solution can be contrasted with the lowest order UV solution by expanding the UV fixed point equation as˙ φ − φ = 0 and then solving to obtain φ(y) UV = φ UV exp( y). (B.13) One could also obtain this solution from expanding the correct leading-order result to lowest order in φ 2 UV . We will refer to this as the lowest order UV solution for φ, which is the solution for the Golberger-Wise field in the UV commonly found in the literature. While this solution is sufficient for the hardwall scenario as φ is always O(1), it is insufficient for a UV region that allows moderately larger values of φ since in that case φ 2 U V is not a small dimensionless quantity, and starts to dominate over the UV value of the 5D cosmological constant term. A comparison between the lowest order solution and our UV solution, along with some exact numerical results are displayed in Figure 11.
We further note that the solution for the scalar field in the finite temperature scenario where there is a black hole horizon also allows an approximate analytic solution. The near horizon limit of the φ equation of motion produces the boundary conditioṅ  Figure 11: Plots for φ and its derivative. Green is the lowest order UV solution, φ 0 e y , red is the analytic UV solution we are presenting, and blue are various full numerical results for the solution parameterized by ∆ φ IR ≡ √ 12 − φ (y IR ). The agreement in the deep UV is good for all cases, as expected, but the lowest order UV solution disagrees severely. Our approximate UV solution outlined above continues to track the unstable UV fixed point up until the IR condensate region develops. but this is the same relation employed as an approximation in the UV region of the zero temperature solution. This means that in order to accommodate a horizon, φ never moves far from the approximate analytic solution studied above. A comparison between φ in the presence of a horizon function and with no horizon function are plotted in figure 12.  Figure 12: A comparison between the φ solution in the presence of a horizon and then with no horizon (numerical and analytical). Since the boundary condition for φ at the horizon is the equation of motion itself to leading order, a horizon can be placed at any point along this trajectory which nominally effects φ numerically, and only at next-to-leading-order analytically.