We study gravitational wave (GW) production in strongly supercooled cosmological phase transitions, taking particular care of models featuring a complex scalar field with a U$(1)$ symmetric potential. We perform lattice simulations of two-bubble collisions to properly model the scalar field gradients, and compute the GW spectrum sourced by them using thin-wall approximation in many-bubble simulations. We find that in the U$(1)$ symmetric case the low-frequency spectrum is $\propto\omega$ whereas for a real scalar field it is $\propto\omega^3$. In both cases the spectrum decays as $\omega^{-2}$ at high frequencies.

It is commonly assumed that the energy density of the Universe was dominated by radiation between reheating after inflation and the onset of matter domination 54,000 years later. While the abundance of light elements indicates that the Universe was radiation dominated during Big Bang Nucleosynthesis (BBN), there is scant evidence that the Universe was radiation dominated prior to BBN. It is therefore possible that the cosmological history was more complicated, with deviations from the standard radiation domination during the earliest epochs. Indeed, several interesting proposals regarding various topics such as the generation of dark matter, matter-antimatter asymmetry, gravitational waves, primordial black holes, or microhalos during a nonstandard expansion phase have been recently made. In this paper, we review various possible causes and consequences of deviations from radiation domination in the early Universe - taking place either before or after BBN - and the constraints on them, as they have been discussed in the literature during the recent years.

Atom interferometers (AIs) on earth and in space offer good capabilities for measuring gravitational waves (GWs) in the mid-frequency deciHz band, complementing the sensitivities of the LIGO/Virgo and LISA experiments and enabling probes of possible modifications of the general relativity predictions for GW propagation. We illustrate these capabilities using the projected sensitivities of the AION (terrestrial) and AEDGE (space-based) AI projects, showing that AION could improve the present LIGO/Virgo direct limit on the graviton mass by a factor $\sim 40$ to $\simeq 10^{-24}\,$eV, and AEDGE could improve the limit by another order of magnitude. AION and AEDGE will also have greater sensitivity than LIGO to some scenarios for Lorentz violation.

We study strongly supercooled cosmological phase transitions. We perform numerical lattice simulations of two-bubble collisions and demonstrate that, depending on the scalar potential, in the collision the field can either bounce to a false vacuum or remain oscillating around the true vacuum. We study if these cases can be distinguished from their gravitational wave signals and discuss the possibility of black hole formation in the bubble collisions.

We discuss the observability of circular polarisation of the stochastic gravitational-wave background (SGWB) generated by helical turbulence following a first-order cosmological phase transition, using a model that incorporates the effects of both direct and inverse energy cascades. We explore the strength of the gravitational-wave signal and the dependence of its polarisation on the helicity fraction, $\zeta_*$, the strength of the transition, $\alpha$, the bubble size, $R_*$, and the temperature, $T_*$, at which the transition finishes. We calculate the prospective signal-to-noise ratios of the SGWB strength and polarisation signals in the LISA experiment, exploring the parameter space in a way that is minimally sensitive to the underlying particle physics model. We find that discovery of SGWB polarisation is generally more challenging than measuring the total SGWB signal, but would be possible for appropriately strong transitions with large bubble sizes and a substantial polarisation fraction.

We outline the experimental concept and key scientific capabilities of AION (Atom Interferometer Observatory and Network), a proposed UK-based experimental programme using cold strontium atoms to search for ultra-light dark matter, to explore gravitational waves in the mid-frequency range between the peak sensitivities of the LISA and LIGO/Virgo/ KAGRA/INDIGO/Einstein Telescope/Cosmic Explorer experiments, and to probe other frontiers in fundamental physics. AION would complement other planned searches for dark matter, as well as probe mergers involving intermediate mass black holes and explore early universe cosmology. AION would share many technical features with the MAGIS experimental programme in the US, and synergies would flow from operating AION in a network with this experiment, as well as with other atom interferometer experiments such as MIGA, ZAIGA and ELGAR. Operating AION in a network with other gravitational wave detectors such as LIGO, Virgo and LISA would also offer many synergies.

