Reheating is a process where the energy density of a dominant component of the universe other than radiation, such as a matter component, is transferred into radiation. It is usually assumed that the temperature of the universe decreases due to cosmic expansion even during the reheating process, in which case the maximal temperature of the universe is much higher than the reheat temperature. We point out that the temperature of the universe during reheating may in fact increase in well-motivated scenarios. We derive the necessary conditions for the temperature to increase during reheating and discuss concrete examples involving a scalar field. We comment on implications for particle physics and cosmology due to an increasing temperature during reheating.

Detection of electron recoils by dark matter (DM) may reveal the structure of the dark sector. We consider a scenario where a heavier DM particle inelastically scatters off an electron and is converted into a lighter DM particle. A small mass difference between the two DM particles is transferred into electron recoil energy. We investigate the DM-electron interaction mediated by a massive dark photon and evaluate the inelastic DM scattering rate, taking account of the atomic structure. It is found that the scattering rate is significantly enhanced because of the small mass splitting, which allows for a small momentum transfer matched with the size of the electron wave function. We show that there exists a viable parameter space which explains the excess of electron recoil events around 2-3 keV recently reported by the XENON1T experiment.

We investigate the possibility that the Peccei-Quinn phase transition occurs at a temperature far below the symmetry breaking scale. Low phase transition temperatures are typical in supersymmetric theories, where symmetry breaking fields have small masses. We find that QCD axions are abundantly produced just after the phase transition. The observed dark matter abundance is reproduced even if the decay constant is much lower than $10^{11}$ GeV. The produced axions tend to be warm. For some range of the decay constant, the effect of the predicted warmness on structure formation can be confirmed by future observations of 21 cm lines. A portion of parameter space requires a mixing between the Peccei-Quinn symmetry breaking field and the Standard Model Higgs, and predicts an observable rate of rare Kaon decays.

We propose a baryogenenesis mechanism that uses a rotating condensate of a Peccei-Quinn (PQ) symmetry breaking field and the dimension-five operator that gives Majorana neutrino masses. The rotation induces charge asymmetries for the Higgs boson and for lepton chirality through sphaleron processes and Yukawa interactions. The dimension-five interaction transfers these asymmetries to the lepton asymmetry, which in turn is transferred into the baryon asymmetry through the electroweak sphaleron process. QCD axion dark matter can be simultaneously produced by dynamics of the same PQ field via kinetic misalignment or parametric resonance, favoring an axion decay constant $f_a \lesssim 10^{10}$ GeV, or by conventional misalignment and contributions from strings and domain walls with $f_a \sim 10^{11}$ GeV. The size of the baryon asymmetry is tied to the mass of the PQ field. In simple supersymmetric theories, it is independent of UV parameters and predicts the supersymmtry breaking mass scale to be $\mathcal{O}(10-10^4)$ TeV, depending on the masses of the neutrinos and whether the condensate is thermalized during a radiation or matter dominated era. We also construct a theory where TeV scale supersymmetry is possible. Parametric resonance may give warm axions, and the radial component of the PQ field may give signals in rare kaon decays from mixing with the Higgs and in dark radiation.

Adding an axion-like particle (ALP) to the Standard Model, with a field velocity in the early universe, simultaneously explains the observed baryon and dark matter densities. This requires one or more couplings between the ALP and photons, nucleons, and/or electrons that are predicted as functions of the ALP mass. These predictions arise because the ratio of dark matter to baryon densities is independent of the ALP field velocity, allowing a correlation between the ALP mass, $m_a$, and decay constant, $f_a$. The predicted couplings are orders of magnitude larger than those for the QCD axion and for dark matter from the conventional ALP misalignment mechanism. As a result, this scheme, ALP cogenesis, is within reach of future experimental ALP searches from the lab and stellar objects, and for dark matter.

$SU(2)_L \times SU(2)_R$ gauge symmetry requires three right-handed neutrinos ($ N _i $), one of which, $N_1$, can be sufficiently stable to be dark matter. In the early universe, $ W _R $ exchange with the Standard Model thermal bath keeps the right-handed neutrinos in thermal equilibrium at high temperatures. $N_1$ can make up all of dark matter if they freeze-out while relativistic and are mildly diluted by subsequent decays of a long-lived and heavier right-handed neutrino, $N_2$. We systematically study this parameter space, constraining the symmetry breaking scale of $SU(2)_R$ and the mass of $N_1$ to a triangle in the $(v_R,M_1)$ plane, with $v_R = (10^6 - 3 \times 10^{12})$ GeV and $M_1 = (2\, {\rm keV} - 1 \, {\rm MeV)}$. Much of this triangle can be probed by signals of warm dark matter, especially if leptogenesis from $N_2$ decay yields the observed baryon asymmetry. The minimal value of $v_R$ is increased to $10^8 \, {\rm GeV}$ for doublet breaking of $SU(2)_R$, and further to $10^9 \, {\rm GeV}$ if leptogenesis occurs via $N_2$ decay, while the upper bound on $M_1$ is reduced to 100 keV. In addition, there is a component of hot $N_1$ dark matter resulting from the late decay of $N_2 \rightarrow N_1 \ell^+ \ell^-$ that can be probed by future cosmic microwave background observations. Interestingly, the range of $v_R$ allows both precision gauge coupling unification and the Higgs Parity understanding of the vanishing of the Standard Model Higgs quartic at scale $v_R$. Finally, we study freeze-in production of $N_1$ dark matter via the $W_R$ interaction, which allows a much wider range of $(v_R,M_1)$.

