We discuss the footprint of evaporation of primordial black holes (PBHs) on stochastic gravitational waves (GWs) induced by scalar perturbations. We consider the case where PBHs once dominated the Universe but eventually evaporated before the big bang nucleosynthesis. The reheating through the PBH evaporation could end with a sudden change in the equation of state of the Universe compared to the conventional reheating caused by particle decay. We show that this "sudden reheating" by the PBH evaporation enhances the induced GWs, whose amount depends on the length of the PBH-dominated era and the width of the PBH mass function. We explore the possibility to constrain the primordial abundance of the evaporating PBHs by observing the induced GWs. We find that the abundance parameter $\beta \gtrsim 10^{-5} \text{ - }10^{-8}$ for $\mathcal{O}(10^3 \text{ - } 10^5) \, \text{g}$ PBHs can be constrained by future GW observations if the width of the mass function is smaller than about a hundredth of the mass.

Axion-like particles can source the baryon asymmetry of our Universe through spontaneous baryogenesis. Here we clarify that this is a generic outcome for essentially any coupling of an axion-like particle to the Standard Model, requiring only a non-zero velocity of the classical axion field while baryon or lepton number violating interactions are present in thermal bath. In particular, coupling the axions only to gluons is sufficient to generate a baryon asymmetry in the presence of electroweak sphalerons or the Weinberg operator. Deriving the transport equation for an arbitrary set of couplings of the axion-like particle, we provide a general framework in which these results can be obtained immediately. If all the operators involved are efficient, it suffices to solve an algebraic equation to obtain the final asymmetries. Otherwise one needs to solve a simple set of differential equations. This formalism clarifies some theoretical subtleties such as redundancies in the axion coupling to the Standard Model particles associated with a field rotation. We demonstrate how our formalism automatically evades potential pitfalls in the calculation of the final baryon asymmetry.

We point out that a model with scalar fields with a large nonminimal coupling to the Ricci scalar, such as Higgs inflation, can be regarded as a nonlinear sigma model (NLSM). The $\sigma$-meson, which is induced by quantum corrections, is identified as the scalaron in this model. Our understanding provides a novel alternative picture for the emergence of the scalaron, which was previously studied in the Jordan frame by computing the running of the term quadratic in the Ricci scalar. We demonstrate that quantum corrections inevitably induce other operators on top of this NLSM, which give rise to the $\sigma$-meson. We confirm that the $\sigma$-meson indeed corresponds to the scalaron. With the help of the $\sigma$-meson, the NLSM is UV-completed to be a linear sigma model (LSM). We show that the LSM only involves renormalizable interactions and hence its perturbative unitarity holds up to the Planck scale unless it hits a Landau pole, which is in agreement with the renormalizability of quadratic gravity.

Particle production in strong electromagnetic fields is a recurring theme in solid state physics, heavy ion collisions, early universe cosmology and formal quantum field theory. In this paper we discuss the Dirac equation in a background of parallel electric and magnetic fields. We review the Schwinger particle production rate, clarify the emergence of the chiral anomaly equation and compute the induced current of charged fermions. We distinguish the contributions from non-perturbative particle production, from the running of the gauge coupling constant and from non-linearities in the effective QED Lagrangian, and clarify how these contributions arise within a single framework. We apply these results to axion inflation. A Chern-Simons coupling between the pseudoscalar particle driving cosmic inflaton and an abelian gauge group induces a dual production of gauge fields and charged fermions. We show that the resulting scalar and gravitational wave power spectra strongly depend on the fermion mass.

The thermal decoupling description of dark matter (DM) and co-annihilating partners is reconsidered. If DM is realized at around the TeV-mass region or above, even the heaviest electroweak force carriers could act as long-range forces, leading to the existence of meta-stable DM bound states. The formation and subsequent decay of the latter further deplete the relic density during the freeze-out process on top of the Sommerfeld enhancement, allowing for larger DM masses. While so far the bound-state formation was described via the emission of an on-shell mediator ($W^{\pm}$, $Z$, $H$, $g$, photon or exotic), we point out that this particular process does not have to be the dominant scattering-bound state conversion channel in general. If the mediator is coupled in a direct way to any relativistic species present in the Early Universe, the bound-state formation can efficiently occur through particle scattering, where a mediator is exchanged virtually. To demonstrate that such a virtually stimulated conversion process can dominate the on-shell emission even for all temperatures, we analyze a simplified model where DM is coupled to only one relativistic species in the primordial plasma through an electroweak-scale mediator. We find that the bound-state formation cross section via particle scattering can exceed the on-shell emission by up to several orders of magnitude.

