We propose a scenario in which a strong Peccei-Quinn (PQ) symmetry breaking in the early universe results in large inhomogeneities of the initial QCD axion field value, leading to the formation of very dense axion bubbles. Some of the axion bubbles subsequently collapse into primordial black holes (PBHs). The spatially homogeneous part of the QCD axion explains dark matter of the universe, while the PBHs arising from the axion bubbles can explain the LIGO events or the seed of supermassive black holes. Interestingly, the mass of PBH is determined by the axion decay constant; for $f_a = 10^{17} (10^{16})$ GeV, the PBH mass is heavier than about $10 (10^4) M_\odot$. In addition, axion miniclusters are also formed from the axion bubbles more abundantly than PBHs, and their masses are expected to be heavier than in the usual scenario based on the spontaneous breaking of the PQ symmetry after inflation.

We revisit the ALP miracle scenario where the inflaton and dark matter are unified by a single axion-like particle (ALP). We first extend our previous analysis on the inflaton dynamics to identify the whole viable parameter space consistent with the CMB observation. Then, we evaluate the relic density of the ALP dark matter by incorporating uncertainties of the model-dependent couplings to the weak gauge bosons as well as the dissipation effect. The preferred ranges of the ALP mass and coupling to photons are found to be $0.01\lesssim m_\phi \lesssim 1$\,eV and $g_{\phi \gamma \gamma} = {\cal O}(10^{-11})$\,GeV$^{-1}$, which slightly depend on these uncertainties. Interestingly, the preferred regions are within reach of future solar axion helioscope experiments, IAXO and TASTE, and laser-based stim

Recently, an anomalous excess was found in the electronic recoil data collected at the XENON1T experiment. The excess may be explained by an axion-like particle (ALP) with mass of a few keV and a coupling to electron of $g_{ae} \sim 10^{-13}$, if the ALP constitutes all or some fraction of local dark matter (DM). In order to satisfy the X-ray constraint, the ALP coupling to photons must be significantly suppressed compared to that to electrons. This strongly suggests that the ALP has no anomalous couplings to photons, i.e., there is no U(1)$_{\rm PQ}$-U(1)$_{\rm em}$-U(1)$_{\rm em}$ anomaly. We show that such anomaly-free ALP DM predicts an X-ray line signal with a definite strength through the operator arising from threshold corrections. The strength of the signal, however, falls slightly short of the projected sensitivity of the ATHENA X-ray observatory unless the ALP mass is heavier than about (4-5) keV. The abundance of ALP DM can be explained by the misalignment mechanism, or by thermal production if it constitutes a part of DM. In particular, we find that the anomalous excess reported by the XENON1T experiment as well as the stellar cooling anomalies from white dwarfs and red giants can be explained simultaneously better when the ALP constitutes about 10% of DM. As concrete models, we revisit the leptophilic anomaly-free ALP DM considered in [Phys.Lett.B 734 (2014) 178-182] as well as an ALP model based on 2HDM.

We study the stochastic axion dark matter scenario in the axion landscape, where one of the axions is light and stable and therefore explains dark matter. If the axion mass at the potential is a typical value of the curvature along the direction, the potential can be well approximated by a quadratic mass term. On the other hand, if the axion mass happens to be suppressed in the vicinity of the minimum, the potential may be approximated by a quartic potential plus a suppressed quadratic one, for which the initial angle, and thus the axion abundance, can be significantly suppressed compared to the quadratic case. We delineate the viable parameter region by taking account of various observational constraints, and find that a broader range of the inflation scale is allowed. Also, if the curvature of the potential is suppressed over a certain range of the potential, the onset of coherent oscillations can be delayed. Then, the axion dark matter with a small decay constant is possible. We also discuss the {\pi}nflation mechanism to realize the hilltop initial condition in the stochastic axion scenario.

The abundance of the QCD axion is known to be suppressed if the Hubble parameter during inflation, $H_{\rm inf}$, is lower than the QCD scale, and if the inflation lasts sufficiently long. We show that the tight upper bound on the inflation scale can be significantly relaxed if the eternal old inflation is driven by the standard-model Higgs field trapped in a false vacuum at large field values. Specifically, $H_{\rm inf}$ can be larger than $100\,$GeV if the false vacuum is located above the intermediate scale. We also discuss the slow-roll inflation after the tunneling from the false vacuum to the electroweak vacuum.

We propose a class of single-field, slow-roll inflation models in which a typical number of $e$-folds can be extremely large. The key point is to introduce a very shallow local minimum near the top of the potential in a hilltop inflation model. In particular, a typical number of $e$-folds is enhanced if classical behavior dominates around the local minimum such that the inflaton probability distribution is drifted to the local minimum as a whole. After the inflaton escapes from the local minimum due to the stochastic dynamics, the ordinary slow-roll inflation follows and it can generate the primordial density perturbation consistent with observation. Interestingly, our scenario inherits the advantages of the old and new inflation: the typical $e$-folds can be extremely large as in the old inflation, and slow-roll inflation naturally follows after the stochastic regime as in the new inflation. In our numerical example, the typical number of $e$-folds can be as large as $10^{10^{10}}$, which is large enough for various light scalars such the QCD axion to reach the Bunch-Davies distribution.

We show that the initial misalignment angle of the QCD axion (or axion-like particles) can be set very close to $\pi$, if the QCD axion has a mixing with another heavy axion which induces the phase shift $\approx \pi$ after inflation. In the simplest case, the heavy axion plays the role of the inflaton, and we call such inflation as "$\pi$nflation." The basic idea was first proposed by Daido and the present authors in Ref. [1702.03284] in 2017 and more recently discussed in Ref. [1903.00462]. We show that the QCD axion with a decay constant $f_a \gtrsim 3 \times 10^9\,$ GeV can explain dark matter by the $\pi$nflation mechanism. A large fraction of the parameter region has an overlap with the projected sensitivity of ORGAN, MADMAX, TOORAD and IAXO. We also study implications for the effective neutrino species and isocurvature perturbations. The $\pi$nflation can provide an initial condition for the hilltop inflation in the axion landscape, and in a certain set-up, a chain of the hilltop inflation may take place.

Hidden monopole is a plausible dark matter candidate due to its stability, but its direct experimental search is extremely difficult due to feeble interactions with the standard model particles in the minimal form. Then, we introduce an axion, $a$, connecting the hidden monopole and the standard model particles and examine the current limits and future prospects of direct dark matter searches and beam-dump experiments. We find two parameter regions around $m_a = {\cal O}(10)$ MeV, $f_a = {\cal O}(10^{5})$ GeV and $m_a = {\cal O}(100)$ MeV, $f_a = {\cal O}(10^{4})$ GeV where monopole dark matter and the axion are respectively within the reach of the future experiments such as PICO-500 and SHiP. We also note that the hidden photons mainly produced by the axion decay contribute to dark radiation with $\Delta N_{\rm eff} \simeq 0.6$ which can relax the $H_0$ tension.

We study a hilltop inflation model where an axion-like particle (ALP) plays the role of the inflaton. We find that, for a broad class of potentials, the decay constant and the mass at the potential minimum satisfy the relation, $m_\phi \sim 10^{-6} f$, to explain the CMB normalization. The ALP is necessarily coupled to the standard model particles for successful reheating. The ALP with the above relation can be searched at beam dump experiments, e.g., the SHiP experiment, if the inflation scale is sufficiently low. In this case, the ALP decays through the interactions that led to the reheating of the Universe. In other words, the Big Bang may be probed at ground-based experiments.

We derive an anthropic bound on the extra neutrino species, $\Delta N_{\rm eff}$, based on the observation that a positive $\Delta N_{\rm eff}$ suppresses the growth of matter fluctuations due to the prolonged radiation dominated era, which reduces the fraction of matter that collapses into galaxies, hence, the number of observers. We vary $\Delta N_{\rm eff}$ and the positive cosmological constant while fixing the other cosmological parameters. We then show that the probability of finding ourselves in a universe satisfying the current bound is of order a few percents for a flat prior distribution. If $\Delta N_{\rm eff}$ is found to be close to the current upper bound or the value suggested by the $H_0$ tension, the anthropic explanation is not very unlikely. On the other hand, if the upper bound on $\Delta N_{\rm eff}$ is significantly improved by future observations, such simple anthropic consideration does not explain the small value of $\Delta N_{\rm eff}$. We also study simple models where dark radiation consists of relativistic particles produced by heavy scalar decays, and show that the prior probability distribution sensitively depends on the number of the particle species.

We find a class of solutions for a homogeneous and isotropic universe in which the initially expanding universe stops expanding, experiences contraction, and then expands again (the "bounce"), in the framework of Einstein gravity with a real scalar field without violating the null energy condition nor encountering any singularities. Two essential ingredients for the bouncing universe are the positive spatial curvature and the scalar potential which becomes flatter at large field values. Depending on the initial condition, either the positive curvature or the negative potential stops the cosmic expansion and begins the contraction phase. The flat potential plays a crucial role in triggering the bounce. After the bounce, the flat potential naturally allows the universe to enter the slow-roll inflation regime, thereby making the bouncing universe compatible with observations. If the e-folding of the subsequent inflation is just enough, a positive spatial curvature may be found in the future observations. Our scenario nicely fits with the creation of the universe from nothing, which leads to the homogeneous and isotropic universe with positive curvature. As a variant of the mechanism, we also find solutions representing a cyclic universe.

We show that the cosmological abundance of string axions is much smaller than naive estimates if the Hubble scale of inflation, $H_{\rm inf}$, is sufficiently low (but can still be much higher than the axion masses) and if the inflation lasts sufficiently long. The reason is that the initial misalignment angles of the string axions follow the Bunch-Davies distribution peaked at the potential minima. As a result, the cosmological moduli problem induced by the string axions can be significantly relaxed by low-scale inflation, and astrophysical and cosmological bounds are satisfied over a wide range of the mass without any fine-tuning of the initial misalignment angles. Specifically, the axion with its decay constant $f_\phi = 10^{16}$\,GeV satisfies the bounds over $10^{-18}{\rm \, eV} \lesssim m_\phi \lesssim 10{\rm\,TeV}$ for $H_{\rm inf} \lesssim 10{\rm\,keV}- 10^{6}$\,{\rm GeV}. We also discuss cases with multiple axions and the QCD axion.

We study if eternal inflation is realized while satisfying the recently proposed string Swampland criteria concerning the range of scalar field excursion, $|\Delta \phi| < \mathcal{D} \cdot M_{\rm P}$, and the potential gradient, $|\nabla V| > c \cdot V/M_{\rm P}$, where $\mathcal{D}$ and $c$ are constants of order unity, and $M_{\rm P}$ is the reduced Planck mass. We find that only the eternal inflation of chaotic type is possible for $c \sim {\cal O}(0.01)$ and $1/\mathcal{D} \sim {\cal O}(0.01)$, and that the Hubble parameter during the eternal inflation is parametrically close to the Planck scale, and is in the range of $2 \pi c \lesssim H_{\rm inf}/M_{\rm P} < 1/\sqrt{3}$.

We point out that the recently proposed Swampland conjecture on the potential gradient can lead to isocurvature perturbations of dark energy, if the quintessence field acquires large quantum fluctuations during high-scale inflation preferred by the conjecture. Also, if the quintessence field is coupled to a dark sector that contains dark matter, isocurvature perturbation of dark matter is similarly induced. Both isocurvature perturbations can be suppressed if the quintessence potential allows a tracker solution in the early Universe. We find that a vector field of mass $\lesssim {\cal O}(1)$ meV is an excellent dark matter candidate in this context, not only because the right abundance is known to be produced by quantum fluctuations during high-scale inflation without running afoul of isocurvature bounds, but also because its coupling to the quintessence does not spoil the flatness of the potential.