Axion cold dark matter from standard misalignment typically requires a decay constant $f_a~\gtrsim~10^{11}$ GeV. Kinetic misalignment and parametric resonance easily allow lower values of $f_a$ when the radial Peccei-Quinn (PQ) symmetry breaking field takes large initial values. Here, we consider the effects of inflation on kinetic misalignment and parametric resonance. We assume that the initial PQ field value is determined by quantum fluctuations, and is set by the Hubble parameter during inflation, $H_I$, and the PQ field mass. PQ field oscillations begin before or after the completion of reheating after inflation at a temperature $T_R$. We determine the range of $f_a$ and the inflationary parameters $(H_I, T_R)$ consistent with axion dark matter for a quartic potential for the PQ field. We find that $4\times 10^8$ GeV $< f_a < 10^{11}$ GeV can consistently produce axion dark matter. A significant portion of the allowed parameter space predicts rare kaon decays, $K_L \rightarrow (\pi^0 + \rm{missing \; energy})$, and/or suppression of structure formation on small scales.

We propose a mechanism called axiogenesis where the cosmological excess of baryons over antibaryons is generated from the rotation of the QCD axion. The Peccei-Quinn (PQ) symmetry may be explicitly broken in the early universe, inducing the rotation of a PQ charged scalar field. The rotation corresponds to the asymmetry of the PQ charge, which is converted into the baryon asymmetry via QCD and electroweak sphaleron transitions. In the concrete model we explore, interesting phenomenology arises due to the prediction of a small decay constant and the connections with new physics at the LHC and future colliders and with axion dark matter.

We consider a cosmological scenario where a relativistic particle and a stable massive particle are simultaneously produced from the decay of a late-decaying particle after Big-Bang Nucleosynthesis but before matter-radiation equality. The relativistic and massive particles behave as dark radiation and warm dark matter, respectively. Due to a common origin, the warmness and abundances are closely related. We refer to the models that lead to such a scenario as Common Origin of Warm and Relativistic Decay Products (COWaRD). We show that COWaRD predicts a correlation between the amount of dark radiation and suppression of the large scale structure, which can be tested in future precision cosmology observations. We demonstrate that COWaRD is realized, as an example, in a class of supersymmetric axion models and that future observations by the next generation Cosmic Microwave Background, Large Scalar Structure, and 21-cm surveys can reveal the structure of the theory.

In the conventional misalignment mechanism, the axion field has a constant initial field value in the early universe and later begins to oscillate. We present an alternative scenario where the axion field has a nonzero initial velocity, allowing an axion decay constant much below the conventional prediction from axion dark matter. This axion velocity can be generated from explicit breaking of the axion shift symmetry in the early universe, which may occur as this symmetry is approximate.

An exact parity replicates the Standard Model giving a Mirror Standard Model, SM $\leftrightarrow$ SM$'$. This "Higgs Parity" and the mirror electroweak symmetry are spontaneously broken by the mirror Higgs, $\left\langle H'\right\rangle = v' \gg \left\langle H\right\rangle$, yielding the Standard Model Higgs as a Pseudo-Nambu-Goldstone Boson of an approximate $SU(4)$ symmetry, with a quartic coupling $\lambda_{SM}(v') \sim 10^{-3}$. Mirror electromagnetism is unbroken and dark matter is composed of $e'$ and $\bar{e}'$. Direct detection may be possible via the kinetic mixing portal, and in unified theories this rate is correlated with the proton decay rate. With a high reheat temperature after inflation, the $e'$ dark matter abundance is determined by freeze-out followed by dilution from decays of mirror neutrinos, $\nu' \rightarrow \ell H$. Remarkably, this requires $v' \sim (10^8 - 10^{10})$ GeV, consistent with the Higgs mass, and a Standard Model neutrino mass of $(10^{-2} - 10^{-1})$ eV, consistent with observed neutrino masses. The mirror QCD sector exhibits a first order phase transition producing gravitational waves that may be detected by future observations. Mirror glueballs decay to mirror photons giving dark radiation with $\Delta N_{\rm eff} \sim 0.03 - 0.4$. With a low reheat temperature after inflation, the $e'$ dark matter abundance is determined by freeze-in from the SM sector by either the Higgs or kinetic mixing portal.