The coupling of an axion-like particle driving inflation to the Standard Model particle content through a Chern-Simons term generically sources a dual production of massless helical gauge fields and chiral fermions. We demonstrate that the interplay of these two components results in a highly predictive baryogenesis model, which requires no further ingredients beyond the Standard Model. If the helicity stored in the hyper magnetic field and the effective chemical potential induced by the chiral fermion production are large enough to avoid magnetic diffusion from the thermal plasma but small enough to sufficiently delay the chiral plasma instability, then the non-vanishing helicity survives until the electroweak phase transition and sources a net baryon asymmetry which is in excellent agreement with the observed value. If any of these two conditions is violated, the final baryon asymmetry vanishes. The observed baryon asymmetry can be reproduced if the energy scale of inflation is around $H_\text{inf} \sim 10^{10}$ - $10^{12}$ GeV with a moderate dependence on inflation model parameters.

We study the production of chiral fermions in a background of a strong non-abelian gauge field with a non-vanishing Chern-Pontryagin density. We discuss both pair production analogous to the Schwinger effect as well as asymmetric production through the chiral anomaly, sourced by the Chern-Pontryagin density. In abelian gauge theories one may nicely understand these processes by considering that the fermion dispersion relation forms discrete Landau levels. Here we extend this analysis to a non-abelian gauge theory, considering an intrinsically non-abelian isotropic and homogeneous SU(2) gauge field background with a non-vanishing Chern-Pontryagin density. We show that the asymmetric fermion production, together with a non-trivial vacuum contribution, correctly reproduces the chiral anomaly. This indicates that the usual vacuum subtraction scheme, imposing normal ordering, fails in this case. As a concrete example of this gauge field background, we consider chromo-natural inflation. Applying our analysis to this particular model, we compute the backreaction of the generated fermions on the gauge field background. This backreaction receives contributions both from the vacuum through a Coleman-Weinberg-type correction and from the fermion excitations through an induced current.

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 study the dual production of helical Abelian gauge fields and chiral fermions through the Chern-Simons (CS) coupling with a pseudo-scalar inflaton in the presence of a chiral anomaly. Through the CS term, the motion of the inflaton induces a tachyonic instability for one of the two helicities of the gauge field. We show that the resulting helical gauge field necessarily leads to the production of chiral fermions by deforming their Fermi sphere into discrete Landau levels. The population of the lowest Landau level leads to a chiral asymmetry as inferred from the chiral anomaly, while the higher levels are populated symmetrically through pair production. From the backreaction of the fermions on the gauge field production we derive a conservative but stringent upper bound on the magnitude of the gauge fields. Consequently, we find that the scalar perturbations sourced by these helical gauge fields, responsible for enhanced structure formation on small scales, get reduced significantly. We also discuss the fate of the primordial chiral asymmetry and of the helical gauge fields after inflation, and show that the instability in the chiral plasma tends to erase these primordial asymmetries. This result may impact scenarios where the baryon asymmetry of the Universe is connected to primordial magnetic fields.