We revisit the adiabatic conversion between the QCD axion and axion-like particle (ALP) at level crossing, which can occur in the early universe as a result of the existence of a hypothetical mass mixing. This is similar to the Mikheyev-Smirnov-Wolfenstein effect in neutrino oscillations. After refining the conditions for the adiabatic conversion to occur, we focus on a scenario where the ALP produced by the adiabatic conversion of the QCD axion explains the observed dark matter abundance. Interestingly, we find that the ALP decay constant can be much smaller than the ordinary case in which the ALP is produced by the realignment mechanism. As a consequence, the ALP-photon coupling is enhanced by a few orders of magnitude, which is advantageous for the future ALP and axion-search experiments using the ALP-photon coupling.

We propose a scenario where the spontaneous breakdown of the Peccei-Quinn symmetry leads to monopole production. Both the axion and the monopole contribute to dark matter, and the Witten effect on the axion mass is a built-in feature. In the KSVZ-type axion model, seemingly different vacua are actually connected by the hidden gauge symmetry, which makes the axionic string unstable and separate into two Alice strings. The Alice string is attached to a single domain wall due to the QCD instanton effect, solving the domain-wall problem. This is in the same spirit of the Lazarides-Shafi mechanism, although the discrete Peccei-Quinn symmetry is not embedded into the center of the original gauge symmetry. In the DFSZ-type axion model, the domain-wall problem is avoided by the Witten effect. If the Peccei-Quinn symmetry is explicitly broken by a small amount, monopoles acquire a tiny electric charge and become mini-charged dark matter. Interestingly, the quality of the Peccei-Quinn symmetry is closely tied to darkness of dark matter.

We study the cosmological evolution of the QCD axion coupled to hidden photons. For a moderately strong coupling, the motion of the axion field leads to an explosive production of hidden photons by tachyonic instability. We use lattice simulations to evaluate the cosmological abundance of the QCD axion. In doing so, we incorporate the backreaction of the produced hidden photons on the axion dynamics, which becomes significant in the non-linear regime. We find that the axion abundance is suppressed by at most ${\cal O}(10^{2})$ for the decay constant $f_a = 10^{16}$ GeV, compared to the case without the coupling. For a sufficiently large coupling, the motion of the QCD axion becomes strongly damped, and as a result, the axion abundance is enhanced. Our results show that the cosmological upper bound on the axion decay constant can be relaxed by a few hundred for a certain range of the coupling to hidden photons.

We show that the axion coupling to photons can be enhanced in simple models with a single Peccei-Quinn field, if the gauge coupling unification is realized by a large kinetic mixing $\chi = {\cal O}(0.1)$ between hypercharge and unbroken hidden U(1)$_H$. The key observation is that the U(1)$_H$ gauge coupling should be rather strong to induce such large kinetic mixing, leading to enhanced contributions of hidden matter fields to the electromagnetic anomaly. We find that the axion-photon coupling is enhanced by about a factor of 10-100 with respect to the GUT-axion models with $E/N = 8/3$.

We calculate a dipole-dipole potential between fermions mediated by a light pseudoscalar, axion, paying a particular attention to the overall sign. While the sign of the potential is physical and important for experiments to discover or constrain the axion coupling to fermions, there is often a sign error in the literature. The purpose of this short note is to clarify the sign issue of the axion-mediated dipole-dipole potential. As a by-product, we find a sign change of the dipole-dipole potenital due to the different spin of the mediating particle.

The string theory predicts many light fields called moduli and axions, which cause a cosmological problem due to the overproduction of their coherent oscillation after inflation. One of the prominent solutions is an adiabatic suppression mechanism, which, however, is non-trivial to achieve in the case of axions because it necessitates a large effective mass term which decreases as a function of time. The purpose of this paper is twofold. First, we provide an analytic method to calculate the cosmological abundance of coherent oscillation in a general situation under the adiabatic suppression mechanism. Secondly, we apply our method to some concrete examples, including the one where a string axion acquires a large effective mass due to the Witten effect in the presence of hidden monopoles.

We propose a scenario where both inflation and dark matter are described by a single axion-like particle (ALP) in a unified manner. In a class of the minimal axion hilltop inflation, the effective masses at the maximum and mimimum of the potential have equal magnitude but opposite sign, so that the ALP inflaton is light both during inflation and in the true vacuum. After inflation, most of the ALPs decay and evaporate into plasma through a coupling to photons, and the remaining ones become dark matter. We find that the observed CMB and matter power spectrum as well as the dark matter abundance point to an ALP of mass $m_\phi = {\cal O}(0.01)$ eV and the axion-photon coupling $g_{\phi \gamma \gamma} ={\cal O}(10^{-11})$GeV$^{-1}$: the ALP miracle. The suggested parameter region is within the reach of the next generation axion helioscope, IAXO. Furthermore, thermalized ALPs contribute to hot dark matter and its abundance is given in terms of the effective number of extra neutrino species, $\Delta N_{\rm eff} \simeq 0.03$, which can be tested by the future CMB experiments. We also discuss a case with multiple ALPs, where the coupling to photons can be enhanced in the early Universe by an order of magnitude or more, which enlarges the parameter space for the ALP miracle. The heavy ALP plays a role of the waterfall field in hybrid inflation, and reheats the Universe, and it can be searched for in various experiments such as SHiP.

If there are a plethora of axions in nature, they may have a complicated potential and create an axion landscape. We study a possibility that one of the axions is so light that it is cosmologically stable, explaining the observed dark matter density. In particular we focus on a case in which two (or more) shift-symmetry breaking terms conspire to make the axion sufficiently light at the potential minimum. In this case the axion has a flat-bottomed potential. In contrast to the case in which a single cosine term dominates the potential, the axion abundance as well as its isocurvature perturbations are significantly suppressed. This allows an axion with a rather large mass to serve as dark matter without fine-tuning of the initial misalignment, and further makes higher-scale inflation to be consistent with the scenario.

We evaluate the impact of domain-wall annihilation on the currently ongoing and planned gravitational wave experiments, including a case in which domain walls experience a frictional force due to interactions with the ambient plasma. We show the sensitivity reach in terms of physical parameters, namely, the wall tension and the annihilation temperature. We find that a Higgs portal scalar, which stabilizes the Higgs potential at high energy scales, can form domain walls whose annihilation produces a large amount of gravitational waves within the reach of the advanced LIGO experiment (O5). Domain wall annihilation can also generate baryon asymmetry if the scalar is coupled to either SU(2)$_L$ gauge fields or the $(B-L)$ current. This is a variant of spontaneous baryogenesis, but it naturally avoids the isocurvature constraint due to the scaling behavior of the domain-wall evolution. We delineate the parameter space where the domain-wall baryogenesis works successfully and discuss its implications for the gravitational wave experiments.

A QCD axion with a time-dependent decay constant has been known to be able to accommodate high-scale inflation without producing topological defects or too large isocurvature perturbations on CMB scales. We point out that a dynamical decay constant also has the effect of enhancing the small-scale axion isocurvature perturbations. The enhanced axion perturbations can even exceed the periodicity of the axion potential, and thus lead to the formation of axionic domain walls. Unlike the well-studied axionic walls, the walls produced from the enhanced perturbations are not bounded by cosmic strings, and thus would overclose the universe independently of the number of degenerate vacua along the axion potential.

We study a Higgs inflation model with a running kinetic term, taking account of the renormalization group evolution of relevant coupling constants. Specifically we study two types of the running kinetic Higgs inflation, where the inflaton potential is given by the quadratic or linear term potential in a frame where the Higgs field is canonically normalized. We solve the renormalization group equations at two-loop level and calculate the scalar spectral index and the tensor-to-scalar ratio. We find that, even if the renormalization group effects are included, the quadratic inflation is ruled out by the CMB observations, while the linear one is still allowed.

We study the formation and evolution of topological defects in an aligned axion model with multiple Peccei-Quinn scalars, where the QCD axion is realized by a certain combination of the axions with decay constants much smaller than the conventional Peccei-Quinn breaking scale. When the underlying U(1) symmetries are spontaneously broken, the aligned structure in the axion field space exhibits itself as a complicated string-wall network in the real space. We find that the string-wall network likely survives until the QCD phase transition if the number of the Peccei-Quinn scalars is greater than two. The string-wall system collapses during the QCD phase transition, producing a significant amount of gravitational waves in the nano-Hz range at present. The typical decay constant is constrained to be below O(100) TeV by the pulsar timing observations, and the constraint will be improved by a factor of 2 in the future SKA observations.

We show that the seesaw mechanism as well as leptogenesis are natural outcomes of a viable chaotic inflation in supergravity. The inflation model contains two superfields, the inflaton and stabilizer fields, which, being singlets under the standard model gauge symmetry, naturally couple to the lepton and Higgs doublets. The inflaton decays into leptons and Higgs fields, and the reheating temperature is predicted to be of O(10^{13}) GeV, for which thermal leptogenesis is possible. On the other hand, gravitinos are copiously produced, and various solutions to the gravitino problem are discussed. We also argue that, if the shift symmetry of the inflaton is explicitly broken down to a discrete one, neutrino Yukawa couplings are periodic in the inflaton field, and masses of leptons and Higgs do not blow up even if the inflaton takes super-Planckian field values. The inflaton potential is given by a sum of sinusoidal functions with different height and periodicity, the so-called multi-natural inflation. We show that the predicted scalar spectral index and tensor-to-scalar ratio lie in the region favored by the Planck data.

We study the Witten effect of hidden monopoles on the QCD axion dynamics, and show that its abundance as well as isocurvature perturbations can be significantly suppressed if there is a sufficient amount of hidden monopoles. When the hidden monopoles make up a significant fraction of dark matter, the Witten effect suppresses the abundance of axion with the decay constant smaller than $10^{12}$ GeV. The cosmological domain wall problem of the QCD axion can also be avoided, relaxing the upper bound on the decay constant when the Peccei-Quinn symmetry is spontaneously broken after inflation.

We consider QCD axion models where the Peccei-Quinn symmetry is badly broken by a larger amount in the past than in the present, in order to avoid the axion isocurvature problem. Specifically we study supersymmetric axion models where the Peccei-Quinn symmetry is dynamically broken by either hidden gauge interactions or the $SU(3)_c$ strong interactions whose dynamical scales are temporarily enhanced by the dynamics of flat directions. The former scenario predicts a large amount of self-interacting dark radiation as the hidden gauge symmetry is weakly coupled in the present Universe. We also show that the observed amount of baryon asymmetry can be generated by the QCD axion dynamics via spontaneous baryogenesis. We briefly comment on the case in which the PQ symmetry is broken by a non-minimal coupling to gravity.

The simplest two-field completion of natural inflation has a regime in which both fields are active and in which its predictions are within the Planck 1-$\sigma$ confidence contour. We show this for the original model of natural inflation, in which inflation is achieved through the explicit breaking of a U(1) symmetry. We consider the case in which the mass coming from explicit breaking of this symmetry is comparable to that from spontaneous breaking, which we show is consistent with a hierarchy between the corresponding energy scales. While both masses are comparable when the observable modes left the horizon, the mass hierarchy is restored in the last e-foldings of inflation, rendering the predictions consistent with the isocurvature bounds. For completeness, we also study the predictions for the case in which there is a large hierarchy of masses and an initial period of inflation driven by the (heavy) radial field.

