In the context of string theory, several conjectural conditions have been proposed for low energy effective field theories not to be in swampland, the UV-incomplete class. The recent ones represented by the de Sitter and trans-Planckian censorship conjectures in particular seem to conflict with the inflation paradigm of the early universe. We first point out that scenarios where inflation is repeated several times (multi-phase inflation) can be easily compatible with these conjectures. In other words, we relax the constraint on the single inflation for the large scale perturbations to only continue at least around 10 e-folds. In this context, we then investigate if a spectator field can be a source of the almost scale-invariant primordial perturbations on the large scale. As a consequence of such an isocurvature contribution, the resultant perturbations exhibit the non-vanishing non-Gaussianity in general. Also the perturbation amplitude on smaller scales can be completely different from that on the large scale due to the multiplicity of inflationary phases. These signatures will be a smoking gun of this scenario by the future observations.

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 investigate a possibility of primordial black hole (PBH) formation with a hierarchical mass spectrum in multiple phases of inflation. As an example, we find that one can simultaneously realize a mass spectrum which has recently attracted a lot of attention: stellar-mass PBHs ($\sim\mathcal{O}(10)M_\odot$) as a possible source of binary black holes detected by LIGO/Virgo collaboration, asteroid-mass ($\sim\mathcal{O}(10^{-12})M_\odot$) as a main component of dark matter, and earth-mass ($\sim\mathcal{O}(10^{-5})M_\odot$) as a source of ultrashort-timescale events in Optical Gravitational Lensing Experiment microlensing data. The recent refined de Sitter swampland conjecture may support such a multi-phase inflationary scenario with hierarchical mass PBHs as a transition signal of each inflationary phase.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