We study I-balls/oscillons, which are long-lived, quasi-periodic, and spatially localized solutions in real scalar field theories. Contrary to the case of Q-balls, there is no evident conserved charge that stabilizes the localized configuration. Nevertheless, in many classical numerical simulations, it has been shown that they are extremely long-lived. In this paper, we clarify the reason for the longevity, and show how the exponential separation of time scales emerges dynamically. Those solutions are time-periodic with a typical frequency of a mass scale of a scalar field. This observation implies that they can be understood by the effective theory after integrating out relativistic modes. We find that the resulting effective theory has an approximate global U(1) symmetry reflecting an approximate number conservation in the non-relativistic regime. As a result, the profile of those solutions is obtained via the bounce method, just like Q-balls, as long as the breaking of the U(1) symmetry is small enough. We then discuss the decay processes of the I-ball/oscillon by the breaking of the U(1) symmetry, namely the production of relativistic modes via number violating processes. We show that the imaginary part is exponentially suppressed, which explains the extraordinary longevity of I-ball/oscillon. In addition, we find that there are some attractor behaviors during the evolution of I-ball/oscillon that further enhance the lifetime. The validity of our effective theory is confirmed by classical numerical simulations. Our formalism may also be useful to study condensates of ultra light bosonic dark matter, such as fuzzy dark matter, and axion stars, for instance.

Traditional computations of the dark matter (DM) relic abundance, for models where attractive self-interactions are mediated by light force-carriers and bound states exist, rely on the solution of a coupled system of classical on-shell Boltzmann equations. This idealized description misses important thermal effects caused by the tight coupling among force-carriers and other charged relativistic species. We develop for the first time a comprehensive ab-initio derivation for the description of DM long-range interactions in the presence of a hot and dense plasma background directly from non-equilibrium quantum field theory. Most importantly, the scattering and bound states get strongly mixed in the thermal plasma environment, representing a characteristic difference from a pure vacuum theory computation. The main result of this work is a novel differential equation for the DM number density, written down in a form which is manifestly independent under the choice of what one would interpret as a bound or a scattering state at finite temperature. The collision term, unifying the description of annihilation and bound state decay, turns out to have in general a non-quadratic dependence on the DM number density. This generalizes the form of the conventional Lee-Weinberg equation which is typically adopted to describe the freeze-out process. We prove that our general number density equation is consistent with previous literature results under certain limits.

Coherent oscillation of axions or axion-like particles may give rise to long-lived clumps, called axion stars, because of the attractive gravitational force or its self-interaction. Such a kind of configuration has been extensively studied in the context of oscillons without the effect of gravity, and its stability can be understood by an approximate conservation of particle number in a non-relativistic effective field theory (EFT). We extend this analysis to the case with gravity to discuss the longevity of axion stars and clarify the EFT expansion scheme in terms of gradient energy and Newton's constant. Our EFT is useful to calculate the axion star configuration and its classical lifetime without any ad hoc assumption. In addition, we derive a simple stability condition against small perturbations. Finally, we discuss the consistency of other non-relativistic effective field theories proposed in the literature.

We revise primordial black hole (PBH) formation in the axion-like curvaton model and investigate whether PBHs formed in this model can be the origin of the gravtitational wave (GW) signals detected by the Advanced LIGO. In this model, small-scale curvature perturbations with large amplitude are generated, which is essential for PBH formation. On the other hand, large curvature perturbations also become a source of primordial GWs by their second-order effects. Severe constraints are imposed on such GWs by pulsar timing array (PTA) experiments. We also check the consistency of the model with these constraints. In this analysis, it is important to take into account the effect of non-Gaussianity, which is generated easily in the curvaton model. We see that, if there are non-Gaussianities, the fixed amount of PBHs can be produced with a smaller amplitude of the primordial power spectrum.

Primordial Black Hole (PBH) is one of the leading non-particle candidates for dark matter (DM). Although several observations severely constrain the amount of PBHs, it is recently pointed out that there is an uncertainty on the microlensing constraints below $\sim 10^{-10} M_\odot$ which was ignored originally but may weaken the constraints significantly. In this paper, facing this uncertainty, we investigate the possibility that PBHs can make up all DM in a broad mass spectrum. Moreover, we propose a concrete inflation model which can simultaneously produce PBHs for all DM in a broad mass spectrum around $\mathcal O(10^{-13}) M_\odot$ and PBHs for LIGO events in a sharp mass spectrum at $\mathcal O(10) M_\odot$.