Scalar fields in the minimal supersymmetric standard model may have large field values during inflation. Because of approximate global symmetry, it is plausible that the phase directions of them are nearly massless during inflation and obtain quantum fluctuations, which may be the origin of the cosmic perturbations. If perturbations are produced through Q-ball formation, baryon asymmetry and dark matter can be consistently generated. Significant baryon and dark matter isocurvature perturbations are produced, but they are predicted to nearly compensate each other. The lepton asymmetry is much larger than the baryon asymmetry. The scenario predicts local non-Gaussianity of $f_{\rm NL} = 5/3$. Implication to the mass spectrum of supersymmetric particles is discussed.

We study interactions of cosmological relics, $X$, of mass $m$ and electric charge $qe$ in the galaxy, including thermalization with the interstellar medium, diffusion through inhomogeneous magnetic fields and Fermi acceleration by supernova shock waves. We find that for $m \mathop{}_{\textstyle \sim}^{\textstyle <} 10^{10} q \; ~{\rm GeV}$, there is a large flux of accelerated $X$ in the disk today, with a momentum distribution $\propto 1/p^{2.5}$ extending to $(\beta p)_{max} \sim 5 \times10^4 q \; ~{\rm GeV}$. Even though acceleration in supernova shocks is efficient, ejecting $X$ from the galaxy, $X$ are continually replenished by diffusion into the disk from the halo or confinement region. For $m \mathop{}_{\textstyle \sim}^{\textstyle >} 10^{10} q \; ~{\rm GeV}$, $X$ cannot be accelerated above the escape velocity within the lifetime of the shock. The accelerated $X$ form a component of cosmic rays that can easily reach underground detectors, as well as deposit energies above thresholds, enhancing signals in various experiments. We find that nuclear/electron recoil experiments place very stringent bounds on $X$ at low $q$; for example, $X$ as dark matter is excluded for $q$ above $10^{-9}$ for any $m$. For larger $q$, stringent bounds on the fraction of dark matter that can be $X$ are set by Cherenkov and ionization detectors. Nevertheless, very small $q$ is highly motivated by the kinetic mixing portal, and we identify regions of $(m,q)$ that can be probed by future experiments.

The Affleck-Dine mechanism creates the baryon asymmetry of the universe from an angular motion of a complex scalar field. The mechanism is efficient and allows for a low reheating temperature of the universe if the scalar field has a flat potential along its radial direction. We consider a possibility where the scalar field is a pseudo-Nambu-Goldstone boson arising from spontaneous breaking of an approximate global symmetry, so that the flatness of the potential is naturally understood. We construct a concrete realization of the idea based on chiral symmetry breaking in an SU(N) gauge theory. The Peccei-Quinn mechanism can be consistently incorporated into the model. We also comment on a possible impact of the model on early universe physics.

Several theoretical motivations point to ultralight QCD axions with large decay constants $f_a \simeq \mathcal{O}(10^{16}-10^{17})$ GeV, to which experimental proposals are dedicated. This regime is known to face the problem of overproduction of axion dark matter from the misalignment mechanism unless the misalignment angle $\theta_{\rm mis}$ is as small as $\mathcal{O}(10^{-3}-10^{-4})$, which is generally considered a fine-tuning problem. We investigate a dynamical explanation for a small $\theta_{\rm mis}$. The axion mass arises from strong dynamics and may be sufficiently enhanced by early dynamics so as to overcome Hubble friction and drive the field value to the bottom of the potential long before the QCD phase transition. Together with an approximate CP symmetry in the theory, this minimum is very closely related to today's value and thus $\theta_{\rm mis}$ can automatically be well under unity. Owing to such efficient relaxation, the isocurvature perturbations are essentially damped. As an existence proof, using supersymmetric theories we illustrate that the Higgs coupling with the inflaton energy can successfully achieve this axion damping in a consistent inflationary cosmology.

The QCD axion serves as a well-motivated dark matter candidate and the misalignment mechanism is known to reproduce the observed abundance with a decay constant $f_a \simeq \mathcal{O}(10^{12})$ GeV for a misalignment angle $\theta_{\rm mis} \simeq \mathcal{O}(1)$. While $f_a \ll 10^{12}$ GeV is of great experimental interest, the misalignment mechanism requires the axion to be very close to the hilltop, i.e. $\theta_{\rm mis} \simeq \pi$. This particular choice of $\theta_{\rm mis}$ has been understood as fine-tuning the initial condition. We offer a dynamical explanation for $\theta_{\rm mis} \simeq \pi$ in a class of models. The axion dynamically relaxes to the minimum of the potential by virtue of an enhanced mass in the early universe. This minimum is subsequently converted to a hilltop because the CP phase of the theory shifts by $\pi$ when one contribution becomes subdominant to another with an opposite sign. We demonstrate explicit and viable examples in supersymmetric models where the higher dimensional Higgs coupling with the inflaton naturally achieves both criteria. Associated phenomenology includes a strikingly sharp prediction of $3 \times 10^9$ GeV $\lesssim f_a \lesssim 10^{10}$ GeV and the absence of isocurvature perturbation.