We propose a variant scenario of spontaneous baryogenesis from asymmetric inflaton based on current-current interactions between the inflaton and matter fields with a non-zero B-L charge. When the inflaton starts to oscillate around the minimum after inflation, it may lead to excitation of a CP-odd component, which induces an effective chemical potential for the B-L number through the current-current interactions. We study concrete inflation models and show that the spontaneous baryogenesis scenario can be naturally implemented in the chaotic inflation in supergravity.

We study a level crossing between the QCD axion and an axion-like particle, focusing on the recently found phenomenon, the axion roulette, where the axion-like particle runs along the potential, passing through many crests and troughs, until it gets trapped in one of the potential minima. We perform detailed numerical calculations to determine the parameter space where the axion roulette takes place, and as a result domain walls are likely formed. The domain wall network without cosmic strings is practically stable, and it is nothing but a cosmological disaster. In a certain case, one can make domain walls unstable and decay quickly by introducing an energy bias without spoiling the Peccei-Quinn solution to the strong CP problem.

We propose a possibility that the inflaton potential is significantly modified after inflation due to heavy field dynamics. During inflation there may be a heavy scalar field stabilized at a value deviated from the low-energy minimum. As the heavy field moves to the low-energy minimum, the inflaton potential could be significantly modified. In extreme cases, the inflaton potential vanishes and the inflaton becomes almost massless at some time after inflation. Such transition of the inflaton potential has interesting implications for primordial density perturbations, reheating, creation of unwanted relics, dark radiation, and experimental search for light degrees of freedom. To be concrete, we consider a chaotic inflation in supergravity where the inflaton mass parameter is promoted to a modulus field, finding that the inflaton becomes stable after the transition and contributes to dark matter. Another example is the new inflation by the MSSM Higgs field which acquires a large expectation value just after inflation, but it returns to the origin after the transition and settles down at the electroweak vacuum. Interestingly, the smallness of the electroweak scale compared to the Planck scale is directly related to the flatness of the inflaton potential.

We propose a new scenario of baryogenesis, in which annihilation of axion domain walls generates a sizable baryon asymmetry. Successful baryogenesis is possible for a wide range of the axion mass and decay constant, $m \simeq 10^8 -10^{13}$ GeV and $f \simeq 10^{13} - 10^{16}$ GeV. Baryonic isocurvature perturbations are significantly suppressed in our model, in contrast to various spontaneous baryogenesis scenarios in the slow-roll regime. In particular, the axion domain wall baryogenesis is consistent with high-scale inflation which generates a large tensor-to-scalar ratio within the reach of future CMB B-mode experiments. We also discuss the gravitational waves produced by the domain wall annihilation and its implications for the future gravitational wave experiments.

Real scalar fields are known to fragment into spatially localized and long-lived solitons called oscillons or $I$-balls. We prove the adiabatic invariance of the oscillons/$I$-balls for a potential that allows periodic motion even in the presence of non-negligible spatial gradient energy. We show that such potential is uniquely determined to be the quadratic one with a logarithmic correction, for which the oscillons/$I$-balls are absolutely stable. For slightly different forms of the scalar potential dominated by the quadratic one, the oscillons/$I$-balls are only quasi-stable, because the adiabatic charge is only approximately conserved. We check the conservation of the adiabatic charge of the $I$-balls in numerical simulation by slowly varying the coefficient of logarithmic corrections. This unambiguously shows that the longevity of oscillons/$I$-balls is due to the adiabatic invariance.

We point out that domain wall formation is a more common phenomenon in the Axiverse than previously thought. Level crossing could take place if there is a mixing between axions, and if some of the axions acquire a non-zero mass through non-perturbative effects as the corresponding gauge interactions become strong. The axion potential changes significantly during the level crossing, which affects the axion dynamics in various ways. We find that, if there is a mild hierarchy in the decay constants, the axion starts to run along the valley of the potential, passing through many crests and troughs, until it gets trapped in one of the minima; the {\it axion roulette}. The axion dynamics exhibits a chaotic behavior during the oscillations, and which minimum the axion is finally stabilized is highly sensitive to the initial misalignment angle. Therefore, the axion roulette is considered to be accompanied by domain wall formation. The cosmological domain wall problem can be avoided by introducing a small bias between the vacua. We discuss cosmological implications of the domain wall annihilation for baryogenesis and future gravitational wave experiments.

The effective potential for the Standard Model Higgs field allows two quasi-degenerate vacua; one is our vacuum at the electroweak scale, while the other is at a much higher scale. The latter minimum may be at a scale much smaller than the Planck scale, if the potential is lifted by new physics. This gives rise to a possibility of domain wall formation after inflation. If the high-scale minimum is a local minimum, domain walls are unstable and disappear through violent annihilation processes, producing a significant amount of gravitational waves. We estimate the amount of gravitational waves produced from unstable domain walls in the Higgs potential and discuss detectability with future experiments.

Multiple axions form a landscape in the presence of various shift symmetry breaking terms. Eternal inflation populates the axion landscape, continuously creating new universes by bubble nucleation. Slow-roll inflation takes place after the tunneling event, if a very flat direction with a super-Planckian decay constant arises due to the alignment mechanism. We study the vacuum structure as well as possible inflationary dynamics in the axion landscape scenario, and find that the inflaton dynamics is given by either natural or multi-natural inflation. In the limit of large decay constant, it is approximated by the quadratic chaotic inflation, which however is disfavored if there is a pressure toward shorter duration of inflation. Therefore, if the spectral index and the tensor-to-scalar ratio turn out to be different from the quadratic chaotic inflation, there might be observable traces of the bubble nucleation. Also, the existence of small modulations to the inflaton potential is a common feature in the axion landscape, which generates a sizable and almost constant running of the scalar spectral index over CMB scales. Non-Gaussianity of equilateral type can also be generated if some of the axions are coupled to massless gauge fields.

We propose an extension of natural inflation, where the inflaton potential is a general periodic function. Specifically, we study elliptic inflation where the inflaton potential is given by Jacobi elliptic functions, Jacobi theta functions or the Dedekind eta function, which appear in gauge and Yukawa couplings in the string theories compactified on toroidal backgrounds. We show that in the first two cases the predicted values of the spectral index and the tensor-to-scalar ratio interpolate from natural inflation to exponential inflation such as $R^2$- and Higgs inflation and brane inflation, where the spectral index asymptotes to $n_s = 1-2/N \simeq 0.967$ for the e-folding number $N = 60$. We also show that a model with the Dedekind eta function gives a sizable running of the spectral index due to modulations in the inflaton potential. Such elliptic inflation can be thought of as a specific realization of multi-natural inflation, where the inflaton potential consists of multiple sinusoidal functions. We also discuss examples in string theory where Jacobi theta functions and the Dedekind eta function appear in the inflaton potential.

We study in detail MSW-like resonant conversions of QCD axions into hidden axions, including cases where the adiabaticity condition is only marginally satisfied, and where anharmonic effects are non-negligible. When the resonant conversion is efficient, the QCD axion abundance is suppressed by the hidden and QCD axion mass ratio. We find that, when the resonant conversion is incomplete due to a weak violation of the adiabaticity, the CDM isocurvature perturbations can be significantly suppressed, while non-Gaussianity of the isocurvature perturbations generically remain unsuppressed. The isocurvature bounds on the inflation scale can therefore be relaxed by the partial resonant conversion of the QCD axions into hidden axions.

The measured values of the Higgs and top masses and of the strong gauge coupling constant point to the near-criticality of the Standard Model, where two vacua at the electroweak and Planck scales are quasi-degenerate. We argue that the criticality is required by the occurrence of an eternal topological inflation induced by the Higgs potential. The role of this inflation is to continuously create sufficiently flat and homogeneous Universe, providing the necessary initial condition for the subsequent slow-roll inflation that generates the density perturbations of the right magnitude. While the condition for the topological Higgs inflation is only marginally satisfied in the Standard Model, it can be readily satisfied if one introduces the right-handed neutrinos and/or the non-minimal coupling to gravity; currently unknown quantum gravity corrections to the potential may also help. We also discuss the $B-L$ Higgs inflation as a possible origin of the observed density perturbations. Its necessary initial condition, the restored $B-L$ symmetry, can be naturally realized by the preceding topological Higgs inflation.

We consider a possibility that one of the flat directions in the minimal supersymmetric standard model plays the role of the inflaton field and realizes large-field inflation. This is achieved by introducing a generalized shift symmetry on the flat direction, which enables us to control the inflaton potential over large field values. After inflation, higher order terms allowed by the generalized shift symmetry automatically cause a helical motion of the field to create the baryon number of the universe, while baryonic isocurvature fluctuations are suppressed.

The BICEP2 experiment determined the Hubble parameter during inflation to be about $10^{14}$ GeV. Such high inflation scale is in tension with the QCD axion dark matter if the Peccei-Quinn (PQ) symmetry remains broken during and after inflation, because too large axion isocurvature perturbations would be generated. The axion isocurvature perturbations can be suppressed if the axion acquires a sufficiently heavy mass during inflation. We show that this is realized if the PQ symmetry is explicitly broken down to a discrete symmetry and if the breaking is enhanced during inflation. We also show that, even when the PQ symmetry becomes spontaneously broken after inflation, such a temporarily enhanced PQ symmetry breaking relaxes the constraint on the axion decay constant.

Hidden axions may be coupled to the standard model particles through a kinetic or mass mixing with QCD axion. We study a scenario in which a hidden axion constitutes a part of or the whole of dark matter and decays into photons through the mixing, explaining the 3.5 keV X-ray line signal. Interestingly, the required long lifetime of the hidden axion dark matter can be realized for the QCD axion decay constant at an intermediate scale, if the mixing is sufficiently small. In such a two component dark matter scenario, the primordial density perturbations of the hidden axion can be highly non-Gaussian, leading to a possible dispersion in the X-ray line strength from various galaxy clusters and near-by galaxies. We also discuss how the parallel and orthogonal alignment of two axions affects their couplings to gauge fields.

We revisit a polynomial chaotic inflation model in supergravity which we proposed soon after the Planck first data release. Recently some issues have been raised in Ref.[12], concerning the validity of our polynomial chaotic inflation model. We study the inflaton dynamics in detail, and confirm that the inflaton potential is very well approximated by a polynomial potential for the parameters of our interest in any practical sense, and in particular, the spectral index and the tensor-to-scalar ratio can be estimated by single-field approximation. This justifies our analysis of the polynomial chaotic inflation in supergravity.

We propose a landscape of many axions, where the axion potential receives various contributions from shift symmetry breaking effects. We show that the existence of the axion with a super-Planckian decay constant is very common in the axion landscape for a wide range of numbers of axions and shift symmetry breaking terms, because of the accidental alignment of axions. The effective inflation model is either natural or multi-natural inflation in the axion landscape, depending on the number of axions and the shift symmetry breaking terms. The tension between BICEP2 and Planck could be due to small modulations to the inflaton potential or steepening of the potential along the heavy axions after the tunneling. The total duration of the slow-roll inflation our universe experienced is not significantly larger than $60$ if the typical height of the axion potentials is of order $(10^{16-17}{\rm \,GeV})^4$.

We revisit large field inflation models with modulations in light of the recent discovery of the primordial B-mode polarization by the BICEP2 experiment, which, when combined with the Planck + WP + highL data, gives a strong hint for additional suppression of the CMB temperature fluctuations at small scales. Such a suppression can be explained by a running spectral index. In fact, it was pointed out by two of the present authors (TK and FT) that the existence of both tensor mode perturbations and a sizable running of the spectral index is a natural outcome of large inflation models with modulations. We find that this holds also in the recently proposed multi-natural inflation, in which the inflaton potential consists of multiple sinusoidal functions and therefore the modulations are a built-in feature.

