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.

The QCD axion and an axion-like particle (ALP) are compelling candidates of dark matter. For the QCD axion, it is known that when the Peccei-Quinn (PQ) symmetry is spontaneously broken after inflation, the large initial fluctuation can lead to axion clump formation. On the other hand, when the symmetry is already broken during inflation, it has been believed that the axion clump formation does not occur due to the small amplitude of the initial axion fluctuation. We revisit this prevailing understanding, considering both the QCD axion and an ALP. We find that for the QCD axion, the clump formation does not occur even if we consider an extremely fine-tuned initial condition. Meanwhile, it turns out that for an ALP which allows a more general potential form, the clump formation can take place through the tachyonic instability or/and the resonance instability, considering a multiple cosine potential.

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.

It was recently shown that a coherent oscillation of an axion can cause an efficient parametric resonance, leading to a prominent emission of the gravitational waves (GWs). In this paper, conducting the Floquet analysis, we investigate the parametric resonance instability, which potentially triggers the emission of the GWs from axions. Such a resonance instability takes place, when the time evolution of the background field significantly deviates from the harmonic oscillation. Therefore, the resonance instability cannot be described by the Mathieu equation, whose stability/instability chart is well known. In this paper, introducing an explicitly calculable parameter $\tilde{q}$, which can be used to classify different types of the parametric resonance described by the general Hill's equation, we investigate the stability/instability chart for the general Hill's equation. This can also apply to the case where the background oscillation is anharmonic. We show that the flapping resonance instability, which takes place for $\tilde{q}=O(1)$, typically leads to the most significant growth of the inhomogeneous modes among the self-resonance instability. We also investigate whether the flapping resonance takes place for the cosine potential or not.

We show that the number of primordial black holes (PBHs) which is originated from primordial density perturbations with moderately-tilted power spectrum fluctuates following the log-normal distribution, while it follows the Poisson distribution if the spectrum is steeply blue. The log-normal, as well as the Poisson, fluctuation of the PBH number behaves as an isocurvature mode and affects the matter power spectrum and the halo mass function in a different way from those for the Poisson case. The future 21cm observation can potentially put a stronger constraint on the PBH fraction than the current one in a wide mass range, $10^{-5}M_\odot$--$10M_\odot$.

Axions predicted in string theory may have a scalar potential which has a much shallower potential region than the conventional cosine potential. We first show that axions which were located at such shallow potential regions generically undergo prominent resonance instabilities: the well-known narrow resonance and/or the flapping resonance, which has not been well investigated. We also study non-linear dynamics of axions caused by these resonance instabilities based on lattice simulation. We find that string axions in various mass ranges generate gravitational waves (GWs) with peaks at various frequencies determined by the mass scales, dubbed the GW forest. This may allow us to explore string axiverse through future multi-frequency GW observations. We also investigate GWs produced by the axion which accounts for present dark matter component.

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 revisit the generation of dark matter isocurvature perturbations in the curvaton model in greater detail, both analytically and numerically. As concrete examples, we investigate the cases of thermally decoupled dark matter and axionic dark matter. We show that the radiation produced by the decay of the curvaton, which has not been taken into account in previous analytical studies, can significantly affect the amplitude of isocurvature perturbations. In particular, we find that they are drastically suppressed even when the dark matter freeze-out (or the onset of the axion oscillations for axionic dark matter) occurs before the curvaton decays, provided the freeze-out takes place deep in the curvaton-dominated Universe. As a consequence, we show that the current observational isocurvature constraints on the curvaton parameters are not as severe as usually thought.

We show that the stochastic evolution of an interacting system of the Higgs and a spectator scalar field naturally gives rise to an enhanced probability of settling down at the electroweak vacuum at the end of inflation. Subsequent destabilization due to parametric resonance between the Higgs and the spectator field can be avoided in a wide parameter range. We further argue that the spectator field can play the role of dark matter.

We discuss early structure formation of small scales sourced by primordial black holes (PBHs) which constitute a small part of present cold dark matter component. We calculate the mass function and power spectrum of haloes originated from the Poisson fluctuations of PBH number and show that the number of small haloes is significantly modified in the presence of PBHs even if their fraction accounts for only $10^{-4}$--$10^{-3}$ of total dark matter abundance. We then compute the subsequent 21cm signature from those haloes. We find that PBHs can provide major contributions at high redshifts within the detectability of future experiments such as Square Kilometer Array, and provide a forecast constraint on the PBH fraction.

