Abstract
We estimate the gravitational wave spectra generated in strongly supercooled phase transitions by bubble collisions and fluid motion. We derive analytically in the thin-wall approximation the efficiency factor that determines the share of the energy released in the transition between the scalar field and the fluid. We perform numerical simulations including the efficiency factor as a function of bubble radius separately for all points on the bubble surfaces to take into account their different collision times. We find that the efficiency factor does not significantly change the gravitational wave spectra and show that the result can be approximated by multiplying the spectrum obtained without the efficiency factor by its value at the radius \(R_{\textrm{eff}} \simeq 5/\beta \), where \(\beta \) is the approximate inverse duration of the transition. We also provide updated fits for the gravitational wave spectra produced in strongly supercooled transitions from both bubble collisions and fluid motion depending on the behaviour of the sources after the collision.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
The first observations of gravitational waves (GWs) by LIGO/Virgo signified the beginning of a new era in astrophysics and cosmology. While up to now all observed events were produced by compact object binaries [1,2,3,4], this new messenger brings hope also for detection of primordial signals in the form of stochastic GW backgrounds. Given the tremendous advancements in sensitivity that are expected throughout a broad frequency spectrum with the upcoming experiments [5,6,7,8,9,10,11,12,13,14,15], the prospects for probing the early Universe processes are great even though the compact object binaries that will contribute to the stochastic GW background make the detection of its primordial components more difficult [16]. Interestingly, the recent pulsar timing observations [17,18,19,20] feature a common spectrum-process which could be an early indication of the upcoming first detection of a stochastic GW background, potentially of primordial origin [21,22,23,24,25,26,27,28,29,30,31,32,33].
Many high-energy processes, including phase transitions [34, 35], cosmic strings [36] and inflation [37], occurring in the early Universe may generate a detectable stochastic GW background. In this paper we focus on cosmological first-order phase transitions featured in various particle physics models. They are intensive processes where bubbles of the new phase nucleate, expand and eventually convert the whole Universe in the true vacuum phase [38]. Interactions between the expanding bubble walls and the surrounding fluid cause motion and inhomogeneities in the fluid, and both the collisions of the bubble walls and the motion of fluid inhomogeneities source GWs [39, 40]. The resulting GW spectra from these components have been extensively studied with numerical and semi-analytical methods (see e.g. [41,42,43,44,45,46,47] for recent progress). These studies indicate that different sources active during the transition can produce different GW spectra.
In order to determine the GW spectrum generated in a phase transition in a given particle physics model, we need to estimate how much each of the GW sources contributes to the final GW spectrum. The vacuum energy released in the transition is split between the gradient energy of the scalar field bubble wall and motion in the fluid. How the total released energy is split depends on the strength of the interactions between the wall and the particles in the fluid, and on the strength of the transition.
In strongly supercooled phase transitions it is possible that the interactions of the bubble wall with fluid do not stop the wall from accelerating before it collides with other bubbles. In this case most of the released energy is in the bubble walls and the bubble collisions give the dominant contribution to the GW spectrum. This can happen in particular in quasi-conformal models [48,49,50,51,52,53,54,55,56,57,58,59,60,61]. If the interactions instead are sufficiently strong, the bubble wall reaches a terminal velocity before the collisions and majority of the released energy goes into fluid motion. This is the typical case in extensions of the Standard Model featuring polynomial scalar potentials [62,63,64,65,66,67,68,69,70,71,72].
In this paper we derive analytically an efficiency factor that determines how large is the contribution from each of the GW sources. We perform numerical simulations of the phase transition, describing both of the GW sources, bubble walls and fluid motion, in the thin-wall limit, to show how the efficiency factor affects the final GW spectrum. Moreover, we derive analytically the probability density function for the radius at which a given point on the surface of a bubble collides with another bubble and verify the results against our numerical simulations. Finally we also provide updated fits to the spectral shapes of the GW signals that can be produced by all sources active in very strong phase transitions.
2 Energy budget
We estimate how the energy released in the bubble expansion is shared between the scalar field gradients and the fluid motion in strongly supercooled phase transitions by studying the bubble expansion under the influence of pressure terms caused by the interactions of the wall with the ambient fluid. We perform the computation consistently in the thin-wall limit, which gives a good description of the system if the bubble reaches ultra-relativistic velocities. The following analysis improves earlier approximations used in the literature [58, 71, 73].
In the thin-wall limit the energy carried by the bubble walls can be modeled using a simple analytical prescription. This assumes that the bubble walls are spherical shells with a certain surface energy density and the interactions of the walls with the ambient fluid are local. In this limit, neglecting the expansion of the Universe,Footnote 1 the evolution of the bubble radius R can be described by the equation of motion [74]
that arises from energy conservation of the coupled system of fluid and the scalar field bubble. Here \(\Delta P(\dot{R})\) denotes the pressure difference across the bubble wall and \(\sigma \) is the surface tension of the wall. The latter is defined through the scalar potential V as [38]
where \(\varphi _c>0\) denotes the field value at which the potential energy is the same as in the false vacuum that lies at the origin, \(V(\varphi _c) = V(0)\).
In terms of the Lorentz factor of the bubble wall, \(\gamma = 1/\sqrt{1-\dot{R}^2}\), the equation of motion (1) is given by
The bubble nucleates at rest, \(\gamma =1\), \(\textrm{d}\gamma /\textrm{d}R = 0\), with an initial radius \(R_0\). By Eq. (3) we can relate the wall tension to the initial radius as \(R_0 = 2\sigma /\Delta P_0\), where \(\Delta P_0 \equiv \Delta P(\gamma =1)\). For a constant pressure difference, \(\Delta P = \Delta P_0\), the solution of the equation of motion is
The total pressure difference across the bubble wall, accounting for \(1\rightarrow 1\) scatterings and \(1\rightarrow N\) splittings at the bubble wall, is given by
where \(\Delta V\) denotes the potential energy difference between the minima. The pressure arising from \(1\rightarrow 1\) scatterings quickly reaches a constant value in the relativistic limit [74, 75]. Subsequently, the \(\gamma \) dependence of the total pressure difference arises only from \(\Delta P_{1\rightarrow N}\), for which we consider two forms. The first one, suggested in [76,77,78] is linear in the Lorentz factor \(\Delta P_{1\rightarrow N} = {\tilde{\Delta }} P_{1\rightarrow N} \gamma \), and the second one, suggested in [79, 80], is quadratic in the Lorentz factor, \(\Delta P_{1\rightarrow N} = {\tilde{\Delta }} P_{1\rightarrow N} \gamma ^2\) In both cases \({\tilde{\Delta }} P_{1\rightarrow N}\) is a constant.
By plugging \(\Delta P(\gamma )\) into (3), we can easily solve \(\gamma \) as a function of R. When \(\Delta P_{1\rightarrow N} \ll \Delta V - \Delta P_{1\rightarrow 1}\) the solution can be approximated by Eq. (4). Assuming in addition that \(R\gg R_0\), the Lorentz factor grows linearly with the radius, \(\gamma \approx 2R/(3R_0)\). Eventually, as the bubble wall accelerates, \(\gamma \) becomes large enough for the \(1\rightarrow N\) splittings to be important, \(\Delta P_{1\rightarrow N} \sim \Delta V - \Delta P_{1\rightarrow 1}\), after which it asymptotically reaches the value
where \(c=1,2\) depending on the scaling of the \(1\rightarrow N\) pressure, \(\Delta P_{1\rightarrow N} \propto \gamma ^c\). The change from the linear growth to asymptotically constant behaviour occurs when the radius reaches
The solution \(\gamma (R)\), in the limit \(\gamma _{\textrm{eq}}\gg 1\), can be approximated by a simple broken power-law
In Fig. 1 we show the full solution \(\gamma (R)\) for different values of \(\gamma _{\textrm{eq}}\) for both \(\Delta P_{1\rightarrow N}\propto \gamma \) (solid) and \(\Delta P_{1\rightarrow N}\propto \gamma ^2\) (dashed) together with the above approximation shown by the dotted black curves. In both cases the transition from linear growth, \(\gamma \propto R\) to the constant value \(\gamma \approx \gamma _\textrm{eq}\) is quite fast, and the difference between the two cases is small. The main effect of the scaling of \(\Delta P_{1\rightarrow N}\) is that it changes \(\gamma _{\textrm{eq}}\) and \(R_{\textrm{eq}}\).
We define the efficiency factor \(\kappa \) as the fraction of the total released energy within a unit solid angle that goes into the bubble wall energy,
The rest of the released energy, \(1-\kappa (R)\), goes into fluid motion. This is a good approximation for strongly supercooled transitions. In weaker transitions one also needs to keep in mind that some of the energy going into the fluid will be lost on its heating which will reduce the overall GW signal from the fluid motion [71, 81].
Using the approximation (8), we can express the efficiency factor as
where
is a constant, \({\mathcal {K}} < 1\). The parameters \(\alpha \) and \(\alpha _\infty \) are defined by scaling with the radiation energy density \(\rho _R\) as \(\alpha = \Delta V/\rho _R\) and \(\alpha _\infty = \Delta P_{1\rightarrow 1}/\rho _R\) (see [71] for more details). Typically for strongly supercooled transitions \({\mathcal {K}} \approx 1\). As shown in the right panel of Fig. 1, the efficiency factor remains constant at \(R\ll R_{\textrm{eq}}\) and decreases as \(\kappa \propto 1/R\) at \(R\gg R_{\textrm{eq}}\). In the same way as for \(\gamma (R)\), the difference between the cases \(\Delta P_{1\rightarrow N}\propto \gamma \) and \(\Delta P_{1\rightarrow N}\propto \gamma ^2\) is small.
3 Bubble nucleation and collisions
Soon after the bubble has nucleated, we can neglect its initial radius, and, if the friction terms are sufficiently small (\(\gamma _{\textrm{eq}}\gg 1\)), we can approximate that the bubble radius grows as \(R = t-t_n\), where \(t_n\) denotes the nucleation time of the bubble. Moreover, assuming that the bubbles are much smaller than the Hubble horizon, we can neglect the expansion of the Universe. The expected number of bubbles reaching a given point is then given by
where \(\Gamma (t)\) denotes the bubble nucleation rate per unit time and volume, and the probability that the given point still is in the false vacuum at time t is
Let us consider a bubble nucleation rate \(\Gamma (t) = C e^{A(t)}\). Around the time \(t_*\) when the transition proceeds, we can expand A(t) to get \(\Gamma (t) = C e^{A(t_*) + \beta (t-t_*)} = \Gamma _0 e^{\beta t}\), where \(\beta \equiv \textrm{d}\ln \Gamma /\textrm{d}t |_{t=t_*}\) and \(\Gamma _0 \equiv C e^{A(t_*) - \beta t_*}\). As the transition is not an instantaneous process, the choice of \(t_*\) includes some ambiguity. It is convenient to choose \(t_*\) by requiring that \(P(t_*) = 1/e\), which gives \(\Gamma _0 = \beta ^4/(8\pi )\), and \(N(t) = e^{\beta t}\).
Next, let us consider a point on the surface of a bubble that nucleated at time \(t_n\). If the point is still in the false vacuum when the radius of the bubble is R, then it has stayed in the false vacuum for the whole time \(0 \le t-t_n < R\). The probability for this is \(P(t_n + R)\). So, the probability that a bubble nucleated within time \(t_n< t < t_n + \textrm{d}t_n\) in a volume V, and that a point on its surface is still in the false vacuum at radius R, is given by \(\textrm{d}t_n \,V \,\Gamma (t_n) P(t_n + R)\). By integrating this over the nucleation time \(t_n\) we get the probability density function for the radius at which a bubble surface element collides with the surface of another bubble,
which we normalize to unity, \(\int \textrm{d}R_c \,p(R_c) = 1\). For the exponential bubble nucleation rate, \(\Gamma (t) \propto e^{\beta t}\), this gives (independently of the prefactor \(\Gamma _0\))Footnote 2
The above result provides a good cross-check for the numerical simulations that we will use for the GW computation. In Fig. 2 the solid curve and the gray band indicate the mean and variance of the \(R_c\) distribution obtained from 90 simulations with simulation volume \((16/\beta )^3\) and each of the simulation including at least 70 bubbles. In these simulations we nucleate thin-wall bubbles according to the exponential rate inside a cubic box with periodic boundary conditions, evolve them according to \(R = t-t_n\), discretise the bubble surfaces, and find the radius at which each of the points on the bubble surface collide with a wall of another bubble using the cosine rule. We label the bubbles with index j and denote the position vectors of the bubble centers by \({\vec {x}}_j\). Consider a point defined by the angles \(\theta \) and \(\phi \) on the surface of the bubble \(j=j'\). The radius at which that point collides with a surface of another bubble is given by
where the minimum is taken over all bubbles, \(d_{j}^2 \equiv |{\vec {x}}_j - {\vec {x}}_{j'}|^2\) is the distance between the bubble nucleation centers, \(\Delta t_{j} \equiv t_{n,j} - t_{n,j'}\) is the time between their nucleation, and \(\theta _{j}\) is the angle between the vector \({\vec {x}}_j - {\vec {x}}_{j'}\) and the vector corresponding to the angles \(\theta \) and \(\phi \). As shown in Fig. 2, the simulation result agrees well with the analytical result (15).
A widely used approximation for the bubble radius upon collision comes from the bubble number density \(n_{\textrm{bubbles}} = \int \textrm{d}t_n \Gamma (t_n) P(t_n)\) which leads to \( R_* = n_{\textrm{bubbles}}^{-1/3} = (8\pi )^{1/3}/\beta \). From \(p(R_c)\) we can calculate moments of the bubble radius when a bubble surface element collides with the surface of another bubble, \(\langle R_c^n\rangle = \int \textrm{d}R_c\, R_c^n p(R_c)\). For the exponential bubble nucleation rate this gives
Given that the released energy scales with the radius to the third power, this leads to a different estimate of the average bubble radius \(\langle R_c^3 \rangle ^{1/3} = 6^{1/3}/\beta \) more appropriate for computation of the GW spectrum.
4 Gravitational waves
The energy released in the bubble expansion is divided between the bubble wall and the fluid shell that follows right behind the wall. Both the bubble walls and the fluid shells source GWs. We model these sources in the thin-wall limit and calculate the GW spectrum accounting for the efficiency factor \(\kappa (R)\) for the bubble collisions and \(1-\kappa (R)\) for the fluid motion. The modeling of the fluid motion in the thin-wall limit is based on the assumption that the released energy going to fluid motion is strongly localized in a thin shell. Before collision this fluid shell is right behind the bubble wall and after the collision it propagates to the same direction as before the collision,however , depending on how strong the interaction are between the fluid and the scalar field, it’s velocity can slow down to the speed of sound.
We calculate the GW spectrum as e.g. in Refs. [43, 44] assuming that, as in the previous section, the bubble nucleation follows exponential rate per unit volume, \(\Gamma \propto e^{\beta t}\). Each of the contributions (\(l = \) bubbles, fluid) to the resulting energy spectrum of GWs can be expressed as
where
encodes the spectral shape of the signal. The integral is over the wavevector \(\vec {k}\) directions, and the integrand is \(\propto V_s/\beta ^5\) if the volume \(V_s\) over which we average the GW energy spectrum is sufficiently big.
Using the thin-wall limit, the functions \(C_{l,+}\) and \(C_{l,\times }\) in the direction \({\hat{k}} = (0,0,1)\), can be expressed as
The sum runs over all the bubbles nucleated in the volume \(V_s\), \(t_{n,j}\) is the time of nucleation of the bubble j, \(z_j\) is the z coordinate of its center, and \(R_j = v_l (t-t_{n,j})\), where \(v_l\) is the bubble wall/fluid shell velocity, denotes the radius of the bubble/fluid shell j at time t. For the bubble walls we use \(v_{\textrm{bubbles}}=1\) both before and after the collision, whereas for the fluid shells we use \(v_{\textrm{fluid}}=1\) before the collision and after the collision we consider two cases: \(v_{\textrm{fluid}} = 1\) and \(v_{\textrm{fluid}} = c_s = 1/\sqrt{3}\). The former is appropriate for very strong transitions [82], whereas the latter is realized for weaker transitions [45]. The functions \(g_{+,\times }\) read \(g_+(\phi ) = \cos (2\phi )\) and \(g_\times (\phi ) = \sin (2\phi )\).
The function \(f_l(R)\) encodes the scaling of the GW source [44]. For the bubble collisions contribution, we follow the results of lattice simulations [43, 44], which showed that the maximum of the stress-energy tensor scales as \(T_{rr} \propto R^{-\xi }\) after the collision. The power \(\xi \) in general depends on the underlying particle physics model. In particular, it was shown in [43] that breaking of a global symmetry corresponds to \(\xi =2\) while in [44] it was shown that in models where the phase transition breaks a gauge symmetry correspond to \(\xi =3\). Accounting also for the efficiency factor \(\kappa \), the \(f_l\) function for bubble collisions is given by
where \(R_c\) denotes the bubble radius at the moment of collision, \(t=t_c\). In contrast with Refs. [43, 44], where \(R_c\) was determined numerically by the bisection method, we find \(R_c\) using Eq. (16).
Also for the fluid motion we assume that the maximum of the stress-energy tensor scales as \(R^{-\xi }\) after the collision. The function \(f_l\) for fluid motion then reads
In the perfect fluid description, that assumes the fluid to remain in local equilibrium at all times, the transverse-traceless part of the stress energy tensor of the fluid reads \(T_{ij} = \gamma ^2 v_i v_j w\), where \(\vec {v}\) is the fluid velocity and w is its enthalpy density. By the interactions of the fluid with the wall, an overdense fluid shell with radial velocity \(v_r > 0\) builds up around the bubble wall. If the wall reaches a terminal velocity, the fluid shell settles into a self-similar profile [81]. We expect that this shell continues to propagate after the bubble wall collides with the wall of another bubble without changing its shape significantly in the collision. This behaviour was first observed for weaker transitions, \(\alpha < {\mathcal {O}}(0.1)\), in lattice simulations [83, 84]. For our case of very strong transitions, \(\alpha \ge {\mathcal {O}}(10)\), we used a simplified simulation involving only the fluid and assuming extra symmetry of the system to retain only one spatial dimension as in [45].Footnote 3 We begun with simulating collisions of two planar shells and verified that they are not significantly modified and instead simply propagate onward. We next simulated the evolution of spherically symmetric fluid shells after the collision. Figure 3 shows an illustrative example of our results. We find that the maximum of the rr component of the stress energy tensor scales as \(T_{rr} \propto R^{-3}\). This matches to the same scaling found in [45] for weak transitions and motivates us to consider \(\xi =3\) for the scaling of the fluid related GW source. For comparison, we consider also \(\xi = 2\) which corresponds to the bulk flow model [86].
In principle, after the GW generation through relativistic fluid shells finishes, one would expect to enter the period of sound waves [41, 83, 84, 87] and perhaps also turbulence [88,89,90,91]. However, in the very strong transitions, that are our primary interest, we expect that the fluid will remain in the relativistic shells for a long time after the transition (at least until the shell radius has grown by \({\mathcal {O}}(1)\) factor). The energy carried by the inhomogeneities has then significantly diluted once the sound wave and turbulence periods begin and, therefore, we expect the main contribution on the GW spectrum in very strong transitions to arise from the relativistic fluid shells or the scalar field bubbles themselves. Thus, in the present analysis we neglect the periods of sound waves and turbulence.
For certain simple forms of the \(f_l\) function the time integral in Eq. (20) can be done analytically, which makes the simulation significantly faster. In particular, it can be done analytically if \(f_l\) is a broken power-law with integer powers. We consider the form (10) for the efficiency factor \(\kappa \) with \(c=1\). Strictly speaking our results then hold for the case that \(\Delta P_{1\rightarrow N} \propto \gamma \). However, since the difference in \(\kappa (R)\) for \(c=1\) and \(c=2\) is very small, our results give a good approximation also of the case \(\Delta P_{1\rightarrow N} \propto \gamma ^2\). The pressure \(\Delta P_{1\rightarrow N}\) mainly just determines the asymptotic radius \(R_{\textrm{eq}}\) through Eqs. (7) and (6). In our simulations \(R_{\textrm{eq}}\) is an input parameter, and we perform the numerical simulations for several values of \(R_{\textrm{eq}}\). We also assume that \({\mathcal {K}}\approx 1\), which typically holds for strongly supercooled transitions, so that
5 Results
We perform 90 simulations with simulation volume \((16/\beta )^3\), each including at least 70 bubbles, for a range of \(R_{\textrm{eq}}\) values including both signals due to bubble walls and the surrounding fluid in each of the cases described in the previous section. From the simulations we compute the spectral shape function (19). In each case we fit the data combined from the 90 simulations with a broken power-law spectrum of the form
where \(a,b>0\) determine the low and high frequency power-law tails of the spectrum, \(c>0\) the width of the transition between these power-laws, while \(f_p\) and A the peak frequency and amplitude of the spectrum respectively. The resulting GW spectra are shown in Fig. 4 with the solid curves. The color coding indicates different values of \(R_{\textrm{eq}}\).
For the solid curves in Fig. 4 the efficiency factor is directly included in the simulation as in Eqs. (21) and (22). A commonly used approximation for the effect of the efficiency factor on the GW spectrum is to multiply the spectra obtained for bubble collisions and fluid motion without any efficiency factor by \(\kappa (R_{\textrm{eff}})^2\) and by \((1-\kappa (R_{\textrm{eff}}))^2\), respectively. To check this, we have computed the amplitude of the GW spectrum in each case relative to the \(R_{\textrm{eq}}\) case that gives the largest amplitude and fitted \(R_{\textrm{eff}}\). The data points and resulting fits for all cases are shown in Fig. 5 and the corresponding fitted values of \(R_{\textrm{eff}}\) in the last line of Table 1. We find that the effect of the efficiency factor is almost independent of the behaviour of the GW source after the collisions. In all cases our results give \(R_{\textrm{eff}} \simeq 5/\beta \), showing that the often used approximation with \(R_\textrm{eff} = (8\pi )^{1/3}/\beta \approx 2.9/\beta \) slightly underestimates \(R_{\textrm{eff}}\). Moreover, the results of applying the fitted efficiency factor are shown in Fig. 4 by the dashed curves. For these curves we have used the mean values given in Table 1 that are obtained by averaging over the fits with different equilibrium radius, except for the amplitude for which we use the strongest signal for each source. We see that the dashed curves agree very well with the fully numerical results shown by the solid curves. This shows that the efficiency factor does not change the shape of the GW spectrum but gives only an overall suppression factor.
To summarize, we have shown that for very strong transitions, \(\alpha \gg \alpha _\infty \), the GW spectrum from bubble collisions and from fluid motion, accounting for the distribution of energy between these sources, is given by
where the efficiency factor is given by Eq. (23). The fitted values of the parameters \(A, a, b, c, f_p\) and \(R_\textrm{eff}\) are given in Table 1. For weaker transitions, \(\alpha \lesssim \alpha _\infty \), also the prefactor \({\mathcal {K}}\), given in Eq. (11), needs to be accounted for, as well as the suppression arising from heating of the fluid around the bubble wall [81]. In the limit of large wall velocity appropriate for strong transitions this reduction takes a simple form [71]
where \(\alpha _{\textrm{eff}} = [1-\kappa (R_{\textrm{eff}})]\alpha \).
The GW spectrum today can be obtained from (25) by accounting for the scaling of the amplitude and frequency with the scale factor [43]:Footnote 4
where h denotes the dimensionless Hubble constant, \(h = 0.674\) [93], and \(h_*\) the inverse Hubble time at the transition redshifted to its value today
Here \(T_*\) denotes the temperature after the transition (including reheating) and \(g_*\) the effective number of relativistic degrees of freedom at that temperature. At scales larger than the horizon scale at the time of the transition the source is not coherent and, consequently, in standard radiation domination the slope of the spectrum changes to \(\Omega _{\textrm{GW}}\propto f^3\) for \(f < h_*\) [94, 95].Footnote 5
6 Conclusions
In this paper we have revisited the energy budget of strong first-order phase transitions to verify its impact on the produced gravitational wave spectra. We have gone beyond the current state-of-art by including the efficiency as a function of radius of the bubble accounting for the collision time of each point on the bubble surface. We have utilised numerical simulations randomly nucleating bubbles in a three dimensional box with periodic boundaries and used these to compute the GW spectra. This has allowed us to confirm that a simplified treatment of simply scaling entire spectra with an efficiency factor computed at some characteristic radius is accurate as the spectral shapes do not change due to the efficiency factor significantly. We did, however, find that in order to accurately describe the results the characteristic radius used in the simplified calculation should be around \(R_{\textrm{eff}}\approx 5/\beta \) rather than the usually employed average bubble separation \(R_*= (8\pi )^\frac{1}{3} /\beta \approx 2.9 /\beta \).
In each simulation we have also took into account the scaling of the GW sources after the collision in order to provide new fits for the resulting spectra from strongly supercooled transitions. The results are shown in Table 1 and, starting from strongest transitions, include bubble collision spectra for both \(T_{rr}\propto R^{-3}\), appropriate for gauge symmetry breaking, and \(T_{rr}\propto R^{-2}\), appropriate for global symmetry breaking. Going towards slightly weaker transitions, we have provided the spectrum generated by fluid motion with the scaling \(T_{rr}\propto R^{-3}\) and assuming the fluid remains in the form of relativistic shocks \(v_{\textrm{fluid}}=1\) after the transition. For transitions which are not extremely strong, we have show results closer to the sound wave picture in which the velocity of the fluid quickly relaxes to the speed of sound \(v_{\textrm{fluid}}=c_s\), again assuming the scaling \(T_{rr}\propto R^{-3}\). Finally, for illustration, we also provide fluid spectra assuming the scaling \(T_{rr}\propto R^{-2}\).
Taking into account that for very relativistic walls the fluid profiles are extremely peaked, we have thus show that the final GW spectrum will be indistinguishable from an even stronger transition where bubble collisions would be the main source. Only for weaker transitions where the hydrodynamical effects change the propagation speed of the fluid shells, the spectrum diverges from the spectrum arising from bubble collisions.
Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors’ comment: All relevant data is contained in the manuscript.]
Notes
The assumption of neglecting the expansion of the Universe is in the end related to the nucleation rate of bubbles, and it is valid if the bubble radius at collision moment is much smaller than the Hubble radius, \(\langle R_c \rangle \ll 1/H\), which translates to the requirement that \(\beta /H \gg 1\) (see Sect. 3).
This agrees with the distribution on the bubble lifetime derived in Ref. [41]. Our result can be generalized to wall velocities \(v_w < 1\) simply by dividing \(R_c\) by \(v_w\) in the exponent.
We use the integration scheme devised in [85] in order to reduce the numerical diffusion. The parameters used in the simulations are \(\Delta x=10^{-4}R_c\) for the step in radius and a tenth of that in time. We also verified reducing the step sizes by an order of magnitude does not modify the results.
Here for simplicity while red-shifting we assumed radiation dominated expansion from the transition time up to the matter-radiation equality. For a review of alternative scenarios and their impact on the spectra see Ref. [92].
The low frequency slope is also changed by possible modifications of the expansion rate [96,97,98] although the only scenario in which the signal is not diminished is when the modification in question is itself caused by the transition for instance through slow decay of the scalar field leading to a period of matter domination [58].
References
B.P. Abbott et al. (LIGO Scientific, Virgo), Phys. Rev. X 9, 031040 (2019). https://doi.org/10.1103/PhysRevX.9.031040. arXiv:1811.12907 [astro-ph.HE]
R. Abbott et al. (LIGO Scientific, Virgo), Phys. Rev. X 11, 021053 (2021). https://doi.org/10.1103/PhysRevX.11.021053. arXiv:2010.14527 [gr-qc]
R. Abbott et al. ( LIGO Scientific, VIRGO), (2021). arXiv:2108.01045 [gr-qc]
R. Abbott et al. (LIGO Scientific, VIRGO, KAGRA) (2021). arXiv:2111.03606 [gr-qc]
M. Punturo et al., Class. Quantum Gravity 27, 194002 (2010). https://doi.org/10.1088/0264-9381/27/19/194002
S. Hild et al., Class. Quantum Gravity 28, 094013 (2011). https://doi.org/10.1088/0264-9381/28/9/094013. arXiv:1012.0908 [gr-qc]
G. Janssen et al., PoS AASKA14, 037 (2015). https://doi.org/10.22323/1.215.0037. arXiv:1501.00127 [astro-ph.IM]
P.W. Graham, J.M. Hogan, M.A. Kasevich, S. Rajendran, Phys. Rev. D 94, 104022 (2016). https://doi.org/10.1103/PhysRevD.94.104022. arXiv:1606.01860 [physics.atom-ph]
P. Amaro-Seoane et al. (LISA) (2017). arXiv:1702.00786 [astro-ph.IM]
P.W. Graham, J.M. Hogan, M.A. Kasevich, S. Rajendran, R.W. Romani (MAGIS) (2017). arXiv:1711.02225 [astro-ph.IM]
L. Badurina et al., JCAP 05, 011 (2020). https://doi.org/10.1088/1475-7516/2020/05/011. arXiv:1911.11755 [astro-ph.CO]
Y.A. El-Neaj et al. (AEDGE), EPJ Quantum Technol. 7, 6 (2020). https://doi.org/10.1140/epjqt/s40507-020-0080-0. arXiv:1908.00802 [gr-qc]
A. Bertoldi et al., Exp. Astron. 51, 1417 (2021). https://doi.org/10.1007/s10686-021-09701-3
I. Alonso et al., EPJ Quant. Technol. 9(1), 30 (2022). https://doi.org/10.1140/epjqt/s40507-022-00147-w. arXiv:2201.07789 [astro-ph.IM]
L. Badurina, O. Buchmueller, J. Ellis, M. Lewicki, C. McCabe, V. Vaskonen, Philos. Trans. A. Math. Phys. Eng. Sci. 380, 20210060 (2021). https://doi.org/10.1098/rsta.2021.0060. arXiv:2108.02468 [gr-qc]
M. Lewicki, V. Vaskonen (2021). arXiv:2111.05847 [astro-ph.CO]
Z. Arzoumanian et al. (NANOGrav), Astrophys. J. Lett. 905, L34 (2020). https://doi.org/10.3847/2041-8213/abd401. arXiv:2009.04496 [astro-ph.HE]
B. Goncharov et al., Astrophys. J. Lett. 917, L19 (2021). https://doi.org/10.3847/2041-8213/ac17f4. arXiv:2107.12112 [astro-ph.HE]
S. Chen et al., Mon. Not. R. Astron. Soc. 508, 4970 (2021). https://doi.org/10.1093/mnras/stab2833. arXiv:2110.13184 [astro-ph.HE]
J. Antoniadis et al., Mon. Not. R. Astron. Soc. 510, 4873 (2022). https://doi.org/10.1093/mnras/stab3418. arXiv:2201.03980 [astro-ph.HE]
J. Ellis, M. Lewicki, Phys. Rev. Lett. 126, 041304 (2021). https://doi.org/10.1103/PhysRevLett.126.041304. arXiv:2009.06555 [astro-ph.CO]
S. Blasi, V. Brdar, K. Schmitz, Phys. Rev. Lett. 126, 041305 (2021). https://doi.org/10.1103/PhysRevLett.126.041305. arXiv:2009.06607 [astro-ph.CO]
V. Vaskonen, H. Veermäe, Phys. Rev. Lett. 126, 051303 (2021). https://doi.org/10.1103/PhysRevLett.126.051303. arXiv:2009.07832 [astro-ph.CO]
V. De Luca, G. Franciolini, A. Riotto, Phys. Rev. Lett. 126, 041303 (2021). https://doi.org/10.1103/PhysRevLett.126.041303. arXiv:2009.08268 [astro-ph.CO]
Y. Nakai, M. Suzuki, F. Takahashi, M. Yamada, Phys. Lett. B 816, 136238 (2021). https://doi.org/10.1016/j.physletb.2021.136238. arXiv:2009.09754 [astro-ph.CO]
W. Ratzinger, P. Schwaller, SciPost Phys. 10, 047 (2021). https://doi.org/10.21468/SciPostPhys.10.2.047. arXiv:2009.11875 [astro-ph.CO]
K. Kohri, T. Terada, Phys. Lett. B 813, 136040 (2021). https://doi.org/10.1016/j.physletb.2020.136040. arXiv:2009.11853 [astro-ph.CO]
S. Vagnozzi, Mon. Not. R. Astron. Soc. 502, L11 (2021). https://doi.org/10.1093/mnrasl/slaa203. arXiv:2009.13432 [astro-ph.CO]
A. Neronov, A. Roper Pol, C. Caprini, D. Semikoz, Phys. Rev. D 103, 041302 (2021). arXiv:2009.14174 [astro-ph.CO]
J.J. Blanco-Pillado, K.D. Olum, J.M. Wachter, Phys. Rev. D 103, 103512 (2021). https://doi.org/10.1103/PhysRevD.103.103512. arXiv:2102.08194 [astro-ph.CO]
D. Wang (2022). arXiv:2201.09295 [astro-ph.CO]
A. Roper Pol, C. Caprini, A. Neronov, D. Semikoz, Phys. Rev. D 105, 123502 (2022). https://doi.org/10.1103/PhysRevD.105.123502. arXiv:2201.05630 [astro-ph.CO]
R.Z. Ferreira, A. Notari, O. Pujolas, F. Rompineve (2022). arXiv:2204.04228 [astro-ph.CO]
C. Caprini et al., JCAP 04, 001 (2016). https://doi.org/10.1088/1475-7516/2016/04/001. arXiv:1512.06239 [astro-ph.CO]
C. Caprini et al., JCAP 03, 024 (2020). https://doi.org/10.1088/1475-7516/2020/03/024. arXiv:1910.13125 [astro-ph.CO]
P. Auclair et al., JCAP 04, 034 (2020). https://doi.org/10.1088/1475-7516/2020/04/034. arXiv:1909.00819 [astro-ph.CO]
N. Bartolo et al., JCAP 12, 026 (2016). https://doi.org/10.1088/1475-7516/2016/12/026. arXiv:1610.06481 [astro-ph.CO]
S.R. Coleman, Phys. Rev. D 15, 2929 (1977). https://doi.org/10.1103/PhysRevD.16.1248. [Erratum: Phys. Rev. D 16, 1248 (1977)]
A. Kosowsky, M.S. Turner, Phys. Rev. D 47, 4372 (1993). https://doi.org/10.1103/PhysRevD.47.4372. arXiv:astro-ph/9211004
M. Kamionkowski, A. Kosowsky, M.S. Turner, Phys. Rev. D 49, 2837 (1994). https://doi.org/10.1103/PhysRevD.49.2837. arXiv:astro-ph/9310044
M. Hindmarsh, M. Hijazi, JCAP 12, 062 (2019). https://doi.org/10.1088/1475-7516/2019/12/062. arXiv:1909.10040 [astro-ph.CO]
D. Cutting, M. Hindmarsh, D.J. Weir, Phys. Rev. Lett. 125, 021302 (2020). https://doi.org/10.1103/PhysRevLett.125.021302. arXiv:1906.00480 [hep-ph]
M. Lewicki, V. Vaskonen, Eur. Phys. J. C 80, 1003 (2020). https://doi.org/10.1140/epjc/s10052-020-08589-1. arXiv:2007.04967 [astro-ph.CO]
M. Lewicki, V. Vaskonen, Eur. Phys. J. C 81, 437 (2021). https://doi.org/10.1140/epjc/s10052-021-09232-3. arXiv:2012.07826 [astro-ph.CO]. [Erratum: Eur. Phys. J. C 81, 1077 (2021)]
R. Jinno, T. Konstandin, H. Rubira, JCAP 04, 014 (2021). https://doi.org/10.1088/1475-7516/2021/04/014. arXiv:2010.00971 [astro-ph.CO]
J. Dahl, M. Hindmarsh, K. Rummukainen, D.J. Weir, Phys. Rev. D 106, 063511 (2022). https://doi.org/10.1103/PhysRevD.106.063511. arXiv:2112.12013 [gr-qc]
D. Cutting, E. Vilhonen, D.J. Weir, Phys. Rev. D 106(10), 103524 (2022). https://doi.org/10.1103/PhysRevD.106.103524. arXiv:2204.03396 [astro-ph.CO]
R. Jinno, M. Takimoto, Phys. Rev. D 95, 015020 (2017). https://doi.org/10.1103/PhysRevD.95.015020. arXiv:1604.05035 [hep-ph]
S. Iso, P.D. Serpico, K. Shimada, Phys. Rev. Lett. 119, 141301 (2017). https://doi.org/10.1103/PhysRevLett.119.141301. arXiv:1704.04955 [hep-ph]
L. Marzola, A. Racioppi, V. Vaskonen, Eur. Phys. J. C 77, 484 (2017). https://doi.org/10.1140/epjc/s10052-017-4996-1. arXiv:1704.01034 [hep-ph]
T. Prokopec, J. Rezacek, B. Świeżewska, JCAP 02, 009 (2019). https://doi.org/10.1088/1475-7516/2019/02/009. arXiv:1809.11129 [hep-ph]
C. Marzo, L. Marzola, V. Vaskonen, Eur. Phys. J. C 79, 601 (2019). https://doi.org/10.1140/epjc/s10052-019-7076-x. arXiv:1811.11169 [hep-ph]
P. Baratella, A. Pomarol, F. Rompineve, JHEP 03, 100 (2019). https://doi.org/10.1007/JHEP03(2019)100. arXiv:1812.06996 [hep-ph]
B. Von Harling, A. Pomarol, O. Pujolàs, F. Rompineve, JHEP 04, 195 (2020). https://doi.org/10.1007/JHEP04(2020)195. arXiv:1912.07587 [hep-ph]
M. Aoki, J. Kubo, JCAP 04, 001 (2020). https://doi.org/10.1088/1475-7516/2020/04/001. arXiv:1910.05025 [hep-ph]
L. Delle Rose, G. Panico, M. Redi, A. Tesi, JHEP 04, 025 (2020). https://doi.org/10.1007/JHEP04(2020)025. arXiv:1912.06139 [hep-ph]
X. Wang, F.P. Huang, X. Zhang, JCAP 05, 045 (2020). https://doi.org/10.1088/1475-7516/2020/05/045. arXiv:2003.08892 [hep-ph]
J. Ellis, M. Lewicki, V. Vaskonen, JCAP 11, 020 (2020). https://doi.org/10.1088/1475-7516/2020/11/020. arXiv:2007.15586 [astro-ph.CO]
I. Baldes, Y. Gouttenoire, F. Sala, JHEP 04, 278 (2021). https://doi.org/10.1007/JHEP04(2021)278. arXiv:2007.08440 [hep-ph]
I. Baldes, Y. Gouttenoire, F. Sala, G. Servant, JHEP 07, 084 (2022). https://doi.org/10.1007/JHEP07(2022)084. arXiv:2110.13926 [hep-ph]
M. Lewicki, O. Pujolàs, V. Vaskonen, Eur. Phys. J. C 81, 857 (2021). https://doi.org/10.1140/epjc/s10052-021-09669-6. arXiv:2106.09706 [astro-ph.CO]
C. Grojean, G. Servant, Phys. Rev. D 75, 043507 (2007). https://doi.org/10.1103/PhysRevD.75.043507. arXiv:hep-ph/0607107
G.C. Dorsch, S.J. Huber, J.M. No, Phys. Rev. Lett. 113, 121801 (2014). https://doi.org/10.1103/PhysRevLett.113.121801. arXiv:1403.5583 [hep-ph]
P. Huang, A.J. Long, L.-T. Wang, Phys. Rev. D 94, 075008 (2016). https://doi.org/10.1103/PhysRevD.94.075008. arXiv:1608.06619 [hep-ph]
M. Artymowski, M. Lewicki, J.D. Wells, JHEP 03, 066 (2017). https://doi.org/10.1007/JHEP03(2017)066. arXiv:1609.07143 [hep-ph]
V. Vaskonen, Phys. Rev. D 95, 123515 (2017). https://doi.org/10.1103/PhysRevD.95.123515. arXiv:1611.02073 [hep-ph]
G.C. Dorsch, S.J. Huber, T. Konstandin, J.M. No, JCAP 05, 052 (2017). https://doi.org/10.1088/1475-7516/2017/05/052. arXiv:1611.05874 [hep-ph]
J. Ellis, M. Lewicki, J.M. No, JCAP 04, 003 (2019). https://doi.org/10.1088/1475-7516/2019/04/003. arXiv:1809.08242 [hep-ph]
A. Beniwal, M. Lewicki, M. White, A.G. Williams, JHEP 02, 183 (2019). https://doi.org/10.1007/JHEP02(2019)183. arXiv:1810.02380 [hep-ph]
M. Fairbairn, E. Hardy, A. Wickens, JHEP 07, 044 (2019). https://doi.org/10.1007/JHEP07(2019)044. arXiv:1901.11038 [hep-ph]
J. Ellis, M. Lewicki, J.M. No, V. Vaskonen, JCAP 06, 024 (2019). https://doi.org/10.1088/1475-7516/2019/06/024. arXiv:1903.09642 [hep-ph]
M. Lewicki, M. Merchand, M. Zych, JHEP 02, 017 (2022). https://doi.org/10.1007/JHEP02(2022)017. arXiv:2111.02393 [astro-ph.CO]
R.-G. Cai, S.-J. Wang, JCAP 03, 096 (2021). https://doi.org/10.1088/1475-7516/2021/03/096. arXiv:2011.11451 [astro-ph.CO]
M. Lewicki, V. Vaskonen, H. Veermäe, Phys. Rev. D 106(10), 103501 (2022). https://doi.org/10.1103/PhysRevD.106.103501. arXiv:2205.05667 [astro-ph.CO]
D. Bodeker, G.D. Moore, JCAP 05, 009 (2009). https://doi.org/10.1088/1475-7516/2009/05/009. arXiv:0903.4099 [hep-ph]
D. Bodeker, G.D. Moore, JCAP 05, 025 (2017). https://doi.org/10.1088/1475-7516/2017/05/025. arXiv:1703.08215 [hep-ph]
A. Azatov, M. Vanvlasselaer, JCAP 01, 058 (2021). https://doi.org/10.1088/1475-7516/2021/01/058. arXiv:2010.02590 [hep-ph]
Y. Gouttenoire, R. Jinno, F. Sala, JHEP 05, 004 (2022). https://doi.org/10.1007/JHEP05(2022)004. arXiv:2112.07686 [hep-ph]
S. Höche, J. Kozaczuk, A.J. Long, J. Turner, Y. Wang, JCAP 03, 009 (2021). https://doi.org/10.1088/1475-7516/2021/03/009. arXiv:2007.10343 [hep-ph]
M. Barroso Mancha, T. Prokopec, B. Swiezewska, JHEP 01, 070 (2021). https://doi.org/10.1007/JHEP01(2021)070. arXiv:2005.10875 [hep-th]
J.R. Espinosa, T. Konstandin, J.M. No, G. Servant, JCAP 06, 028 (2010). https://doi.org/10.1088/1475-7516/2010/06/028. arXiv:1004.4187 [hep-ph]
R. Jinno, H. Seong, M. Takimoto, C.M. Um, JCAP 10, 033 (2019). https://doi.org/10.1088/1475-7516/2019/10/033. arXiv:1905.00899 [astro-ph.CO]
M. Hindmarsh, S.J. Huber, K. Rummukainen, D.J. Weir, Phys. Rev. D 96, 103520 (2017). https://doi.org/10.1103/PhysRevD.96.103520. arXiv:1704.05871 [astro-ph.CO]. [Erratum: Phys. Rev. D 101, 089902 (2020)]
M. Hindmarsh, S.J. Huber, K. Rummukainen, D.J. Weir, Phys. Rev. D 92, 123009 (2015). https://doi.org/10.1103/PhysRevD.92.123009. arXiv:1504.03291 [astro-ph.CO]
A. Kurganov, E. Tadmor, J. Comput. Phys. 160, 241 (2000). https://doi.org/10.1006/jcph.2000.6459
T. Konstandin, JCAP 03, 047 (2018). https://doi.org/10.1088/1475-7516/2018/03/047. arXiv:1712.06869 [astro-ph.CO]
M. Hindmarsh, Phys. Rev. Lett. 120, 071301 (2018). https://doi.org/10.1103/PhysRevLett.120.071301. arXiv:1608.04735 [astro-ph.CO]
A. Roper Pol, S. Mandal, A. Brandenburg, T. Kahniashvili, A. Kosowsky, Phys. Rev. D 102, 083512 (2020). https://doi.org/10.1103/PhysRevD.102.083512. arXiv:1903.08585 [astro-ph.CO]
T. Kahniashvili, A. Brandenburg, G. Gogoberidze, S. Mandal, A. Roper Pol, Phys. Rev. Res. 3, 013193 (2021). https://doi.org/10.1103/PhysRevResearch.3.013193. arXiv:2011.05556 [astro-ph.CO]
A. Roper Pol, S. Mandal, A. Brandenburg, T. Kahniashvili, JCAP 04, 019 (2022). https://doi.org/10.1088/1475-7516/2022/04/019. arXiv:2107.05356 [gr-qc]
P. Auclair, C. Caprini, D. Cutting, M. Hindmarsh, K. Rummukainen, D.A. Steer, D.J. Weir, JCAP 09, 029 (2022). https://doi.org/10.1088/1475-7516/2022/09/029. arXiv:2205.02588 [astro-ph.CO]
R. Allahverdi et al. (2020). https://doi.org/10.21105/astro.2006.16182. arXiv:2006.16182 [astro-ph.CO]
N. Aghanim et al. (Planck), Astron. Astrophys. 641, A6 (2020). https://doi.org/10.1051/0004-6361/201833910. arXiv:1807.06209 [astro-ph.CO]. [Erratum: Astron. Astrophys. 652, C4 (2021)]
C. Caprini, R. Durrer, T. Konstandin, G. Servant, Phys. Rev. D 79, 083519 (2009). https://doi.org/10.1103/PhysRevD.79.083519. arXiv:0901.1661 [astro-ph.CO]
R.-G. Cai, S. Pi, M. Sasaki, Phys. Rev. D 102, 083528 (2020). https://doi.org/10.1103/PhysRevD.102.083528. arXiv:1909.13728 [astro-ph.CO]
G. Barenboim, W.-I. Park, Phys. Lett. B 759, 430 (2016). https://doi.org/10.1016/j.physletb.2016.06.009. arXiv:1605.03781 [astro-ph.CO]
A. Hook, G. Marques-Tavares, D. Racco, JHEP 02, 117 (2021). https://doi.org/10.1007/JHEP02(2021)117. arXiv:2010.03568 [hep-ph]
Y. Gouttenoire, G. Servant, P. Simakachorn (2021). arXiv:2111.01150 [hep-ph]
Acknowledgements
This work was supported by the Spanish MINECO grants FPA2017-88915-P and SEV-2016-0588, the Spanish MICINN grants IJC2019-041533-I and PID2020-115845GB-I00/AEI/10.13039/501100011033, the grant 2017-SGR-1069 from the Generalitat de Catalunya, the Polish National Science Center grant 2018/31/D/ST2/02048, and the Polish National Agency for Academic Exchange within Polish Returns Programme under agreement PPN/PPO/2020/1/00013/U/00001. IFAE is partially funded by the CERCA program of the Generalitat de Catalunya.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Funded by SCOAP3. SCOAP3 supports the goals of the International Year of Basic Sciences for Sustainable Development.
About this article
Cite this article
Lewicki, M., Vaskonen, V. Gravitational waves from bubble collisions and fluid motion in strongly supercooled phase transitions. Eur. Phys. J. C 83, 109 (2023). https://doi.org/10.1140/epjc/s10052-023-11241-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-023-11241-3