We argue that there may arise a relatively flat inflaton potential over super-Planckian field values with an approximate shift symmetry, if the coefficients of the kinetic terms for many singlet scalars are subject to a certain random distribution. The inflaton potential generically contains various shift-symmetry breaking terms, leading to a possibly large deviation of the predicted values of the spectral index and tensor-to-scalar ratio from those of the simple quadratic chaotic inflation. We revisit a polynomial chaotic inflation in supergravity as such.

We show that the standard model Higgs field can realize the quadratic chaotic inflation, if the kinetic term is significantly modified at large field values. This is a simple realization of the so-called running kinetic inflation. The point is that the Higgs field respects an approximate shift symmetry at high energy scale. The tensor-to-scalar ratio is predicted to be $r \simeq 0.13 - 0.16$, which nicely explains the primordial B-mode polarization, $r=0.20^{+0.07}_{-0.05}$, recently discovered by the BICEP2 experiment. In particular, allowing small modulations induced by the shift symmetry breaking, the negative running spectral index can also be induced. The reheating temperature is expected to be so high that successful thermal leptogenesis is possible. The suppressed quartic coupling of the Higgs field at high energy scales may be related to the Higgs chaotic inflation.

If inflation scale is high, light scalars acquire large quantum fluctuations during inflation.If sufficiently long-lived, they will give rise to CDM isocurvature perturbations, which are highly constrained by the Planck data. Focusing on string axions as such light scalars, we show that thermal inflation can provide a sufficiently large entropy production to dilute the CDM isocurvature perturbations. Importantly, efficient dilution is possible for the string axions, because effectively no secondary coherent oscillations are induced at the end of thermal inflation, in contrast to the moduli fields. We also study the viability of the axion dark matter with mass of about 7 keV as the origin of the 3.5 keV X-ray line excess, in the presence of large entropy production.

We propose a multi-natural inflation model in which the single-field inflaton potential consists of two or more sinusoidal potentials that are comparable in size, but have different periodicity with a possible non-zero relative phase. The model is versatile enough to realize both large-field and small-field inflation. We show that, in a model with two sinusoidal potentials, the predicted values of the spectral index and tensor-to-scalar ratio lie within the $1\sigma$ region of the Planck data. In particular, there is no lower bound on the decay constants in contrast to the original natural inflation. We also show that, in a certain limit, multi-natural inflation can be approximated by a hilltop quartic inflation model.

We propose a scenario in which extra light fermions in an $E_6$-inspired U(1) extension of the standard model constitute the dark matter, as a simple variation of our model for dark radiation presented in 2010. Interestingly, for the light fermions of mass about 7 keV, they radiatively decay into active neutrinos and photons with a lifetime in the range of $10^{27}-10^{28}$ seconds, which naturally explains the recently discovered 3.5 keV X-ray line signal.

We show that the recently proposed multi-natural inflation can be realized within the framework of 4D ${\cal N}=1$ supergravity. The inflaton potential mainly consists of two sinusoidal potentials that are comparable in size, but have different periodicity with a possible non-zero relative phase. For a sub-Planckian decay constant, the multi-natural inflation model is reduced to axion hilltop inflation. We show that, taking into account the effect of the relative phase, the spectral index can be increased to give a better fit to the Planck results, with respect to the hilltop quartic inflation. We also consider a possible UV completion based on a string-inspired model. Interestingly, the Hubble parameter during inflation is necessarily smaller than the gravitino mass, avoiding possible moduli destabilization. Reheating processes as well as non-thermal leptogenesis are also discussed.

The most naive interpretation of the BICEP2 data is the chaotic inflation by an inflaton with a quadratic potential. When combined with supersymmetry, we argue that the inflaton plays the role of right-handed scalar neutrino based on rather general considerations. The framework suggests that the right-handed sneutrino tunneled from a false vacuum in a landscape to our vacuum with a small negative curvature and suppressed scalar perturbations at large scales.

Gravitinos are generically produced by inflaton decays, which place tight constraints on inflation models as well as supersymmetry breaking scale. We revisit the gravitino production from decays of the inflaton and the supersymmetry breaking field, based on a chaotic inflation model suggested by the recent BICEP2 result. We study cosmological constraints on thermally and non-thermally produced gravitinos for a wide range of the gravitino mass, and show that there are only three allowed regions of the gravitino mass: $m_{3/2}\lesssim 16$eV, $m_{3/2}\simeq 10$--$1000$TeV and $m_{3/2}\gtrsim 10^{13}$GeV.

We revisit the recently proposed multi-natural inflation and its realization in supergravity in light of the BICEP2 results. Multi-natural inflation is a single-field inflation model where the inflaton potential consists of multiple sinusoidal functions, and it is known that a sizable running spectral index can be generated, which relaxes the tension between the BICEP2 and the Planck results. In this paper we show that multi-natural inflation can accommodate a wide range of values of $(n_s, r)$, including the spectral index close to or even above unity. This will be be favored if the tension is resolved by other sources such as dark radiation, hot dark matter, or non-zero neutrino mass. We also discuss the implications for the implementation in string theory.

The recent unidentified 3.5 keV X-ray line signal can be explained by decaying moduli dark matter with a cutoff scale one order of magnitude smaller than the Planck scale. We show that such modulus field with the low cutoff scale follows a time-dependent potential minimum and its abundance is reduced by the adiabatic suppression mechanism. As a result the modulus abundance can naturally be consistent with the observed dark matter abundance without any fine-tuning of the initial oscillation amplitude.

We propose a scenario where the saxion dominates the energy density of the Universe and reheats the standard model sector via the dilatonic coupling, while its axionic partner contributes to dark matter decaying into photons via the same operator in supersymmetry. Interestingly, for the axion mass $m_a \simeq 7$keV and the decay constant $f_a \simeq 10^{14-15}$GeV, the recently discovered X-ray line at $3.5$keV in the XMM Newton X-ray observatory data can be explained. We discuss various cosmological aspects of the $7$keV axion dark matter such as the production of axion dark matter, the saxion decay process, hot dark matter and isocurvature constraints on the axion dark matter, and the possible baryogenesis scenarios.

The recently discovered X-ray line at about $3.5\,$keV can be explained by sterile neutrino dark matter with mass, $m_s \simeq 7\,$keV, and the mixing, $\sin^2 2\theta \sim 10^{-10}$. Such sterile neutrino is more long-lived than estimated based on the seesaw formula, which strongly suggests an extra flavor structure in the seesaw sector. We show that one can explain both the small mass and the longevity based on the split flavor mechanism where the breaking of flavor symmetry is tied to the breaking of the $B-L$ symmetry. In a supersymmetric case we find that the $7\,$keV sterile neutrino implies the gravitino mass about $100\,$TeV.

We propose a chaotic inflation model in which the lightest right-handed sneutrino serves as the inflaton and the predicted values of the spectral index and tensor-to-scalar ratio are consistent with the Planck data. Interestingly, the observed magnitude of primordial density perturbations is naturally explained by the inflaton mass of order $10^{13}$GeV, which is close to the right-handed neutrino mass scale suggested by the seesaw mechanism and the neutrino oscillation experiments. We find that the agreement of the two scales becomes even better in the neutrino mass anarchy. We show that the inflation model can be embedded into supergravity and discuss thermal history of the Universe after inflation such as non-thermal leptogenesis by the right-handed sneutrino decays and the modulus dynamics.

We revisit spectator field models including curvaton and modulated reheating scenarios, specifically focusing on their viability in the new Planck era, based on the derived expression for the spectral index in general spectator field models. Importantly, the recent Planck observations give strong preference to a red-tilted power spectrum, while the spectator field models tend to predict a scale-invariant one. This implies that, during inflation, either (i) the Hubble parameter varies significantly as in chaotic inflation, or (ii) a scalar potential for the spectator field has a relatively large negative curvature. Combined with the tight constraint on the non-Gaussianity, the Planck data provides us with rich implications for various spectator field models.

The presence of a hot dark matter component has been hinted at 3 sigma by a combination of the results from different cosmological observations. We examine a possibility that pseudo Nambu-Goldstone bosons account for both hot and cold dark matter components. We show that the QCD axions can do the job for the axion decay constant f_a < O(10^10) GeV, if they are produced by the saxion decay and the domain wall annihilation. We also investigate the cases of thermal QCD axions, pseudo Nambu-Goldstone bosons coupled to the standard model sector through the Higgs portal, and axions produced by modulus decay.

We revisit the isocurvature density perturbations induced by quantum fluctuations of the axion field by extending a recently developed analytic method and approximations to a time-dependent scalar potential, which enables us to follow the evolution of the axion until it starts to oscillate. We find that, as the initial misalignment angle approaches the hilltop of the potential, the isocurvature perturbations become significantly enhanced, while the non-Gaussianity parameter increases slowly but surely. As a result, the isocurvature constraint on the inflation scale is tightened as Hinf < O(100) GeV for the axion decay constant fa < 10^10 GeV, near the smaller end of the axion dark matter window. We also derive useful formulae for the power spectrum and non-Gaussianity of the isocurvature perturbations.

Sterile neutrino dark matter of mass O(1-10) keV decays into an active neutrino and an X-ray photon, and the non-observation of the corresponding X-ray line requires the sterile neutrino to be more long-lived than estimated based on the seesaw formula : the longevity problem. We show that, if one or more of the B-L Higgs fields are charged under a flavor symmetry (or discrete R symmetry), the split mass spectrum for the right-handed neutrinos as well as the required longevity is naturally realized. We provide several examples in which the predicted the X-ray flux is just below the current bound.

Axions seem ubiquitous in string theories and some of them may be produced non-thermally by heavy scalar decays, contributing to dark radiation. We study various cosmological effects of photons produced from the axionic dark radiation through axion-photon conversion in the presence of primordial magnetic fields, and derive tight constraints on the combination of the axion-photon coupling and the primordial magnetic field.

We revisit a single-field new inflation model based on a discrete R symmetry. Interestingly, the inflaton dynamics naturally leads to a heavy gravitino of mass m_3/2 = O(1-100) TeV, which is consistent with the standard-model like Higgs boson of mass mh = 126 GeV. However, the predicted spectral index n_s = 0.94 is in tension with the Planck result, n_s = 0.9603+-0.073. We show that the spectral index can be increased by allowing a small constant term in the superpotential during inflation. The required size of the constant is close to the largest allowed value for successful inflation, and it may be a result of a pressure toward larger values in the landscape. Alternatively, such constant term may arise in association with supersymmetry breaking required to cancel the negative cosmological constant from the inflaton sector.

We consider a simple class of models where dark radiation has self-interactions and therefore does not free stream. Such dark radiation has no anisotropic stress (or viscosity), leaving a distinct signature on the CMB angular power spectrum. Specifically we study a possibility that hidden gauge bosons and/or chiral fermions account for the excess of the effective number of neutrino species. They have gauge interactions and remain light due to the unbroken hidden gauge symmetry, leading to \Delta N_{\rm eff} \simeq 0.29 in some case.

We present a general polynomial chaotic inflation model in supergravity, for which the predicted spectral index and tensor-to-scalar ratio can lie within the 1 sigma region allowed by the Planck results. Most importantly, the predicted tensor-to-scalar ratio is large enough to be probed in the on-going and future B-mode experiments. We study the inflaton dynamics and the subsequent reheating process in a couple of specific examples. The non-thermal gravitino production from the inflaton decay can be suppressed in a case with a discrete Z_2 symmetry. We find that the reheating temperature can be naturally as high as O(10^{9-10}) GeV, sufficient for baryon asymmetry generation through (non-)thermal leptogenesis.

