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.

We study initial conditions for inflation in scenarios where the inflaton potential has a plateau shape. Such models are those most favored by Planck data and can be obtained in a large number of model classes. As a representative example, we consider Higgs inflation with and without an $R^2$ term in the context of Palatini gravity. We show that inflation with a large number of e-folds generically occurs in a large part of the parameter space without any fine-tuning of parameters even when the scale of inflation and the inflaton field value during inflation are much smaller than the Planck scale. We discuss consequences for detection of primordial gravitational waves and spectral tilt of curvature perturbations, as well as the recently proposed "Trans-Planckian Censorship" conjecture.

Gravitational waves (GW) produced in the early Universe contribute to the number of relativistic degrees of freedom, $N_{\rm eff}$, during Big Bang Nucleosynthesis (BBN). By using the constraints on $N_{\rm eff}$, we present a new bound on how much the observable Universe could have expanded during cosmic inflation. The new bound is roughly four orders of magnitude more stringent than previous bounds. We also discuss the sensitivities of the current and planned GW observatories such as LIGO and LISA, and show that the constraints they could impose are always less stringent than the BBN bound.

We present an introduction to cosmic inflation in the context of Palatini gravity, which is an interesting alternative to the usual metric theory of gravity. In the latter case only the metric $g_{\mu\nu}$ determines the geometry of space-time, whereas in the former case both the metric and the space-time connection $\Gamma^\lambda_{\mu\nu}$ are a priori independent variables - a choice which can lead to a theory of gravity different from the metric one. In scenarios where the field(s) responsible for cosmic inflation are coupled non-minimally to gravity or the gravitational sector is otherwise extended, assumptions of the underlying gravitational degrees of freedom can have a big impact on the observational consequences of inflation. We demonstrate this explicitly by reviewing several interesting and well-motivated scenarios including Higgs inflation, $R^2$ inflation, and $\xi$-attractor models. We also discuss some prospects for future research and argue why $r=10^{-3}$ is a particularly important goal for future missions that search for signatures of primordial gravitational waves.

We show that a non-minimal coupling to gravity can not only make some inflationary models consistent with cosmological data, similar to the case of Higgs inflation, but can also invoke slow-roll violation to realize graceful exit from inflation. In particular, this is the case in models where a destabilizing mechanism that ends inflation should be assumed when the model is minimally coupled to gravity. As explicit examples, we consider the power-law and inverse monomial inflation models with a non-minimal coupling to gravity. While these models are excluded in the minimally coupled case, we show that they can become viable again in non-minimally coupled scenarios. We also argue that in most scenarios we consider reheating can be naturally realized via gravitational particle production but that this depends on the underlying theory of gravity in a non-trivial way.

If the inflationary phase lasted longer than the minimal period, the length scales observed today originate from modes that were smaller than the Planck length during inflation. It was recently argued that this "trans-Planckian problem" can never arise in a consistent string theory framework, which places a stringent constraint on the energy scale of inflation, $V^{1/4}\lesssim 10^9$ GeV. In this paper, we show that this requirement corresponds to a very small Hubble scale during inflation, $H_{\rm inf}\lesssim 1$ GeV, and therefore has serious consequences on scenarios where the dark matter density was generated by amplification of quantum fluctuations during inflation. We also present a class of inflationary models which both satisfy the above limit for the scale of inflation and are in perfect agreement with observational data.

Dark matter (DM) may have its origin in a pre-Big Bang epoch, the cosmic inflation. Here, we consider for the first time a broad class of scenarios where a massive free scalar field unavoidably reaches an equilibrium between its classical and quantum dynamics in a characteristic time scale during inflation and sources the DM density. The study gives the abundance and perturbation spectrum of any DM component sourced by the scalar field. We show that this class of scenarios generically predicts enhanced structure formation, allowing one to test models where DM interacts with matter only gravitationally.