The stochastic approach aims at describing the long-wavelength part of quantum fields during inflation by a classical stochastic theory. It is usually formulated in terms of Langevin equations, giving rise to a Fokker-Planck equation for the probability distribution function of the fields, and possibly their momenta. The link between these two descriptions is ambiguous in general, as it depends on an implicit discretisation procedure, the two prominent ones being the It\^o and Stratonovich prescriptions. Here we show that the requirement of general covariance under field redefinitions is verified only in the latter case, however at the expense of introducing spurious `frame' dependences. This stochastic anomaly disappears when there is only one source of stochasticity, like in slow-roll single-field inflation, but manifests itself when taking into account the full phase space, or in the presence of multiple fields. Despite these difficulties, we use physical arguments to write down a covariant Fokker-Planck equation that describes the diffusion of light scalar fields in non-linear sigma models in the overdamped limit. We apply it to test scalar fields in de Sitter space and show that some statistical properties of a class of two-field models with derivative interactions can be reproduced by using a correspondence with a single-field model endowed with an effective potential. We also present explicit results in a simple extension of the single-field $\lambda \phi^4$ theory to a hyperbolic field space geometry. The difficulties we describe seem to be the stochastic counterparts of the notoriously difficult problem of maintaining general covariance in quantum theories, and the related choices of operator ordering and path-integral constructions. Our work thus opens new avenues of research at the crossroad between cosmology, statistical physics, and quantum field theory.

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

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

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

We show that a detectable tensor-to-scalar ratio $(r\ge 10^{-3})$ on the CMB scale can be generated even during extremely low energy inflation which saturates the BBN bound $\rho_{\rm inf}\approx (30 {\rm MeV})^4$. The source of the gravitational waves is not quantum fluctuations of graviton but those of $SU(2)$ gauge fields, energetically supported by coupled axion fields. The curvature perturbation, the backreaction effect and the validity of perturbative treatment are carefully checked. Our result indicates that measuring $r$ alone does not immediately fix the inflationary energy scale.

The prospects of future galaxy surveys for non-Gaussianity measurements call for the development of robust techniques for computing the bispectrum of primordial cosmological perturbations. In this paper, we propose a novel approach to the calculation of the squeezed bispectrum in multiple-field inflation. With use of the $\delta N$ formalism, our framework sheds new light on the recently pointed out difference between the squeezed bispectrum for global observers and that for local observers, while allowing one to calculate both. For local observers in particular, the squeezed bispectrum is found to vanish in single-field inflation. Furthermore, our framework allows one to go beyond the near-equilateral ("small hierarchy") limit, and to automatically include intrinsic non-Gaussianities that do not need to be calculated separately. The explicit computational programme of our method is given and illustrated with a few examples.

We propose a novel scenario to produce abundant primordial black holes (PBHs) in new inflation which is a second phase of a double inflation in the supergravity frame work. In our model, some preinflation phase before the new inflation is assumed and it would be responsible for the primordial curvature perturbations on the cosmic microwave background scale, while the new inflation produces only the small scale perturbations. Our new inflation model has linear, quadratic, and cubic terms in its potential and PBH production corresponds with its flat inflection point. The linear term can be interpreted to come from a supersymmetry-breaking sector, and with this assumption, the vanishing cosmological constant condition after inflation and the flatness condition for the inflection point can be consistently satisfied.

We revisit the constraints on the small scale density perturbations ($10^4\,\mathrm{Mpc}^{-1}\lesssim k \lesssim10^5\,\mathrm{Mpc}^{-1}$) from the modification of the freeze-out value of the neutron-proton ratio at big-bang nucleosynthesis era. Around the freeze-out temperature $T\sim 0.5\,\mathrm{MeV}$, the universe can be divided into several local patches which have different temperatures since any perturbation which enters the horizon after the neutrino decoupling has not diffused yet. Taking account of this situation, we calculate the freeze-out value in detail. We find that the small scale perturbations decrease the n-p ratio in contrast to previous works. With use of the latest observed $^4$He abundance, we obtain the constraint on the power spectrum of the curvature perturbations as $\Delta^2_\mathcal{R}\lesssim 0.018$ on $10^4\,\mathrm{Mpc}^{-1}\lesssim k \lesssim 10^5\,\mathrm{Mpc}^{-1}$.

We studied the possibility whether the massive primordial black holes (PBHs) surviving today can be produced in hybrid inflation. Though it is of great interest since such PBHs can be the candidate for dark matter or seeds of the supermassive black holes in galaxies, there have not been quantitatively complete works yet because of the non-perturbative behavior around the critical point of hybrid inflation. Therefore, combining the stochastic and $\delta N$ formalism, we numerically calculated the curvature perturbations in a non-perturbative way and found, without any specific assumption of the types of hybrid inflation, PBHs are rather overproduced when the waterfall phase of hybrid inflation continues so long that the PBH scale is well enlarged and the corresponding PBH mass becomes sizable enough.

Primordial black holes (PBHs) are theoretical black holes which may be formed during the radiation dominant era and, basically, caused by the gravitational collapse of radiational overdensities. It has been well known that in the context of the structure formation in our Universe such collapsed objects, e.g., halos/galaxies, could be considered as bias tracers of underlying matter fluctuations and the halo/galaxy bias has been studied well. Employing a peak-background split picture which is known to be a useful tool to discuss the halo bias, we consider the large scale clustering behavior of the PBH and propose an almost mass-independent constraint to the scenario that dark matters (DMs) consist of PBHs. We consider the case where the statistics of the primordial curvature perturbations is almost Gaussian, but with small local-type non-Gaussianity. If PBHs account for the DM abundance, such a large scale clustering of PBHs behaves as nothing but the matter isocurvature perturbation and constrained strictly by the observations of cosmic microwave backgrounds (CMB). From this constraint, we show that, in the case a certain single field causes both CMB temperature perturbations and PBH formations, the PBH-DM scenario is excluded even with quite small local-type non-Gaussianity, $|f_\mathrm{NL}|\sim\mathcal{O}(0.01)$, while we give the constraints to parameters in the case where the source field of PBHs is different from CMB perturbations.

There has been a growing evidence for the existence of magnetic fields in the extra-galactic regions, while the attempt to associate their origin with the inflationary epoch alone has been found extremely challenging. We therefore take into account the consistent post-inflationary evolution of the magnetic fields that are originated from vacuum fluctuations during inflation. In the model of our interest, the electromagnetic (EM) field is coupled to a pseudo-scalar inflaton $\phi$ through the characteristic term $\phi F\tilde F$, breaking the conformal invariance. This interaction dynamically breaks the parity and enables a continuous production of only one of the polarization states of the EM field through tachyonic instability. The produced magnetic fields are thus helical. We find that the dominant contribution to the observed magnetic fields in this model comes from the modes that leave the horizon near the end of inflation, further enhanced by the tachyonic instability right after the end of inflation. The EM field is subsequently amplified by parametric resonance during the period of inflaton oscillation. Once the thermal plasma is formed (reheating), the produced helical magnetic fields undergo a turbulent process called inverse cascade, which shifts their peak correlation scales from smaller to larger scales. We consistently take all these effects into account within the regime where the perturbation of $\phi$ is negligible and obtain $B_{\rm eff} \sim 10^{-19}$G, indicating the necessity of additional mechanisms to accommodate the observations.

We calculate the CMB $\mu$-distortion and the angular power spectrum of its cross-correlation with the temperature anisotropy in the presence of the non-Gaussian neutrino isocurvature density (NID) mode. While the pure Gaussian NID perturbations give merely subdominant contribution to $<\mu>$ and vanishing $< \mu T>$, the latter quantity can be large enough to be detected in the future when the NID perturbations $\mathcal S(\mathbf x)$ are proportional to the square of a Gaussian field $g(\mathbf x)$, i.e. $\mathcal S({\mathbf x})\propto g^2({\mathbf x})$. In particular, large $< \mu T>$ can be realized since Gaussian-squared perturbations can yield a relatively large bispectrum, satisfying the constraints from the power spectrum of CMB anisotropies, i.e. $\mathcal{P}_\mathcal{SS}(k_0) \sim\mathcal{P}_g^2(k_0)\lesssim10^{-10}$ at $k_0=0.05$ Mpc$^{-1}$. We also forecast constraints from the CMB temperature and E-mode polarisation bispectra, and show that $\mathcal{P}_g(k_0)\lesssim10^{-5}$ would be allowed from Planck data. We find that $< \mu >$ and $|l(l+1)C^{\mu T}_l|$ can respectively be as large as $10^{-9}$ and $10^{-14}$ with uncorrelated scale-invariant NID perturbations for $\mathcal{P}_g(k_0)=10^{-5}$. When the spectrum of the Gaussian field is blue-tilted (with spectral index $n_g \simeq 1.5$), $< \mu T>$ can be enhanced by an order of magnitude.

In our previous paper, we have proposed a new algorithm to calculate the power spectrum of the curvature perturbations generated in inflationary universe with use of the stochastic approach. Since this algorithm does not need the perturbative expansion with respect to the inflaton fields on super-horizon scale, it works even in highly stochastic cases. For example, when the curvature perturbations are very large or the non-Gaussianities of the curvature perturbations are sizable, the perturbative expansion may break down but our algorithm enables to calculate the curvature perturbations. We apply it to two well-known inflation models, chaotic and hybrid inflation, in this paper. Especially for hybrid inflation, while the potential is very flat around the critical point and the standard perturbative computation is problematic, we successfully calculate the curvature perturbations.

We propose a new approach for calculating the curvature perturbations produced during inflation in the stochastic formalism. In our formalism, the fluctuations of the e-foldings are directly calculated without perturbatively expanding the inflaton field and they are connected to the curvature perturbations by the $\delta N$ formalism. The result automatically includes the contributions of the higher order perturbations because we solve the equation of motion non-perturbatively. In this paper, we analytically prove that our result (the power spectrum and the nonlinearity parameter) is consistent with the standard result in single field slow-roll inflation. We also describe the algorithm for numerical calculations of the curvature perturbations in more general inflation models.