We study the dynamics of cosmological phase transitions in the case of small velocities of bubble walls, $v_w<0.1$. We discuss the conditions in which this scenario arises in a physical model, and we compute the development of the phase transition. We consider different kinds of approximations and refinements for relevant aspects of the dynamics, such as the dependence of the wall velocity on hydrodynamics, the distribution of the latent heat, and the variation of the nucleation rate. Although in this case the common simplifications of a constant wall velocity and an exponential nucleation rate break down due to reheating, we show that a delta-function rate and a velocity which depends linearly on the temperature give a good description of the dynamics and allow to solve the evolution analytically. We also consider a Gaussian nucleation rate, which gives a more precise result for the bubble size distribution. We discuss the implications for the computation of cosmic remnants.

Strongly first-order phase transitions, i.e., those with a large order parameter, are characterized by a considerable supercooling and high velocities of phase transition fronts. A very strong phase transition may have important cosmological consequences due to the departures from equilibrium caused in the plasma. In general, there is a limit to the strength, since the metastability of the old phase may prevent the transition to complete. Near this limit, the bubble nucleation rate achieves a maximum and thus departs from the widely assumed behavior in which it grows exponentially with time. We study the dynamics of this kind of phase transitions. We show that in some cases a gaussian approximation for the nucleation rate is more suitable, and in such a case we solve analytically the evolution of the phase transition. We compare the gaussian and exponential approximations with realistic cases and we determine their ranges of validity. We also discuss the implications for cosmic remnants such as gravitational waves.

We investigate the production of a stochastic background of gravitational waves in the electroweak phase transition. We consider extensions of the Standard Model which can give very strongly first-order phase transitions, such that the transition fronts either propagate as detonations or run away. To compute the bubble wall velocity, we estimate the friction with the plasma and take into account the hydrodynamics. We track the development of the phase transition up to the percolation time, and we calculate the gravitational wave spectrum generated by bubble collisions, magnetohydrodynamic turbulence, and sound waves. For the kinds of models we consider, we find parameter regions for which the gravitational waves are potentially observable at the planned space-based interferometer eLISA. In such cases, the signal from sound waves is generally dominant, while that from bubble collisions is the least significant of them. Since the sound waves and turbulence mechanisms are diminished for runaway walls, the models with the best prospects of detection at eLISA are those which do not have such solutions. In particular, we find that heavy extra bosons provide stronger gravitational wave signals than tree-level terms.

In cosmological first-order phase transitions, gravitational waves are generated by the collisions of bubble walls and by the bulk motions caused in the fluid. A sizeable signal may result from fast-moving walls. In this work we study the hydrodynamics associated to the fastest propagation modes, namely, ultra-relativistic detonations and runaway solutions. We compute the energy injected by the phase transition into the fluid and the energy which accumulates in the bubble walls. We provide analytic approximations and fits as functions of the net force acting on the wall, which can be readily evaluated for specific models. We also study the back-reaction of hydrodynamics on the wall motion, and we discuss the extrapolation of the friction force away from the ultra-relativistic limit. We use these results to estimate the gravitational wave signal from detonations and runaway walls.

The subsonic expansion of bubbles in a strongly first-order electroweak phase transition is a convenient scenario for electroweak baryogenesis. For most extensions of the Standard Model, stationary subsonic solutions (i.e., deflagrations) exist for the propagation of phase transition fronts. However, deflagrations are known to be hydrodynamically unstable for wall velocities below a certain critical value. We calculate this critical velocity for several extensions of the Standard Model and compare with an estimation of the wall velocity. In general, we find a region in parameter space which gives stable deflagrations as well as favorable conditions for electroweak baryogenesis.

The growth of bubbles in cosmological first-order phase transitions involves nontrivial hydrodynamics. For that reason, the study of the propagation of phase transition fronts often requires several approximations. A frequently used approximation consists in describing the two phases as being composed only of radiation and vacuum energy (the so-called bag equation of state). We show that, in realistic models, the speed of sound in the low-temperature phase is generally smaller than that of radiation, and we study the hydrodynamics in such a situation. We find in particular that a new kind of hydrodynamical solution may be possible, which does not arise in the bag model. We obtain analytic results for the efficiency of the transfer of latent heat to bulk motions of the plasma, as a function of the speed of sound in each phase.