LIGO-Virgo collaboration has found black holes as heavy as $M \sim 30M_\odot$ through the detections of the gravitational waves emitted during their mergers. Primordial black holes (PBHs) produced by inflation could be an origin of such events. While it is tempting to presume that these PBHs constitute all Dark Matter (DM), there exists a number of constraints for PBHs with $\mathcal{O} (10) M_\odot$ which contradict with the idea of PBHs as all DM. Also, it is known that weakly interacting massive particle (WIMP) that is a common DM candidate is almost impossible to coexist with PBHs. These observations motivate us to pursue another candidate of DM. In this paper, we assume that the string axion solving the strong CP problem makes up all DM, and discuss the coexistence of string axion DM and inflationary PBHs for LIGO events.

We study the stability of the electroweak vacuum in low-scale inflation models whose Hubble parameter is much smaller than the instability scale of the Higgs potential. In general, couplings between the inflaton and Higgs are present, and hence we study effects of these couplings during and after inflation. We derive constraints on the couplings between the inflaton and Higgs by requiring that they do not lead to catastrophic electroweak vacuum decay, in particular, via resonant production of the Higgs particles.

Primordial black holes (PBHs) are one of the candidates to explain the gravitational wave (GW) signals observed by the LIGO detectors. Among several phenomena in the early Universe, cosmic inflation is a major example to generate PBHs from large primordial density perturbations. In this paper, we discuss the possibility to interpret the observed GW events as mergers of PBHs which are produced by cosmic inflation. The primordial curvature perturbation should be large enough to produce a sizable amount of PBHs and thus we have several other probes to test this scenario. We point out that the current pulsar timing array (PTA) experiments already put severe constraints on GWs generated via the second-order effects, and that the observation of the cosmic microwave background (CMB) puts severe restriction on its $\mu$ distortion. In particular, it is found that the scalar power spectrum should have a very sharp peak at $k \sim 10^{6}$ Mpc$^{-1}$ to fulfill the required abundance of PBHs while evading constraints from the PTA experiments together with the $\mu$ distortion. We propose a mechanism which can realize such a sharp peak. In the future, simple inflation models that generate PBHs via almost Gaussian fluctuations could be probed/excluded.

Following a new microlensing constraint on primordial black holes (PBHs) with $\sim10^{20}$--$10^{28}\,\mathrm{g}$~[1], we revisit the idea of PBH as all Dark Matter (DM). We have shown that the updated observational constraints suggest the viable mass function for PBHs as all DM to have a peak at $\simeq 10^{20}\,\mathrm{g}$ with a small width $\sigma \lesssim 0.1$, by imposing observational constraints on an extended mass function in a proper way. We have also provided an inflation model that successfully generates PBHs as all DM fulfilling this requirement.

We study non-thermal gravitino production in the minimal supergravity inflation. In this minimal model utilizing orthogonal nilpotent superfields, the particle spectrum includes only graviton, gravitino, inflaton, and goldstino. We find that a substantial fraction of the cosmic energy density can be transferred to the longitudinal gravitino due to non-trivial change of its sound speed. This implies either a breakdown of the effective theory during reheating or a serious gravitino problem.

We study particle production at the preheating era in inflation models with nonminimal coupling $\xi \phi^2R$ and quartic potential $\lambda \phi^4/4$ for several cases: real scalar inflaton, complex scalar inflaton and Abelian Higgs inflaton. We point out that the preheating proceeds much more violently than previously thought. If the inflaton is a complex scalar, the phase degree of freedom is violently produced at the first stage of preheating. If the inflaton is a Higgs field, the longitudinal gauge boson production is similarly violent. This is caused by a spike-like feature in the time dependence of the inflaton field, which may be understood as a consequence of the short time scale during which the effective potential or kinetic term changes suddenly. The produced particles typically have very high momenta $k \lesssim \sqrt{\lambda}M_\text{P}$. The production might be so strong that almost all the energy of the inflaton is carried away within one oscillation for $\xi^2\lambda \gtrsim {\mathcal O}(100)$. This may partly change the conventional understandings of the (p)reheating after inflation with the nonminimal coupling to gravity such as Higgs inflation. We also discuss the possibility of unitarity violation at the preheating stage.

