We present a very minimal model for baryogenesis by a dark first-order phase transition. It employs a new dark $SU(2)_{D}$ gauge group with two doublet Higgs bosons, two lepton doublets, and two singlets. The singlets act as a neutrino portal that transfer the generated asymmetry to the Standard Model. The model predicts $\Delta N_\text{eff} = 0.09-0.13$ detectable by future experiments as well as possible signals from exotic decays of the Higgs and $Z$ bosons and stochastic gravitational waves.

We study the energy budget of a first-order cosmological phase transition, which is an important factor in the prediction of the resulting gravitational wave spectrum. Formerly, this analysis was based mostly on simplified models as for example the bag equation of state. Here, we present a model-independent approach that is exact up to the temperature dependence of the speed of sound in the broken phase. We find that the only relevant quantities that enter in the hydrodynamic analysis are the speed of sound in the broken phase and a linear combination of the energy and pressure differences between the two phases which we call pseudotrace (normalized to the enthalpy in the broken phase). The pseudotrace quantifies the strength of the phase transition and yields the conventional trace of the energy-momentum tensor for a relativistic plasma (with speed of sound squared of one third). We study this approach in several realistic models of the phase transition and also provide a code snippet that can be used to determine the efficiency coefficient for a given phase transition strength and speed of sound. It turns out that our approach is accurate to the percent level for moderately strong phase transitions, while former approaches give at best the right order of magnitude.

We study the Schr\"odinger-Poisson (SP) method in the context of cosmological large-scale structure formation in an expanding background. In the limit $\hbar \to 0$, the SP technique can be viewed as an effective method to sample the phase space distribution of cold dark matter that remains valid on non-linear scales. We present results for the 2D and 3D matter correlation function and power spectrum at length scales corresponding to the baryon acoustic oscillation (BAO) peak. We discuss systematic effects of the SP method applied to cold dark matter and explore how they depend on the simulation parameters. In particular, we identify a combination of simulation parameters that controls the scale-independent loss of power observed at low redshifts, and discuss the scale relevant to this effect.

We study the power spectrum of dark matter density fluctuations in the framework of the Effective Field Theory of Large Scale Structures (EFTofLSS) up to three loop orders. In principle, several counter-terms may be needed to handle the short-distance sensitivity in perturbation theory. However, we show that a small number of extra coefficients are sufficient to match numerical simulations with percent accuracy when a generic renormalization prescription is implemented (allowing for running of the individual counter-terms). We show that the level of accuracy increases with respect to the two loop results, up to $k \simeq 0.4\, h$Mpc$^{-1}$ at redshift $z=0$, although the overall improvement is somewhat marginal. At the same time, we argue there is evidence that the behavior of the loop expansion in the EFTofLSS is typical of an asymptotic series, already on the brink of its maximum predictive power (at $z=0$). Hence, the inclusion of higher orders will likely deteriorate the matching to data, even at moderate values of $k$. Part of the reason for this behavior is due to large contributions to the (renormalized) power spectrum at three loop order from mildly non-linear scales, even after the UV counter-terms are included. In conclusion, the EFTofLSS to three loop orders provides the best approximation to the (deterministic part of the) power spectrum in the weakly non-linear regime at $z=0$, and higher loops are not expected to improve our level of accuracy.

We introduce a model for matters-genesis in which both the baryonic and dark matter asymmetries originate from a first-order phase transition in a dark sector with an $SU(3) \times SU(2) \times U(1)$ gauge group and minimal matter content. In the simplest scenario, we predict that dark matter is a dark neutron with mass either $m_n = 1.33$ GeV or $m_n = 1.58$ GeV. Alternatively, dark matter may be comprised of equal numbers of dark protons and pions. This model, in either scenario, is highly discoverable through both dark matter direct detection and dark photon search experiments. The strong dark matter self interactions may ameliorate small-scale structure problems, while the strongly first-order phase transition may be confirmed at future gravitational wave observatories.