In a cosmological first-order phase transition, bubbles of the stable phase nucleate and expand in the supercooled metastable phase. In many cases, the growth of bubbles reaches a stationary state, with bubble walls propagating as detonations or deflagrations. However, these hydrodynamical solutions may be unstable under corrugation of the interface. Such instability may drastically alter some of the cosmological consequences of the phase transition. Here, we study the hydrodynamical stability of deflagration fronts. We improve upon previous studies by making a more careful and detailed analysis. In particular, we take into account the fact that the equation of motion for the phase interface depends separately on the temperature and fluid velocity on each side of the wall. Fluid variables on each side of the wall are similar for weakly first-order phase transitions, but differ significantly for stronger phase transitions. As a consequence, we find that, for large enough supercooling, any subsonic wall velocity becomes unstable. Moreover, as the velocity approaches the speed of sound, perturbations become unstable on all wavelengths. For smaller supercooling and small wall velocities, our results agree with those of previous works. Essentially, perturbations on large wavelengths are unstable, unless the wall velocity is higher than a critical value. We also find a previously unobserved range of marginally unstable wavelengths. We analyze the dynamical relevance of the instabilities, and we estimate the characteristic time and length scales associated to their growth. We discuss the implications for the electroweak phase transition and its cosmological consequences.

The steady state propagation of a phase transition front is classified, according to hydrodynamics, as a deflagration or a detonation, depending on its velocity with respect to the fluid. These propagation modes are further divided into three types, namely, weak, Jouguet, and strong solutions, according to their disturbance of the fluid. However, some of these hydrodynamic modes will not be realized in a phase transition. One particular cause is the presence of instabilities. In this work we study the linear stability of weak detonations, which are generally believed to be stable. After discussing in detail the weak detonation solution, we consider small perturbations of the interface and the fluid configuration. When the balance between the driving and friction forces is taken into account, it turns out that there are actually two different kinds of weak detonations, which behave very differently as functions of the parameters. We show that the branch of stronger weak detonations are unstable, except very close to the Jouguet point, where our approach breaks down.

In cosmological first-order phase transitions, the microscopic interaction of the phase transition fronts with non-equilibrium plasma particles manifests itself macroscopically as friction forces. In general, it is a nontrivial problem to compute these forces, and only two limits have been studied, namely, that of very slow walls and, more recently, ultra-relativistic walls which run away. In this paper we consider ultra-relativistic velocities and show that stationary solutions still exist when the parameters allow the existence of runaway walls. Hence, we discuss the necessary and sufficient conditions for the fronts to actually run away. We also propose a phenomenological model for the friction, which interpolates between the non-relativistic and ultra-relativistic values. Thus, the friction depends on two friction coefficients which can be calculated for specific models. We then study the velocity of phase transition fronts as a function of the friction parameters, the thermodynamic parameters, and the amount of supercooling.

We study the generation of gravitational waves in the electroweak phase transition. We consider a few extensions of the Standard Model, namely, the addition of scalar singlets, the minimal supersymmetric extension, and the addition of TeV fermions. For each model we consider the complete dynamics of the phase transition. In particular, we estimate the friction force acting on bubble walls, and we take into account the fact that they can propagate either as detonations or as deflagrations preceded by shock fronts, or they can run away. We compute the peak frequency and peak intensity of the gravitational radiation generated by bubble collisions and turbulence. We discuss the detectability by proposed spaceborne detectors. For the models we considered, runaway walls require significant fine tuning of the parameters, and the gravitational wave signal from bubble collisions is generally much weaker than that from turbulence. Although the predicted signal is in most cases rather low for the sensitivity of LISA, models with strongly coupled extra scalars reach this sensitivity for frequencies $f\sim 10^{-4}\,\mathrm{Hz}$, and give intensities as high as $h^2\Omega_{\mathrm{GW}}\sim 10^{-8}$.

We consider the motion of planar phase-transition fronts in first-order phase transitions of the Universe. We find the steady state wall velocity as a function of a friction coefficient and thermodynamical parameters, taking into account the different hydrodynamic modes of propagation. We obtain analytical approximations for the velocity by using the thin wall approximation and the bag equation of state. We compare our results to those of numerical calculations and discuss the range of validity of the approximations. We analyze the structure of the stationary solutions. Multiple solutions may exist for a given set of parameters, even after discarding non-physical ones. We discuss which of these will be realized in the phase transition as the stationary wall velocity. Finally, we discuss on the saturation of the friction at ultra-relativistic velocities and the existence of runaway solutions.