We consider scenarios where the inflaton field decays dominantly to a hidden dark matter (DM) sector. By studying the typical behavior of the Standard Model (SM) Higgs field during inflation, we derive a relation between the primordial tensor-to-scalar ratio $r$ and amplitude of the residual DM isocurvature perturbations $\beta$ which is typically generated if the DM is thermally decoupled from the SM sector. We consider different expansion histories and find that if the Universe was radiation- or matter-dominated after inflation, a future discovery of primordial DM isocurvature will rule out all simple scenarios of this type because generating observable $\beta$ from the Higgs is not possible without violating the bounds on $r$. Seen another way, the Higgs field is generically not a threat to models where both the inflaton and DM reside in a decoupled sector. However, this is not necessarily the case for an early kination-dominated epoch, as then the Higgs can source sizeable $\beta$. We also discuss why the Higgs cannot source the observed curvature perturbation at large scales in any of the above cases but how the field can still be the dominant source of curvature perturbations at small scales.

In order to accommodate the QCD axion as the dark matter (DM) in a model in which the Peccei-Quinn (PQ) symmetry is broken before the end of inflation, a relatively low scale of inflation has to be invoked in order to avoid bounds from DM isocurvature fluctuations, $H_* \lesssim \mathcal{O}(10^9)\,$GeV. We construct a simple model in which the Standard Model Higgs field is non-minimally coupled to Palatini gravity and acts as the inflaton, leading to a scale of inflation $H_* \sim 10^8\,$GeV. When the energy scale at which the PQ symmetry breaks is much larger than the scale of inflation, we find that in this scenario the required axion mass for which the axion constitutes all DM is $m_0 \lesssim 0.05{\rm \,\mu eV}$ for a quartic Higgs self-coupling $\lambda_\phi = 0.1$, which correspond to the PQ breaking scale $v_\sigma \gtrsim 10^{14}\,$GeV and tensor-to-scalar ratio $r \sim 10^{-12}$. Future experiments sensitive to the relevant QCD axion mass scale can therefore shed light on the physics of the Universe before the end of inflation.

We consider light scalar fields during inflation and show how the stochastic spectral expansion method can be used to calculate two-point correlation functions of an arbitrary local function of the field in de Sitter space. In particular, we use this approach for a massive scalar field with quartic self-interactions to calculate the fluctuation spectrum of the density contrast and compare it to other approximations. We find that neither Gaussian nor linear approximations accurately reproduce the power spectrum, and in fact always overestimate it. For example, for a scalar field with only a quartic term in the potential, $V=\lambda\phi^4/4$, we find a blue spectrum with spectral index $n-1=0.579\sqrt{\lambda}$.

We study models of inflation where the scalar field $\phi$ that drives inflation is coupled non-minimally to gravity via $\xi \phi^2 R$, or where the gravity sector is enlarged by an $R^2$ term. We consider the original Higgs inflation, Starobinsky inflation, and two different versions of a scenario where the inflaton is a scalar field other than the Higgs, and discuss if they can be distinguished from each other by measuring the tensor-to-scalar ratio and runnings of the spectral index of primordial curvature perturbations, on top of the amplitude and spectral index of the perturbations. We consider both metric and Palatini theories of gravity, showing how detailed studies of non-minimally coupled models can help to identify the inflaton field and how they may provide for a way to also distinguish between different theories of gravity in the present context.

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.

We present the first end-to-end nonperturbative analysis of the gravitational wave power spectrum from a thermal first-order electroweak phase transition (EWPT), using the framework of dimensionally reduced effective field theory and pre-existing nonperturbative simulation results. We are able to show that a first-order EWPT in any beyond the Standard Model (BSM) scenario that can be described by a Standard Model-like effective theory at long distances will produce gravitational wave signatures too weak to be observed at existing and planned detectors. This implies that colliders are likely to provide the best chance of exploring the phase structure of such theories, while transitions strong enough to be detected at gravitational wave experiments require either previously neglected higher-dimension operators or light BSM fields to be included in the dimensionally reduced effective theory and therefore necessitate dedicated nonperturbative studies. As a concrete application, we analyze the real singlet-extended Standard Model and identify regions of parameter space with single-step first-order transitions, comparing our findings to those obtained using a fully perturbative method. We discuss the prospects for exploring the electroweak phase diagram in this model at collider and gravitational wave experiments in light of our nonperturbative results.