Our electroweak vacuum may be metastable in light of the current experimental data of the Higgs/top quark mass. If this is really the case, high-scale inflation models require a stabilization mechanism of our vacuum during inflation. A possible candidate is the Higgs-inflaton/-curvature coupling because it induces an additional mass term to the Higgs during the slow roll regime. However, after the inflation, the additional mass term oscillates, and it can potentially destabilize our electroweak vacuum via production of large Higgs fluctuations during the inflaton oscillation era. In this paper, we study whether or not the Higgs-inflaton/-curvature coupling can save our vacuum by properly taking account of Higgs production during the preheating stage. We put upper bounds on the Higgs-inflaton/-curvature coupling, and discuss possible dynamics that might relax them.

We revisit the nonthermal gravitino production at the (p)reheating stage after inflation. Particular attention is paid to large field inflation models with a $\mathbb{Z}_2$ symmetry, for which the previous perturbative analysis is inapplicable; and inflation models with a stabilizer superfield, which have not been studied non-perturbatively. It is found that in single-superfield inflation models (without the stabilizer field), nonthermal production of the transverse gravitino can be cosmologically problematic while the abundance of the longitudinal gravitino is small enough. In multi-superfield inflation models (with the stabilizer field), production of the transverse and longitudinal gravitinos is significantly suppressed, and they are cosmologically harmless. We also clarify the relation between the background field method used in the preheating context and the standard perturbative decay method to estimate the gravitino abundance.

We revisit the compatibility between the chaotic inflation, which provides a natural solution to the initial condition problem, and the metastable electroweak vacuum, which is suggested by the results of LHC and the current mass measurements of top quark and Higgs boson. It is known that the chaotic inflation poses a threat to the stability of the electroweak vacuum because it easily generates large Higgs fluctuations during inflation or preheating and triggers the catastrophic vacuum decay. In this paper, we propose a simple cosmological solution in which the vacuum is stabilized during chaotic inflation, preheating and after that. This simple solution naturally predicts the formation of primordial black holes. We find interesting parameter regions where the present dark matter density is provided by them. Also, the thermal leptogenesis can be accommodated in our scenario.

We study production of light particles due to oscillation of the Hubble parameter or the scale factor. Any coherently oscillating scalar field, irrespective of its energy fraction in the universe, imprints such an oscillating feature on them. Not only the Einstein gravity but extended gravity models, such as models with non-minimal (derivative) coupling to gravity and $f(R)$ gravity, lead to oscillation of the scale factor. We present a convenient way to estimate the gravitational particle production rate in these circumstances. Cosmological implications of gravitational particle production, such as dark matter/radiation and moduli problem, are discussed. For example, if the theory is described solely by the standard model plus the Peccei-Quinn sector, the Starobinsky $R^2$ inflation may lead to observable amount of axion dark radiation.

It is known that the present electroweak vacuum is likely to be metastable and it may lead to a serious instability during/after inflation. We propose a simple solution to the problem of vacuum instability during/after inflation. If there is a moduli field which has Planck-suppressed interactions with the standard model fields, the Higgs quartic coupling in the early universe naturally takes a different value from the present one. A slight change of the quartic coupling in the early universe makes the Higgs potential absolutely stable and hence we are free from the vacuum instability during/after inflation.

We investigate the thermalization process of the Universe after inflation to determine the evolution of the effective temperature. The time scale of thermalization is found to be so long that it delays the evolution of the effective temperature, and the resulting maximal temperature of the Universe can be significantly lower than the one obtained in the literature. Our results clarify the finite density corrections to the effective potential of a scalar field and also processes of heavy particle production. In particular, we find that the maximum temperature of the Universe may be at most electroweak scale if the reheating temperature is as low as ${\cal O} (1)$ MeV, which implies that the electroweak symmetry may be marginally restored. In addition, it is noticeable that the dark matter may not be produced from thermal plasma in such a low reheating scenario, since the maximum temperature can be smaller than the conventional estimation by five orders of magnitude. We also give implications to the Peccei-Quinn mechanism and the Affleck-Dine baryogenesis.