The cosmological remnants of a first-order phase transition generally depend on the perturbations that the walls of expanding bubbles originate in the plasma. Several of the formation mechanisms occur when bubbles collide and lose their spherical symmetry. However, spherical bubbles are often considered in the literature, in particular for the calculation of gravitational waves. We study the steady state motion of bubble walls for different bubble symmetries. Using the bag equation of state, we discuss the propagation of phase transition fronts as detonations and subsonic or supersonic deflagrations. We consider the cases of spherical, cylindrical and planar walls, and compare the energy transferred to bulk motions of the relativistic fluid. We find that the different wall geometries give similar perturbations of the plasma. For the case of planar walls, we obtain analytical expressions for the kinetic energy in the bulk motions. As an application, we discuss the generation of gravitational waves.

We study the velocity of bubble walls in the electroweak phase transition. For several extensions of the Standard Model, we estimate the friction and calculate the wall velocity, taking into account the hydrodynamics. We find that deflagrations are generally more likely than detonations. Nevertheless, for models with extra bosons, which give a strongly first-order phase transition, the deflagration velocity is in general quite high, $0.1\lesssim v_w\lesssim 0.6$. Therefore, such phase transitions may produce an important signal of gravitational waves. On the other hand, models with extra fermions which are strongly coupled to the Higgs boson may provide a strongly first-order phase transition and small velocities, $10^{-2}\lesssim v_w\lesssim 10^{-1}$, as required by electroweak baryogenesis.

We study the steady state motion of bubble walls in cosmological phase transitions. Taking into account the boundary and continuity conditions for the fluid variables, we calculate numerically the wall velocity as a function of the nucleation temperature, the latent heat, and a friction parameter. We determine regions in the space of these parameters in which detonations and/or deflagrations are allowed. In order to apply the results to a physical case, we calculate these quantities in a specific model, which consists of an extension of the Standard Model with singlet scalar fields. We also obtain analytic approximations for the wall velocity, both in the case of deflagrations and of detonations.

The walls of bubbles in a first-order phase transition can propagate either as detonations, with a velocity larger than the speed of sound, or deflagrations, which are subsonic. We calculate the gravitational radiation that is produced by turbulence during a phase transition which develops via deflagration bubbles. We take into account the fact that a deflagration wall is preceded by a shock front which distributes the latent heat throughout space and influences other bubbles. We show that turbulence can induce peak values of $\Omega_{GW}$ as high as $\sim 10^{-9}$. We discuss the possibility of detecting at LISA gravitational waves produced in the electroweak phase transition with wall velocities $v_w\lesssim 10^{-1}$, which favor electroweak baryogenesis.

Cosmological phase transitions are predicted by Particle Physics models, and have a variety of important cosmological consequences, which depend strongly on the dynamics of the transition. In this work we investigate in detail the general features of the development of a first-order phase transition. We find thermodynamical constraints on some quantities that determine the dynamics, namely, the latent heat, the radiation energy density and the false-vacuum energy density. Using a simple model with a Higgs field, we study numerically the amount and duration of supercooling and the subsequent reheating and phase coexistence. We analyze the dependence of the dynamics on the different parameters of the model, namely, the energy scale, the number of degrees of freedom and the couplings of the scalar field with bosons and fermions. We also inspect the implications for the cosmological outcomes of the phase transition.

We suggest that the current acceleration of the universe may be explained by the vacuum energy of a hidden sector which is stuck in a state of equilibrium between phases. The phases are associated to a late-time first-order phase transition, where phase coexistence originates at a temperature $T_c \sim 10^{-3}eV$ and lasts until temperature falls below $T\sim 10^{-4}eV$. During phase coexistence, the energy density has an effective cosmological constant component with the observed magnitude. This scenario does not require supercooling and may arise naturally in realistic models.

It is known that baryon number inhomogeneities may arise as a consequence of electroweak baryogenesis. Their geometry, size, and amplitude depend on the parameters that characterize the baryogenesis mechanism, as well as on those that determine the phase transition dynamics. We investigate this parametric dependance. We show that in the case of the minimal supersymmetric standard model, the geometry of the inhomogeneities most probably consists of spherical regions of high density surrounded by low-density walls, in contrast to the case of the minimal standard model. In this supersymmetric extension we find that density contrasts of up to a factor of 100 may arise. This amplitude increases for higher values of the latent heat or lower values of the bubble wall tension, and can be significantly larger in different extensions of the standard model. Such inhomogeneities may thus affect the dynamics of the subsequent quark-hadron phase transition.