We investigate the potential for observing gravitational waves from cosmological phase transitions with LISA in light of recent theoretical and experimental developments. Our analysis is based on current state-of-the-art simulations of sound waves in the cosmic fluid after the phase transition completes. We discuss the various sources of gravitational radiation, the underlying parameters describing the phase transition and a variety of viable particle physics models in this context, clarifying common misconceptions that appear in the literature and identifying open questions requiring future study. We also present a web-based tool, PTPlot, that allows users to obtain up-to-date detection prospects for a given set of phase transition parameters at LISA.

We discuss vacuum stability in Froggatt-Nielsen (FN) models. One concern in FN models is that for large flavon VEVs the running of the quartic Higgs coupling is enhanced what might lead to a more severe instability compared to the Standard Model (SM). We study this issue using the renormalization-group improved scalar potential. Another issue is that the mixing between the Higgs and the flavon can potentially destabilize the potential. However, taking current bounds on the flavon phenomenology into account, we find that both effects do not lead to an instability that is more severe than in the SM.

We study scalar bubble collisions in first-order phase transitions focusing on the relativistic limit. We propose 'trapping equation' which describes the wall behavior after collision, and test it with numerical simulations in several setups. We also examine the energy dynamics after collision and discuss its implications to gravitational wave production.

We perform simulations in a simple model that aims to mimic the hydrodynamic evolution of a relativistic fluid during a cosmological first-order phase transitions. The observable we are concerned with is hereby the spectrum of gravitational radiation produced by colliding fluid shells. We present simple parameterizations of our results as functions of the wall velocity, the duration of the phase transition and the latent heat. We also improve on previous results in the envelope approximation and compare with hydrodynamic simulations.

The Lyman-$\alpha$ forest is a valuable probe of dark matter models featuring a scale-dependent suppression of the power spectrum as compared to $\Lambda$CDM. In this work, we present a new estimator of the Lyman-$\alpha$ flux power spectrum that does not rely on hydrodynamical simulations. Our framework is characterized by nuisance parameters that encapsulate the complex physics of the intergalactic medium and sensitivity to highly non-linear small-scale modes. After validating the approach based on high-resolution hydrodynamical simulations for $\Lambda$CDM, we derive conservative constraints on interacting dark matter models from BOSS Lyman-$\alpha$ data on large scales, k<0.02(km/s)^(-1), with the relevant nuisance parameters left free in the model fit. The estimator yields lower bounds on the mass of cannibal dark matter, where freeze-out occurs through 3-to-2 annihilation, in the MeV range. Furthermore, we find that models of dark matter interacting with dark radiation, which have been argued to address the $H_0$ and $\sigma_8$ tensions, are compatible with BOSS Lyman-$\alpha$ data.

We perform a quantitative comparison between N-body simulations and the Schr\"odinger-Poisson system in 1+1 dimensions. In particular, we study halo formation with different initial conditions. We observe the convergence of various observables in the Planck constant h-bar and also test virialization. We discuss the generation of higher order cumulants of the particle distribution function which demonstrates that the Schr\"odinger-Poisson equations should not be perceived as a generalization of the dust model with quantum pressure but rather as one way of sampling the phase space of the Vlasov-Poisson system -- just as N-body simulations. Finally, we quantitatively recover the scaling behavior of the halo density profile from N-body simulations.

We show that simple Two Higgs Doublet models still provide a viable explanation for the matter-antimatter asymmetry of the Universe via electroweak baryogenesis, even after taking into account the recent order-of-magnitude improvement on the electron-EDM experimental bound by the ACME Collaboration. Moreover we show that, in the region of parameter space where baryogenesis is possible, the gravitational wave spectrum generated at the end of the electroweak phase transition is within the sensitivity reach of the future space-based interferometer LISA.