It has recently been suggested that the Standard Model Higgs boson could act as the inflaton while minimally coupled to gravity - given that the gravity sector is extended with an $\alpha R^2$ term and the underlying theory of gravity is of Palatini, rather than metric, type. In this paper, we revisit the idea and correct some shortcomings in earlier studies. We find that in this setup the Higgs can indeed act as the inflaton and that the tree-level predictions of the model for the spectral index and the tensor-to-scalar ratio are $n_s\simeq 0.941$, $r\simeq 0.3/(1+10^{-8}\alpha)$, respectively, for a typical number of e-folds, $N=50$, between horizon exit of the pivot scale $k=0.05\, {\rm Mpc}^{-1}$ and the end of inflation. Even though the tensor-to-scalar ratio is suppressed compared to the usual minimally coupled case and can be made compatible with data for large enough $\alpha$, the result for $n_s$ is in severe tension with the Planck results. We briefly discuss extensions of the model.

If cosmic inflation was driven by an electrically neutral scalar field stable on cosmological time scales, the field necessarily constitutes all or part of dark matter (DM). We study this possibility in a scenario where the inflaton field $s$ resides in a hidden sector, which is coupled to the Standard Model sector through the Higgs portal $\lambda_{hs} s^2\mathcal{H}^\dagger\mathcal{H}$ and non-minimally to gravity via $\xi_s s^2 R$. We study scenarios where the field $s$ first drives inflation, then reheats the Universe, and later constitutes all DM. We consider two benchmark scenarios where the DM abundance is generated either by production during reheating or via non-thermal freeze-in. In both cases, we take into account all production channels relevant for DM in the mass range from keV to PeV scale. On the inflationary side, we compare the dynamics and the relevant observables in two different but well-motivated theories of gravity (metric and Palatini), discuss multifield effects in case both fields ($s$ and $h$) were dynamical during inflation, and take into account the non-perturbative nature of particle production during reheating. We find that, depending on the initial conditions for inflation, couplings and the DM mass, the scenario works well especially for large DM masses, $10^2$ GeV$\lesssim m_{s}\lesssim 10^6$ GeV, although there are also small observationally allowed windows at the keV and MeV scales. We discuss how the model can be tested through astrophysical observations.

We study production of self-interacting dark matter (DM) during an early matter-dominated phase. As a benchmark scenario, we consider a model where the DM consists of singlet scalar particles coupled to the visible Standard Model (SM) sector via the Higgs portal. We consider scenarios where the initial DM abundance is set by either the usual thermal freeze-out or an alternative freeze-in mechanism, where DM was never in thermal equilibrium with the SM sector. For the first time, we take the effect of self-interactions within the hidden sector into account in determining the DM abundance, reminiscent to the Strongly Interacting Massive Particle (SIMP) scenario. In all cases, the number density of DM may change considerably compared to the standard radiation-dominated case, having important observational and experimental ramifications.

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 observed dark matter abundance in the Universe can be fully accounted for by a minimally coupled spectator scalar field that was light during inflation and has sufficiently strong self-coupling. In this scenario, dark matter was produced during inflation by amplification of quantum fluctuations of the spectator field. The self-interaction of the field suppresses its fluctuations on large scales, and therefore avoids isocurvature constraints. The scenario does not require any fine-tuning of parameters. In the simplest case of a single real scalar field, the mass of the dark matter particle would be in the range $1~{\rm GeV}\lesssim m\lesssim 10^8~{\rm GeV}$, depending on the scale of inflation, and the lower bound for the quartic self-coupling is $\lambda\gtrsim 0.45$.

