Gravitational-wave detectors on earth have detected gravitational waves from merging compact objects in the local Universe. In future we will detect gravitational waves from higher-redshift sources, which trace the high-redshift structure formation history. That is, by observing high-redshift gravitational-wave events we will be able to probe structure formation history. This will provide additional insight into the early Universe when primordial fluctuations are generated and also into the nature of dark matter.

We explore the possibility of using a gigaparsec-scale local void to reconcile the Hubble tension. Such a gigaparsec-scale void can be produced by multi-stream inflation where different parts of the observable universe follow different inflationary trajectories. The impact of such a void for cosmological observations is studied, especially those involving supernovae, Baryon Acoustic Oscillations (BAO) and the kinetic Sunyaev-Zel'dovich (kSZ) effect. As a benchmark model, a 1.7Gpc scale with boundary width 0.7Gpc and density contrast -0.14 may ease the Hubble tension.

We revisit the effects of an early matter-dominated era on gravitational waves induced by scalar perturbations. We carefully take into account the evolution of the gravitational potential, the source of these induced gravitational waves, during a gradual transition from an early matter-dominated era to the radiation-dominated era, where the transition timescale is comparable to the Hubble time at that time. Realizations of such a gradual transition include the standard perturbative reheating with a constant decay rate. Contrary to previous works, we find that the presence of an early matter-dominated era does not necessarily enhance the induced gravitational waves due to the decay of the gravitational potential around the transition from an early matter-dominated era to the radiation-dominated era.

We study gravitational waves induced from the primordial scalar perturbations at second order around the reheating of the Universe. We consider reheating scenarios in which a transition from an early matter dominated era to the radiation dominated era completes within a timescale much shorter than the Hubble time at that time. We find that an enhanced production of induced gravitational waves occurs just after the reheating transition because of fast oscillations of scalar modes well inside the Hubble horizon. This enhancement mechanism just after an early matter-dominated era is much more efficient than a previously known enhancement mechanism during an early matter era, and we show that the induced gravitational waves could be detectable by future observations if the reheating temperature $T_{\text{R}}$ is in the range $T_\text{R} \lesssim 7\times 10^{-2}$GeV or $20 \, \text{GeV} \lesssim T_\text{R} \lesssim 2 \times 10^7 \, \text{GeV}$. This is the case even if the scalar perturbations on small scales are not enhanced relative to those on large scales, probed by the observations of the cosmic microwave background.

The possibility that primordial black hole binary mergers of stellar mass can explain the signals detected by the gravitational-wave interferometers has attracted much attention. In this scenario, primordial black holes can comprise only part of the entire dark matter, say, of order 0.1 %. This implies that most of the dark matter is accounted for by a different component, such as Weakly Interacting Massive Particles. We point out that in this situation, very compact dark matter minihalos, composed of the dominant component of the dark matter, are likely to be formed abundantly in the early Universe, with their formation redshift and abundance depending on primordial non-Gaussianity. They may be detected in future experiments via pulsar observations.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

We study the gravitational production of super-Hubble-mass dark matter in the very early universe. We first review the simplest scenario where dark matter is produced mainly during slow roll inflation. Then we move on to consider the cases where dark matter is produced during the transition period between inflation and the subsequent cosmological evolution. The limits of smooth and sudden transitions are studied, respectively. The relic abundances and the cosmological collider signals are calculated.

We show that the effect of initial non-Gaussian clustering can significantly enhance the event rate for primordial black hole (PBH) coalescence. The impact of such clustering is studied in a specific scenario of multi-stream inflation. Initial clustering enables the possibility of detecting massive PBH coalescence by space-based gravitational wave interferometers such as LISA and DECIGO/BBO. The parameter regime for the ground-based detectors to detect PBH coalescence is also extended.

Compared to primordial perturbations on large scales, roughly larger than $1$ megaparsec, those on smaller scales are not severely constrained. We revisit the issue of probing small-scale primordial perturbations using gravitational waves (GWs), based on the fact that, when large-amplitude primordial perturbations on small scales exist, GWs with relatively large amplitudes are induced at second order in scalar perturbations, and these induced GWs can be probed by both existing and planned gravitational-wave projects. We use accurate methods to calculate these induced GWs and take into account sensitivities of different experiments to induced GWs carefully, to report existing and expected limits on the small-scale primordial spectrum.