The measured (central) values of the Higgs and top quark masses indicate that the Standard Model (SM) effective potential develops an instability at high field values. The scale of this instability, determined as the Higgs field value at which the potential drops below the electroweak minimum, is about $10^{11}$ GeV. However, such a scale is unphysical as it is not gauge-invariant and suffers from a gauge-fixing uncertainty of up to two orders of magnitude. Subjecting our system, the SM, to several probes of the instability (adding higher order operators to the potential; letting the vacuum decay through critical bubbles; heating up the system to very high temperature; inflating it) and asking in each case physical questions, we are able to provide several gauge-invariant scales related with the Higgs potential instability.

We study the general NMSSM with an emphasis on the parameter regions with a very strong first-order electroweak phase transition (EWPT). In the presence of heavy fields coupled to the Higgs sector, the analysis can be problematic due to the existence of sizable radiative corrections. In this paper we propose a subtraction scheme that helps to circumvent this problem. For simplicity we focus on a parameter region that is by construction hidden from the current collider searches. The analysis proves that (at least) in the identified parameter region the EWPT can be very strong and striking gravitational wave signals can be produced. The corresponding gravitational stochastic background can potentially be detected at the planned space-based gravitational wave observatory eLISA, depending on the specific experiment design that will be approved.

We investigate the potential for the eLISA space-based interferometer to detect the stochastic gravitational wave background produced by strong first-order cosmological phase transitions. We discuss the resulting contributions from bubble collisions, magnetohydrodynamic turbulence, and sound waves to the stochastic background, and estimate the total corresponding signal predicted in gravitational waves. The projected sensitivity of eLISA to cosmological phase transitions is computed in a model-independent way for various detector designs and configurations. By applying these results to several specific models, we demonstrate that eLISA is able to probe many well-motivated scenarios beyond the Standard Model of particle physics predicting strong first-order cosmological phase transitions in the early Universe.