We consider first-order phase transitions of the Universe in the radiation-dominated era. We argue that in general the velocity of interfaces is non-relativistic due to the interaction with the plasma and the release of latent heat. We study the general evolution of such slow phase transitions, which comprise essentially a short reheating stage and a longer phase equilibrium stage. We perform a completely analytical description of both stages. Some rough approximations are needed for the first stage, due to the non-trivial relations between the quantities that determine the variation of temperature with time. The second stage, instead, is considerably simplified by the fact that it develops at a constant temperature, close to the critical one. Indeed, in this case the equations can be solved exactly, including back-reaction on the expansion of the Universe. This treatment also applies to phase transitions mediated by impurities. We also investigate the relations between the different parameters that govern the characteristics of the phase transition and its cosmological consequences, and discuss the dependence of these parameters with the particle content of the theory.

Nucleation of branes by a four-form field has recently been considered in string motivated scenarios for the neutralization of the cosmological constant. An interesting question in this context is whether the nucleation of stacks of coincident branes is possible, and if so, at what rate does it proceed. Feng et al. have suggested that, at high ambient de Sitter temperature, the rate may be strongly enhanced, due to large degeneracy factors associated with the number of light species living on the worldsheet. This might facilitate the quick relaxation from a large effective cosmological constant down to the observed value. Here, we analyse this possibility in some detail. In four dimensions, and after the moduli are stabilized, branes interact via repulsive long range forces. Because of that, the Coleman-de Luccia (CdL) instanton for coincident brane nucleation may not exist, unless there is some short range interaction which keeps the branes together. If the CdL instanton exists, we find that the degeneracy factor depends only mildly on the ambient de Sitter temperature, and does not switch off even in the case of tunneling from flat space. This would result in catastrophic decay of the present vacuum. If, on the contrary, the CdL instanton does not exist, coindident brane nucleation may still proceed through a "static" instanton, representing pair creation of critical bubbles -- a process somewhat analogous to thermal activation in flat space. In that case, the branes may stick together due to thermal symmetry restoration, and the pair creation rate depends exponentially on the ambient de Sitter temperature, switching off sharply as the temperature approaches zero. Such static instanton may be well suited for the "saltatory" relaxation scenario proposed by Feng et al.

The latent heat released during the expansion of bubbles in the electroweak phase transition reheats the plasma and causes the bubble growth to slow down. This decrease of the bubble wall velocity affects the result of electroweak baryogenesis. Since the efficiency of baryogenesis peaks for a wall velocity $\sim 10^{-2}$, the resulting baryon asymmetry can either be enhanced or suppressed, depending on the initial value of the wall velocity. We calculate the evolution of the phase transition taking into account the release of latent heat. We find that, although in the SM the baryon production is enhanced by this effect, in the MSSM it causes a suppression to the final baryon asymmetry.

We investigate the evolution of the electroweak phase transition, using a one-Higgs effective potential that can be regarded as an approximation for the Minimal Supersymmetric Standard Model. The phase transition occurs in a small interval around a temperature T_t below the critical one. We calculate this temperature as a function of the parameters of the potential and of a damping coefficient related to the viscosity of the plasma. The parameters that are relevant for baryogenesis, such as the velocity and thickness of the walls of bubbles and the value of the Higgs field inside them, change significantly in the range of temperatures where the first-order phase transition can occur. However, we find that in the likely interval for T_t there is no significant variation of these parameters. Furthermore, the temperature T_t is in general not far below the temperature at which bubbles begin to nucleate.

The broken-symmetry electroweak vacuum is destabilized in the presence of a magnetic field stronger than a critical value. Such magnetic field may be generated in the phase transition and restore the symmetry inside the bubbles. A numerical calculation indicates that the first-order phase transition is delayed but may be completed for a sufficient low value of the Higgs mass unless the magnetic field is extremely high.

We discuss the stability of an electroweak string with axions in its core, which give to the configuration a quasi-topological property, and compare it with other modifications using instantons in the thin wall approximation.

We describe different electroweak strings with axionic content, including non-topological configurations calculated numerically, and show their possible influence on baryogenesis indicating that they may constitute a mechanism competitive to that of bubble nucleation with two Higgs-doublets.