The twin Higgs mechanism is a solution to the little hierarchy problem in which the top partner is neutral under the Standard Model (SM) gauge group. The simplest mirror twin Higgs (MTH) model -- where a $\mathbf{Z}_2$ symmetry copies each SM particle -- has too many relativistic degrees of freedom to be consistent with cosmological observations. We demonstrate that MTH models can have an observationally viable cosmology if the twin mass spectrum leads to twin neutrino decoupling before the SM and twin QCD phase transitions. Our solution requires the twin photon to have a mass of $\sim 20$ MeV and kinetically mix with the SM photon to mediate entropy transfer from the twin sector to the SM. This twin photon can be robustly discovered or excluded by future experiments. Additionally, the residual twin degrees of freedom present in the early Universe in this scenario would be detectable by future observations of the cosmic microwave background.

Vector bosons heavier than $10^{-22}$ eV can be viable dark matter candidates with distinctive experimental signatures. Ultralight dark matter generally requires a non-thermal origin to achieve the observed density, while still behaving like a pressureless fluid at late times. We show that such a production mechanism naturally occurs for vectors whose mass originates from a dark Higgs. If the dark Higgs has a large field value after inflation, the energy in the Higgs field can be efficiently transferred to vectors through parametric resonance. Computing the resulting abundance and spectra requires careful treatment of the transverse and longitudinal components, whose dynamics are governed by distinct differential equations. We study these equations in detail and find that the mass of the vector may be as low as $10 ^{ - 18 }$ eV, while making up the dominant dark matter abundance. This opens up a wide mass range of vector dark matter as cosmologically viable, further motivating their experimental search.

We revisit the non-thermal dark matter (DM) production during the thermalization and reheating era after inflation. The decay of inflaton produces high-energy particles that are thermalized to complete the reheating of the Universe. Before the thermalization is completed, DM can be produced from a collision between the high-energy particles and/or the ambient plasma. We calculate the DM abundance produced from these processes for the case where the cross section of the DM production is proportional to the $n$-th power of the center of mass energy. We find that the collision between the high-energy particles is almost always dominant for $n \gtrsim 4$ while it is subdominant for $n \lesssim 2$. The production from the ambient plasma is dominant when $n \lesssim 3$ and the reheating temperature is of the order of or larger than the DM mass. The production from a collision between the high-energy particle and the ambient plasma is important for $n \lesssim 2$ and the reheating temperature is much lower than the DM mass.

We investigate if the observed small and nearly scale-invariant primordial cosmic perturbation, i.e. the perturbation amplitude $P_\zeta\sim10^{-9}$ and the spectral index $n_s \simeq 0.965$, is typical in the landscape of vacua after imposing anthropic selections on them. We consider the situation where the universe begins from a metastable vacuum driving a precedent inflation, a curvature-dominated open universe is created by tunneling, and the curvature energy is inflated away by new inflation. We argue that the initial inflaton field value is homogeneous but typically non-zero because of the quantum fluctuation of long wavelength modes created during the precedent inflation, and only the universe which accidentally has a small inflaton field value is anthropically selected. We show that this bias, together with certain distributions of inflation model parameters that are physically well-motivated, makes the observed small and nearly scale-invariant spectrum typical.

The QCD axion is a good dark matter candidate. The observed dark matter abundance can arise from misalignment or defect mechanisms, which generically require an axion decay constant $f_a \sim \mathcal{O}(10^{11})$ GeV (or higher). We introduce a new cosmological origin for axion dark matter, parametric resonance from oscillations of the Peccei-Quinn symmetry breaking field, that requires $f_a \sim (10^8 -10^{11})$ GeV. The axions may be warm enough to give deviations from cold dark matter in Large Scale Structure.

The domain wall problem of the Peccei-Quinn mechanism can be solved if the Peccei-Quinn symmetry is explicitly broken by a small amount. Domain walls decay into axions, which may account for dark matter of the universe. This scheme is however strongly constrained by overproduction of axions unless the phase of the explicit breaking term is tuned. We investigate the case where the universe is matter-dominated around the temperature of the MeV scale and domain walls decay during this matter dominated epoch. We show how the viable parameter space is expanded.

The large hierarchy between the Planck scale and the weak scale can be explained by the dynamical breaking of supersymmetry in strongly coupled gauge theories. Similarly, the hierarchy between the Planck scale and the energy scale of inflation may also originate from strong dynamics, which dynamically generate the inflaton potential. We present a model of the hidden sector which unifies these two ideas, i.e., in which the scales of inflation and supersymmetry breaking are provided by the dynamics of the same gauge group. The resultant inflation model is chaotic inflation with a fractional power-law potential in accord with the upper bound on the tensor-to-scalar ratio. The supersymmetry breaking scale can be much smaller than the inflation scale, so that the solution to the large hierarchy problem of the weak scale remains intact. As an intrinsic feature of our model, we find that the sgoldstino, which might disturb the inflationary dynamics, is automatically stabilized during inflation by dynamically generated corrections in the strongly coupled sector. This renders our model a field-theoretical realization of what is sometimes referred to as sgoldstino-less inflation.