We propose a chaotic inflation model in supergravity based on polynomial interactions of the inflaton. Specifically we study the chaotic inflation model with quadratic, cubic and quartic couplings in the scalar potential and show that the predicted scalar spectral index and tensor-to-scalar ratio can lie within the 1 sigma region allowed by the Planck results.

We point out that the cosmological moduli problem is not necessarily resolved even if the modulus mass is heavier than O(10)TeV, contrary to the common wisdom. The point is that, in many scenarios where the lightest moduli fields are stabilized by supersymmetry breaking effects, those moduli fields tend to mainly decay into almost massless axions, whose abundance is tightly constrained by the recent Planck results. We study the moduli-induced axion problem in concrete examples, and discuss possible solutions. The problem and its solutions are widely applicable to decays of heavy scalar fields which dominate the energy density of the Universe, for instance, the reheating of the inflaton.

We evaluate the effective mass of a scalar field phi coupled to thermal plasma through Planck-suppressed interactions. We find it useful to rescale the coupled fields so that all the phi-dependences are absorbed into the yukawa and gauge couplings, which allows us to read off the leading order contributions to the effective mass \tilde m_{phi} from the 2-loop free energy calculated with the rescaled couplings. We give an analytical expression for \tilde m_{phi} at a sufficiently high temperature in the case where phi is coupled to the MSSM chiral superfields through non-minimal Kahler potential. We find that \tilde m_{phi}^2 is about 10^{-3} H^2 \sim 10^{-2} H^2 for typical parameter sets, where H is the Hubble expansion rate in the radiation-dominated era.

We study non-Gaussianity of density perturbations generated by an axionic curvaton, focusing on the case that the curvaton sits near the hilltop of the potential during inflation. Such hilltop curvatons can generate a red-tilted density perturbation spectrum without invoking large-field inflation. We show that, even when the curvaton dominates the Universe, the non-Gaussianity parameter fNL is positive and mildly increases towards the hilltop of the curvaton potential, and that fNL = O(10) is a general and robust prediction of such hilltop axionic curvatons. In particular, we find that the non-Gaussianity parameter is bounded as fNL <~ 30 - 40 for a range of the scalar spectral index, ns = 0.94 - 0.99, and that fNL = 20 - 40 is realized for the curvaton mass m_\sigma = 10 - 10^6 GeV and the decay constant f = 10^{12} - 10^{17} GeV. One of the plausible candidates for the axionic curvaton is an imaginary component of a modulus field with mass of order 10 - 100 TeV and decay constant of 10^{16} - 10^{17} GeV. We also discuss extreme cases where the curvaton drives a second inflation and find that fNL is typically smaller compared to non-inflating cases.

The baryon-dark matter coincidence is a long-standing issue. Interestingly, the recent observations suggest the presence of dark radiation, which, if confirmed, would pose another coincidence problem of why the density of dark radiation is comparable to that of photons. These striking coincidences may be traced back to the dark sector with particle contents and interactions that are quite similar, if not identical, to the standard model: a dark parallel world. It naturally solves the coincidence problems of dark matter and dark radiation, and predicts a sterile neutrino(s) with mass of ${\cal O}(0.1 - 1)$\,eV, as well as self-interacting dark matter made of the counterpart of ordinary baryons. We find a robust prediction for the relation between the abundance of dark radiation and the sterile neutrino, which can serve as the smoking-gun evidence of the dark parallel world.

We show that a sizable running spectral index suggested by the recent SPT data can be explained in the axionic curvaton model with a potential that consists of two sinusoidal contributions of different height and period. We find that the running spectral index is generically given by d ns/dlnk ~ (2pi/dN)(n_s - 1), where dN is the e-folds during one period of modulations. In the string axiverse, axions naturally acquire a mass from multiple contributions, and one of the axions may be responsible for the density perturbations with a sizable running spectral index via the curvaton mechanism. We note that the axionic curvaton model with modulations can also accommodate the red-tilted spectrum with a negligible running, without relying on large-field inflation.

We argue that dark radiation is naturally generated from the decay of the overall volume modulus in the LARGE volume scenario. We consider both sequestered and non-sequestered cases, and find that the axionic superpartner of the modulus is produced by the modulus decay and it can account for the dark radiation suggested by observations, while the modulus decay through the Giudice-Masiero term gives the dominant contribution to the total decay rate. In the sequestered case, the lightest supersymmetric particles produced by the modulus decay can naturally account for the observed dark matter density. In the non-sequestered case, on the other hand, the supersymmetric particles are not produced by the modulus decay, since the soft masses are of order the heavy gravitino mass. The QCD axion will then be a plausible dark matter candidate.

It is known that gravitinos are non-thermally produced in inflaton decay processes, which excludes many inflation models for a wide range of the gravitino mass. We find that the constraints from the gravitino overproduction can be greatly relaxed if the supersymmetry breaking field is much lighter than the inflaton, and if the dynamical scale of the supersymmetry breaking is higher than the inflaton mass. In particular, we show that many inflation models then become consistent with the pure gravity mediation with O(100)TeV gravitino which naturally explains the recently observed Higgs boson mass of about 125GeV.

In hybrid inflation, the inflaton generically has a tadpole due to gravitational effects in supergravity, which significantly changes the inflaton dynamics in high-scale supersymmetry. We point out that the tadpole can be cancelled if there is a supersymmetry breaking singlet with gravitational couplings, and in particular, the cancellation is automatic in no-scale supergravity. We consider the LARGE volume scenario as a concrete example and discuss the compatibility between the hybrid inflation and the moduli stabilization. We also point out that the dark radiation generated by the overall volume modulus decay naturally relaxes a tension between the observed spectral index and the prediction of the hybrid inflation.

We consider the cosmological moduli problem in the context of high-scale supersymmetry breaking suggested by the recent discovery of the standard-model like Higgs boson. In order to solve the notorious moduli-induced gravitino problem, we focus on the LARGE volume scenario, in which the modulus decay into gravitinos can be kinematically forbidden. We then consider the Affleck-Dine mechanism with or without an enhanced coupling with the inflaton, taking account of possible Q-ball formation. We show that the baryon asymmetry of the present Universe can be generated by the Affleck-Dine mechanism in LARGE volume scenario, solving the moduli and gravitino problems. We also find that the overall volume modulus decays into a pair of the axionic superpartners, which contribute to the extra relativistic degrees of freedom.

Gauge mediation provides us with a complete picture of supersymmetry breaking and its mediation within the effective field theories, and thus allows us to discuss consistencies with low-energy particle physics as well as cosmological observations. We study in detail the cosmological evolution of the pseudo-modulus field in the supersymmetry breaking sector and also the production of the gravitinos in the early Universe in a simple (but a complete) model of gauge mediation. Under fairly reasonable assumptions, it is found that there exists a parameter region where dark matter of the Universe is explained by both thermally and non-thermally produced gravitinos, while the baryon asymmetry of the Universe is generated through the thermal leptogenesis.

We show that the cosmological moduli problem is solved, without relying on huge late-time entropy production, if the universal cutoff scale of the theory is a few orders of magnitude smaller than the Planck scale. We obtain a general estimate of the modulus abundance in terms of the inflationary scale and the reheating temperature, and find in particular that the reheating temperature can be high enough for the non-thermal leptogenesis to work.

We propose a new inflation model in which a gauge singlet inflaton turns into the Higgs condensate after inflation. The inflationary path is characterized by a moduli space of supersymmetric vacua spanned by the inflaton and Higgs field. The inflation energy scale is related to the soft supersymmetry breaking, and the Hubble parameter during inflation is smaller than the gravitino mass. The initial condition for the successful inflation is naturally realized by the pre-inflation in which the Higgs plays a role of the waterfall field.

We show that heavy supersymmetric particles around O(100) TeV to O(1) PeV naturally appear in new inflation in which the Higgs boson responsible for the breaking of U(1)B-L plays the role of inflaton. Most important, the supersymmetric breaking scale is bounded above by the inflationary dynamics, in order to suppress the Coleman-Weinberg potential which would otherwise spoil the slow-roll inflation. Our scenario has rich phenomenological and cosmological implications: the Higgs boson mass at around 125 GeV can be easily explained, non-thermal leptogenesis works automatically, the gravitino production from inflaton decay is suppressed, the dark matter is either the lightest neutralino or the QCD axion, and the upper bound on the inflation scale for the modulus stabilization can be marginally satisfied.

We show that the modulated reheating mechanism can naturally account for dark radiation, whose existence is hinted by recent observations of the cosmic microwave background radiation and the primordial Helium abundance. In this mechanism, the inflaton decay rate depends on a light modulus which acquires almost scale-invariant quantum fluctuations during inflation. We find that the light modulus is generically produced by the inflaton decay and therefore a prime candidate for the dark radiation. Interestingly, an almost scale-invariant power spectrum predicted in the modulated reheating mechanism gives a better fit to the observation in the presence of the extra radiation. We discuss the production mechanism of the light modulus in detail taking account of its associated isocurvature fluctuations. We also consider a case where the modulus becomes the dominant component of dark matter.

We show that the standard-model Higgs boson mass m_h is correlated with the spectral index of density perturbation n_s in the inflation scenario with the inflaton being identified with the B-L Higgs boson. The Higgs boson mass ranges from m_h ~ 120GeV to 140GeV for n_s ~ 0.95 - 0.96. In particular, as n_s approaches to 0.96, the Higgs mass is predicted to be in the range of 125GeV to 140GeV in the case of relatively light gauginos, and 120GeV to 135GeV in the case where all SUSY particle masses are of the same order. This will be tested soon by the LHC experiment and the Planck satellite. The relation is due to the PeV-scale supersymmetry required by the inflationary dynamics. We also comment on the cosmological implications of our scenario such as non-thermal leptogenesis and dark matter.

We investigate density perturbations sourced by a curvaton with a generic energy potential. The key feature of a curvaton potential which deviates from a quadratic is that the curvaton experiences a non-uniform onset of its oscillation. This sources additional contributions to the resulting density perturbations, and we especially find that the non-Gaussianity parameter f_NL can become large with either sign no matter whether the curvaton dominates or subdominates the universe when it decays. Such non-quadratic curvaton potentials are required in order to produce a red-tilted density perturbation spectrum (without invoking large-field inflation), and are also motivated by explicit curvaton models based on microscopic physics. We further apply our generic results to the case where the curvaton is a pseudo-Nambu-Goldstone (NG) boson of a broken U(1) symmetry, and show that the resulting density perturbations are strongly enhanced towards the hilltop region of the energy potential, accompanied by a mild increase of the non-Gaussianity. Such hilltop NG curvatons can produce observationally suggested density perturbations under wide ranges of inflation/reheating scales, and further predict the non-Gaussianity of the density perturbations to lie within the range 10 <~ f_NL <~ 30.

We investigate an inflation model with the inflaton being identified with a Higgs boson responsible for the breaking of U(1)B-L symmetry. We show that supersymmetry must remain a good symmetry at scales one order of magnitude below the inflation scale, in order for the inflation model to solve the horizon and flatness problems, as well as to account for the observed density perturbation. The upper bound on the soft supersymmetry breaking mass lies between 1TeV and 10^3TeV. Interestingly, our finding opens up a possibility that universes with the low-scale supersymmetry are realized by the inflationary selection. Our inflation model has rich implications; non-thermal leptogenesis naturally works, and the gravitino and moduli problems as well as the moduli destabilization problem can be solved or ameliorated; the standard-model higgs boson receives a sizable radiative correction if the supersymmertry breaking takes a value on the high side ~10^3TeV.

