Many models of physics beyond the Standard Model predict a strong first-order phase transition (SFOPT) in the early Universe that leads to observable gravitational waves (GWs). In this paper, we propose a novel method for presenting and comparing the GW signals that are predicted by different models. Our approach is based on the observation that the GW signal has an approximately model-independent spectral shape. This allows us to represent it solely in terms of a finite number of observables, that is, a set of peak amplitudes and peak frequencies. As an example, we consider the GW signal in the real-scalar-singlet extension of the Standard Model (xSM). We construct the signal region of the xSM in the space of observables and show how it will be probed by future space-borne interferometers. Our analysis results in sensitivity plots that are reminiscent of similar plots that are typically shown for dark-matter direct-detection experiments, but which are novel in the context of GWs from a SFOPT. These plots set the stage for a systematic model comparison, the exploration of underlying model-parameter dependencies, and the construction of distribution functions in the space of observables. In our plots, the experimental sensitivities of future searches for a stochastic GW signal are indicated by peak-integrated sensitivity curves. A detailed discussion of these curves, including fit functions, is contained in a companion paper [2002.04615]. The data and code that we used in our analysis can be downloaded from Zenodo [https://doi.org/10.5281/zenodo.3699415].

We consider inflation within a model framework where the Higgs boson arises as a pseudo-Goldstone boson associated with the breaking of a global symmetry at a scale significantly larger than the electroweak one. We show that in such a model the scalar self-couplings can be parametrically suppressed and, consequently, the non-minimal couplings to gravity can be of order one or less, while the inflationary predictions of the model remain compatible with the precision cosmological observations. Furthermore, in the model we study, the existence of the electroweak scale is entirely due to the inflaton field. Our model therefore suggests that inflation and low energy particle phenomenology may be more entwined than assumed so far.

We investigate an extension of the Standard Model containing two Higgs doublets and a singlet scalar field (2HDSM). We show that the model can have a strongly first-order phase transition and give rise to the observed baryon asymmetry of the Universe, consistent with all experimental constraints. In particular, the constraints from the electron and neutron electric dipole moments are less constraining here than in pure two-Higgs-doublet model (2HDM). The two-step, first-order transition in 2HDSM, induced by the singlet field, may lead to strong supercooling and low nucleation temperatures in comparison with the critical temperature, $T_n \ll T_c$, which can significantly alter the usual phase-transition pattern in 2HD models with $T_n \approx T_c$. Furthermore, the singlet field can be the dark matter particle. However, in models with a strong first-order transition its abundance is typically but a thousandth of the observed dark matter abundance.

Motivated by the possibility to explain dark matter abundance and strong electroweak phase transition, we consider simple extensions of the Standard Model containing singlet fields coupled with the Standard Model via a scalar portal. Concretely, we consider a basic portal model consisting of a singlet scalar with $Z_2$ symmetry and a model containing a singlet fermion connected with the Standard Model fields via a singlet scalar portal. We perform a Monte Carlo analysis of the parameter space of each model, and we find that in both cases the dark matter abundance can be produced either via freeze-out or freeze-in mechanisms, but only in the latter model one can obtain also a strong electroweak phase transition required by the successful electroweak baryogenesis. We impose the direct search limits and consider systematically the possibility that the model produces only a subdominant portion of the dark matter abundance. We also study the renormalization group evolution of the couplings of the model to determine if the scalar sector of the model remains stable and perturbative up to high scales. With explicit examples of benchmark values of the couplings at weak scale, we show that this is possible. Models of this type are further motivated by the possibility that the excursions of the Higgs field at the end of inflation are large and could directly probe the instability region of the Standard Model.