We derive a non-perturbative equation for the large scale structure power spectrum of long-wavelength modes. Thereby, we use an operator product expansion together with relations between the three-point function and power spectrum in the soft limit. The resulting equation encodes the coupling to ultraviolet (UV) modes in two time-dependent coefficients, which may be obtained from response functions to (anisotropic) parameters, such as spatial curvature, in a modified cosmology. We argue that both depend weakly on fluctuations deep in the UV. As a byproduct, this implies that the renormalized leading order coefficient(s) in the effective field theory (EFT) of large scale structures receive most of their contribution from modes close to the non-linear scale. Consequently, the UV dependence found in explicit computations within standard perturbation theory stems mostly from counter-term(s). We confront a simplified version of our non-perturbative equation against existent numerical simulations, and find good agreement within the expected uncertainties. Our approach can in principle be used to precisely infer the relevance of the leading order EFT coefficient(s) using small volume simulations in an `anisotropic separate universe' framework. Our results suggest that the importance of these coefficient(s) is a $\sim 10 \%$ effect, and plausibly smaller.

We study soft limits of correlation functions for the density and velocity fields in the theory of structure formation. First, we re-derive the (resummed) consistency conditions at unequal times using the eikonal approximation. These are solely based on symmetry arguments and are therefore universal. Then, we explore the existence of equal-time relations in the soft limit which, on the other hand, depend on the interplay between soft and hard modes. We scrutinize two approaches in the literature: The time-flow formalism, and a background method where the soft mode is absorbed into a locally curved cosmology. The latter has been recently used to set up (angular averaged) `equal-time consistency relations'. We explicitly demonstrate that the time-flow relations and `equal-time consistency conditions' are only fulfilled at the linear level, and fail at next-to-leading order for an Einstein de-Sitter universe. While applied to the velocities both proposals break down beyond leading order, we find that the `equal-time consistency conditions' quantitatively approximates the perturbative results for the density contrast. Thus, we generalize the background method to properly incorporate the effect of curvature in the density and velocity fluctuations on short scales, and discuss the reasons behind this discrepancy. We conclude with a few comments on practical implementations and future directions.

We compute non-linear corrections to the matter power spectrum taking the time- and scale-dependent free-streaming length of neutrinos into account. We adopt a hybrid scheme that matches the full Boltzmann hierarchy to an effective two-fluid description at an intermediate redshift. The non-linearities in the neutrino component are taken into account by using an extension of the time-flow framework. We point out that this remedies a spurious behaviour that occurs when neglecting non-linear terms for neutrinos. This behaviour is related to how efficiently short modes decouple from long modes and can be traced back to the violation of momentum conservation if neutrinos are treated linearly. Furthermore, we compare our results at next to leading order to various other methods and quantify the accuracy of the fluid description. Due to the correct decoupling behaviour of short modes, the two-fluid scheme is a suitable starting point to compute higher orders in perturbations or for resummation methods.

By means of a relativistic microscopic approach we calculate the expansion velocity of bubbles generated during a first-order electroweak phase transition. In particular, we use the gradient expansion of the Kadanoff-Baym equations to set up the fluid system. This turns out to be equivalent to the one found in the semi-classical approach in the non-relativistic limit. Finally, by including hydrodynamic deflagration effects and solving the Higgs equations of motion in the fluid, we determine velocity and thickness of the bubble walls. Our findings are compared with phenomenological models of wall velocities. As illustrative examples, we apply these results to three theories providing first-order phase transitions with a particle content in the thermal plasma that resembles the Standard Model.

We analyze the dark matter power spectrum at three-loop order in standard perturbation theory of large scale structure. We observe that at late times the loop expansion does not converge even for large scales (small momenta) well within the linear regime, but exhibits properties compatible with an asymptotic series. We propose a technique to restore the convergence in the limit of small momentum, and use it to obtain a perturbative expansion with improved convergence for momenta in the range where baryonic acoustic oscillations are present. Our results are compared with data from N-body simulations at different redshifts, and we find good agreement within this range.

We discuss the convergence of cosmological perturbation theory. We prove that the polynomial enhancement of the non-linear corrections expected from the effects of soft modes is absent in equal-time correlators like the power or bispectrum. We first show this at leading order by resumming the most important corrections of soft modes to an arbitrary skeleton of hard fluctuations. We derive the same result in the eikonal approximation, which also allows us to show the absence of enhancement at any order. We complement the proof by an explicit calculation of the power spectrum at two-loop order, and by further numerical checks at higher orders. Using these insights, we argue that the modification of the power spectrum from soft modes corresponds at most to logarithmic corrections. Finally, we discuss the asymptotic behavior in the large and small momentum regimes and identify the expansion parameter pertinent to non-linear corrections.

We address electroweak baryogenesis in the context of composite Higgs models, pointing out that modifications to the Higgs and top quark sectors can play an important role in generating the baryon asymmetry. Our main observation is that composite Higgs models that include a light, gauge singlet scalar in the spectrum [as in the model based on the symmetry breaking pattern SO(6)/SO(5)], provide all necessary ingredients for viable baryogenesis. In particular, the singlet leads to a strongly first-order electroweak phase transition and introduces new sources of CP violation in dimension-five operators involving the top quark. We discuss the amount of baryon asymmetry produced and the experimental constraints on the model.

It is well known that the electroweak phase transition (EWPhT) in extensions of the Standard Model with one real scalar singlet can be first-order for realistic values of the Higgs mass. We revisit this scenario with the most general renormalizable scalar potential systematically identifying all regions in parameter space that develop, due to tree-level dynamics, a potential barrier at the critical temperature that is strong enough to avoid sphaleron wash-out of the baryon asymmetry. Such strong EWPhTs allow for a simple mean-field approximation and an analytic treatment of the free-energy that leads to very good theoretical control and understanding of the different mechanisms that can make the transition strong. We identify a new realization of such mechanism, based on a flat direction developing at the critical temperature, which could operate in other models. Finally, we discuss in detail some special cases of the model performing a numerical calculation of the one-loop free-energy that improves over the mean-field approximation and confirms the analytical expectations.

We discuss a hydrodynamic obstruction to bubble wall acceleration during a cosmological first-order phase transition. The obstruction results from the heating of the plasma in the compression wave in front of the phase transition boundary. We provide a simple criterion for the occurrence of the obstruction at subsonic bubble wall velocity in terms of the critical temperature, the phase transition temperature, and the latent heat of the model under consideration. The criterion serves as a sufficient condition of subsonic bubble wall velocities as required by electroweak baryogenesis.

The study of the hydrodynamics of bubble growth in first-order phase transitions is very relevant for electroweak baryogenesis, as the baryon asymmetry depends sensitively on the bubble wall velocity, and also for predicting the size of the gravity wave signal resulting from bubble collisions, which depends on both the bubble wall velocity and the plasma fluid velocity. We perform such study in different bubble expansion regimes, namely deflagrations, detonations, hybrids (steady states) and runaway solutions (accelerating wall), without relying on a specific particle physics model. We compute the efficiency of the transfer of vacuum energy to the bubble wall and the plasma in all regimes. We clarify the condition determining the runaway regime and stress that in most models of strong first-order phase transitions this will modify expectations for the gravity wave signal. Indeed, in this case, most of the kinetic energy is concentrated in the wall and almost no turbulent fluid motions are expected since the surrounding fluid is kept mostly at rest.

We discuss the question if the upcoming generation of collider and low-energy experiments can successfully probe the nature of the electroweak phase transition. In particular, we are interested in phase transitions strong enough for electroweak baryogenesis or even for a production of gravitational radiation observable by the Big Bang Observer. As an explicit example, we present an analysis in a singlet extension of the Standard Model. We focus on the region in parameter space where the model develops no significant deviation in its low energy phenomenology from the Standard Model. Nevertheless, this class of models can develop a very strong phase transition.

In this paper we discuss some general aspects of the gravitational wave background arising from post-inflationary short-lasting cosmological events such as phase transitions. We concentrate on the physics which determines the shape and the peak frequency of the gravitational wave spectrum. We then apply our general findings to the case of bubble collisions during a first order phase transition and compare different results in the recent literature.

Using the worldline method, we derive an effective action of the bosonic sector of the Standard Model by integrating out the fermionic degrees of freedom. The CP violation stemming from the complex phase in the CKM matrix gives rise to CP-violating operators in the one-loop effective action in the next-to-leading order of a gradient expansion. We calculate the prefactor of the appropriate operators and give general estimates of CP violation in the bosonic sector of the Standard Model. In particular, we show that the effective CP violation for weak gauge fields is not suppressed by the Yukawa couplings of the light quarks and is much larger than the bound given by the Jarlskog determinant.

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.

We reexamine the production of gravitational waves by bubble collisions during a first-order phase transition. The spectrum of the gravitational radiation is determined by numerical simulations using the "envelope approximation". We find that the spectrum rises as f^3.0 for small frequencies and decreases as f^-1.0 for high frequencies. Thus, the fall-off at high frequencies is significantly slower than previously stated in the literature. This result has direct impact on detection prospects for gravity waves originating from a strong first-order electroweak phase transition at space-based interferometers, such as LISA or BBO. In addition, we observe a slight dependence of the peak frequency on the bubble wall velocity.

During a strongly first-order phase transition gravitational waves are produced by bubble collisions and turbulent plasma motion. We analyze the relevant characteristics of the electroweak phase transition in the nMSSM to determine the generated gravitational wave signal. Additionally, we comment on correlations between the production of gravitational waves and baryogenesis. We conclude that the gravitational wave relic density in this model is generically too small to be detected in the near future by the LISA experiment. We also consider the case of a "Standard Model" with dimension-six Higgs potential, which leads to a slightly stronger signal of gravitational waves.

We analyze the left-right symmetric type I+II seesaw mechanism, where an eight-fold degeneracy among the mass matrices of heavy right-handed neutrinos M_R is known to exist. Using the stability property of the solutions and their ability to lead to successful baryogenesis via leptogenesis as additional criteria, we discriminate among these eight solutions and partially lift their eight-fold degeneracy. In particular, we find that viable leptogenesis is generically possible for four out of the eight solutions.