We consider a novel cosmological scenario in which a curvaton is long-lived and plays the role of cold dark matter (CDM) in the presence of a short, secondary inflation. Non-trivial evolution of the large scale cosmological perturbation in the curvaton scenario can affect the duration of the short term inflation, resulting in the inhomogeneous end of inflation. Non-linear parameters of the curvature perturbation are predicted to be fNL ~ 5/4 and gNL ~ 0. The curvaton abundance can be well diluted by the short-term inflation and accordingly, it does not have to decay into the Standard Model particles. Then the curvaton can account for the present CDM with the isocurvature perturbation being sufficiently suppressed because both the adiabatic and CDM isocurvature perturbations have the same origin. As an explicit example, we consider the thermal inflation scenario and a string axion as a candidate for this curvaton-dark matter. We further discuss possibilities to identify the curvaton-dark matter with the QCD axion.

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

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.

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.

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 study the impact of a temperature-dependent curvaton decay rate on the primordial curvature perturbation generated in the curvaton scenario. Using the familiar sudden decay approximation, we obtain an analytical expression for the curvature perturbation after the decay of the curvaton. We then investigate numerically the evolution of the background and of the perturbations during the decay. We first show that the instantaneous transfer coefficient, related to the curvaton energy fraction at the decay, can be extended into a more general parameter, which depends on the net transfer of the curvaton energy into radiation energy or, equivalently, on the total entropy ratio after the complete curvaton decay. We then compute the curvature perturbation and compare this result with the sudden decay approximation prediction.

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 investigate the gravitational wave background induced by the first order scalar perturbations in the curvaton models. We consider the quadratic and axion-like curvaton potential which can generate the blue-tilted power spectrum of curvature perturbations on small scales and derive the maximal amount of gravitational wave background today. We find the power spectrum of the induced gravitational wave background has a characteristic peak at the frequency corresponding to the scale reentering the horizon at the curvaton decay, in the case where the curvaton does not dominate the energy density of the Universe. We also find the enhancement of the amount of the gravitational waves in the case where the curvaton dominates the energy density of the Universe. Such induced gravitational waves would be detectable by the future space-based gravitational wave detectors or pulsar timing observations.

It is known that supersymmetric hybrid inflation model may require severe tunings on the initial condition for large gravitino mass of order 100 - 1000 TeV due to the constant term in the superpotential. We propose a modified hybrid inflation model, where the constant term is suppressed during inflation and generated after inflation by replacing a constant term with dynamical field. In this modified model, successful inflation consistent with large gravitino mass takes place without severe tunings on the initial condition. Constraint from cosmic strings is also relaxed.

We argue that the existence of the cold dark matter is explained by primordial black holes.We show that a significant number of primordial black holes can be formed in an axion-like curvaton model, in which the highly blue-tilted power spectrum of primordial curvature perturbations is achieved.It is found that the produced black holes with masses $\sim 10^{20} -10^{38} \mathrm{g}$ account for the present cold dark matter.We also argue the possibility of forming the primordial black holes with mass $\sim 10^5 M_{\odot}$ as seeds of the supermassive black holes.

We show that the smooth hybrid inflation is naturally realized in a framework of supersymmetric axion model. Identifying the Peccei-Quinn scalar fields as a part of the infaton sector, successful inflation takes place reproducing the amplitude and spectral index of the curvature perturbation observed by WMAP. A relatively large axion isocurvature perturbation and its non-Gaussianity are predicted in our model. The saxion coherent oscillation has a large amplitude and dominates the Universe. The subsequent decay of the saxion produces huge amount of entropy, which dilutes unwanted relics. Winos, the lightest supersymmetric particles in this scenario, are produced non-thermally in the decay and account for dark matter.

We show that the hybrid inflation is naturally realized in the framework of a supersymmetric axion model, which is consistent with the WMAP observation if the Peccei-Quinn symmetry breaking scale is around 10^{15}GeV. By solving the post inflationary scalar dynamics, it is found that the scalar partner of the axion, saxion, oscillates with large amplitude and its decay produces a huge entropy and dilutes the axion. As a result, the axion coherent oscillation can be the dominant component of the dark matter in the Universe. Cosmological gravitino and axino problems are solved.

We show that a supersymmetric axion model naturally induces a hybrid inflation with the waterfall field identified as a Peccei-Quinn scalar. The Peccei-Quinn scale is predicted to be around 10^{15}GeV for reproducing the large-scale density perturbation of the Universe. After the built-in late-time entropy-production process, the axion becomes a dark matter candidate. Several cosmological implications are discussed.