We derive a lower bound on the merger rate of primordial black hole (PBH) binaries by estimating the maximal fraction of binaries that were perturbed between formation in the early Universe and merger, and computing a conservative merger rate of perturbed binaries. This implies robust constraints on the PBH abundance in the range $1-100 M_\odot$. We further show that LIGO/Virgo design sensitivity has the potential to reach the PBH mass range of $10^{-2}-10^3 M_\odot$. The constraint from the merger rate of perturbed binaries is stronger if PBHs are initially spatially clustered.

We propose in this White Paper a concept for a space experiment using cold atoms to search for ultra-light dark matter, and to detect gravitational waves in the frequency range between the most sensitive ranges of LISA and the terrestrial LIGO/Virgo/KAGRA/INDIGO experiments. This interdisciplinary experiment, called Atomic Experiment for Dark Matter and Gravity Exploration (AEDGE), will also complement other planned searches for dark matter, and exploit synergies with other gravitational wave detectors. We give examples of the extended range of sensitivity to ultra-light dark matter offered by AEDGE, and how its gravitational-wave measurements could explore the assembly of super-massive black holes, first-order phase transitions in the early universe and cosmic strings. AEDGE will be based upon technologies now being developed for terrestrial experiments using cold atoms, and will benefit from the space experience obtained with, e.g., LISA and cold atom experiments in microgravity. This paper is based on a submission (v1) in response to the Call for White Papers for the Voyage 2050 long-term plan in the ESA Science Programme. ESA limited the number of White Paper authors to 30. However, in this version (v2) we have welcomed as supporting authors participants in the Workshop on Atomic Experiments for Dark Matter and Gravity Exploration held at CERN: ({\tt https://indico.cern.ch/event/830432/}), as well as other interested scientists, and have incorporated additional material.

We study the generation of intergalactic magnetic fields in two models for first-order phase transitions in the early Universe that have been studied previously in connection with the generation of gravitational waves (GWs): the Standard Model supplemented by an $|H|^6$ operator (SM+$H^6$) and a classically scale-invariant model with an extra gauged U(1) $B - L$ symmetry (SM$_{B-L}$). We consider contributions to magnetic fields generated by bubble collisions and by turbulence in the primordial plasma, and we consider the hypotheses that helicity is seeded in the gauge field or kinetically. We study the conditions under which the intergalactic magnetic fields generated may be larger than the lower bounds from blazar observations, and correlate them with the observability of GWs and possible collider signatures. In the SM+$H^6$ model bubble collisions alone cannot yield large enough magnetic fields, whereas turbulence may do so. In the SM$_{B-L}$ model bubble collisions and turbulence may both yield magnetic fields above the blazar bound unless the B$-$L gauge boson is very heavy. In both models there may be observable GW and collider signatures if sufficiently large magnetic fields are generated.

We consider a real scalar singlet field which provides a strong first-order electroweak phase transition via its coupling to the Higgs boson, and gives a $CP$ violating contribution on the top quark mass via a dimension-6 operator. We study the correlation between the baryon-to-entropy ratio produced by electroweak baryogenesis, and the gravitational wave signal from the electroweak phase transition. We show that future gravitational wave experiments can test, in particular, the region of the model parameter space where the observed baryon-to-entropy ratio can be obtained even if the new physics scale, which is explicit in the dimension-6 operator, is high.

We present a novel mechanism for the production of primordial black holes (PBHs). The mechanism is based on a period of thermal inflation followed by fast-roll inflation due to tachyonic mass of order the Hubble scale. Large perturbations are generated at the end of the thermal inflation as the thermal inflaton potential turns from convex to concave. These perturbations can lead to copious production of PBHs when the relevant scales re-enter horizon. We show that such PBHs can naturally account for the observed dark matter in the Universe when the mass of the thermal inflaton is about $10^6\,$GeV and its coupling to the thermal bath preexisting the late inflation is of order unity. We consider also the possibility of forming the seeds of the supermassive black holes. In this case we find that the mass of the thermal inflaton is about $1\,$GeV, but its couplings have to be very small, $\sim 10^{-7}$. Finally we study a concrete realisation of our mechanism through a running mass model.

Within classically conformal models, the spontaneous breaking of scale invariance is usually associated to a strong first order phase transition that results in a gravitational wave background within the reach of future space-based interferometers. In this paper we study the case of the classically conformal gauged B-L model, analysing the impact of this minimal extension of the Standard Model on the dynamics of the electroweak symmetry breaking and derive its gravitational wave signature. Particular attention is paid to the problem of vacuum stability and to the role of the QCD phase transition, which we prove responsible for concluding the symmetry breaking transition in part of the considered parameter space. Finally, we calculate the gravitational wave signal emitted in the process, finding that a large part of the parameter space of the model can be probed by LISA.

We derive efficiency factors for the production of gravitational waves through bubble collisions and plasma-related sources in strong phase transitions, and find the conditions under which the bubble collisions can contribute significantly to the signal. We use lattice simulations to clarify the dependence of the colliding bubbles on their initial state. We illustrate our findings in two examples, the Standard Model with an extra $|H|^6$ interaction and a classically scale-invariant $U(1)_{\rm B-L}$ extension of the Standard Model. The contribution to the GW spectrum from bubble collisions is found to be negligible in the $|H|^6$ model, whereas it can play an important role in parts of the parameter space in the scale-invariant $U(1)_{\rm B-L}$ model. In both cases the sound-wave period is much shorter than a Hubble time, suggesting a significant amplification of the turbulence-sourced signal. We find, however, that the peak of the plasma-sourced spectrum is still produced by sound waves with the slower-falling turbulence contribution becoming important off-peak.

We revisit the cosmological and astrophysical constraints on the fraction of the dark matter in primordial black holes (PBHs) with an extended mass function. We consider a variety of mass functions, all of which are described by three parameters: a characteristic mass and width and a dark matter fraction. Various observations then impose constraints on the dark matter fraction as a function of the first two parameters. We show how these constraints relate to those for a monochromatic mass function, demonstrating that they usually become more stringent in the extended case than the monochromatic one. Considering only the well-established bounds, and neglecting the ones that depend on additional astrophysical assumptions, we find that there are three mass windows, around $4\times 10^{-17}M_\odot,$ $2\times 10^{-14}M_\odot$ and $25-100M_\odot$, where PBHs can constitute all dark matter. However, if one includes all the bounds, PBHs can only constitute of order $10\%$ of the dark matter.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

The abundance of primordial black holes (PBHs) in the mass range $0.1 - 10^3 M_\odot$ can potentially be tested by gravitational wave observations due to the large merger rate of PBH binaries formed in the early universe. To put the estimates of the latter on a firmer footing, we first derive analytical PBH merger rate for general PBH mass functions while imposing a minimal initial comoving distance between the binary and the PBH nearest to it, in order to pick only initial configurations where the binary would not get disrupted. We then study the formation and evolution of PBH binaries before recombination by performing N-body simulations. We find that the analytical estimate based on the tidally perturbed 2-body system strongly overestimates the present merger rate when PBHs comprise all dark matter, as most initial binaries are disrupted by the surrounding PBHs. This is mostly due to the formation of compact N-body systems at matter-radiation equality. However, if PBHs make up a small fraction of the dark matter, $f_{\rm PBH} \lesssim 10\%$, these estimates become more reliable. In that case, the merger rate observed by LIGO imposes the strongest constraint on the PBH abundance in the mass range $2 - 160 M_\odot$. Finally, we argue that, even if most initial PBH binaries are perturbed, the present BH-BH merger rate of binaries formed in the early universe is larger than $\mathcal{O}(10)\,{\rm Gpc}^{-3} {\rm yr}^{-1}\, f_{\rm PBH}^3$

We study production of dark matter (DM) in models with a non-standard expansion history. We consider both freeze-out and freeze-in mechanisms for producing the observed DM abundance in a model where the DM consists of scalar singlet particles coupled to the Standard Model sector via the Higgs portal. We show that a non-standard expansion phase can lead to a significant change in the DM abundance and therefore to observational ramifications. For example, for DM freeze-in the required portal coupling can be much larger, whereas for DM freeze-out much smaller values become allowed. We evaluate the relevant constraints and discuss prospects for direct detection of such DM.

The EDGES experiment has recently measured an anomalous global 21-cm spectrum due to hydrogen absorptions at redshifts of about $z\sim 17$. Model independently, the unusually low temperature of baryons probed by this observable sets strong constraints on any physical process that transfers energy into the baryonic environment at such redshifts. Here we make use of the 21-cm spectrum to derive bounds on the energy injection due to a possible population of ${\cal O}(1-100) M_\odot$ primordial black holes, which induce a wide spectrum of radiation during the accretion of the surrounding gas. After calculating the total radiative intensity of a primordial black hole population, we estimate the amount of heat and ionisations produced in the baryonic gas and compute the resulting thermal history of the Universe with a modified version of RECFAST code. Finally, by imposing that the temperature of the gas at $z\sim 17$ does not exceed the indications of EDGES, we constrain the possible abundance of primordial black holes. Depending on uncertainties related to the accretion model, we find that ${\cal O}(10) M_\odot$ primordial black holes can only contribute to a fraction $f_{\rm PBH}<(1-10^{-3})$ of the total dark matter abundance.

We study possible effects of a dark matter (DM) core on the maximum mass of a neutron star (NS), on the mass-radius relation and on the NS tidal deformability parameter $\Lambda$. We show that all these quantities would in general be reduced in the presence of a DM core. In particular, our calculations indicate that the presence of a DM core with a mass fraction $\sim 5\%$ could affect significantly the interpretation of these NS data as constraints on the nuclear equation of state (EOS), potentially excluding some EOS models on the basis of the measured mass of PSR J0348+0432, while allowing other EOS models to become consistent with the LIGO/Virgo upper limit on $\Lambda$. Specific scenarios for generating such DM cores are explored in an Appendix.

The radiation emitted by horizonless exotic compact objects (ECOs), such as wormholes, 2-2-holes, fuzzballs, gravastars, boson stars, collapsed polymers, superspinars etc., is expected to be strongly suppressed when compared to the radiation of black holes. If large primordial curvature fluctuations collapse into such objects instead of black holes, they do not evaporate or evaporate much slower than black holes and could thus constitute all of the dark matter with masses below $M < 10^{-16}M_\odot.$ We reevaluate the relevant experimental constraints for light ECOs in this mass range and show that very large new parameter space down to ECO masses $M\sim 10\,{\rm TeV}$ opens up for light primordial dark matter. A new dedicated experimental program is needed to test this mass range of primordial dark matter.

We study the production of primordial black hole (PBH) binaries and the resulting merger rate, accounting for an extended PBH mass function and the possibility of a clustered spatial distribution. Under the hypothesis that the gravitational wave events observed by LIGO were caused by PBH mergers, we show that it is possible to satisfy all present constraints on the PBH abundance, and find the viable parameter range for the lognormal PBH mass function. The non-observation of gravitational wave background allows us to derive constraints on the fraction of dark matter in PBHs, which are stronger than any other current constraint in the PBH mass range $0.5-30M_\odot$. We show that the predicted gravitational wave background can be observed by the coming runs of LIGO, and non-observation would indicate that the observed events are not of primordial origin. As the PBH mergers convert matter into radiation, they may have interesting cosmological implications, for example, in the context of relieving the tension between the high and low redshift measurements of the Hubble constant. However, we find that these effects are negligible as, after recombination, no more that $1\%$ of DM can be converted into gravitational waves.

The recent observation of the 21-cm absorption signal by the EDGES experiment provides a new observational window into the dynamics of the young Universe. Based on this result, we constrain the properties of ALPs and excited dark matter via the energy injections into the gaseous medium that these models produce. In particular, we derive bounds that outperform the present constraints for energy injections in the 10.2 eV to 50 keV range by analysing the intensity of Lyman-$\alpha$ photons and the electron emissivity produced by ionising radiation. Our results also show that once the standard temperature of the soft radiation background is assumed, the explanations of the observed 3.5 keV photon signal within the considered dark matter model are excluded at 95\% confidence level.

Motivated by the recent detection of the gravitational wave signal emitted by a binary neutron star merger, we analyse the possible impact of dark matter on such signals. We show that dark matter cores in merging neutron stars may yield an observable supplementary peak in the gravitational wave power spectral density following the merger, which could be distinguished from the features produced by the neutron components.

We study models of quartic inflation where the inflaton field $\phi$ is coupled non-minimally to gravity, $\xi \phi^2 R$, and perform a study of quantum corrections in curved space-time at one-loop level. We specifically focus on comparing results between the metric and Palatini theories of gravity. Transformation from the Jordan to the Einstein frame gives different results for the two formulations and by using an effective field theory expansion we derive the appropriate $\beta$-functions and the renormalisation group improved effective potentials in curved space for both cases in the Einstein frame. In particular, we show that in both formalisms the Einstein frame depends on the order of perturbation theory but that the flatness of the potential is unaltered by quantum corrections.

The recently claimed anomaly in the measurement of the 21 cm hydrogen absorption signal by EDGES at $z\sim 17$, if cosmological, requires the existence of new physics. The possible attempts to resolve the anomaly rely on either (i) cooling the hydrogen gas via new dark matter-hydrogen interactions or (ii) modifying the soft photon background beyond the standard CMB one, as possibly suggested also by the ARCADE~2 excess. We argue that solutions belonging to the first class are generally in tension with cosmological dark matter probes once simple dark sector models are considered. Therefore, we propose soft photon emission by light dark matter as a natural solution to the 21 cm anomaly, studying a few realizations of this scenario. We find that the signal singles out a photophilic dark matter candidate characterised by an enhanced collective decay mechanism, such as axion mini-clusters.

We consider an extension of the Standard Model with a singlet sector consisting of a real (pseudo)scalar and a Dirac fermion coupled with the Standard Model only via the scalar portal. We assume that the portal coupling is weak enough for the singlet sector not to thermalize with the Standard Model allowing the production of singlet particles via the freeze-in mechanism. If the singlet sector interacts with itself sufficiently strongly, it may thermalize within itself, resulting in dark matter abundance determined by the freeze-out mechanism operating within the singlet sector. We investigate this scenario in detail. In particular, we show that requiring the absence of inflationary isocurvature fluctuations provides lower bounds on the magnitude of the dark sector self-interactions and in parts of the parameter space favors sufficiently large self-couplings, supported also by the features observed in the small-scale structure formation.

We present an overview of scenarios where the observed Dark Matter (DM) abundance consists of Feebly Interacting Massive Particles (FIMPs), produced non-thermally by the so-called freeze-in mechanism. In contrast to the usual freeze-out scenario, frozen-in FIMP DM interacts very weakly with the particles in the visible sector and never attained thermal equilibrium with the baryon-photon fluid in the early Universe. Instead of being determined by its annihilation strength, the DM abundance depends on the decay and annihilation strengths of particles in equilibrium with the baryon-photon fluid, as well as couplings in the DM sector. This makes frozen-in DM very difficult but not impossible to test. In this review, we present the freeze-in mechanism and its variations considered in the literature (dark freeze-out and reannihilation), compare them to the standard DM freeze-out scenario, discuss several aspects of model building, and pay particular attention to observational properties and general testability of such feebly interacting DM.

We study production of primordial black holes (PBHs) during an early matter-dominated phase. As a source of perturbations, we consider either the inflaton field with a running spectral index or a spectator field that has a blue spectrum and thus provides a significant contribution to the PBH production at small scales. First, we identify the region of the parameter space where a significant fraction of the observed dark matter can be produced, taking into account all current PBH constraints. Then, we present constraints on the amplitude and spectral index of the spectator field as a function of the reheating temperature. We also derive constraints on the running of the inflaton spectral index, ${\rm d}n/{\rm d}{\rm ln}k \lesssim -0.002$, which are comparable to those from the Planck satellite for a scenario where the spectator field is absent.

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.

We consider a fermionic dark matter candidate interacting via a scalar mediator coupled with the Standard Model through a Higgs portal. We consider general setting including both scalar and pseudoscalar interactions between the scalar and fermion, and illustrate the relevant features for dark matter abundance, direct search limits and collider constraints. The case where dark matter has a self-interaction strength $\left\langle \sigma_V \right\rangle/m_\psi \sim 0.1-1 \,\mathrm{cm}^2/\mathrm{g}$ is strongly constrained, in particular by the Big Bang Nucleosynthesis. We show that these constraints can be alleviated by introducing a new light sterile neutrino $N$. The allowed region for the extended model consists of a triangle at $10 (\sin \theta)^{-2/5} \,{\rm MeV} \lesssim m_N \lesssim 1\,{\rm GeV}$.

We show that decoupled hidden sectors can have observational consequences. As a representative model example, we study dark matter production in the Higgs portal model with one real singlet scalar $s$ coupled to the Standard Model Higgs via $\lambda_{\rm hs}\Phi^\dagger\Phi s^2$ and demonstrate how the combination of non-observation of cosmological isocurvature perturbations and astrophysical limits on dark matter self-interactions imply stringent bounds on the magnitude of the scalar self-coupling $\lambda_{\rm s}s^4$. For example, for dark matter mass $m_{\rm s}=10$ MeV and Hubble scale during cosmic inflation $H_*=10^{12}$ GeV, we find $10^{-4}\lesssim \lambda_{\rm s}\lesssim 0.2$.

We consider a scenario where the inflaton decays to a hidden sector thermally decoupled from the visible Standard Model sector. A tiny portal coupling between the hidden and the visible sectors later heats the visible sector so that the Standard Model degrees of freedom come to dominate the energy density of the Universe before Big Bang Nucleosynthesis. We find that this scenario is viable, although obtaining the correct dark matter abundance and retaining successful Big Bang Nucleosynthesis is not obvious. We also show that the isocurvature perturbations constituted by a primordial Higgs condensate are not problematic for the viability of the scenario.

We study a singlet scalar extension of the Standard Model. The singlet scalar is coupled non-minimally to gravity and assumed to drive inflation, and also couple sufficiently strongly with the SM Higgs field in order to provide for a strong first order electroweak phase transition. Requiring the model to describe inflation successfully, be compatible with the LHC data, and yield a strong first order electroweak phase transition, we identify the regions of the parameter space where the model is viable. We also include a singlet fermion with scalar coupling to the singlet scalar to probe the sensitivity of the constraints on additional degrees of freedom and their couplings in the singlet sector. We also comment on the general feasibility of these fields to act as dark matter.

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.

We consider portal models which are ultraweakly coupled with the Standard Model, and confront them with observational constraints on dark matter abundance and isocurvature perturbations. We assume the hidden sector to contain a real singlet scalar $s$ and a sterile neutrino $\psi$ coupled to $s$ via a pseudoscalar Yukawa term. During inflation, a primordial condensate consisting of the singlet scalar $s$ is generated, and its contribution to the isocurvature perturbations is imprinted onto the dark matter abundance. We compute the total dark matter abundance including the contributions from condensate decay and nonthermal production from the Standard Model sector. We then use the Planck limit on isocurvature perturbations to derive a novel constraint connecting dark matter mass and the singlet self coupling with the scale of inflation: $m_{\rm DM}/{\rm GeV}\lesssim 0.2\lambda_{\rm s}^{\scriptscriptstyle 3/8} \left(H_*/10^{\scriptscriptstyle 11}{\rm GeV}\right)^{\scriptscriptstyle -3/2}$. This constraint is relevant in most portal models ultraweakly coupled with the Standard Model and containing light singlet scalar fields.

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.