Phase transition and gravitational wave phenomenology of scalar conformal extensions of the Standard Model

TL;DR

This paper studies how scalar conformal extensions of the Standard Model can produce strong first-order phase transitions and detectable gravitational wave signals, emphasizing the importance of model complexity.

Contribution

It analyzes phase transition dynamics and gravitational wave signals in scalar conformal extensions, highlighting the limitations of minimal models and the potential of next-to-minimal models.

Findings

Minimal models struggle with phase transition dynamics due to thermal corrections.

Next-to-minimal models achieve correct electroweak symmetry breaking.

Next-to-minimal models predict strong gravitational wave signals.

Abstract

Thermal corrections in classically conformal models typically induce a strong first-order electroweak phase transition, thereby resulting in a stochastic gravitational wave background that could be detectable at gravitational wave observatories. After reviewing the basics of classically conformal scenarios, in this paper we investigate the phase transition dynamics in a thermal environment and the related gravitational wave phenomenology within the framework of scalar conformal extensions of the Standard Model. We find that minimal extensions involving only one additional scalar field struggle to reproduce the correct phase transition dynamics once thermal corrections are accounted for. Next-to-minimal models, instead, yield the desired electroweak symmetry breaking and typically result in a very strong gravitational wave signal.

Figures

Tables

Table 1. Table 1: The values of the parameters adopted in the plot of fig.  4 .

${\lambda }_p$	$v_{\varphi }/$GeV	$T_n/$GeV	$T_{\ast }/$GeV	$T_c/$GeV	$\beta /H_{\ast }$
$1$	${10}^4$	$5.68$	$493$	$1940$	$23.5$
$2$	${10}^4$	$512$	$1200$	$2990$	$70.1$
$1$	${10}^9$	$6.76\times {10}^4$	$4.88\times {10}^7$	$1.94\times {10}^8$	$13.6$
$2$	${10}^9$	$2.11\times {10}^7$	$9.00\times {10}^7$	$2.99\times {10}^8$	$32.5$

Equations

$$V=\sum _{i,j,k,l=1}^n{\lambda }_{ijkl}{\varphi }_i{\varphi }_j{\varphi }_k{\varphi }_l+\sum _{k=1}^n\frac{g_k{M}_k^4}{64{\pi }^2}\log \left(\frac{M_k^2}{{\mu }^2}\right)+\delta V,$$

$$V$$

$$A=\sum _{k=1}^n\frac{g_k{W}_k^2}{64{\pi }^2}\log W_k^2,B=\sum _{k=1}^n\frac{g_k{W}_k^4}{32{\pi }^2}.$$

$$V=B{\varphi }^4\left(\log \left(\frac{\varphi }{v_{\varphi }}\right)-\frac{1}{4}\right),$$

$$V_T=\sum _{k=1}^n{J}_T(M_k)+\sum _{k\in \mathrm{bosons}}{D}_T(M_k,{\Pi }_k),$$

$$J_T(M_k)=g_{k}T\int \frac{{\mathrm{d}}^{3}p}{(2\pi {)}^3}\log (1\mp e^{-E/T}),$$

$$D_T(M_k,{\Pi }_k)=\frac{g_{k}T}{12\pi }\left(M_k^3-(M_k^2+{\Pi }_k(T){)}^{\frac{3}{2}}\right),$$

$$J_T(M_k\ll T)$$

