We undertake a careful analysis of stochastic gravitational wave production from cosmological phase transitions in an expanding universe, studying both a standard radiation as well as a matter dominated history. We analyze in detail the dynamics of the phase transition, including the false vacuum fraction, bubble lifetime distribution, bubble number density, mean bubble separation, etc., for an expanding universe. We also study the full set of differential equations governing the evolution of plasma and the scalar field during the phase transition and generalize results obtained in Minkowski spacetime. In particular, we generalize the sound shell model to the expanding universe and determine the velocity field power spectrum. This ultimately provides an accurate calculation of the gravitational wave spectrum seen today for the dominant source of sound waves. For the amplitude of the gravitational wave spectrum visible today, we find a suppression factor arising from the finite lifetime of the sound waves and compare with the commonly used result in the literature, which corresponds to the asymptotic value of our suppression factor. We point out that the asymptotic value is only applicable for a very long lifetime of the sound waves, which is highly unlikely due to the onset of shocks, turbulence and other damping processes. We also point out that features of the gravitational wave spectral form may hold the tantalizing possibility of distinguishing between different expansion histories using phase transitions.

We present a dedicated complementarity study of gravitational wave and collider measurements of the simplest extension of the Higgs sector: the singlet scalar augmented Standard Model. We study the following issues: (i) the electroweak phase transition patterns admitted by the model, and the proportion of parameter space for each pattern; (ii) the regions of parameter space that give detectable gravitational waves at future space-based detectors; and (iii) the current and future collider measurements of di-Higgs production, as well as searches for a heavy weak diboson resonance, and how these searches interplay with regions of parameter space that exhibit strong gravitational wave signals. We carefully investigate the behavior of the normalized energy released during the phase transition as a function of the model parameters, address subtle issues pertaining to the bubble wall velocity, and provide a description of different fluid velocity profiles. On the collider side, we identify the subset of points that are most promising in terms of di-Higgs and weak diboson production studies while also giving detectable signals at LISA, setting the stage for future benchmark points that can be used by both communities.

Models of particle physics that feature phase transitions typically provide predictions for stochastic gravitational wave signals at future detectors and such predictions are used to delineate portions of the model parameter space that can be constrained. The question is: how precise are such predictions? Uncertainties enter in the calculation of the macroscopic thermal parameters and the dynamics of the phase transition itself. We calculate such uncertainties with increasing levels of sophistication in treating the phase transition dynamics. Currently, the highest level of diligence corresponds to careful treatments of the source lifetime; mean bubble separation; going beyond the bag model approximation in solving the hydrodynamics equations and explicitly calculating the fraction of energy in the fluid from these equations rather than using a fit; and including fits for the energy lost to vorticity modes and reheating effects. The lowest level of diligence incorporates none of these effects. We compute the percolation and nucleation temperatures, the mean bubble separation, the fluid velocity, and ultimately the gravitational wave spectrum corresponding to the level of highest diligence for three explicit examples: SMEFT, a dark sector Higgs model, and the real singlet-extended Standard Model (xSM). In each model, we contrast different levels of diligence in the calculation and find that the difference in the final predicted signal can be several orders of magnitude. Our results indicate that calculating the gravitational wave spectrum for particle physics models and deducing precise constraints on the parameter space of such models continues to remain very much a work in progress and warrants care.

No abstract

We study the properties of Bose-Einstein Condensate (BEC) systems consisting of two scalars, focusing on both the case where the BEC is stellar scale as well as the case when it is galactic scale. After studying the stability of such systems and making contact with existing single scalar limits, we undertake a numerical study of the two interacting scalars using Einstein-Klein-Gordon (EKG) equations, including both non-gravitational self-interactions and interactions between the species. We show that the presence of extra scalars and possible interactions between them can leave unique imprints on the BEC system mass profile, especially when the system transitions from being dominated by one scalar to being dominated by the other. At stellar scales (nonlinear regime), we observe that a repulsive interaction between the two scalars of the type +ϕ1 2 ϕ2 2 can stabilize the BEC system and support it up to high compactness, a phenomenon only known to exist in the +ϕ4 system. We provide simple analytic understanding of this behavior and point out that it can lead to interesting gravitational wave signals at LIGO-Virgo. At galactic scales, on the other hand, we show that two-scalar BECs can address the scaling problem that arises when one uses ultralight dark matter mass profiles to fit observed galactic core mass profiles. In the end, we construct a particle model of two ultralight scalars with the repulsive ϕ1 2 ϕ2 2 interaction using collective symmetry breaking. We develop a fast numerical code that utilizes the relaxation method to solve the EKG system, which can be easily generalized to multiple scalars.