It is shown that a positive non-minimal coupling of the Higgs field to gravity can solve the two problems in inflation models in which postinflationary universe is dominated by an energy with stiff equation of state such as a kination, namely, overproduction of gravitons in gravitational reheating scenario, and overproduction of curvature perturbation from Higgs condensation. Furthermore, we argue that the non-minimal coupling parameter can be constrained more stringently with the progress in observations of large-scale structure and cosmic microwave background.

High energy collisions of particles may have created tiny black holes in the early Universe, which might leave stable remnants instead of fully evaporating as a result of Hawking radiation. If the reheating temperature was sufficiently close to the fundamental gravity scale, which can be different from the usual Planck scale depending of the presence and properties of spatial extra-dimensions, the formation rate could have been sufficiently high and hence such remnants could account for the entire cold dark matter of the Universe.

Recently, the formation of primordial black holes (PBHs) from the collapse of primordial fluctuations has received much attention. The abundance of PBHs formed during radiation domination is sensitive to the tail of the probability distribution of primordial fluctuations. We quantify the level of fine-tuning due to this sensitivity. For example, if the main source of dark matter is PBHs with mass $10^{-12}M_\odot$, then anthropic reasoning suggests that the dark matter to baryon ratio should range between 1 and 300. For this to happen, the root-mean-square amplitude of the curvature perturbation has to be fine-tuned within a $7.1\%$ range. As another example, if the recently detected gravitational-wave events are to be explained by PBHs, the corresponding degree of fine-tuning is $3.8\%$. We also find, however, that these fine-tunings can be relaxed if the primordial fluctuations are highly non-Gaussian, or if the PBHs are formed during an early-matter-dominated phase. We also note that no fine-tuning is needed for the scenario of a reheating of the universe by evaporated PBHs with Planck-mass relics left to serve as dark matter.

If primordial black holes (PBHs) form directly from inhomogeneities in the early Universe, then the number in the mass range $10^5 -10^{12}M_{\odot}$ is severely constrained by upper limits to the $\mu$ distortion in the cosmic microwave background (CMB). This is because inhomogeneities on these scales will be dissipated by Silk damping in the redshift interval $5\times 10^4\lesssim z\lesssim2\times 10^6$. If the primordial fluctuations on a given mass scale have a Gaussian distribution and PBHs form on the high-$\sigma$ tail, as in the simplest scenarios, then the $\mu$ constraints exclude PBHs in this mass range from playing any interesting cosmological role. Only if the fluctuations are highly non-Gaussian, or form through some mechanism unrelated to the primordial fluctuations, can this conclusion be obviated.

We revisit constraints on small-scale primordial power from annihilation signals from dark matter minihalos. Using gamma rays and neutrinos from extragalactic minihalos and assuming the delta-function primordial spectrum, we show the dependence of the constraints on annihilation modes, the mass of dark matter, and the annihilation cross section. We report both conservative constraints by assuming minihalos are fully destructed when becoming part of halos originating from the standard almost-scale invariant primordial spectrum, and optimistic constraints by neglecting destruction.

The small-scale crisis, discrepancies between observations and N-body simulations, may imply suppressed matter fluctuations on subgalactic distance scales. Such a suppression could be caused by some early-universe mechanism (e.g., broken scale invariance during inflation), leading to a modification of the primordial power spectrum at the onset of the radiation-domination era. Alternatively, it may be due to nontrivial dark-matter properties (e.g., new dark-matter interactions or warm dark matter) that affect the matter power spectrum at late times, during radiation domination, after the perturbations re-enter the horizon. We show that early- and late-time suppression mechanisms can be distinguished by measurement of the $\mu$ distortion to the frequency spectrum of the cosmic microwave background. This is because the $\mu$ distortion is suppressed, if the power suppression is primordial, relative to the value expected from the dissipation of standard nearly scale-invariant fluctuations. We emphasize that the standard prediction of the $\mu$ distortion remains unchanged in late-time scenarios even if the dark-matter effects occur before or during the era (redshifts $5\times 10^4 \lesssim z \lesssim 2\times 10^6$) at which $\mu$ distortions are generated.