$$c_k=\{\begin{array}{ll}1& \text{ (bosons)}\\ -1/2& \text{ (fermions)}\end{array}.$$

$$V_{eff.}=V+V_T=B{\varphi }^4\left(\log \left(\frac{\varphi }{v_{\varphi }}\right)-\frac{1}{4}\right)+C{T}^2{\varphi }^2,$$

$$V_T=J_T(M_{\sigma })+D_T(M_{\sigma },{\Pi }_{\sigma }),$$

$$\Gamma (T)\simeq T^4{\left(\frac{S_3}{2\pi T}\right)}^{\frac{3}{2}}\exp \left(-\frac{S_3}{T}\right),$$

$$S_3=4\pi \int r^2\mathrm{d}r\left(\frac{1}{2}{\left(\frac{\mathrm{d}\varphi }{\mathrm{d}r}\right)}^2+V_{eff.}(\varphi ,T)\right)$$

$$\frac{{\mathrm{d}}^2\varphi }{\mathrm{d}{r}^2}+\frac{2}{r}\frac{\mathrm{d}\varphi }{\mathrm{d}r}=\frac{\mathrm{d}{V}_{eff.}}{\mathrm{d}\varphi }$$

$$\frac{4\pi }{3}\frac{\Gamma (T_n)}{H(T_n{)}^4}\simeq 1.$$

$$\alpha (T)=\frac{1}{{\rho }_{\gamma }}\left(\Delta V_{eff.}-\frac{T}{4}\Delta \left(\frac{\mathrm{d}{V}_{eff.}}{\mathrm{d}T}\right)\right),$$

$$\alpha (T)\simeq \frac{V(0)-V(v_{\varphi })}{{\rho }_{\gamma }(T)}\simeq 6\times 10^{-4}\frac{{\lambda }_p^2{v}_{\varphi }^4}{g_{\ast }{T}^4},$$

$${\Omega }_{\mathrm{gw}}{h}^2=\frac{4.9\times 10^{-6}{\left(\frac{f}{f_{\mathrm{env}}}\right)}^{2.8}}{1+2.8{\left(\frac{f}{f_{\mathrm{env}}}\right)}^{3.8}}{\left(\frac{H_{\ast }}{\beta }\right)}^2{\left(\frac{100}{g_{\ast }}\right)}^{\frac{1}{3}}.$$

$$\frac{f_{\mathrm{env}}}{\mathrm{Hz}}=3.5\times 10^{-6}\left(\frac{\beta }{H_{\ast }}\right)\left(\frac{T_{\ast }}{100\mathrm{GeV}}\right){\left(\frac{g_{\ast }}{100}\right)}^{\frac{1}{6}}.$$

$$\frac{\beta }{H_{\ast }}=T_n{\frac{\mathrm{d}}{\mathrm{d}T}\frac{S_3}{T}|}_{T=T_n}.$$

$$V(H,s)={\lambda }_h(H^{†}H{)}^2+\frac{{\lambda }_{hs}}{2}(H^{†}H)s^2+\frac{{\lambda }_s}{4}{s}^4,$$

$$\{\begin{array}{l}h=\varphi \cos \theta \\ s=\varphi \sin \theta \end{array}.$$

$${\tan }^2{\theta }^{\ast }=\frac{2{\lambda }_h-{\lambda }_{hs}}{2{\lambda }_s-{\lambda }_{hs}},$$

$${\lambda }_{hs}^2-4{\lambda }_h{\lambda }_s=0,$$

$$B=\frac{1}{32{\pi }^2}\left(W_{\sigma }^4+3W_Z^4+6W_W^4-12W_t^4\right),$$

$$\frac{W_Z^2}{{\cos }^2{\theta }^{\ast }}=\frac{g_L^2+g_Y^2}{4},\frac{W_W^2}{{\cos }^2{\theta }^{\ast }}=\frac{g_L^2}{4},\frac{W_t^2}{{\cos }^2{\theta }^{\ast }}=\frac{y_t^2}{2}.$$

$$\{\begin{array}{l}\varphi =h\cos {\theta }^{\ast }+s\sin {\theta }^{\ast }\\ \sigma =-h\sin {\theta }^{\ast }+s\cos {\theta }^{\ast }\end{array}$$

$$M_{\varphi }^2=4B{v}_{\varphi }^2,M_{\sigma }^2=W_{\sigma }^2{v}_{\varphi }^2,$$

$${\lambda }_s\simeq \frac{M_h^2{v}_h^2}{2v_{\varphi }^4},{\lambda }_{hs}\simeq -\frac{M_h^2}{v_{\varphi }^2},{\sin }^2\theta \simeq 1-\frac{v_h^2}{v_{\varphi }^2}.$$

$$V(H,s,s^{\prime })=$$

Full text

I Introduction

The LIGO collaboration has recently observed for the first time the direct effects of gravitational waves on matter Abbott et al. (2016a), marking the beginning of gravitational wave astronomy. This new experimental field pursues the fascinating possibility that important information about the evolution of the Universe could be encoded in a gravitational wave background. In light of this, as the boundaries of the observed gravitational wave spectrum are to be considerably extended by forthcoming space-based interferometers such as LISA Baker et al. (2007), we can expect that gravitational wave astronomy will soon produce new important observables and benchmarks for models of particle physics and gravitation. In this regard, the possible occurrence of phase transitions in the evolution of our Universe is one of the matters within the reach of current and next-generation gravitational waves experiments.

The dynamics of phase transitions have been investigated in connection to a number of different topics, ranging from the problem of the baryon asymmetry of the Universe No (2011); Huang et al. (2016a); Chala et al. (2016); Katz and Riotto (2016); Artymowski et al. (2017); Vaskonen (2017); Dorsch et al. (2017); Beniwal et al. (2017); Kobakhidze et al. (2017); Saeedhosini and Tofighi (2017) to the natural appearance of such phenomena in several high energy completions of the Standard Model (SM) Schwaller (2015); Kakizaki et al. (2015); Jinno et al. (2016); Huber et al. (2016); Leitao and Megevand (2016); Jaeckel et al. (2016); Dev and Mazumdar (2016); Jinno and Takimoto (2017); Huang et al. (2016b); Hashino et al. (2017); Kubo and Yamada (2016); Balazs et al. (2017); Huang and Zhang (2017); Baldes (2017); Chao et al. (2017); Huang and Zhang (2017); Lewis and Sullivan (2017); Tsumura et al. (2017). Phase transitions which directly or indirectly result in the generation of the electroweak scale have also been analysed within the framework of classically conformal (or scale invariant) models Espinosa and Quiros (2007); Espinosa et al. (2008); Konstandin and Servant (2011a, b); Servant (2014); Fuyuto and Senaha (2015); Sannino and Virkajärvi (2015), often in connection to other open problems of contemporary physics such as the origins of Dark Matter (DM) and the Inflationary dynamics Cooper and Venturi (1981); Espinosa et al. (2008); Foot et al. (2010); Ishiwata (2012); Okada and Orikasa (2012); Heikinheimo et al. (2014); Hambye and Strumia (2013); Farzinnia et al. (2013); Khoze (2013); Gabrielli et al. (2014); Allison et al. (2014, 2015); Kannike et al. (2014, 2015); Karam and Tamvakis (2015); Kannike et al. (2016); Wang et al. (2016); Marzola et al. (2016); Kannike et al. (2017a); Karam and Tamvakis (2016); Kannike et al. (2017b); Marzola and Racioppi (2016). An important feature of classically conformal models is that the phase transition associated to the insurgence of the electroweak scale is generally of the first order and very strong. It is then plausible that the related dynamics result in sizeable gravitational wave signals, which would be nowadays encoded in a stochastic background.

Attracted by this possibility, in this paper we focus on the classically conformal extensions of the SM expanding on previous analysis Espinosa and Quiros (2007); Espinosa et al. (2008); Sannino and Virkajärvi (2015) of the topic, accounting for the impact of thermal corrections on the dynamics of phase transition and extending the phenomenology of these models to cover their possible gravitational wave signal. In more detail, we demonstrate in a general way that scalar conformal models have indeed the capability to give rise to gravitational wave signals that current and future dedicated observations could detect. By applying the developed formalism we then find that the simplest scalar conformal extension of the SM is, at best, strongly constrained by the phenomenology associated to the phase transition once thermal corrections and perturbativity arguments are taken into account. The next-to-minimal models that rely on two new scalar fields, instead, bypass the shortcomings affecting the minimal scenario once the phase transition dynamics is relegated to early epochs in the evolution of the Universe. In this case, we find that the scenarios generally give rise to very strong gravitational wave signals detectable at the current and next-generation gravitational wave interferometers.

The paper is organised as follows: in section II we review the computation of the phase transition temperature and of the emitted gravitational wave spectrum for a classically conformal theory with two scalar fields. The results of this simple model are discussed in relation to the expected sensitivities of aLIGO and LISA observatories. In section III we apply our formalism to a conformal scalar singlet extensions of the SM, extending the work in Espinosa et al. (2008) to scenarios where the Higgs vacuum expectation value (VEV) is induced by the new scalar via a portal coupling. The next-to-minimal models are addressed in section IV, whereas our conclusions are presented in section V.

II Gravitational wave signatures of conformal models

We briefly review here the key steps in the computation of the potential at finite temperature for a conformal scalar field model. At zero temperature, the one-loop scalar potential of $n$ scalar fields $\phi_{j}$, $j\in\{1,\dots,n\}$, is given by the Colemann-Weinberg result Coleman and Weinberg (1973)

[TABLE]

where $\delta V$ contains the counterterms, $\mu$ is the renormalisation scale and $M_{k}$ and $g_{k}$ are, respectively, the field dependent tree level mass and the number or intrinsic degrees or freedom of the particle $k$. Notice that in our convention $g_{k}$ assumes positive values for bosons and negative ones for fermions.

Consider now a direction in the scalar field space defined by $\phi=\sum_{j=1}^{n}a_{j}\phi_{j}$, where $\sum_{j=1}^{n}a_{j}^{2}=1$. Along this direction the tree level mass of the scalar field $\phi$ can be written as $M_{k}=W_{k}\phi$, where $W_{k}$ depends only on adimensional couplings. The potential in eq. (1) along the direction $\phi$ is then written as Sannino and Virkajärvi (2015)

[TABLE]

where

[TABLE]

We require that the tree level potential is flat along the direction $\phi$, $\lambda_{\phi}=0$, and set the counterterm $\delta\lambda$ via the renormalisation condition ${\rm d}^{4}V/{\rm d}\phi^{4}|_{\phi=e^{-11/6}v_{\phi}}=0$. Here $v_{\phi}$ is the VEV of $\phi$ induced in the spontaneous breaking of the symmetry via the Coleman-Weinberg mechanism. In this way the scalar potential along $\phi$ finally reads

[TABLE]

and for the mass of $\phi$ we have $M_{\phi}^{2}=4Bv_{\phi}^{2}$.

The one-loop finite temperature corrections to the above scalar potential are given by Laine and Vuorinen (2016)

[TABLE]

where the thermal integral $J_{T}$ is specified by

[TABLE]

and the upper (lower) sign is for bosons (fermions). The contribution from the re-summed daisy diagrams instead amounts to

[TABLE]

and depends on the Debye mass $\Pi_{k}$ of the boson $k$. Notice that for $T\gg M_{k}$ the thermal integral can be approximated as

[TABLE]

where

[TABLE]

In this regime, the one-loop finite temperature effective potential along the $\phi$ direction is then given by

[TABLE]

where we defined $C=\sum_{k}c_{k}g_{k}W_{k}^{2}/12\geq 0$. Notice that the second derivative of the scalar potential at $\phi=0$ matches $CT^{2}$, so the thermal potential has necessarily a local minimum at $\phi=0$ independently of the specifics of the underlying conformal model.

For sake of definiteness, consider now a minimal scenario comprising $n=2$ scalar fields, $\phi$ and $\sigma$, with a $\mathbb{Z}_{2}$ symmetry that bars linear terms in the potential. Suppose also that the scalar field $\phi$, which lies along the flat direction of the tree-level potential, give a mass $M_{\sigma}^{2}=\lambda_{p}\phi^{2}/2$ to the field $\sigma$ via a (positive) portal coupling $\lambda_{p}$, so that $B=\lambda_{p}^{2}/128\pi^{2}$. The one-loop thermal part of the potential is given by

[TABLE]

where $\Pi_{\sigma}=\lambda_{p}T^{2}/3$ and $g_{\sigma}=1$. We find that using the approximation $V_{T}=M_{\sigma}^{2}T^{2}/12$ leads to underestimating the reference phase transition temperatures ($T_{c}$ and $T_{n}$, as explained later) by a factor of less than $10$ but does not alter the qualitative discussion of the example at hand. The evolution of the total potential $V_{T}$ along the $\phi$ direction as a function of temperature is shown in fig. 1. The temperature at which the minima at $\phi=0$ and $\phi\neq 0$ are degenerate is the critical temperature $T_{c}$.

The phase transition

As made clear from fig. 1, for temperatures of the thermal bath large enough to lift the minimum of the potential, the field111Or, rather, its expectation value which obeys the classical equations of motions. is drawn towards the origin, in a way that $v_{\phi}=0$ and possible symmetries are restored. The field is stuck at this point until $T<T_{c}$, when the potential develops again a new global minimum characterised by a non-zero VEV. As mentioned before, however, thermal corrections in conformal models necessarily result in a potential barrier which separates the origin, a local minimum of the potential, from the minimum corresponding to the true vacuum of the theory. Because such potential barrier disappears only for $T=0$, the phase transition from $v_{\phi}=0$ to $v_{\phi}\neq 0$ in conformal models is always of first order for any finite temperature. We remark that this is a model independent result which applies to all classically conformal scenarios, including the conformal extensions of the Standard Model. Consequently to the abrupt transition, the phase transition proceeds via nucleation and consequent expansion of bubbles inside of which the field is in the broken phase of the theory. The bubble nucleation rate per unit of time and volume is given by Linde (1983)

[TABLE]

where

[TABLE]

is the three-dimensional Euclidean action for an O(3)-symmetric bubble222We remark that the action is minimised by an O(3)-symmetric solution rather than an O(4) one. This is expected whenever the potential barrier arises purely from thermal effects.. The largest contribution to the bubble nucleation rate arises from the path (in the field space) which minimizes $S_{3}$, obtained by solving the equation

[TABLE]

with boundary conditions ${\rm d}\phi/{\rm d}r=0$ at $r=0$, and $\phi\to 0$ at $r\to\infty$. The bubble nucleation temperature, $T_{n}$, is then defined as the temperature at which the probability of producing at least one bubble per horizon volume in Hubble time approaches the unity:

[TABLE]

We plot in fig. 2 the bubble nucleation temperature as a function of the portal coupling $\lambda_{p}$ and the VEV $v_{\phi}$. Requiring that the phase transition occurs above the big bang nucleosynthesis (BBN) temperature, $T_{\rm BBN}\simeq 4$ MeV, excludes values of the portal coupling \lambda_{p}\mathrel{\raise 1.29167pt\hbox{<\kern-7.5pt\lower 4.30554pt\hbox{\sim}}}0.5. Notice also that for \lambda_{p}\mathrel{\raise 1.29167pt\hbox{<\kern-7.5pt\lower 4.30554pt\hbox{\sim}}}2 the bubble nucleation temperature is much lower than the critical temperature $T_{c}$, implying a large amount of supercooling. The transition is then very strong and we can consequently expect a sizeable gravitational wave signal.

Gravitational wave production

The ratio of the vacuum energy released during the phase transition to the energy density of the radiation bath at a temperature $T$ is given by Espinosa et al. (2010)

[TABLE]

where we indicated with $\Delta X$ the difference $X(0)-X(v_{\phi})$ for a quantity $X$. As shown in fig. 3, the value of $\alpha$ for $T\ll T_{c}$ can be very well approximated by

[TABLE]

where $g_{*}$ is the effective number of relativistic degrees of freedom in the thermal plasma. Given that $T_{c}\propto v_{\phi}$, the curves shown in fig. 3 do not depend on the value of $v_{\phi}$.

The thin dashed lines in fig. 2 denote instead the configurations of the model for which $\alpha(T)\simeq 1$, with increasingly larger values falling on the right hand side of the lines. We see that for small \lambda_{p}\mathrel{\raise 1.29167pt\hbox{<\kern-7.5pt\lower 4.30554pt\hbox{\sim}}}2 the vacuum energy released in the transition is much larger than the energy density in the radiation bath, $\alpha\gg 1$, as suggested by the large hierarchy between critical and nucleation temperature. In this regime we then expect that in classically conformal models the dynamics of phase transition result in a substantial reheating of the system, consequent to the scalar field tunnelling through the thermal potential barrier Megevand and Ramirez (2017).

For these values of the parameters, it can be shown that plasma effects do not play an important role in the bubble expansion Caprini et al. (2016). Then, as the bubble wall velocity approaches the speed of light, the gravitational wave signal arises purely from the scalar field contribution Kamionkowski et al. (1994); Caprini et al. (2016), yielding:

[TABLE]

Here $H_{*}$ is the value of the Hubble parameter at $T=T_{*}$, corresponding to the temperature of the radiation bath after the phase transition, whereas $f_{\rm env}$ is the redshifted peak frequency of the spectrum as measurable today,

[TABLE]

The parameter $\beta$ describes instead the duration of the transition,

[TABLE]

In fig. 4 we compare the sensitivity of current and future generation of gravitational waves observatories with the gravitational wave spectrum obtained in the considered general model for the benchmark points given in table 1. As we can see, these experiments have the capability to detect the gravitational echoes of the phase transition in classically conformal models and the results we obtain demonstrate the potential impact of gravitational wave phenomenology on these scenarios.

III Minimal conformal extension of the Standard Model

We apply now the formalism we exemplified for our two-scalars model to scenarios previously considered in the literature, starting with the minimal conformal scalar extension of the SM Allison et al. (2014, 2015); Sannino and Virkajärvi (2015); Marzola and Racioppi (2016). The tree level potential is given in this case by333The Lagrangian we consider admits a $\mathbb{Z}_{2}$ symmetry that, being preserved by the symmetry breaking, could lead to the formation of domain walls. The issue is avoided in UV completion of the proposed scenario which explicitly break the $\mathbb{Z}_{2}$ symmetry or extend it to a continuous gauge group.

[TABLE]

where $H$ is the SM Higgs doublet and $s$ is a real scalar which transforms as a singlet under the gauge symmetry of the SM. Notice that successful electroweak symmetry breaking requires $\lambda_{hs}<0$.

In order to find the direction in the field space where the minimum of the one-loop potential lies, we rewrite the physical Higgs field (in unitary gauge), $h$, and $s$ in polar coordinates $(\phi,\theta)$:

[TABLE]

Then, we solve for the angle $\theta=\theta^{*}$ such that ${\rm d}V/{\rm d}\theta|_{\theta^{*}}=0$, obtaining

[TABLE]

and by imposing the condition

[TABLE]

the tree-level potential is flat along the $\theta=\theta^{*}$ direction. We remark that the choice of couplings encoded in eq. (24) is meant to guarantee the existence of a direction in the potential, corresponding to $\theta=\theta^{*}$, along which quantum corrections dominate over the tree-level contribution. In this way the Coleman-Weinberg mechanism is successfully implemented and the potential acquires the form in eq. (4) along such direction. The $B$ parameter is given here by

[TABLE]

where $W_{\sigma}^{2}=-\lambda_{hs}$ is the contribution arising from the tree-level mass of the eigenstate $\sigma$, perpendicular to $\phi$. The remaining quantities depend instead on the SM parameters as follows:

[TABLE]

The expressions for the mass eigenstates in terms of the original fields $h$ and $s$ are given by

[TABLE]

whereas the associated masses are

[TABLE]

with $v_{\phi}^{2}=v_{h}^{2}+v_{s}^{2}$.

Scenario I: $v_{s}<v_{h}$

With the above formalism at hand we now investigate the scenario where $\phi$ corresponds to the scalar particle observed at the LHC Aad et al. (2012); Chatrchyan et al. (2012). The mixing angle $\theta^{*}$ is constrained by the LHC Higgs phenomenology to the range $\cos\theta^{*}>0.85$ Alanne et al. (2014), which can be recast through the relation $\tan^{2}\theta^{*}=v_{s}^{2}/v_{h}^{2}$ as an upper bound on the VEV of $s$: $v_{s}<152$ GeV. We compute $\lambda_{h}$ from $M_{\phi}^{2}=M_{h}^{2}=(126\,{\rm GeV})^{2}$, whereas $\lambda_{s}$ and $\lambda_{hs}$ are determined through eq. (23) and  (24) by using $\tan^{2}\theta^{*}=v_{s}^{2}/v_{h}^{2}$. As for the remaining parameters, we take the Higgs VEV at $v_{h}=246$ GeV and the following values for the SM couplings: $g_{L}=0.648$, $g_{Y}=0.359$, and $y_{t}=0.951$. The portal coupling corresponding to $\lambda_{p}$ in the example of section II is given here by $\lambda_{p}=\sqrt{2}\,4\pi M_{\phi}/v_{\phi}$.

Fig. 5 shows the parameters of the model as a function of $v_{\phi}$. We see that in the region allowed by the LHC $\lambda_{s}$ and $\lambda_{hs}$ assume very large values: $\lambda_{s}\gtrsim 5$ and $|\lambda_{hs}|\gtrsim 4$. As a consequence, the scenario is strongly impaired by the presence of a Landau pole at relatively low energies. At best, for the smallest allowed magnitudes of the couplings $\lambda_{s}\simeq 5$, $\lambda_{hs}\simeq-4$ and $\lambda_{h}\simeq 0.8$, we estimate with 1-loop RGEs that the Landau pole appears at a scale $\Lambda\sim$ 1 TeV. The consistency of the scenario then imposes the presence of new physics below such scale, clashing with the null results of current LHC searches.

Scenario II: $v_{h}<v_{s}$

We analyse next the complementary case in which the detected Higgs boson corresponds to $\sigma$. The mixing angle imposed by the LHC Higgs phenomenology is then large, $\sin\theta^{*}>0.85$, and consequently $v_{s}>397$ GeV. For $v_{h}\ll v_{s}$ we obtain $\lambda_{h}\simeq M_{h}^{2}/(2v^{2})=0.131$ for the Higgs quartic coupling, while the $s$ quartic coupling, the portal coupling and the mixing angle are given by:

[TABLE]

The portal coupling along the $\phi$ direction matches here $\lambda_{p}=2M_{h}^{2}/v_{\phi}^{2}$ and, in the region consistent with the LHC data, is typically small: $\lambda_{p}\lesssim 0.14$. We plot in fig. 6 the values obtained for this parameter against the mixing angle. By comparing these results to the ones obtained in section II, we conclude that this scenario is excluded for reheating temperatures high enough to restore the electroweak symmetry via thermal effects. The region of the parameter space associated to the scenario allowed by the LHC constraints leads in fact to a bubble nucleation temperature that violates the lower bound posed by the BBN temperature.

IV Next-to-minimal model

We showed in the previous section that the minimal conformal extension of the SM is strongly disfavoured by the electroweak phase transition phenomenology and perturbativity arguments. Here we consider instead the next-to-minimal scenario Foot et al. (2010); Ishiwata (2012); Farzinnia et al. (2013); Khoze (2013); Gabrielli et al. (2014), where two new real scalar fields $s$ and $s^{\prime}$, both singlets under the gauge symmetry of the SM, couple to the Higgs boson. As customary in this framework, we assign a $\mathbb{Z}_{2}$ symmetry to exclude terms containing odd powers of the new fields in the Lagrangian of the model. The tree level scalar potential then reads

[TABLE]

In this analysis we limit ourselves to the case $\lambda_{ss^{\prime}},\lambda_{s^{\prime}}>0$, so that $v_{\sigma}\equiv 0$, and set $\lambda_{hs^{\prime}}=0$ for simplicity. As we will see, this simplified scenario is sufficient to show the potential impact of gravitational wave experiments.

In fact, although the presence of a new scalar cannot prevent the Landau pole from appearing at relatively low scale in the setup of the scenario I considered above, the extra portal coupling allows us to recover a correct phase transition dynamics in the complementary case of scenario II. Supposing that $v_{s}\gg v_{h}$, the dynamics of this next-to-minimal scenario differs from that of the former case only by the definition of $B$, which includes here the term arising from the portal coupling between $s^{\prime}$ and $s$:

[TABLE]

The BBN constraint that impairs the scenario II can then be overcome provided that $B$ is dominated by the new contribution above.

By identifying $\lambda_{p}\simeq\lambda_{ss^{\prime}}$ we can discuss the phenomenology of the model by using the results obtained in section II. In this regard, we plot in fig. 7 the parameters of the model as a function of the VEV $v_{\phi}$, which as usual lies along the flat direction of the tree-level potential. As we can see, by considering large values of $v_{\phi}$ we can allow for lower values of $\lambda_{p}$ (i.e. of $\lambda_{ss^{\prime}}$), so that the parameters of the model retain perturbativity up to the Grand Unified Theory scale (for $\lambda_{ss^{\prime}}\lesssim 1.0$) or even Planck scale (in this case $\lambda_{ss^{\prime}}\lesssim 0.9$). The BBN constraint, instead, is straightforwardly satisfied whenever $\lambda_{ss^{\prime}}\gtrsim 0.5$ for the reasons previously explained.

As demonstrated in Gabrielli et al. (2014), in absence of effects that break the $\mathbb{Z}_{2}$ symmetry imposed on the Lagrangian, $s^{\prime}$ is a stable particle which can play the role of DM candidate. A first rough estimate of the the relative relic abundance shows that the observed value can be matched through the freeze-out mechanism444For $\lambda_{ss^{\prime}}\sim\mathcal{O}(1)$ the dark matter candidate thermalises in the early Universe provided that $\lambda_{hs}\gtrsim 10^{-7}$. via annihilations of $s^{\prime}$ to the $ss$ final state for a range of values of the involved parameters. A more careful assessment of the DM relic density in the scenario and the detailed gravitational wave phenomenology of the considered next-to-minimal conformal extension of the SM will be presented in a forthcoming analysis.

V Conclusions

In this paper we studied the phase transition dynamics of scalar classically conformal scenarios posing particular attention to their possible gravitational wave signatures. After having reviewed the basis of the framework and showed that it necessarily leads to a first order phase transition at finite temperature, we studied the properties of a general two-scalars model that captures the gist of the conformal extensions of the Standard Model. We find that the phase transition is generally very strong and leads to the production of a stochastic gravitational wave background which can be observed in current and next-generation dedicated experiments.

We then applied the analysis for the minimal classically conformal scalar extension of the Standard Model, where the presence of a new scalar field coupled to the Higgs boson implements a Coleman-Weinberg type of potential. The scenario can be analysed in two limits, depending on the hierarchy between the vacuum expectation values of the involved scalar fields. In the case where the vacuum expectation value of the singlet scalar field is smaller than that of the Higgs field, the quartic couplings of the model are very large and result in a Landau pole at the TeV scale. Although the problem could be solved by invoking the presence of new physics, the null result from the LHC disfavour this possibility.

In the complementary regime, where the vacuum expectation value of the new scalar boson is larger than the Higgs field one, we find that the model retains its perturbativity up to the Planck scale. In spite of that, the smallness of the portal coupling between the two scalars delays the electroweak phase transition to temperatures below the Big Bang Nucleosynthesis one, generally excluding the scenario.

Lastly, we considered a next-to-minimal conformal extension of the Standard Model involving an additional scalar field. In this case we showed that the presence of such particle allows to satisfy the constraint from Big Bang Nucleosynthesis and perturbativity, which impaired the minimal extension. We find that, in the considered regime, the next-to-minimal extension accomplishes a first order electroweak phase transition which gives rise to a sizeable gravitational wave signal, demonstrating the capability of dedicated experiments to explore the scenario.

Acknowledgements

This work was supported by the Estonian Research Council grant IUT23-6, PUT1026 and ERDF Centre of Excellence project TK133.