We study cosmological evolution after inflation in models with non-minimal derivative coupling to gravity. The background dynamics is solved and particle production associated with rapidly oscillating Hubble parameter is studied in detail. In addition, production of gravitons through the non-minimal derivative coupling with the inflaton is studied. We also find that the sound speed squared of the scalar perturbation oscillates between positive and negative values when the non-minimal derivative coupling dominates over the minimal kinetic term. This may lead to an instability of this model. We point out that the particle production rates are the same as those in the Einstein gravity with the minimal kinetic term, if we require the sound speed squared is positive definite.

We point out that the inflaton inevitably couples to all non-conformally coupled matters gravitationally through an oscillation in the Hubble parameter or the cosmic scale factor. It leads to particle production during the inflaton oscillation regime, which is most efficient just after inflation. Moreover, the analysis is extended to the model with non-minimal inflaton couplings to gravity, in which the Hubble parameter oscillates more violently. We apply our results to the graviton production by the inflaton: gravitons are also produced just after inflation, but the non-minimal coupling does not induce inflaton decay into the graviton pair.

We consider background dynamics of generalized Galileon theories in the context of inflation, where gravity and inflaton are non-minimally coupled to each other. In the inflaton oscillation regime, the Hubble parameter and energy density oscillate violently in many cases, in contrast to the Einstein gravity with minimally coupled inflaton. However, we find that there is an adiabatic invariant in the inflaton oscillation regime in any generalized Galileon theory. This adiabatic invariant is useful in estimating the expansion law of the universe and also the particle production rate due to the oscillation of the Hubble parameter.

If a real scalar field is dominated by non-relativistic modes, then it approximately conserves its particle number and obeys an equation that governs a complex scalar field theory with a conserved global U(1) symmetry. From this fact, it is shown that the I-ball (oscillon) can be naturally understood as a projection (e.g., real part) of the non-relativistic Q-ball solution. In particular, we clarify that the stability of the I-ball is guaranteed by the U(1) symmetry in the corresponding complex scalar field theory as long as the non-relativistic condition holds. We also discuss the longevity of I-ball from the perspective of the complex scalar field in terms of U(1) charge violating processes.

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.

We revisit the cosmological aspects of axion models. In the high-scale inflation scenario, the Peccei-Quinn (PQ) symmetry is likely to be restored during/after inflation. If the curvature of the PQ scalar potential at the origin is smaller than its vacuum expectation value; for instance in a class of SUSY axion models, thermal inflation happens before the radial component of the PQ scalar (saxion) relaxes into the global minimum of the potential and the decay of saxion coherent oscillation would produce too much axion dark radiation. In this paper, we study how to avoid the overproduction of axion dark radiation with some concrete examples. We show that, by taking account of the finite-temperature dissipation effect appropriately, the overproduction constraint can be relaxed since the PQ scalar can take part in the thermal plasma again even after the PQ phase transition. We also show that it can be further relaxed owing to the late time decay of another heavy CP-odd scalar, if it is present.

We show that the local type non-Gaussianity in a class of curvaton models is suppressed, i.e. the non-linearity parameters f_{NL} and those related with higher order statistics can be at most O(1), even if the curvaton energy density is subdominant at the decay. This situation is naturally realized in a very simple curvaton potential with quadratic term plus quartic term.

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.

We propose an economical model in which a singlet Z$_2$-odd scalar field accounts for the primordial inflation and the present dark matter abundance simultaneously in the light of recent BICEP2 result. Interestingly, the reheating temperature and the thermal dark matter abundance are closely connected by the same interaction between the singlet scalar and the standard model Higgs. In addition, the reheating temperature turns out to be quite high, $T_\text{R} \gtrsim 10^{12}$ GeV, and hence the thermal leptogenesis is compatible with this model. Therefore, it can be one of the simplest cosmological scenarios.

We revisit the dynamics of the curvaton in detail taking account of effects from thermal environment, effective potential and decay/dissipation rate for general field values and couplings. We also consider the curvature perturbation generated through combinations of various effects: large scale modulation of the oscillation epoch, the effective dissipation rate and the timing at which the equation of state changes. In particular, we find that it tends to be difficult to explain the observed curvature perturbation by the curvaton mechanism without producing too large non-Gaussianity if the curvaton energy density is dissipated through thermal effects. Therefore, the interaction between the curvaton and light elements in the thermal plasma should be suppressed in order for the curvaton to survive the thermal dissipation.