We propose that the stability of dark matter is ensured by a discrete subgroup of the U(1)B-L gauge symmetry, Z_2(B-L). We introduce a set of chiral fermions charged under the U(1)B-L in addition to the right-handed neutrinos, and require the anomaly-cancellation conditions associated with the U(1)B-L gauge symmetry. We find that the possible number of fermions and their charges are tightly constrained, and that non-trivial solutions appear when at least five additional chiral fermions are introduced. The Fermat theorem in the number theory plays an important role in this argument. Focusing on one of the solutions, we show that there is indeed a good candidate for dark matter, whose stability is guaranteed by Z_2(B-L).

We argue that strong dynamics at the Planck scale can solve the cosmological moduli problem. We discuss its implications for inflation models, and find that a certain type of multi-field inflation model is required for this mechanism to work, since otherwise it would lead to the serious eta-problem. Combined with the inflaton-induced gravitino problem, we show that a chaotic inflation with a discrete symmetry naturally avoids both problems. Interestingly, the focus point supersymmetry is predicted when this mechanism is applied to the Polonyi model.

We argue that a large negative running spectral index, if confirmed, might suggest that there are abundant structures in the inflaton potential, which result in a fairly large (both positive and negative) running of the spectral index at all scales. It is shown that the center value of the running spectral index suggested by the recent CMB data can be easily explained by an inflaton potential with superimposed periodic oscillations. In contrast to cases with constant running, the perturbation spectrum is enhanced at small scales, due to the repeated modulations. We mention that such features at small scales may be seen by 21 cm observations in the future.

Recent cosmological observations, such as the measurement of the primordial 4He abundance, CMB, and large scale structure, give preference to the existence of extra radiation component, Delta N_nu > 0. The extra radiation may be accounted for by particles which were in thermal equilibrium and decoupled before the big bang nucleosynthesis. Broadly speaking, there are two possibilities: 1) there are about 10 particles which have very weak couplings to the standard model particles and decoupled much before the QCD phase transition; 2) there is one or a few light particles with a reasonably strong coupling to the plasma and it decouples after the QCD phase transition. Focusing on the latter case, we find that a light chiral fermion is a suitable candidate, which evades astrophysical constraints. Interestingly, such a scenario may be confirmed at the LHC. As a concrete example, we show that such a light fermion naturally appears in the E_6-inspired GUT.

We construct a chaotic inflation model in which the Higgs fields play the role of the inflaton in the standard model as well as in the singlet extension of the supersymmetric standard model. The key idea is to allow a non-canonical kinetic term for the Higgs field. The model is a realization of the recently proposed running kinetic inflation, in which the coefficient of the kinetic term grows as the inflaton field. The inflaton potential depends on the structure of the Higgs kinetic term. For instance, the inflaton potential is proportional to phi^2 and phi^{2/3} in the standard model and NMSSM, respectively. It is also possible to have a flatter inflaton potential.

We study various inflation models in the Jordan frame supergravity with a logarithmic Kahler potential. We find that, in a class of inflation models containing an additional singlet in the superpotential, three types of inflation can be realized: the Higgs-type inflation, power-law inflation, and chaotic inflation with/without a running kinetic term. The former two are possible if the holomorphic function dominates over the non-holomorphic one in the frame function, while the chaotic inflation occurs when both are comparable. Interestingly, the fractional-power potential can be realized by the running kinetic term. We also discuss the implication for the Higgs inflation in supergravity.

The seesaw mechanism in models with extra dimensions is shown to be generically consistent with a broad range of Majorana masses. The resulting democracy of scales implies that the seesaw mechanism can naturally explain the smallness of neutrino masses for an arbitrarily small right-handed neutrino mass. If the scales of the seesaw parameters are split, with two right-handed neutrinos at a high scale and one at a keV scale, one can explain the matter-antimatter asymmetry of the universe, as well as dark matter. The dark matter candidate, a sterile right-handed neutrino with mass of several keV, can account for the observed pulsar velocities and for the recent data from Chandra X-ray Observatory, which suggest the existence of a 5 keV sterile right-handed neutrino.

We study a recently proposed running kinetic inflation model in which the inflaton potential becomes flat due to rapid growth of the kinetic term at large inflaton field values. As concrete examples, we build a variety of chaotic inflation models in supergravity with e.g. quadratic, linear, and fractional-power potentials. The power of the potential generically increases after inflation, and the inflaton is often massless at the potential minimum in the supersymmetric limit, which leads to many interesting phenomena. First, the light inflaton mass greatly relaxes severe thermal and non-thermal gravitino problems. Secondly, the kination epoch is naturally present after inflation, which may enhance the gravity waves. Thirdly, since the inflaton is light, it is likely coupled to the Higgs sector for successful reheating. The inflaton and its superpartner, inflatino, may be produced at the LHC. Interestingly, the inflatino can be dark matter, if it is the lightest supersymmetric particle.

We propose a class of inflation models in which the coefficient of the inflaton kinetic term rapidly changes with energy scale. This may occur especially if the inflaton moves over a long distance during inflation as in the case of large-scale inflation. The peculiar behavior of the kinetic term opens up a new way to construct an inflation model. As a concrete example we construct a linear inflation model in supergravity. It is straightforward to build a chaotic inflation model with a fractional power along the same line. Interestingly, the potential takes a different form after inflation because of the running kinetic term.

We study Q-ball formation in the expanding universe on 1D, 2D and 3D lattice simulations. We obtain detailed Q-ball charge distributions, and find that the distribution is peaked at Q^{3D}_{peak} \simeq 1.9\times 10^{-2} (|\Phi_{in}|/m)^2, which is greater than the existing result by about 60%. Based on the numerical simulations, we discuss how the Q-ball formation proceeds. Also we make a comment on possible deviation of the charge distributions from what was conjectured in the past.

We argue that a primordial black hole is a natural and unique candidate for all dark matter. We show that, in a smooth-hybrid new double inflation model, a right amount of the primordial black holes, with a sharply-defined mass, can be produced at the end of the smooth-hybrid regime, through preheating. We first consider masses < 10^(-7)M_sun which are allowed by all the previous constraints. We next discuss much heavier mass 10^5 M_sun hinted at by entropy, and galactic size evolution, arguments. Effects on the running of the scalar spectral index are computed.

While the excess in cosmic-ray electrons and positrons reported by PAMELA and Fermi may be explained by dark matter decaying primarily into charged leptons, this does not necessarily mean that dark matter should not have any hadronic decay modes. In order to quantify the allowed hadronic activities, we derive constraints on the decay rates of dark matter into WW, ZZ, hh, qqbar and gg using the Fermi and HESS gamma-ray data. We also derive gamma-ray constraints on the leptonic e+e-, mu+mu- and tau+tau- final states. We find that dark matter must decay primarily into mu+mu- or tau+tau- in order to simultaneously explain the reported excess and meet all gamma-ray constraints.

Non-Gaussianity may exist in the CDM isocurvature perturbation. We provide general expressions for the bispectrum and trispectrum of both adiabatic and isocurvature pertubations. We apply our result to the QCD axion case, and found a consistency relation between the coefficients of the bispectrum and trispectrum : tau_{NL}^(iso)~10^3 [f_{NL}^(iso)]^{4/3}, if the axion is dominantly produced by quantum fluctuation. Thus future observations of the trispectrum, as well as the bispectrum, will be important for understanding the origin of the CDM and baryon asymmetry.

We study cosmological implications of an R-axion, a pseudo Nambu-Goldstone boson associated with spontaneous breaking of an U(1)_R symmetry, focusing on its quantum fluctuations generated during inflation. We show that, in the anomaly mediation, the R-axion decays into a pair of gravitinos, which eventually decay into the visible particles producing the neutralino LSP. As a result, the quantum fluctuations of the R-axion are inherited by the cold dark matter isocurvature density perturbation with potentially large non-Gaussianity. The constraints on the inflation scale and the initial misalignment are derived.

We show that the PAMELA anomaly in the positron fraction as well as the ATIC/PPB-BETS excesses in the e^- + e^+ flux are simultaneously explained in our scenario that a hidden U(1)H gauge boson constitutes dark matter of the Universe and decays into the standard-model particles through a kinetic mixing with an U(1)B-L gauge boson. Interestingly, the B-L charge assignment suppresses an antiproton flux in consistent with the PAMELA and BESS experiments, while the hierarchy between the B-L symmetry breaking scale and the weak scale naturally leads to the right lifetime of O(10^26) seconds.

We provide a scenario in which a hidden U(1) gauge boson constitutes dark matter of the Universe and decays into the standard-model particles through a kinetic mixing with an $U(1)_{B-L}$ gauge boson. Interestingly, our model can naturally account for the steep rise in the positron fraction recently reported by PAMELA. Moreover, we find that due to the charge assignment of $U(1)_{B-L}$, only a small amount of antiprotons are produced in the decay, which is also consistent with the PAMELA and other observational data.

We show that, assuming the dispersion relation proposed recently by Horava in the context of quantum gravity, radiation energy density exhibits a peculiar dependence on the scale factor; the radiation energy density decreases proportional to a^{-6}. This simple scaling can have an impact on cosmology. As an example, we show that the resultant baryon asymmetry as well as the stochastic gravity waves can be enhanced. We also discuss current observational constraint on the dispersion relation.

We show that R-parity violating decay of Wino dark matter of mass \sim 3 TeV can naturally account for the flux and spectral shape of the cosmic-ray electrons and positrons observed by the PAMELA and Fermi satellites. To provide a theoretical basis for the scenario, we also present a model that trilinear R-parity breaking appears with a coefficient suppressed by powers of the gravitino mass, which naturally leads to the Wino lifetime of O(10^26) sec.

A new and powerful probe of the origin and evolution of structures in the Universe has emerged and been actively developed over the last decade. In the coming decade, non-Gaussianity, i.e., the study of non-Gaussian contributions to the correlations of cosmological fluctuations, will become an important probe of both the early and the late Universe. Specifically, it will play a leading role in furthering our understanding of two fundamental aspects of cosmology and astrophysics: (i) the physics of the very early universe that created the primordial seeds for large-scale structures, and (ii) the subsequent growth of structures via gravitational instability and gas physics at later times. To date, observations of fluctuations in the Cosmic Microwave Background (CMB) and the Large-Scale Structure of the Universe (LSS) have focused largely on the Gaussian contribution as measured by the two-point correlations (or the power spectrum) of density fluctuations. However, an even greater amount of information is contained in non-Gaussianity and a large discovery space therefore still remains to be explored. Many observational probes can be used to measure non-Gaussianity, including CMB, LSS, gravitational lensing, Lyman-alpha forest, 21-cm fluctuations, and the abundance of rare objects such as clusters of galaxies and high-redshift galaxies. Not only does the study of non-Gaussianity maximize the science return from a plethora of present and future cosmological experiments and observations, but it also carries great potential for important discoveries in the coming decade.

We point out that the long-standing strong CP problem may be resolved by an anthropic argument. The key ideas are: (i) to allow explicit breaking(s) of the Peccei-Quinn symmetry which reduces the strong CP problem to the cosmological constant problem, and (ii) to conjecture that the probability distribution of the vacuum energy has a mild pressure towards higher values. The cosmological problems of the (s)axion with a large Peccei-Quinn scale are absent in our mechanism, since the axion acquires a large mass from the explicit breaking.

We study a scenario that a U(1) hidden gaugino constitutes the dark matter in the Universe and decays into a lepton and slepton pair through a mixing with a U(1)B-L gaugino. We find that the dark-matter decay can account for the recent PAMELA and ATIC anomalies in the cosmic-ray positrons and electrons without an overproduction of antiprotons.