In all supersymmetric theories, gravitinos, with mass suppressed by the Planck scale, are an obvious candidate for dark matter; but if gravitinos ever reached thermal equilibrium, such dark matter is apparently either too abundant or too hot, and is excluded. However, in theories with an axion, a saxion condensate is generated during an early era of cosmological history and its late decay dilutes dark matter. We show that such dilution allows previously thermalized gravitinos to account for the observed dark matter over very wide ranges of gravitino mass, keV < $m_{3/2}$ < TeV, axion decay constant, $10^9$ GeV < $f_a$ < $10^{16}$ GeV, and saxion mass, 10 MeV < $m_s$ < 100 TeV. Constraints on this parameter space are studied from BBN, supersymmetry breaking, gravitino and axino production from freeze-in and saxion decay, and from axion production from both misalignment and parametric resonance mechanisms. Large allowed regions of $(m_{3/2}, f_a, m_s)$ remain, but differ for DFSZ and KSVZ theories. Superpartner production at colliders may lead to events with displaced vertices and kinks, and may contain saxions decaying to $(WW,ZZ,hh), gg, \gamma \gamma$ or a pair of Standard Model fermions. Freeze-in may lead to a sub-dominant warm component of gravitino dark matter, and saxion decay to axions may lead to dark radiation.

We present a simple and natural dark sector model in which dark matter particles arise as composite states of hidden strong dynamics and their stability is ensured by accidental symmetries. The model has only a few free parameters. In particular, the gauge symmetry of the model forbids the masses of dark quarks, and the confinement scale of the dynamics provides the unique mass scale of the model. The gauge group contains an Abelian symmetry $U(1)_D$, which couples the dark and standard model sectors through kinetic mixing. This model, despite its simple structure, has rich and distinctive phenomenology. In the case where the dark pion becomes massive due to $U(1)_D$ quantum corrections, direct and indirect detection experiments can probe thermal relic dark matter which is generically a mixture of the dark pion and the dark baryon, and the Large Hadron Collider can discover the $U(1)_D$ gauge boson. Alternatively, if the dark pion stays light due to a specific $U(1)_D$ charge assignment of the dark quarks, then the dark pion constitutes dark radiation. The signal of this radiation is highly correlated with that of dark baryons in dark matter direct detection.

In a Mirror Twin World with a maximally symmetric Higgs sector the little hierarchy of the Standard Model can be significantly mitigated, perhaps displacing the cutoff scale above the LHC reach. We show that consistency with observations requires that the Z2 parity exchanging the Standard Model with its mirror be broken in the Yukawa couplings. A minimal such effective field theory, with this sole Z2 breaking, can generate the Z2 breaking in the Higgs sector necessary for the Twin Higgs mechanism. The theory has constrained and correlated signals in Higgs decays, direct Dark Matter Detection and Dark Radiation, all within reach of foreseen experiments, over a region of parameter space where the fine-tuning for the electroweak scale is 10-50%. For dark matter, both mirror neutrons and a variety of self-interacting mirror atoms are considered. Neutrino mass signals and the effects of a possible additional Z2 breaking from the vacuum expectation values of B-L breaking fields are also discussed.

This report, based on the Dark Sectors workshop at SLAC in April 2016, summarizes the scientific importance of searches for dark sector dark matter and forces at masses beneath the weak-scale, the status of this broad international field, the important milestones motivating future exploration, and promising experimental opportunities to reach these milestones over the next 5-10 years.

An interesting possibility for dark matter is a scalar particle of mass of order 10 MeV-1 GeV, interacting with a U(1) gauge boson (dark photon) which mixes with the photon. We present a simple and natural model realizing this possibility. The dark matter arises as a composite pseudo Nambu-Goldstone boson (dark pion) in a non-Abelian gauge sector, which also gives a mass to the dark photon. For a fixed non-Abelian gauge group, SU(N), and a U(1) charge of the constituent dark quarks, the model has only three free parameters: the dynamical scale of the non-Abelian gauge theory, the gauge coupling of the dark photon, and the mixing parameter between the dark and standard model photons. In particular, the gauge symmetry of the model does not allow any mass term for the dark quarks, and stability of the dark pion is understood as a result of an accidental global symmetry. The model has a significant parameter space in which thermal relic dark pions comprise all of the dark matter, consistently with all experimental and cosmological constraints. In a corner of the parameter space, the discrepancy of the muon g-2 between experiments and the standard model prediction can also be ameliorated due to a loop contribution of the dark photon. Smoking-gun signatures of the model include a monophoton signal from the e+ e- collision into a photon and a "dark rho meson." Observation of two processes in e+ e- collision, the mode into the dark photon and that into the dark rho meson, would provide strong evidence for the model.

Spectator field models such as the curvaton scenario and the modulated reheating are attractive scenarios for the generation of the cosmic curvature perturbation, as the constraints on inflation models are relaxed. In this paper, we discuss the effect of Hubble induced masses on the dynamics of spectator fields after inflation. We pay particular attention to the Hubble induced mass by the kinetic energy of an oscillating inflaton, which is generically unsuppressed but often overlooked. In the curvaton scenario, the Hubble induced mass relaxes the constraint on the property of the inflaton and the curvaton, such as the reheating temperature and the inflation scale. We comment on the implication of our discussion for baryogenesis in the curvaton scenario. In the modulated reheating, the predictions of models e.g. the non-gaussianity can be considerably altered. Furthermore, we propose a new model of the modulated reheating utilizing the Hubble induced mass which realizes a wide range of the local non-gaussianity parameter.