Primordial black hole (PBH) mergers have been proposed as an explanation for the gravitational wave events detected by the LIGO collaboration. Such PBHs may be formed in the early Universe as a result of the collapse of extremely rare high-sigma peaks of primordial fluctuations on small scales, as long as the amplitude of primordial perturbations on small scales is enhanced significantly relative to the amplitude of perturbations observed on large scales. One consequence of these small-scale perturbations is generation of stochastic gravitational waves that arise at second order in scalar perturbations, mostly before the formation of the PBHs. These induced gravitational waves have been shown, assuming gaussian initial conditions, to be comparable to the current limits from the European Pulsar Timing Array, severely restricting this scenario. We show, however, that models with enhanced fluctuation amplitudes typically involve non-gaussian initial conditions. With such initial conditions, the current limits from pulsar timing can be evaded. The amplitude of the induced gravitational-wave background can be larger or smaller than the stochastic gravitational-wave background from supermassive black hole binaries.

Recently we have proposed a novel method to probe primordial gravitational waves from upper bounds on the abundance of primordial black holes (PBHs). When the amplitude of primordial tensor perturbations generated in the early Universe is fairly large, they induce substantial scalar perturbations due to their second-order effects. If these induced scalar perturbations are too large when they reenter the horizon, then PBHs are overproduced, their abundance exceeding observational upper limits. That is, primordial tensor perturbations on superhorizon scales can be constrained from the absence of PBHs. In our recent paper we have only shown simple estimations of these new constraints, and hence in this paper, we present detailed derivations, solving the Einstein equations for scalar perturbations induced at second order in tensor perturbations. We also derive an approximate formula for the probability density function of induced density perturbations, necessary to relate the abundance of PBHs to the primordial tensor power spectrum, assuming primordial tensor perturbations follow Gaussian distributions. Our new upper bounds from PBHs are compared with other existing bounds obtained from big bang nucleosynthesis, cosmic microwave background, LIGO/Virgo and pulsar timing arrays.

We propose a novel method to probe primordial gravitational waves by means of primordial black holes (PBHs). When the amplitude of primordial tensor perturbations on comoving scales much smaller than those relevant to Cosmic Microwave Background is very large, it induces scalar perturbations due to second-order effects substantially. If the amplitude of resultant scalar perturbations becomes too large, then PBHs are overproduced to a level that is inconsistent with a variety of existing observations constraining their abundance. This leads to upper bounds on the amplitude of initial tensor perturbations on super-horizon scales. These upper bounds from PBHs are compared with other existing bounds.

We construct cosmological long-wavelength solutions without symmetry in general gauge conditions compatible with the long-wavelength scheme. We then specify the relationship among the solutions in different time slicings. Applying this general framework to spherical symmetry, we derive the correspondence relation between long-wavelength solutions in the constant mean curvature slicing with conformally flat spatial coordinates and asymptotic quasihomogeneous solutions in the comoving gauge and compare the numerical results of PBH formation in these two different approaches. To discuss the PBH formation, it is convenient and conventional to use $\tilde{\delta}_{c}$, the value which the averaged density perturbation at threshold in the comoving slicing would take at horizon entry in the lowest-order long-wavelength expansion. We numerically find that within compensated models, the sharper the transition from the overdense region to the FRW universe is, the larger the $\tilde{\delta}_{c}$ becomes. We suggest that, for the equation of state $p=(\Gamma-1)\rho$, we can apply the analytic formulas for the minimum $\tilde{\delta}_{c, {\rm min}}\simeq [3\Gamma/(3\Gamma+2)]\sin^{2}\left[\pi\sqrt{\Gamma-1}/(3\Gamma-2)\right]$ and the maximum $\tilde{\delta}_{c, {\rm max}}\simeq 3\Gamma/(3\Gamma+2)$. As for the threshold peak value of the curvature variable $\psi_{0,c}$, we find that the sharper the transition is, the smaller the $\psi_{0,c}$ becomes. We analytically explain this feature. Using simplified models, we also analytically deduce an environmental effect that $\psi_{0,c}$ can be significantly larger (smaller) if the underlying density perturbation of much longer wavelength is positive (negative).