We show that, if decaying gravitino dark matter is responsible for the PAMELA and ATIC/PPB-BETS anomalies in the cosmic-ray electron and positron fluxes, both a reheating temperature and a gluino mass are constrained from above. In particular, the gluino mass is likely within the reach of LHC, if the observed baryon asymmetry is explained by thermal leptogenesis scenario.

Light scalars may be ubiquitous in nature, and their quantum fluctuations can produce large non-Gaussianity in the cosmic microwave background temperature anisotropy. The non-Gaussianity may be accompanied with a small admixture of isocurvature perturbations, which often have correlations with the curvature perturbations. We present a general method to calculate the non-Gaussianity in the adiabatic and isocurvature perturbations with and without correlations, and see how it works in several explicit examples. We also show that they leave distinct signatures on the bispectrum of the cosmic microwave background temperature fluctuations.

The hidden scalar field, which couples to the visible sector only through Planck-suppressed interactions, is a candidate for dark matter owing to its long lifetime. Decay of such a scalar field offers observational tests of this scenario. We show that decay of the hidden scalar field can explain the observed excess of high-energy positrons/electrons observed by ATIC/PPB-BETS, for a suitable choice of the mass and the vacuum expectation value of the field. We also show that the same choice of the parameters gives the observed dark matter abundance. Such a remarkable coincidence suggests that the Planck-suppressed interactions may be responsible for the observed excess in the cosmic-ray positrons/electrons.

We argue that we may be able to sort out dark matter models in which electrons are generated through the annihilation and/or decay of dark matter, by using a fact that the initial energy spectrum is reflected in the cosmic-ray electron flux observed at the Earth even after propagation through the galactic magnetic field. To illustrate our idea we focus on three representative initial spectra: (i)monochromatic (ii)flat and (iii)double-peak ones. We find that those three cases result in significantly different energy spectra, which may be probed by the Fermi satellite in operation or an up-coming cosmic-ray detector such as CALET.

If dark matter possesses a lepton number, it is natural to expect the dark-matter annihilation and/or decay mainly produces the standard model leptons, while negligible amount of the antiproton is produced. To illustrate such a simple idea, we consider a scenario that a right-handed sneutrino dark matter decays into the standard model particles through tiny R-parity violating interactions. Interestingly enough, charged leptons as well as neutrinos are directly produced, and they can lead to a sharp peak in the predicted positron fraction. Moreover, the decay of the right-handed sneutrino also generates diffuse continuum gamma rays which may account for the excess observed by EGRET, while the primary antiproton flux can be suppressed. Those predictions on the cosmic-ray fluxes of the positrons, gamma rays and antiprotons will be tested by the PAMELA and FGST observatories.

We study a scenario that a hidden gauge boson constitutes the dominant component of dark matter and decays into the standard model particles through a gauge kinetic mixing. Interestingly, gamma rays and positrons produced from the decay of hidden gauge boson can explain both the EGRET excess of diffuse gamma rays and the HEAT anomaly in the positron fraction. The spectra of the gamma rays and the positrons have distinctive features; the absence of line emission of the gamma ray and a sharp peak in the positron fraction. Such features may be observed by the GLAST and PAMELA satellites.

We study non-Gaussianity induced by a pseudo Nambu-Goldstone boson with a cosine-type scalar potential. We focus on how the non-Gaussianity is affected when the pseudo Nambu-Goldstone boson rolls down from near the top of the scalar potential where the deviation from a quadratic potential is large. We find that the resultant non-Gaussianity is similar to that obtained in the quadratic potential, if the pseudo Nambu-Goldstone boson accounts for the curvature perturbation; the non-Gaussianity is enhanced, otherwise.

We study a scenario that large non-Gaussianity arises from the baryon asymmetry of the Universe. There are baryogenesis scenarios containing a light scalar field, which may result in baryonic isocurvature perturbations with some amount of non-Gaussianity. As an explicit example we consider the Affleck-Dine mechanism and show that a flat direction of the supersymmeteric standard model can generate large non-Gaussianity in the curvature perturbations, satisfying the observational constraints on the baryonic isocurvature perturbations. The sign of a non-linearity parameter, f_{NL}, is negative, if the Affleck-Dine mechanism accounts for the observed baryon asymmetry; otherwise it can be either positive or negative.

We develop a formalism to study non-Gaussianity in both curvature and isocurvature perturbations. It is shown that non-Gaussianity in the isocurvature perturbation between dark matter and photons leaves distinct signatures in the CMB temperature fluctuations, which may be confirmed in future experiments, or possibly, even in the currently available observational data. As an explicit example, we consider the QCD axion and show that it can actually induce sizable non-Gaussianity for the inflationary scale, H_{inf} = O(10^9 - 10^{11})GeV.

We point out that a light scalar field fluctuating around a symmetry-enhaced point can generate large non-Gaussianity in density fluctuations. We name such a particle as an "ungaussiton", a scalar field dominantly produced by the quantum fluctuations,generating sizable non-Gaussianity in the density fluctuations. We derive a consistency relation between the bispectrum and the trispectrum, tau_NL = 10^3 f_NL^(4/3), which can be extended to arbitrary high order correlation functions. If such a relation is confirmed by future observations, it will strongly support this mechanism.

If gaugino condensations occur in the early universe, domain walls are produced as a result of the spontaneous breaking of a discrete R symmetry. Those domain walls eventually annihilate with one another, producing the gravitational waves. We show that the gravitational waves can be a probe for measuring the gravitino mass, if the constant term in the superpotential is the relevant source of the discrete R symmetry breaking.

The cosmic abundance of a long-lived charged particle such as a stau is tightly constrained by the catalyzed big bang nucleosynthesis. One of the ways to evade the constraints is to dilute those particles by a huge entropy production. We evaluate the dilution factor in a case that non-relativistic matter dominates the energy density of the universe and decays with large entropy production. We find that large Q balls can do the job, which is naturally produced in the gauge-mediated supersymmetry breaking scenario.

We discuss inflaton decay in supergravity, taking account of the gravitational effects. It is shown that, if the inflaton has a nonzero vacuum expectation value, it generically couples to any matter fields that appear in the superpotential at the tree level, and to any gauge sectors through anomalies in the supergravity. Through these processes, the inflaton generically decays into the supersymmetry breaking sector, producing many gravitinos. The inflaton also directly decays into a pair of the gravitinos. We derive constraints on both inflation models and supersymmetry breaking scenarios for avoiding overproduction of the gravitinos. Furthermore, the inflaton naturally decays into the visible sector via the top Yukawa coupling and SU(3)_C gauge interactions.

We discuss a scenario that gravitinos produced non-thermally by an inflaton decay constitute dark matter in the present universe. We find that this scenario is realized for wide ranges of the inflaton mass and the vacuum expectation value. What is intriguing about this scenario is that the gravitino dark matter can have a relatively large free streaming length at matter-radiation equality, which can be probed by future observation on QSO-galaxy strong lens system.

We present models of decaying particles for increasing the effective number of neutrinos N_nu after big bang nucleosynthesis but before the structure formation begins. We point out that our scenario not only solves the discrepancy between the constraints on N_nu from these two epochs, but also provides a possible answer to deeper inconsistency in the estimation of the matter power spectrum amplitude at small scales, represented by sigma_8, between the WMAP and some small scale matter power measurements such as the Lyman-alpha forest and weak lensing. We consider (a) saxion decay into two axions; (b) gravitino decay into axino and axion; (c) Dirac right-handed sneutrino decay into gravitino and right-handed neutrino.

We derive constraint on the effective number of neutrino species N_nu from the cosmic microwave background power spectrum of the WMAP and galaxy clustering power spectrum of the SDSS luminous red galaxies (LRGs). Using these two latest data sets of CMB and galaxy clustering alone, we obtain the limit 0.9 < N_nu < 8.2 (95% C.L.) for the power-law LambdaCDM flat universe, with no external prior. The lower limit corresponds to the lower bound on the reheating temperature of the universe T_R > 2 MeV.

We propose a retrofitted gravity mediation model which alleviates the gravitino overproduction from decays of an inflaton and a supersymmetry breaking field. In the model, we introduce an approximate U(1) symmetry under which the supersymmetry breaking field is charged, although it is broken by a mass term of messenger fields to generate gaugino masses of order the weak scale. In a low-scale inflation model, we find regions in which the gravitino overproduction problem is avoided.

We study the decay of the inflaton in no-scale supergravity and show that decay due to the gravitational interactions through supergravity effects is highly suppressed relative to the case in minimal supergravity or models with a generic Kahler potential. We also show that decay to gravitinos is suppressed. We demonstrate that decay and sufficient reheating are possible with the introduction of a non-trivial gauge kinetic term. This channel may be dominant in no-scale supergravity, yet yields a re-heating temperature which is low enough to avoid the gravitino problem while high enough for Big Bang Nucleosynthesis and baryogenesis.

We point out that the inflaton spontaneously decays into any gauge bosons and gauginos via the super-Weyl, Kahler and sigma-model anomalies in supergravity, once the inflaton acquires a non-vanishing vacuum expectation value. In particular, in the dynamical supersymmetry breaking scenarios, the inflaton necessarily decays into the supersymmetry breaking sector, if the inflaton mass is larger than the dynamical scale. This generically causes the overproduction of the gravitinos, which severely constrains the inflation models.

We argue that a non-thermal leptogenesis occurs spontaneously, without direct couplings of the inflaton with right-handed neutrinos, in a wide class of high-scale inflation models such as the chaotic and hybrid inflation. It is only a finite vacuum expectation value of the inflaton, or more precisely, a linear term in the Kahler potential, that is a prerequisite for the spontaneous non-thermal leptogenesis. To exemplify how it works, we show that a chaotic inflation model in supergravity naturally produces a right amount of baryon asymmetry via the spontaneous non-thermal leptogenesis. We also discuss the gravitino production from the inflaton.

We point out that supergravity effects enable the inflaton to decay into all matter fields, including the visible and the supersymmetry breaking sectors, once the inflaton acquires a non-vanishing vacuum expectation value. The new decay processes have great impacts on cosmology; the reheating temperature is bounded below; the gravitinos are produced by the inflaton decay in a broad class of the dynamical supersymmetry breaking models. We derive the bounds on the inflaton mass and the vacuum expectation value, which severely constrain high-scale inflations such as the hybrid and chaotic inflation models.

We show that the gravitino-overproduction problem is prevalent among inflation models in supergravity. An inflaton field \phi generically acquires (effective) non-vanishing auxiliary field G_\phi, if the K\"ahler potential is non-minimal. The inflaton field then decays into a pair of the gravitinos. We extensively study the cosmological constraints on G_\phi for a wide range of the gravitino mass. For many inflation models we explicitly estimate G_\phi, and show that the gravitino-overproduction problem severely constrains the inflation models, unless such an interaction as K = \kappa/2 |\phi|^2 z^2 + h.c. is suppressed (here z is the field responsible for the supersymmetry breaking). We find that many of them are already excluded or on the verge of, if \kappa \sim O(1).

A heavy scalar field such as moduli or an inflaton generally mixes with a field responsible for the supersymmetry breaking. We study the scalar decay into the standard model particles and their superpartners, gravitinos, and the supersymmetry breaking sector, particularly paying attention to decay modes that proceed via the mixing between the scalar and the supersymmetry breaking field. The impacts of the new decay processes on cosmological scenarios are also discussed; the modulus field generically produces too much gravitinos, and most of the inflation models tend to result in too high reheating temperature and/or gravitino overproduction.