Motivated by gauge coupling unification and dark matter, we present an extension to the Standard Model where both are achieved by adding an extra new matter multiplet. Such considerations lead to a Grand Unified Theory with very heavy WIMPzilla dark matter, which has mass greater than ~10^7 GeV and must be produced before reheating ends. Naturally, we refer to this scenario as GUTzilla dark matter. Here we present a minimal GUTzilla model, adding a vector-like quark multiplet to the Standard Model. Proton decay constraints require the new multiplet to be both color and electroweak charged, which prompts us to include a new confining SU(3) gauge group that binds the multiplet into a neutral composite dark matter candidate. Current direct detection constraints are evaded due to the large dark matter mass; meanwhile, next-generation direct detection and proton decay experiments will probe much of the parameter space. The relic abundance is strongly dependent on the dynamics of the hidden confining sector, and we show that dark matter production during the epoch of reheating can give the right abundance.

We consider a cosmologically consistent scenario with a heavy Polonyi field. The Polonyi field with a mass of ${\cal O}(100){\rm\,TeV}$ decays before the Big-Bang Nucleosynthesis (BBN) and avoids the severe constraint from the BBN. However, the abundance of the Lightest Supersymmetric Particle (LSP) produced from the decay often exceeds the observed dark matter density. In our scenario, the dark matter density is obtained by the LSP abundance with an aid of entropy production, and baryon asymmetry is generated by the Affleck-Dine mechanism. We show that the observed baryon-to-dark matter ratio of ${\cal O}(0.1\mathchar`-1)$ is naturally explained in sequestering models with a QCD axion.

The Peccei-Quinn mechanism suffers from the problem of the isocurvature perturbations. The isocurvature perturbations are suppressed if the Peccei-Quinn breaking scale is large during inflation. The oscillation of the Peccei-Quinn breaking field after inflation, however, leads to the formation of domain walls due to the parametric resonance effect. In this paper, we discuss the evolution of the Peccei-Quinn breaking field after inflation in detail, and propose a model where the parametric resonance is ineffective and hence domain walls are not formed. We also discuss consistency of our model with supersymmetric theory.

We point out that the prediction of the minimal chaotic inflation model is altered if a scalar field takes a large field value close to the Planck scale during inflation due to a negative Hubble induced mass. In particular, we show that the inflaton potential is effectively suppressed at a large inflaton field value in the presence of such a scalar field. The scalar field may be identified with the standard model Higgs field or flat directions in supersymmetric theory. With such spontaneous suppression, we find that the minimal chaotic inflation model, especially the model with a quadratic potential, is consistent with recent observations of the cosmic microwave background fluctuation without modifying the inflation model itself.

We consider an initial condition problem in a nearly quadratic chaotic inflation model in supergravity. We introduce shift symmetry breaking not only in the superpotential but also in the Kahler potential. In this model the inflaton potential is nearly quadratic for inflaton field values around the Planck scale, but deviates from the quadratic one for larger field values. As a result, the prediction on the tensor-to-scalar ratio can be smaller than that of a purely quadratic model. Due to the shift symmetry breaking in the Kahler potential, the inflaton potential becomes steep for large inflaton field values, which may prevent inflation from naturally taking place in a closed universe. We estimate an upper bound on the magnitude of the shift symmetry breaking so that inflation takes place before a closed universe with a Planck length size collapses, which yields a lower bound on the tensor-to-scalar ratio, $r \mathop{}_{\textstyle \sim}^{\textstyle >} 0.1$.

The observation of the B-mode in the cosmic microwave background radiation combined with the so-called Lyth bound suggests the trans-Planckian variation of the inflaton field during inflation. Such a large variation generates concerns over inflation models in terms of the effective field theory below the Planck scale. If the inflaton resides in a Riemann surface and the inflaton potential is a multivalued function of the inflaton field when it is viewed as a function on a complex plane, the Lyth bound can be satisfied while keeping field values in the effective field theory within the Planck scale. We show that a multivalued inflaton potential can be realized starting from a single-valued Lagrangian of the effective field theory below the Planck scale.

Matter isocurvature perturbations are strictly constrained from cosmic microwave background observations. We study a sneutrino curvaton model where both cold dark matter (CDM)/baryon isocurvature perturbations are generated. In our model, total matter isocurvature perturbations are reduced since the CDM/baryon isocurvature perturbations compensate for each other. We show that this model can not only avoid the stringent observational constraints but also suppress temperature anisotropies on large scales, which leads to improved agreement with observations.

We construct a natural inflation model in supergravity where the inflaton is identified with a modulus field possessing a shift symmetry. The superpotential for the inflaton is generated by meson condensation due to strong dynamics with deformed moduli constraints. In contrast to models based on gaugino condensation, the inflaton potential is generated without $R$-symmetry breaking and hence does not depend on the gravitino mass. Thus, our model is compatible with low scale supersymmetry.

The discovery of the primordial B-mode polarisation by the BICEP2 experiment indicates inflation with a relatively high energy scale. Taking this indication into account, we propose consistent scenarios to account for the observed baryon and dark matter densities in gravity and gauge mediated supersymmetry breaking models. The baryon asymmetry is explained by the Afflck-Dine mechanism, which requires relatively low reheating temperature to avoid a sizable baryonic isocurvature perturbation. The low reheating temperature then requires non-thermal production of dark matter to account for the correct relic density of dark matter. Our scenarios can account for the observations of baryon and dark matter density in gravity and gauge mediation and predict some parameters, including the mass of dark matter.

A model of natural inflation with an effectively trans-Planckian decay constant can be easily achieved by the "phase locking" mechanism while keeping field values in the effective field theory within the Planck scale. We give detailed description of "phase locked" inflation based on this mechanism. We also construct supersymmetric natural inflation based on this mechanism and show that the model is consistent with low scale supersymmetry. We also investigate couplings of the inflaton with the minimal supersymmetric standard model to achieve an appropriate reheating process. Interestingly, in a certain class of models, we find that the inflation scale is related to the mass of the right-handed neutrino in a consistent way with the seesaw mechanism.

Chaotic inflation based on a simple monomial scalar potential, V(phi) ~ phi^p, is an attractive large-field model of inflation capable of generating a sizable tensor-to-scalar ratio r. Therefore, assuming that future CMB observations will confirm the large r value reported by BICEP2, it is important to determine what kind of dynamical mechanism could possibly endow the inflaton field with such a simple effective potential. In this paper, we answer this question in the context of field theory, i.e. in the framework of dynamical chaotic inflation (DCI), where strongly interacting supersymmetric gauge dynamics around the scale of grand unification dynamically generate a fractional power-law potential via the quantum effect of dimensional transmutation. In constructing explicit models, we significantly extend our previous work, as we now consider a large variety of possible underlying gauge dynamics and relax our conditions on the field content of the model. This allows us to realize almost arbitrary rational values for the power p in the inflaton potential. The present paper may hence be regarded as a first step towards a more complete theory of dynamical chaotic inflation.

We investigate compatibility of the pure gravity mediation (or the minimal split supersymmetry) with chaotic inflation models in supergravity. We find that an approximate Z2 parity of the inflaton is useful to suppress gravitino production from thermal bath and to obtain consistent inflation dynamics. We discuss production of the lightest supersymmetric particle through the decay of the inflaton with the approximate Z2 symmetry and find that a large gravitino mass is favored to avoid the overproduction of the lightest supersymmetric particle, while a lower gravitino mass requires tuning of parameters. This may explain why the nature has chosen the gravitino mass of O(100) TeV rather than O(100) GeV.

The BICEP2 collaboration has recently reported a large tensor fluctuation in the cosmic microwave background, which suggests chaotic inflation models. In this letter, we reconsider the chaotic inflation model in the supergravity. We introduce a non-holomorphic shift-symmetry breaking parameter, which we expect to exist in general, and discuss its effect on the inflaton dynamics. We show that the model predicts a sizable deviation from the original chaotic inflation model and the predicted tensor fluctuation can lie between the BICEP2 result and the upper bound given by the Planck experiment with a small shift-symmetry breaking parameter. The model is characterized by only two parameters, which yields predictability and testability in future experiments.

If reheating of the Universe takes place via Planck-suppressed decay, it seems that the thermalization of produced particles might be delayed, since they have large energy/small number densities and number violating large angle scatterings which decrease the momentum of particles by large amount are inefficient correspondingly. In this paper, we study the thermalization of such "under occupied" decay products in detail, following recent developments in understanding the thermalization of non-abelian plasma. Contrary to the above naive expectation, it is shown that in most cases thermalization after/during reheating occurs instantaneously by properly taking account of scatterings with small angles and of particles with small momenta. In particular, the condition for instantaneous thermalization before the completion of reheating is found to be $\alpha^{8/5} \gg (m_\phi / M_{\rm pl}) (M_{\rm pl}^2 \Gamma_\phi / m_\phi^3)^{1/5}$, which is much milder than that obtained in previous works with small angle scatterings taken into account.

We estimate dark matter density for the Universe with a reheating temperature smaller than the mass of dark matter, assuming dark matter to be a weakly interacting massive particle. During the reheating process, an inflaton decays and releases high-energy particles, which are scattered inelastically by the thermal plasma and emit many particles. Dark matters are produced through these inelastic scattering processes and pair creation processes by high-energy particles. We properly take account of the Landau-Pomeranchuk-Migdal effect on inelastic processes and show that the resultant energy density of dark matter is much larger than that estimated in the literature and can be consistent with that observed when the mass of dark matter is larger than $O(100)$ GeV.

The measurement of a large tensor-to-scalar ratio by the BICEP2 experiment, r = 0.20 (-0.05)(+0.07), severely restricts the landscape of viable inflationary models and shifts attention once more towards models featuring large inflaton field values. In this context, chaotic inflation based on a fractional power-law potential that is dynamically generated by the dynamics of a strongly coupled supersymmetric gauge theory appears to be particularly attractive. We revisit this class of inflation models and find that, in the light of the BICEP2 measurement, models with a non-minimal gauge group behind the dynamical model seem to be disfavored, while the model with the simplest group, i.e. SU(2), is consistent with all results. We also discuss how the dynamical model can be distinguished from the standard chaotic inflation model based on a quadratic inflaton potential.

The large value of the tensor-to-scalar ratio in the cosmic microwave background radiation reported by the BICEP2 collaboration gives strong impact on models of supersymmetry (SUSY). The large ratio indicates inflation with a high-energy scale and thus a high reheating temperature in general, and various SUSY models suffer from the serious gravitino and Polonyi problems. In this article, we discuss a class of the high-scale SUSY breaking models which are completely free from those problems. With especially focusing on the dark matter relic abundance, we examine how the BICEP2 result narrows down the parameter space of the models, assuming the simplest chaotic inflation model. We find that the mass of the dark matter is predicted to be less than about 1 TeV thanks to the non-thermal production in the early universe through the decay of abundant gravitinos produced after the reheating process. We also discuss implications in some details to dark matter searches at collider and indirect dark matter detection experiments.

We investigate the consistency of a scenario in which the baryon asymmetry, dark matters, as well as the cosmic density perturbation are generated simultaneously through the evaporation of primordial black holes (PBHs). This scenario can explain the coincidence of the dark matter and the baryon density of the universe, and is free from the isocurvature perturbation problem. We show that this scenario predicts the masses of PBHs, right-handed neutrinos and dark matters, the Hubble scale during inflation, the non-gaussianity and the running of the spectral index. We also discuss the testability of the scenario by detecting high frequency gravitational waves from PBHs.

We present a novel mechanism to generate the cosmic perturbation from evaporation of primordial black holes. A mass of a field is fluctuated if it is given by a vacuum expectation value of a light scalar field because of the quantum fluctuation during inflation. The fluctuated mass causes variations of the evaporation time of the primordial black holes. Therefore provided the primordial black holes dominate the universe when they evaporate, primordial cosmic perturbations are generated. We find that the amplitude of the large scale curvature perturbation generated in this scenario can be consistent with the observed value. Interestingly, our mechanism works even if all fields which are responsible for inflation and the generation of the cosmic perturbation are decoupled from the visible sector except for the gravitational interaction. An implication to the running spectral index is also discussed.

In supersymmetric theories, the $R$ symmetry plays a unique role in suppressing a constant term in the superpotential. In single chiral field models of spontaneous breaking of a discrete $R$ symmetry, an $R$-breaking field can be a good candidate for an inflaton in new inflation models. In this paper, we revisit the compatibility of the single-field $R$-breaking new inflation model with the results of the Planck experiment. As a result, we find that the model predicts a lower limit on the gravitino mass, $m_{3/2}>O(100)$ TeV. This lower limit is consistent with the observed Higgs mass of 126 GeV when the masses of the stops are of order the gravitino mass scale.

The Polonyi field is a necessary ingredient in any viable scenario of gravity mediated supersymmetry breaking. However, it is known that the presence of the Polonyi field leads to several serious cosmological problems, which are collectively referred to as the Polonyi problem. We show that the Polonyi problem can be solved if the Polonyi field couples to a pseudo modulus in the superpotential and this pseudo modulus has a large field expectation value during inflation. To illustrate our idea, we construct an explicit model which can be readily connected to scenarios of gravity mediation. The generation of the mass parameters contained in our model by strong gauge dynamics is also commented on.

Models of chaotic inflation with a fractional power-law potential are not only viable but also testable in the foreseeable future. We show that such models can be realized in simple strongly coupled supersymmetric gauge theories. In these models, the energy scale during inflation is dynamically generated by the dimensional transmutation due to the strong gauge dynamics. Therefore, such models not only explain the origin of the fractional power in the inflationary potential but also provide a reason why the energy scale of inflation is much smaller than the Planck scale.

We propose a curvaton model in which the initial condition of the curvaton oscillation is determined by its attractor behavior during inflation. Assuming a chaotic inflation model, we find that the initial condition determined by the attractor behavior is appropriate to generate a sizable non-Gaussianity contribution to the curvature perturbation, which will be tested in the foreseeable future. Implications on the thermal history of the universe and on particle physics models are also discussed.

The recently observed mass ~ 125 GeV for the Higgs boson suggests a high-energy scale SUSY breaking, above O(10) TeV. It is, however, very puzzling why nature chooses such a high energy scale for the SUSY breaking, if the SUSY is a solution to the hierarchy problem. We show that the pure gravity mediation provides us with a possible solution to this puzzle if the topological inflation is the last inflation in the early universe. We briefly discuss a chaotic inflation model in which a similar solution can be obtained.