We study multifield inflation in scenarios where the fields are coupled non-minimally to gravity via $\xi_I(\phi^I)^n g^{\mu\nu}R_{\mu\nu}$, where $\xi_I$ are coupling constants, $\phi^I$ the fields driving inflation, $g_{\mu\nu}$ the space-time metric, $R_{\mu\nu}$ the Ricci tensor, and $n>0$. We consider the so-called $\alpha$-attractor models in two formulations of gravity: in the usual metric case where $R_{\mu\nu}=R_{\mu\nu}(g_{\mu\nu})$, and in the Palatini formulation where $R_{\mu\nu}$ is an independent variable. As the main result, we show that, regardless of the underlying theory of gravity, the field-space curvature in the Einstein frame has no influence on the inflationary dynamics at the limit of large $\xi_I$, and one effectively retains the single-field case. However, the gravity formulation does play an important role: in the metric case the result means that multifield models approach the single-field $\alpha$-attractor limit, whereas in the Palatini case the attractor behaviour is lost also in the case of multifield inflation. We discuss what this means for distinguishing between different models of inflation.

We study the formation of primordial black holes (PBHs) in the early Universe during a period of slow reheating after inflation. We demonstrate how the PBH formation mechanism may change even before the end of the matter-dominated phase and calculate the expected PBH mass function. We find that there is a threshold for the variance of the density contrast, $\sigma_c \simeq 0.05$, below which the transition occurs even before reheating, with this having important consequences for the PBH mass function. We also show that there is a maximum cut-off for the PBH mass at around $100\,M_{\odot}$, below which the subdominant radiation bath affects PBH production, making the scenario particularly interesting for the recent LIGO observations of black hole mergers.

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.

We show that in the Feebly Interacting Massive Particle (FIMP) model of Dark Matter (DM), one may express the inflationary energy scale $H_*$ as a function of three otherwise unrelated quantities, the DM isocurvature perturbation amplitude, its mass and its self-coupling constant, independently of the tensor-to-scalar ratio. The FIMP model assumes that there exists a real scalar particle that alone constitutes the DM content of the Universe and couples to the Standard Model via a Higgs portal. We consider carefully the various astrophysical, cosmological and model constraints, accounting also for variations in inflationary dynamics and the reheating history, to derive a robust estimate for $H_*$ that is confined to a relatively narrow range. We point out that, within the context of the FIMP DM model, one may thus determine $H_*$ reliably even in the absence of observable tensor perturbations.

We study observational consequences arising from dark matter (DM) of non-thermal origin, produced by dark freeze-out from a hidden sector heat bath. We assume this heat bath was populated by feebly-coupled mediator particles, produced via a Higgs portal interaction with the Standard Model (SM). The dark sector then attained internal equilibrium with a characteristic temperature different from the SM photon temperature. We find that even if the coupling between the DM and the SM sectors is very weak, the scenario allows for indirect observational signals. We show how the expected strength of these signals depends on the temperature of the hidden sector at DM freeze-out.

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 study Quadratic Inflation with the inflaton field $\phi$ coupled non-minimally to the curvature scalar $R$, so that the potential during inflation is of the form $V\propto m^2\phi^2+\xi R\phi^2$. We show that with a suitable choice of the non-minimal coupling strength, $\xi=\mathcal{O}(10^{-3})$, one can resurrect the success of the scenario when compared against the Planck and BICEP2/Keck Array data, and that in the region of the parameter space which is still allowed the model predicts values of the tensor-to-scalar ratio in the range $0.01\leq r < 0.12$, making it possible to either confirm the scenario or rule it out already by the current or near-future experiments, such as BICEP3 or LiteBIRD. However, we show that in this case the near-future observations are unlikely to be able to distinguish between the metric and Palatini formulations of gravity.

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.

We study dark matter production in scenarios where a scale invariant hidden sector interacts with the Standard Model degrees of freedom via a Higgs portal $\lambda \Phi^\dagger\Phi s^2$. If the hidden sector is very weakly coupled to the SM but exhibits strong interactions within its own particle species, the dark matter abundance may arise as a result of a dark freeze-out occurring in the hidden sector. Because of scale invariance, the free parameters in the hidden sector are determined and the dark matter candidate exhibits a 'WIMP miracle of the second kind'. Demonstrating the predictive power of scale invariance, we carry out thorough analysis of dark matter production in several benchmark scenarios where the hidden sector contains either a scalar, fermion (sterile neutrino), or vector dark matter, and discuss the observational consequences of these scenarios.

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 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$.

