Some Cosmological Implications of Hidden Sectors
We discuss some cosmological implications of extensions of the Standard Model with hidden sector scalars coupled to the Higgs boson. We put special emphasis on the conformal case, in which the electroweak symmetry is broken radiatively with a Higgs mass above the experimental limit. Our refined analysis of the electroweak phase transition in this kind of models strengthens the prediction of a strongly first-order phase transition as required by electroweak baryogenesis. We further study gravitational wave production and the possibility of low-scale inflation as well as a viable dark matter candidate.
The Standard Model of particle physics (SM) is nowadays considered to be an effective theory valid only up to a certain physical cutoff scale. Even though there exist a large variety of extensions of the SM, models with a hidden sector have recently attracted some attention. We will consider models with additional scalar fields that might transform non-trivially under a hidden gauge group but which are singlets under the SM gauge group. The only renormalizable interaction of such scalars with the SM occurs via the Higgs sector, which in this case serves as a portal to the hidden sector .
In this paper we are concerned with some of the possible cosmological implications of hidden sector extensions of the SM. This is a continuation of the study of the electroweak breaking and phase transition presented in Ref.  and we will provide some technical details that were omitted there. In addition we will present an analysis of other cosmological implications, namely gravitational wave production and dark matter abundance. We also comment on the possibility of low-scale inflation and present a calculation of the bubble wall velocity in case of a first-order electroweak phase transition. As in Ref.  we pay special attention to the classically conformal case which, for a strong coupling between the hidden sector scalars and the Higgs field, can be consistent with the mass bounds on the Higgs particle.
Ii The Model
ii.1 Zero Temperature Potential
We consider a set of real scalar fields coupled to the SM Higgs doublet with the tree level potential
We assume there are no linear or cubic terms in the hidden-sector scalar fields [this can be enforced by some global symmetry, e.g. ]. Besides, we assume that the squared masses of the hidden scalars are semi-positive definite (), such that this global symmetry remains unbroken and no quartic terms are necessary to stabilize the potential.
In the presence of a background Higgs field, , the one-loop effective potential in Landau gauge and scheme is then given by
The subscript denotes the gauge bosons ( and ), top quark, Higgs boson, Goldstone bosons ( and ) and hidden-sector scalar fields with while for gauge bosons and for fermions and scalars. The -dependent tree level masses are
where and denote the SM gauge couplings and the top quark Yukawa coupling.
As it was mentioned in the introduction the case with classical conformal invariance (i.e. and ) is especially interesting. In this situation all masses are proportional to the Higgs vacuum expectation value (VEV) and no dimensionful parameters enter into the tree level potential. However conformal invariance is broken by loop corrections as can be seen in Eq. (2) by the occurrence of the renormalization scale . In this way a mass scale is introduced via dimensional transmutation ; . Notice that an important difference with respect to the pure -theory is that, in the interesting region of the parameter space, the loop contributions are dominated by the hidden-sector scalar (and top) particles. Hence, it is not mandatory to improve the one-loop potential by renormalization group techniques (unlike in the -theory) and moreover the Goldstone and Higgs one-loop contributions to the potential can be safely neglected ; . In the classically conformal case the correct VEV follows from the minimization condition
where GeV is the observed Higgs VEV of the SM and should be chosen near . The potential then reads as
where is the one-loop Higgs mass given by
One can see that the occurrence of a sizable number of hidden-sector
scalars, rather strongly coupled to the Higgs field, can lead to
a Higgs mass above the LEP bound, even if the theory is classically
Given the fact that the dramatic impact on electroweak symmetry breaking we find is due to a sizable number of scalars somewhat strongly coupled to the Higgs, one might worry about the stability of the results when higher-order corrections to the potential are included. It is straightforward to obtain the dominant two-loop radiative corrections to the Higgs potential (those that depend on the top Yukawa coupling and ) by using standard techniques, as e.g. those used in Ref. . We have found that these two-loop effects never modify the structure of the potential in a qualitative way.
Finally, we would like to comment on the influence of the hidden-sector scalars on the cubic Higgs self-coupling. In Ref.  the claim was made that a strong phase transition often would lead to a deviation of the cubic Higgs coupling from its SM value. Taking only into account the top and hidden-sector scalar contributions, one obtains
which will be correlated with the strength of the phase transition in a later section.
ii.2 Finite Temperature Potential
In order to study the electroweak phase transition of the model, we consider the one-loop potential at finite temperature including the resummed Daisy-diagrams. The corresponding contributions are given by
where the holds for fermions (bosons) and are the thermal masses of the different bosonic species. Neglecting small contributions they read
Besides, in the resummed Daisy diagrams only the longitudinal polarizations of the gauge bosons contribute.
Iii Cosmology of Hidden-Sector Scalars
In this section we discuss cosmological implications of the hidden-sector scalar extensions of the SM. Namely, the electroweak phase transition, low-scale inflation, the bubble wall velocity during a first-order phase transition, gravitational wave production and dark matter are analyzed.
iii.1 Electroweak Phase Transition
In order to study the electroweak phase transition, we determine the so-called bounce solution of the three-dimensional Euclidean action that quantifies the tunneling probability in case of a first-order phase transition ; ; .
At finite temperature the bounce solution is obtained by extremizing the action
(where is the radial distance from the center of the bubble) with solutions obeying the boundary conditions
In addition it is understood that the bounce solution starts (at ) close to the global minimum of the potential (the broken phase of the Higgs vacuum).
The tunneling rate per unit volume and time element is approximately given by 
with , such that the average number of bubble nucleations per Hubble volume is given by
where the Hubble parameter is given by
is the effective number of degrees of freedom and GeV is the Planck mass.
Tunneling becomes in principle possible below the temperature at which the two minima of the potential are degenerate, but for almost degenerate vacua, the tunneling rate is still too small to start the phase transition. We define the temperature at which the phase transition starts by the average occurrence of one bubble per Hubble volume
The first nucleation of bubbles will hence approximately take place when
In order to characterize the end of the phase transition the fraction of space that is covered by bubbles can be used. Neglecting overlapping bubbles this is given by
Here is the velocity of the bubble wall and we define the end of the phase transition by
In order to quantify the strength of the phase transition we determine several quantities. These are evaluated at the end of the phase transition, when most cosmological processes such as baryogenesis and gravitational wave production take place. The first quantity is the ratio between the Higgs VEV and the temperature, . This ratio is important for baryogenesis, since suppression of washout effects by sphalerons  requires in the Standard Model. We do not expect this bound to be much different in the present model, since the sphaleron energy is dominated by the contributions from the gauge field configurations excited in the sphaleron rather than the scalar ones . The second quantity is the duration of the phase transition , which is given by
The last quantity we are interested in is the latent heat
The latent heat is usually normalized to the energy density of the radiation in the plasma, through the dimensionless parameter
The quantities and , as well as the bubble velocity , are the key parameters that govern gravitational wave production (discussed in a later section).
For our numerical examples we take, as in Ref. , a number of scalars with universal couplings to the Higgs, , and no explicit mass terms, . The results for the electroweak phase transition parameters listed above, as functions of and for several values of the Higgs mass (consistent with electroweak breaking conditions), are plotted in Figs. 1 to 3. For small values of , the phase transition is SM-like and therefore it is of second-order or a cross-over. As expected, the phase transition is in general stronger for larger values of and smaller Higgs masses. The latent heat (as described by ) and the strength of the transition [as measured by ] are both quickly increasing with and larger for smaller (see Fig. 1). In the figures we mark the conformal case with a cross and we see that, even in this case, the model shows a first-order phase transition strong enough to allow for electroweak baryogenesis (see Fig. 3). To the right of that conformal point the Higgs potential of the model has a barrier separating the symmetric and broken phases even at . For too large values of this barrier becomes too high and tunneling by thermal fluctuations is not efficient to trigger the electroweak phase transition. Note how the time of the transition, , gets larger and larger with increasing . Eventually no thermal transition will occur beyond a critical point and one would get stuck in the symmetric minimum (see below).
Finally we point out that in the present model a strong first-order phase transition does not necessarily imply a very large deviation of the cubic Higgs coupling from its SM value. Independently of the value of the Higgs mass, a phase transition that is strong enough for the suppression of sphaleron processes, , is possible for deviations of the cubic coupling as small as , as can be seen from Fig. 4.
iii.2 Low-Scale Inflation
Every time a relatively strong first-order phase transition occurs during the evolution of the Universe the plasma undergoes a stage of large over-cooling. This means that the energy finally released as latent heat is large compared to the thermal energy stored in the plasma. In this context it is worthwhile to ask whether during the stage of over-cooling the expansion of the Universe is significantly accelerated due to the dominance of the vacuum energy, i.e. whether inflation occurs. This might be interesting in order to connect the predictions of inflation to low energy physics, but also for more exotic scenarios like cold electroweak baryogenesis ; .
However if inflation takes place at electroweak scales the problem on whether this scenario allows for a graceful exit arises. It is well known that a realistic first-order phase transition cannot proceed arbitrarily slow, since percolation requires the decay rate of the vacuum to become rather large at a certain temperature . This severely constrains the prospects of low-scale inflation in such scenarios. We will now analyze the two possible scenarios.
The first scenario is that inflation ends by thermal tunneling . In this case a substantial amount of inflation is hardly achieved as can be seen as follows. Suppose that nucleation takes place after e-folds of inflation. In this case the nucleation temperature, , will be very low compared to the temperature of degenerate vacua
The Higgs VEV will have to remain stuck in the symmetric phase down to very small temperatures and the energy density will be dominated by the vacuum energy, , rather than by the thermal energy of the plasma. Roughly speaking, the vacuum energy is related to the temperature of degenerate vacua, such that the Hubble parameter is
Imposing that the thermal decay rate at is larger than the Hubble rate we get the condition
In this regime of very small nucleation temperature the three-dimensional action (that increases with temperature) has therefore to be much smaller than the electroweak scale:
This requires the potential barrier at zero temperature to be very small and we are thus led to the parameter region near the conformal case (in the conformal limit the barrier and the three-dimensional tunnel action vanish). In particular, the parameter cannot be much larger than in the conformal case. However we know that near the conformal case the Higgs VEV does not get stuck at the origin: in fact the phase transition occurs already at temperatures of electroweak size, as can be seen in Fig. 2.
The fact that in the conformal case can be understood as follows: The potential difference between the symmetric minimum and the broken minimum is in this case given by
The comparison of Eq. (29) with the thermal contributions to the potential in Eq. (9) shows that the temperature where the two minima are degenerate is of order of the Higgs mass. At the same time, the potential barrier between the minima is absent at zero temperature, arising solely by temperature effects. Therefore one expects that the phase transition takes place with a nucleation temperature of order . In particular, the temperatures and increase with increasing Higgs mass. All this agrees with the numerical results as shown in Fig. 2.
The second scenario is that the minimum at the origin does not decay by thermal fluctuations but rather through vacuum (quantum) tunneling . Instead of Eq. (14), the tunneling rate for vacuum decay is given by
where denotes the four-dimensional tunnel action
In order to assure that the minimum at the origin does not decay by thermal fluctuations, we should be in a parameter region in which there is a sufficiently large barrier at zero . This occurs when is larger than its conformal value that can be read off (for a fixed Higgs mass) from Fig. 1. Numerically, one finds that in this large-barrier regime the quantum tunneling probability is suppressed (), such that quantum tunneling cannot provide a graceful exit.
This result was to be expected, since to obtain stability under thermal fluctuations, a barrier comparable to the difference in vacuum energy is mandatory. In such regime the tunneling actions can be calculated using the thin-wall approximation and one finds the scaling behavior
This shows that, for , quantum tunneling is typically as unlikely as tunneling by thermal fluctuations.
As a conclusion we see that, in general, low-scale inflation with only one field seems to require two amply separate scales as for example in the case discussed in Ref. .
To end this section we show in Fig. 5 the typical behavior of the tunneling actions, and , as functions of the temperature for GeV and two choices of . For , gets eventually below the critical nucleation value (horizontal line), and the electroweak phase transition takes place. For no satisfactory transition would occur as the strong suppression will hinder percolation.
iii.3 Bubble Wall Velocities
During a strongly first-order phase transition the wall velocity of the expanding nucleated bubbles is an important parameter. For example, the standard picture of electroweak baryogenesis is based on the diffusion of charge-asymmetries into the unbroken phase in front of the wall. This effect is strongly suppressed if the wall expansion is supersonic, making electroweak baryogenesis implausible. On the other hand, gravitational wave production requires large wall velocities. Hence, baryogenesis and gravitational wave production at electroweak scales seem to be mutually exclusive .
The wall velocity does not only depend on the thermo-dynamical characteristics of the phase transition, but also on the particle content of the plasma. In particular, bosonic degrees of freedom that are massless before the phase transition but become heavy due to a strong coupling to the Higgs VEV exert a strong friction force on the wall ; ; . In this way the presence of many hidden-sector scalars leads to subsonic wall velocities and the phase transition proceeds by deflagration. Following the arguments of Ref. , in the present case the friction is dominated by the scalars and is given by
where denotes the thickness of the bubble wall during the phase transition and is the squared Debye mass of the hidden-sector scalars.
The expansion of the bubbles is driven by the pressure produced by the latent heat, . If the friction forces are large, , the wall velocity can be estimated to be . Comparison with Eq. (33) shows that this is the case if the phase transition is weak, in the sense that .
In the opposite regime, in which friction effects from the plasma on the bubble wall are negligible, one expects the phase transition to proceed by detonation. In this case the bubble wall velocity can be determined by a self-consistent calculation that leads to supersonic wall velocities . The wall velocity is then approximately given by
where is the latent heat normalized to the energy density of the plasma as given in Eq. (24). In fact, this value for the wall velocity is only an upper bound, since a larger class of detonation solutions is known to exist , but we use nevertheless this formula in the analysis of gravitational wave production.
The results presented in Fig. 1 show that in principle both possibilities can occur in hidden-sector scalar extensions of the SM without significant tuning.
iii.4 Gravitational Waves
Another smoking-gun signal of a cosmological first-order phase transition is gravitational wave (GW) radiation. When the Higgs bubbles nucleate and expand, a portion of the latent heat is transformed into kinetic energy of the Higgs field and also into bulk motion of the plasma that follows the passing bubble wall profile. When the bubbles finally percolate and collide, this energy is partially released into gravitational waves ; ; ; . Surprisingly, the only parameters that enter into the analysis of the gravitational wave radiation by collisions are the latent heat normalized to the radiation energy , the wall velocity of the bubbles and the duration of the phase transition .
In principle, there might be additional mechanisms of GW production, as e.g. turbulence in the plasma  and/or magnetic fields . However, for very strong phase transitions, the peak frequency of the GW spectrum is shifted to lower frequencies and mostly the high frequency part of the GW spectrum lies in the range of best experimental sensitivity of the planned space-based experiments. The contributions from bubble collisions usually dominate for , such that at the frequency of best sensitivity of LISA or BBO, it suffices to consider the contributions from collisions. A more complete discussion of these issues can be found in Refs. ; ; .
In the following we summarize the formulas for GW production by bubble collisions as recently presented in Ref. . The peak frequency is given by
whereas the energy density in GWs amounts to
The efficiency factor indicates the fraction of latent heat that is transformed into bulk motion of the plasma and finally into gravitational waves. It is given by 
The best sensitivity of BBO (LISA) is at mHz ( mHz) expected to be (). Considering that the GW spectrum from collisions scales approximately as for large frequencies , one obtains for BBO a signal to sensitivity ratio
and for LISA a value that is smaller by four orders of magnitude.
Comparison with the parameters of the phase transition in Fig. 1 shows that, in the present model a signal that is detectable by BBO is feasible, if the parameter is rather close to the critical point beyond which no thermal tunneling occurs, requiring a tuning in at the percent level. For example, the parameters
and to a signal to sensitivity ratio of . On the other hand, no observable traces from the electroweak phase transition are expected at LISA in the present model.
iii.5 Dark Matter
In this section we investigate if the new scalar degrees of freedom constitute a viable dark matter candidate. For simplicity we consider only one hidden-sector scalar as the generalization to several scalars is straightforward. Singlet dark matter has already been extensively discussed in the literature (see e.g. Ref. ; ; ). Here we focus on two aspects. First, we discuss if the same scalar species might be responsible both for a strong phase transition and for dark matter (Ref.  addresses the same question in a different extension of the Higgs sector). Secondly, we focus on the classically conformal case, in which the scalar has no explicit mass term.
The scalar has to be stable to constitute a viable dark matter candidate. This is achieved by the choice of the potential in Eq. (1), since the scalars are protected from decay by a symmetry. In particular we assume that this symmetry is not spontaneously broken. Nevertheless the scalars annihilate and the particle density of the scalar obeys the Boltzmann equation 
where the equilibrium distribution is given by
and denotes the Hubble parameter as given in Eq. (16). Rescaling the distribution functions, , one obtains the equation
The contributions to from annihilation to pairs of Higgs-, W-, Z-bosons and SM-fermions are respectively given by
is the decay width of the Higgs particle in the SM and that of the hidden scalars. Finally, for leptons (quarks).
An approximate solution to this equation was given in Ref. . The scalar freeze out temperature is given by
and typically one finds . The final particle density is
At present, , the total energy density in scalars is
where denotes the critical energy density of the Universe at present.
The dependence of on the scalar mass for fixed coupling is plotted in Fig. 6. Notice that we only plotted the dark matter density for scalar masses that are larger than and hence correspond to a positive mass term in the Lagrangian (according to ).
Besides a logarithmic dependence on the freeze-out temperature, the dark matter density scales for large masses as . Notice that for most annihilation channels are enhanced and the scalar contribution to dark matter is suppressed. Finally, the annihilation cross-section drops considerably below the W-boson threshold, , since if the scalar is light it mostly annihilates into bottom/anti-bottom pairs, which is suppressed by the bottom-quark Yukawa coupling. This leads to an increase of the dark matter density below the W-boson threshold. Notice that taking temperature effects into account, one expects that the annihilation cross-section changes less drastically when the scalar mass is varied. In particular the enhancement close to the Higgs mass is expected to be less prominent. Likewise, the drop below the W-boson threshold proceeds in an interval of width .
Therefore, we see that there are two valid regimes of scalar dark matter. The first option is to increase the scalar mass term , while keeping the coupling fixed. However, even in the case of a rather large number of scalars , this requires scalar masses of order TeV and such scalars cannot be responsible for a strong phase transition. Alternatively, the scalar could be rather light, with , and weakly coupled, such that its annihilation is suppressed. Also in this case, the impact of the scalars on the phase transition is small.
Finally, consider a model without an explicit singlet mass term in the Lagrangian. In Fig. 7 the dark matter density is plotted as a function of for and for two different values of the Higgs mass. The predicted dark matter density typically surpasses the observed one below the W-boson threshold. Besides, in the case that the Higgs boson is lighter than two W-bosons, the resonant enhancement in the decay channel can lead to two additional viable values for the parameter that reproduce the observed dark matter density. Again, such weakly coupled scalars cannot increase the strength of the phase transition considerably. In particular, a classically conformal model requires several additional, strongly coupled scalars to surpass the current bounds on the Higgs mass, see Fig. 1.
In conclusion, if extra scalar degrees of freedom are responsible for a strong electroweak first-order phase transition, as well as for dark matter, it seems that either the coupling constants or the mass terms are non-universal. Scalar dark matter requires either a scalar with a rather large mass TeV, or a rather weak coupling . However both types of scalars cannot contribute significantly to the strength of the phase transition. Hence, the existence of both features in a universal scalar framework would require a very large number of scalars, which we estimate to be .
We have studied several cosmological implications of Standard Model extensions with hidden-sector scalars. In particular, we strengthen the results of  finding that in models with a moderate number of hidden-sector scalars, , the electroweak phase transition is generically of first-order as long as the Higgs mass is not much larger than the electroweak scale and the coupling to the hidden sector is substantial, . An interesting feature of the model is that this property persists even if the theory is classically conformal invariant and the electroweak scale is induced by dimensional transmutation. This was already emphasized in Ref. . We find that the phase transition is in a large portion of the parameter space strong enough to suppress the sphaleron process after the phase transition, as required by electroweak baryogenesis. Besides, we find that sizable production of gravitational radiation requires a tuning of the parameters at the percent level.
Besides a strong first-order phase transition, viable electroweak baryogenesis requires sizable CP violation. Electroweak baryogenesis in non-SUSY models typically utilizes a Higgs VEV that has a changing complex phase during the phase transition. One useful ingredient hence seems to be to complexify the present scalars and to allow for scalar VEVs, but still this would not induce a change in the complex phase of the Higgs VEV such that the introduction of a second Higgs doublet seems unavoidable. Alternatively, one can introduce an additional source of CP violation in the quark sector (see e.g. Ref. ) but undoubtedly CP violation arising from the hidden sector would be much more appealing in our model.
Concerning dark matter, we find that the abundance required by the concordance model can be provided by hidden-sector scalars in two different regimes. In the first, the hidden-sector scalars have moderate couplings but large masses TeV. In the second, the hidden-sector scalars are rather light, . In this case, the scalars cannot annihilate into W-bosons, which greatly enhances the dark matter abundance. Notice that this scenario is compatible with scalars that obtain their mass solely by electroweak symmetry breaking. Nevertheless, neither type of scalar can contribute significantly to the strength of the phase transition, such that a viable dark matter candidate cannot substantially improve the prospects of electroweak baryogenesis compared to the SM. Hence, a simultaneous solution of the dark matter and baryogenesis problems of the Standard Model close to electroweak scales either requires a large number of scalars (in which case we found ), or several types of scalars in the hidden sector with non-uniform masses and/or couplings to the Higgs sector.
J.R.E. thanks A. Casas, A. Riotto and G. Servant for discussions. Work supported in part by the European Commission under the European Union through the Marie Curie Research and Training Networks “Quest for Unification” (MRTN-CT-2004-503369) and “UniverseNet” (MRTN-CT-2006-035863); by a Comunidad de Madrid project (P-ESP-00346); and by CICYT, Spain, under contracts FPA 2007-60252 and FPA 2005-02211. T.K. is supported by the EU FP6 Marie Curie Research & Training Network ’UniverseNet’ (MRTN-CT-2006-035863).
- In this case it is obvious that the one-loop Higgs mass differs significantly from its tree level value . Similarly the Goldstone boson is massless at one-loop unlike the tree level result . This discrepancy between the tree-level and one-loop masses does also appear in the non-conformal case. It usually does not constitute a significant complication, since the loop contributions of the Higgs and Goldstone bosons are small.
- The four-dimensional bounce solution can be trusted as long as its radius is smaller than . Above some (small) temperature the three-dimensional bounce solution should be considered instead.
- R. Schabinger and J. D. Wells, “A minimal spontaneously broken hidden sector and its impact on Higgs boson physics at the Large Hadron Collider,” Phys. Rev. D 72 (2005) 093007 [hep-ph/0509209]. B. Patt and F. Wilczek, “Higgs-field portal into hidden sectors,” [hep-ph/0605188]. M. Bowen, Y. Cui and J. D. Wells, “Narrow trans-TeV Higgs bosons and H h h decays: Two LHC search paths for a hidden sector Higgs boson,” JHEP 0703 (2007) 036 [hep-ph/0701035].
- J. R. Espinosa and M. Quirós, “Novel effects in electroweak breaking from a hidden sector,” Phys. Rev. D 76 (2007) 076004 [hep-ph/0701145].
- S. R. Coleman and E. Weinberg, “Radiative Corrections As The Origin Of Spontaneous Symmetry Breaking,” Phys. Rev. D 7 (1973) 1888.
- R. Jackiw, “Functional evaluation of the effective potential,” Phys. Rev. D 9, 1686 (1974).
- See e.g. J. R. Espinosa and R. J. Zhang, “Complete two-loop dominant corrections to the mass of the lightest CP-even Higgs boson in the minimal supersymmetric standard model,” Nucl. Phys. B 586 (2000) 3 [hep-ph/0003246].
- A. Noble and M. Perelstein, “Higgs Self-Coupling as a Probe of Electroweak Phase Transition,” [hep-ph/0711.3018].
- S. R. Coleman, “The Fate Of The False Vacuum. 1. Semiclassical Theory,” Phys. Rev. D 15, 2929 (1977) [Erratum-ibid. D 16, 1248 (1977)].
- C. G. Callan and S. R. Coleman, “The Fate Of The False Vacuum. 2. First Quantum Corrections,” Phys. Rev. D 16, 1762 (1977).
- A. D. Linde, “Fate Of The False Vacuum At Finite Temperature: Theory And Applications,” Phys. Lett. B 100, 37 (1981); “Decay Of The False Vacuum At Finite Temperature,” Nucl. Phys. B 216 (1983) 421 [Erratum-ibid. B 223 (1983) 544].
- G. R. Farrar and M. E. Shaposhnikov, “Baryon Asymmetry Of The Universe In The Standard Electroweak Theory,” Phys. Rev. D 50 (1994) 774 [hep-ph/9305275].
- F. R. Klinkhamer and N. S. Manton, “A Saddle Point Solution In The Weinberg-Salam Theory,” Phys. Rev. D 30 (1984) 2212.
- J. García-Bellido, D. Y. Grigoriev, A. Kusenko and M. E. Shaposhnikov, “Non-equilibrium electroweak baryogenesis from preheating after inflation,” Phys. Rev. D 60, 123504 (1999) [hep-ph/9902449].
- A. Tranberg and J. Smit, “Baryon asymmetry from electroweak tachyonic preheating,” JHEP 0311 (2003) 016 [hep-ph/0310342].
- A. H. Guth and E. J. Weinberg, “Could The Universe Have Recovered From A Slow First Order Phase Transition?,” Nucl. Phys. B 212, 321 (1983).
- G. Nardini, M. Quirós and A. Wulzer, “A Confining Strong First-Order Electroweak Phase Transition,” JHEP 0709 (2007) 077 [hep-ph/0706.3388].
- S. J. Huber and T. Konstandin, “Production of Gravitational Waves in the nMSSM,” JCAP 0805 (2008) 017 [hep-ph/0709.2091].
- G. D. Moore and T. Prokopec, “Bubble wall velocity in a first order electroweak phase transition,” Phys. Rev. Lett. 75, 777 (1995) [hep-ph/9503296].
- G. D. Moore and T. Prokopec, “How fast can the wall move? A Study of the electroweak phase transition dynamics,” Phys. Rev. D 52, 7182 (1995) [hep-ph/9506475].
- G. D. Moore, “Electroweak bubble wall friction: Analytic results,” JHEP 0003, 006 (2000) [hep-ph/0001274].
- P. J. Steinhardt, “Relativistic Detonation Waves And Bubble Growth In False Vacuum Decay,” Phys. Rev. D 25 (1982) 2074.
- M. Laine, “Bubble Growth As A Detonation,” Phys. Rev. D 49 (1994) 3847 [hep-ph/9309242].
- M. Kamionkowski, A. Kosowsky and M. S. Turner, “Gravitational radiation from first order phase transitions,” Phys. Rev. D 49, 2837 (1994) [astro-ph/9310044].
- R. Apreda, M. Maggiore, A. Nicolis and A. Riotto, “Gravitational waves from electroweak phase transitions,” Nucl. Phys. B 631 (2002) 342 [gr-qc/0107033].
- C. Caprini, R. Durrer and G. Servant, “Gravitational wave generation from bubble collisions in first-order phase Phys. Rev. D 77 (2008) 124015 [astro-ph/0711.2593].
- S. J. Huber and T. Konstandin, “Gravitational Wave Production by Collisions: More Bubbles,” JCAP 0809, 022 (2008) [arXiv:0806.1828 [hep-ph]].
- A. Kosowsky, A. Mack and T. Kahniashvili, “Gravitational radiation from cosmological turbulence,” Phys. Rev. D 66 (2002) 024030 [astro-ph/0111483]. A. D. Dolgov, D. Grasso and A. Nicolis, “Relic backgrounds of gravitational waves from cosmic turbulence,” Phys. Rev. D 66, 103505 (2002) [astro-ph/0206461]. G. Gogoberidze, T. Kahniashvili and A. Kosowsky, “The spectrum of gravitational radiation from primordial turbulence,” Phys. Rev. D 76 (2007) 083002 [astro-ph/0705.1733]. T. Kahniashvili, G. Gogoberidze and B. Ratra, “Gravitational Radiation from Primordial Helical MHD Turbulence,” Phys. Rev. Lett. 100 (2008) 231301 [astro-ph/0802.3524]. A. Megevand, “Gravitational waves from deflagration bubbles in first-order phase transitions,” Phys. Rev. D 78 (2008) 084003 [arXiv:0804.0391 [astro-ph]].
- C. Caprini and R. Durrer, “Gravitational waves from stochastic relativistic sources: Primordial turbulence and magnetic fields,” Phys. Rev. D 74 (2006) 063521 [astro-ph/0603476]. E. Fenu and R. Durrer, “Interactions of cosmological gravity waves and magnetic fields,” [astro-ph/0809.1383].
- A. Nicolis, “Relic gravitational waves from colliding bubbles and cosmic turbulence,” Class. Quant. Grav. 21 (2004) L27 [gr-qc/0303084].
- C. Grojean and G. Servant, “Gravitational waves from phase transitions at the electroweak scale and beyond,” Phys. Rev. D 75 (2007) 043507 [hep-ph/0607107].
- V. Silveira and A. Zee, “Scalar Phantoms,” Phys. Lett. B 161 (1985) 136.
- J. McDonald, “Gauge Singlet Scalars as Cold Dark Matter,” Phys. Rev. D 50 (1994) 3637 [hep-ph/0702143].
- C. P. Burgess, M. Pospelov and T. ter Veldhuis, “The minimal model of nonbaryonic dark matter: A singlet scalar,” Nucl. Phys. B 619 (2001) 709 [hep-ph/0011335].
- T. Hambye and M. H. G. Tytgat, “Electroweak Symmetry Breaking induced by Dark Matter,” Phys. Lett. B 659 (2008) 651 [arXiv:0707.0633 [hep-ph]].
- B. W. Lee and S. Weinberg, “Cosmological lower bound on heavy-neutrino masses,” Phys. Rev. Lett. 39 (1977) 165.
- D. Bodeker, L. Fromme, S. J. Huber and M. Seniuch, “The baryon asymmetry in the standard model with a low cut-off,” JHEP 0502, 026 (2005) [arXiv:hep-ph/0412366].