Supermassive black holes and intermediate mass black holes are believed to exist in the Universe. There is no established astrophysical explanation for their origin and considerations have been made in the literature that those massive black holes (MBHs) may be primordial black holes (PBHs), black holes which are formed in the early universe (well before the matter-radiation equality) due to the direct collapse of primordial overdensities. This paper aims at discussing the possibility of excluding the PBH scenario as the origin of the MBHs. We first revisit the constraints on PBHs obtained from the CMB distortion that the seed density perturbation causes. By adopting a recent computation of the CMB distortion sourced by the seed density perturbation and the stronger constraint on the CMB distortion set by the COBE/FIRAS experiment used in the literature, we find that PBHs in the mass range $6\times 10^4~M_\odot \sim 5 \times 10^{13}~M_\odot$ are excluded. Since PBHs lighter than $6 \times 10^4~M_\odot$ are not excluded from the non-observation of the CMB distortion, we propose a new method which can potentially exclude smaller PBHs as well. Based on the observation that large density perturbations required to create PBHs also result in the copious production of ultracompact minihalos (UCMHs), compact dark matter halos formed at around the recombination, we show that weakly interacting massive particles (WIMPs) as dark matter annihilate efficiently inside UCMHs to yield cosmic rays far exceeding the observed flux. Our bound gives severe restriction on the compatibility between the particle physics models for WIMPs and the PBH scenario as the explanation of MBHs.

Primordial black holes (PBHs) are a useful tool in cosmology to probe primordial inhomogeneities on small scales that reenter the Hubble radius during the radiation dominated epoch. In this paper, a phenomenon we call the double formation of PBHs, described below, is explored. Suppose there exists a highly perturbed region which will collapse to form a PBH after the horizon crossing of this region, and farther that this region is superposed on much larger region, which also collapses upon reentry. One then expects the collapse of the central smaller region at the time of the crossing of this region, followed by another collapse of the larger perturbation at the time of its respective crossing. The smaller PBH, formed earlier, should be swallowed in the second collapse leading to a single larger PBH as the final state. This paper reports the first direct numerical confirmation of such double PBH formation. Related to this, we also discuss the effects of high-frequency modes on the formation of PBHs.

We provide a simple but robust bound on the primordial curvature perturbation in the range $10^4 {\rm Mpc}^{-1} < k < 10^5 {\rm Mpc}^{-1}$, which has not been constrained so far unlike low wavenumber modes. Perturbations on these scales dissipate the energy of their acoustic oscillations by the Silk damping after primordial nucleosynthesis but before the redshift $z \sim 2 \times 10^6$ and reheat the photon bath without invoking CMB distortions. This {\it acoustic reheating} results in the decrease of the baryon-photon ratio. By combining independent measurements probing the nucleosynthesis era and around the recombination epoch, we find an upper bound on the amplitude of the curvature perturbation over the above wavenumber range as ${\cal P}_\zeta < 0.06$. Implications for super massive black holes are also discussed.

Primordial black holes (PBHs) are an important tool in cosmology to probe the primordial spectrum of small-scale curvature perturbations that reenter the cosmological horizon during radiation domination epoch. We numerically solve the evolution of spherically symmetric highly perturbed configurations to clarify the criteria of PBHs formation using an extremely wide class of curvature profiles characterized by five parameters, (in contrast to only two parameters used in all previous papers) which specify the curvature profiles not only at the central region but also at the outer boundary of configurations. It is shown that formation or non-formation of PBHs is determined entirely by only two master parameters one of which can be presented as an integral of curvature over initial configurations and the other is presented in terms of the position of the boundary and the edge of the core.

For an arbitrary strong, spherically symmetric super-horizon curvature perturbation, we present analytical solutions of the Einstein equations in terms of asymptotic expansion over the ratio of the Hubble radius to the length-scale of the curvature perturbation under consideration. To obtain this solution we develop a recursive method of quasi-linearization which reduces the problem to a system of coupled ordinary differential equations for the $N$-th order terms in the asymptotic expansion with sources consisting of a non-linear combination of the lower order terms. We use this solution for setting initial conditions for subsequent numerical computations. For an arbitrary precision requirement predetermined by the intended accuracy and stability of the computer code, our analytical solution yields optimal truncated asymptotic expansion which can be used to find the upper limit on the moment of time when the initial conditions expressed in terms of such truncated expansion should be set. Examples of how these truncated (up to eighth order) solutions provide initial conditions with given accuracy for different radial profiles of curvature perturbations are presented.