In real-time lattice simulations of cosmic strings in the Abelian Higgs model, the broken translational invariance introduces lattice artefacts; relativistic strings therefore decelerate and radiate. We introduce two different methods to construct a moving string on the lattice, and study in detail the lattice effects on moving strings. We find that there are two types of lattice artefact: there is an effective maximum speed with which a moving string can be placed on the lattice, and a moving string also slows down, with the deceleration approximately proportional to the exponential of the velocity. To mitigate this, we introduce and study an improved discretisation, based on the tree-level L\"{u}scher-Weisz action, which is found to reduce the deceleration by an order of magnitude, and to increase the string speed limit by an amount equivalent to halving the lattice spacing. The improved algorithm is expected to be very useful for 3D simulations of cosmic strings in the early universe, where one wishes to simulate as large a volume as possible.

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 present a scenario where a $Z_2$-symmetric scalar field $\phi$ first drives cosmic inflation, then reheats the Universe but remains out-of-equilibrium itself, and finally comprises the observed dark matter abundance, produced by particle decays \`{a} la freeze-in mechanism. We work model-independently without specifying the interactions of the scalar field besides its self-interaction coupling, $\lambda\phi^4$, non-minimal coupling to gravity, $\xi\phi^2R$, and coupling to another scalar field, $g\phi^2\sigma^2$. We find the scalar field $\phi$ serves both as the inflaton and a dark matter candidate if $10^{-9}\lesssim \lambda\lesssim g\lesssim 10^{-7}$ and $3 \rm{keV} \lesssim m_{\rm \phi}\lesssim 85 \rm{MeV}$ for $\xi=\mathcal{O}(1)$. Such a small value of the non-minimal coupling is also found to be of the right magnitude to produce the observed curvature perturbation amplitude within the scenario. We also discuss how the model may be distinguished from other inflationary models of the same type by the next generation CMB satellites.

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.

We show that dark matter abundance and the inflationary scale $H$ could be intimately related. Standard Model extensions with Higgs mediated couplings to new physics typically contain extra scalars displaced from vacuum during inflation. If their coupling to Standard Model is weak, they will not thermalize and may easily constitute too much dark matter reminiscent to the moduli problem. As an example we consider Standard Model extended by a $Z_2$ symmetric singlet $s$ coupled to the Standard Model Higgs $\Phi$ via $\lambda \Phi^{\dag}\Phi s^2$. Dark matter relic density is generated non-thermally for $\lambda \lesssim 10^{-7}$. We show that the dark matter yield crucially depends on the inflationary scale. For $H\sim 10^{10}$ GeV we find that the singlet self-coupling and mass should lie in the regime $\lambda_{\rm s}\gtrsim 10^{-9}$ and $m_{\rm s}\lesssim 50$ GeV to avoid dark matter overproduction.

We study the post-inflationary dynamics of the Standard Model Higgs and a real singlet scalar $s$, coupled together through a renormalizable coupling $\lambda_{sh}h^2s^2$, in a $Z_2$ symmetric model that may explain the observed dark matter abundance and/or the origin of baryon asymmetry. The initial values for the Higgs and $s$ condensates are given by inflationary fluctuations, and we follow their dissipation and relaxation to the low energy vacua. We find that both the lowest order perturbative and the non-perturbative decays are blocked by thermal effects and large background fields and that the condensates decay by two-loop thermal effects. Assuming instant reheating at $T=10^{16}$ GeV, the characteristic temperature for the Higgs condensate thermalization is found to be $T_h \sim 10^{14}$ GeV, whereas $s$ thermalizes typically around $T_s \sim 10^{6}$ GeV. By that time, the amplitude of the singlet is driven very close to the vacuum value by the expansion of the universe, unless the portal coupling takes a value $\lambda_{sh}\lesssim 10^{-7}$ and the singlet $s$ never thermalizes. With these values of the coupling, it is possible to slowly produce a sizeable fraction of the observed dark matter abundance via singlet condensate fragmentation and thermal Higgs scattering. Physics also below the electroweak scale can therefore be affected by the non-vacuum initial conditions generated by inflation.