The evolution of a coherently oscillating scalar field with $Z_2$ symmetry is studied in detail. We calculate the dissipation rate of the scalar field based on the closed time path formalism. Consequently, it is shown that the energy density of the coherent oscillation can be efficiently dissipated if the coupling constant is larger than the critical value, even though the scalar particle is stable due to the $Z_2$ symmetry.

We study the expansion law of the universe dominated by the oscillating scalar field with non-minimal derivative coupling to gravity as G^{\mu \nu} \partial_{\mu} \phi \partial_{\nu} \phi. In this system the Hubble parameter oscillates with a frequency of the effective mass of the scalar field, which formerly caused a difficulty in analyzing how the universe expands. We find an analytical solution for power law potentials and interpret the solution in an intuitive way by using a new invariant of the system. As a result, we find marginally accelerated expansion for the quadratic potential and no accelerated expansion for the potential with higher power.

We study in detail the dynamics of a scalar field in thermal bath with symmetry breaking potential. In particular, we focus on the process of trapping of a scalar field at an enhanced symmetry point through the thermal/non-thermal particle production, taking into account the interactions of produced particles with the standard model particles. As an explicit example, we revisit the saxion dynamics with an initial amplitude much larger than the Peccei-Quinn scale and show that the saxion trapping phenomenon happens for the most cases and it often leads to thermal inflation. We also study the saxion dynamics after thermal inflation, and it is shown that thermal dissipation effect on the saxion can relax the axion overproduction problem from the saxion decay.

There often appear coherently oscillating scalar fields in particle physics motivated cosmological scenarios, which may have rich phenomenological consequences. Scalar fields should somehow interact with background thermal bath in order to decay into radiation at an appropriate epoch, but introducing some couplings to the scalar field makes the dynamics complicated. We investigate in detail the dynamics of a coherently oscillating scalar field, which has renormalizable couplings to another field interacting with thermal background. The scalar field dynamics and its resultant abundance are significantly modified by taking account of following effects : (1) thermal correction to the effective potential, (2) dissipation effect on the scalar field in thermal bath, (3) non-perturbative particle production events and (4) formation of non-topological solitons. There appear many time scales depending on the scalar mass, amplitude, couplings and the background temperature, which make the efficiencies of these effects non-trivial.

We show that decaying axino dark matter with R-parity violation can explain the observed excess of the 130GeV gamma-ray line from the Galactic center in the Fermi data. The branching fraction of the axino decay into monochromatic photons can be O(1), and constraints from continuum gamma-rays and the anti-proton flux are ameliorated. The Peccei-Quinn scale of $O(10^{13}-10^{14})$\,GeV and the R-parity violation parameter of $O(10^{-12}-10^{-11})$ are cosmologically favored.

The inflaton must convert its energy into radiation after inflation, which, in a conventional scenario, is caused by the perturbative inflaton decay. This reheating process would be much more complicated in some cases: the decay products obtain masses from an oscillating inflaton and thermal environment, and hence the conventional reheating scenario can be modified. We study in detail processes of particle production from the inflaton, their subsequent thermalization and evolution of inflaton/plasma system by taking dissipation of the inflaton in a hot plasma into account. It is shown that the reheating temperature is significantly affected by these effects.

We consider a scalar field (called $\phi$) which is very weakly coupled to thermal bath, and study the evolution of its number density. We use the Boltzmann equation derived from the Kadanoff-Baym equations, assuming that the degrees of freedom in the thermal bath are well described as "quasi-particles." When the widths of quasi-particles are negligible, the evolution of the number density of $\phi$ is well governed by a simple Boltzmann equation, which contains production rates and distribution functions both evaluated with dispersion relations of quasi-particles with thermal masses. We pay particular attention to the case that dark matter is non-thermally produced by the decay of particles in thermal bath, to which the above mentioned formalism is applicable. When the effects of thermal bath are properly included, the relic abundance of dark matter may change by $O(10-100%)$ compared to the result without taking account of thermal effects.