We examine non-thermal dark matter production from a late-decaying scalar field, with a particular attention on non-renormalizable operators of D=5 through which the scalar field decays into the standard model particles and their superpartners. We show that almost the same number of superparticles as that of particles are generally produced from the decay. To avoid the gravitino overproduction problem, the decay is favored to proceed via interactions with an intermediate cut-off scale M << M_P. This should be contrasted to the conventional scenario using the modulus decay. The bosonic supersymmetry partner of the axion, i.e., saxion, is proposed as a natural candidate for such late-decaying scalar fields. We find that a right amount of the wino/higgsino dark matter with a mass of O(100) GeV is obtained for the saxion mass around the weak scale and axion decay constant, F_A = O(10^{9-12}) GeV.

We investigate the cosmological moduli problem by studying a modulus decay in detail and find that the branching ratio of the gravitino production is generically of O(0.01-1), which causes another cosmological disaster. Consequently, the cosmological moduli problem cannot be solved simply by making the modulus mass heavier than 100 TeV. We also illustrate our results by explicitly calculating the branching ratio into the gravitinos in the mixed modulus--anomaly/KKLT- and racetrack-type models.

Most of the inflation models end up with non-vanishing vacuum expectation values of the inflaton fields \phi in the true vacuum, which induce, in general, nonvanishing auxiliary field G_\phi for the inflaton potential in supergravity. We show that the presence of nonzero G_\phi gives rise to inflaton decay into a pair of the gravitinos and are thereby severely constrained by cosmology especially if the gravitino is unstable and its mass is in a range of O(100) GeV \sim O(10) TeV. For several inflation models, we explicitly calculate the values of G_\phi and find that most of them are excluded or on the verge of being excluded for the gravitino mass in that range. We conclude that an inflation model with vanishing G_\phi, typically realized in a chaotic inflation, is favored in a sense that it naturally avoids the potential gravitino overproduction problem.

We propose a scenario in which dark energy and dark matter are described in a unified manner. The ultralight pseudo-Nambu-Goldstone (pNG) boson, A, naturally explains the observed magnitude of dark energy, while the bosonic supersymmetry partner of the pNG boson, B, can be a dominant component of dark matter. The decay of B into a pair of electron and positron may explain the 511 keV gamma ray from the Galactic Center.

The 511 keV photons from the galactic center can be explained by positrons produced through Q-ball decay. In the scheme of gauge-mediated supersymmetry breaking, large Q balls with lepton charge are necessarily long-lived. In particular, the lifetime can be as long as (or even longer than) the age of the universe. If kinematically allowed, such large Q balls decay into positrons, which eventually annihilate with electrons into 511 keV photons. Our scenario is realized within the minimal supersymmetric standard model in the inflationary universe, which is very plausible.

We study how the oscillations of the neutrinos affect their thermalization process during the reheating period with temperature O(1) MeV in the early universe. We follow the evolution of the neutrino density matrices and investigate how the predictions of big bang nucleosynthesis vary with the reheating temperature. For the reheating temperature of several MeV, we find that including the oscillations makes different predictions, especially for $^4$He abundance. Also, the effects on the lower bound of the reheating temperature from cosmological observations are discussed.

Q-ball is a non-topological soliton whose stability is ensured by global U(1) symmetry. We study a Q-ball which arises in the Affleck-Dine mechanism for baryogenesis and consider its possible instability due to U(1) breaking term ($A$-term) indispensable for successful baryogenesis. It is found that the instability destroys the Q-ball if its growth rate exceeds inverse of the typical relaxation time scale of the Q-ball. However, the instability is not so strong as it obstructs the cosmological formation of the Q-balls.

It is shown that late-time decay of domain walls can dilute unwanted relics such as moduli, if the universe was dominated by frustrated domain walls with tension $\sigma = (1 - 100 TeV)^3$. Since energy density of the frustrated domain walls decreases as slow as the inverse of the scale factor, an overclosure limit on the axion decay constant $f_a$ is also considerably relaxed. In fact $f_a$ can be as large as the Planck scale, which may enable us to naturally implement the QCD axion in the string scheme. Furthermore, in contrast to thermal inflation models, the Affleck-Dine baryogenesis can generate enough asymmetry to explain the present baryon abundance, even in the presence of late-time entropy production.

We discuss cosmic microwave background anisotropies in models with an unconventional primordial power spectrum. In particular, we consider an initial power spectrum with some ``spiky'' corrections. Interestingly, such a primordial power spectrum generates ``wavy'' structure in the CMB angular power spectrum.

Within the standard big bang nucleosynthesis (BBN) and cosmic microwave background (CMB) framework, the baryon density measured by the Wilkinson Microwave Anisotropy Probe (WMAP) or the primordial D abundance is much higher than the one measured by the $^4$He or $^7$Li abundances. To solve the discrepancy, we propose a scenario in which additional baryons appear after BBN. We show that simply adding the baryons can not be a solution but the existence of a large lepton asymmetry before BBN makes the scenario successful. These extra baryons and leptons, in addition to the initial baryons which exist before the BBN, can be all produced from $Q$-balls.

We study dynamics of flat directions in the minimal supersymmetric standard model taking account of its constituent fields. It is found that there exist new instabilities due to the D-term potential and the nature of these instabilities depends on the eccentricity of the orbit. For a roughly circular orbit, it is similar to narrow-band resonance described by the Mathieu equation. For an elliptic orbit, the instabilities exhibit tachyonic nature. In the narrow-band regime, we show that the existence of the instabilities is crucial to the formation of $Q$-balls. The $Q$-ball formation proceeds through quasi-stable state called ``$I$-ball''. The transition from $I$-balls to $Q$-balls is made efficient by the $D$-term instability and during this process some fraction of the charge of the $Q$-ball is emitted. This discovery may revive the scenario which relates the baryon number and dark matter density of the universe. Furthermore, the tachyonic $D$-term instability causes the drastic decay of the flat direction moving in an orbit with relatively large eccentricity. Thus the evolution of flat directions is completely altered by the appearance of this instability.

We consider the possibility of using one of the $D$-flat directions in the minimal supersymmetric standard model (MSSM) as the inflaton. We show that the flat direction consisting of (first generation) left- and right-handed up-squarks as well as the up-type Higgs boson may play the role of the inflaton if dominant part of the up-quark mass is radiatively generated from supersymmetric loop diagrams. We also point out that, if the R-parity violating Yukawa coupling is of $O(10^{-7})$, $R$-odd $D$-flat directions may be another possible candidate of the inflaton. Such inflation models using $D$-flat directions in the MSSM are not only testable with collider experiments but also advantageous to resolve the problem how the inflaton reheats the universe.

We discuss a spontaneous baryogenesis mechanism in flat directions. After identifying the Nambu-Goldstone mode which derivatively couples to the associated $U$(1) current and rotates due to the A-term, we show that spontaneous baryogenesis can be naturally realized in the context of the flat direction. As applications, we discuss two scenarios of baryogenesis in detail. One is baryogenesis in a flat direction with a vanishing $B-L$ charge, especially, with neither baryon nor lepton charge, which was recently proposed by Chiba and the present authors. The other is a baryogenesis scenario compatible with a large lepton asymmetry.

It is shown that hypothetical neutrino-majoron coupling can suppress neutrino flavor oscillations in the early universe, in contrast to the usual weak interaction case. This reopens a window for a noticeable cosmological lepton asymmetry which is forbidden for the large mixing angle solution in the case of standard interactions of neutrinos.

Flat directions in the minimal supersymmetric standard model are known to deform into non-topological solitons, Q-balls, which generally possess both baryon and lepton asymmetries. We investigate how Q-balls evolve if some of the constituent fields of the flat direction decay into light species. It is found that the Q-balls takes a new configuration whose energy per charge slightly increases due to the decay. Specifically, we show that all the stable Q-balls eventually transform into pure B-balls via the decay into neutrinos.

We present a new mechanism of spontaneous baryogenesis. Usually such mechanisms require a derivative coupling between a scalar field and baryon current. In our model, the scalar field instead derivatively couples to a current associated with some global symmetry $U(1)_Q$. In this case, baryogenesis is still possible provided that an interaction exists, which violates the baryon and $U(1)_Q$ symmetries simultaneously. As a concrete example, we discuss baryogenesis in a flat direction with neither baryon nor lepton charge.

We study the curvaton scenario in supersymmetric framework paying particular attention to the fact that scalar fields are inevitably complex in supersymmetric theories. If there are more than one scalar fields associated with the curvaton mechanism, isocurvature (entropy) fluctuations between those fields in general arise, which may significantly affect the properties of the cosmic density fluctuations. We examine several candidates for the curvaton in the supersymmetric framework, such as moduli fields, Affleck-Dine field, $F$- and $D$-flat directions, and right-handed sneutrino. We estimate how the isocurvature fluctuations generated in each case affect the cosmic microwave background angular power spectrum. With the use of the recent observational result of the WMAP, stringent constraints on the models are derived and, in particular, it is seen that large fraction of the parameter space is excluded if the Affleck-Dine field plays the role of the curvaton field. Natural and well-motivated candidates of the curvaton are also listed.

The recent observation of the cosmic microwave background anisotropy by the WMAP confirmed that the lower multipoles are considerably suppressed. From the standpoint of the cosmic variance, it is nothing but a statistical accident. Alternatively, one can attribute the deficit of fluctuation on the large scale to the cosmic history, which might be explained in the context of the inflationary physics. In this paper, we show that it is possible to explain the suppressed lower multipoles in the hybrid new inflation model.

We investigate whether the Affleck-Dine (AD) mechanism works when the contribution of the two-loop thermal correction to the potential is negative in the gauge-mediated supersymmetry breaking models. The AD field is trapped far away from the origin by the negative thermal correction for a long time until the temperature of the universe becomes low enough. The most striking feature is that the Hubble parameter becomes much smaller than the mass scale of the radial component of the AD field, during the trap. Then, the amplitude of the AD field decreases so slowly that the baryon number is not fixed even after the onset of radial oscillation. The resultant baryon asymmetry crucially depends on whether the Hubble parameter, $H$, is larger than the mass scale of the phase component of the AD field, $M_\theta$, at the beginning of oscillation. If $H < M_\theta$ holds, the formation of Q balls plays an essential role to determine the baryon number, which is found to be washed out due to the nonlinear dynamics of Q-ball formation. On the other hand, if $H > M_\theta$ holds, it is found that the dynamics of Q-ball formation does not affect the baryon asymmetry, and that it is possible to generate the right amount of the baryon asymmetry.

We find that there exists a soliton-like solution ``I-ball'' in theories of a real scalar field if the scalar potential satisfies appropriate conditions. Although the I-ball does not have any topological or global U(1) charges, its stability is ensured by the adiabatic invariance for the oscillating field.

We propose a scenario which can explain large lepton asymmetry and small baryon asymmetry simultaneously. Large lepton asymmetry is generated through Affleck-Dine (AD) mechanism and almost all the produced lepton numbers are absorbed into Q-balls (L-balls). If the lifetime of the L-balls is longer than the onset of electroweak phase transition but shorter than the epoch of big bang nucleosynthesis (BBN), the large lepton asymmetry in the L-balls is protected from sphaleron effects. On the other hand, small (negative) lepton numbers are evaporated from the L-balls due to thermal effects, which are converted into the observed small baryon asymmetry by virtue of sphaleron effects. Large and positive lepton asymmetry of electron type is often requested from BBN. In our scenario, choosing an appropriate flat direction in the minimal supersymmetric standard model (MSSM), we can produce positive lepton asymmetry of electron type but totally negative lepton asymmetry.