Ultralight primordial black holes (PBHs) with masses $\lesssim 10^{15}$g and subatomic Schwarzschild radii, produced in the early Universe, are expected to have evaporated by the current cosmic age due to Hawking radiation. Based on this assumption, a number of constraints on the abundance of ultralight PBHs have been made. However, Hawking radiation has thus far not been verified experimentally. It would, therefore, be of interest if constraints on ultralight PBHs could be placed independent of the assumption of Hawking-radiation. In this paper, we explore the possibility of probing these PBHs, within a narrow mass range, using gravitational-wave (GW) data from the two LIGO detectors. The idea is that large primordial curvature perturbations that result in the formation of PBHs, would also generate GWs through non-linear mode couplings. These induced GWs would produce a stochastic background. Specifically, we focus our attention on PBHs of mass range $\sim 10^{13} - 10^{15}$g for which the induced stochastic GW background peak falls in the sensitivity band of LIGO. We find that, assuming a monochromatic PBH mass distribution, the corresponding GW background would be detectable using presently available LIGO data, provided we neglect the existing constraints on the abundance of PBHs, which are based on Hawking radiation. For more realistic (broader) mass distributions, we find that this stochastic background would be detectable in LIGO's third observing run, even after considering the existing constraints on PBH abundance. A non-detection should enable us to constrain the amplitude of primordial curvature perturbations as well as the abundance of ultralight PBHs.

It is known that geometrical optics no longer applies to the gravitational lensing if the wavelength of a propagating wave becomes comparable to or larger than the Schwarzshild radius of a lensing object. We investigate the propagation of gravitational waves in wave optics, particularly focusing on the difference between their arrival time and the arrival time of light. We argue that, contrary to the observation in the previous work, gravitational waves never arrive at an observer earlier than light when both gravitational waves and light are emitted from a same source simultaneously.

Computations of the primordial black hole (PBH) mass function discussed in the literature have conceptual issues. They stem from that the mass function is a differential quantity and the standard criterion of the PBH formation from the seed primordial fluctuations cannot be directly applied to the computation of the differential quantities. We propose a new criterion of the PBH formation which is an addition of one extra condition to the existing one. By doing this, we derive a formal expression of the PBH mass function without introducing any ambiguous interpretations which exist in the previous studies. Once the underlying primordial fluctuations are specified, the PBH mass function can be in principle determined by the new formula. As a demonstration of our formulation, we compute the PBH mass function analytically for the case where the perturbations are Gaussian and the space is one dimension.

We propose a simple but novel cosmological scenario where both the Planck mass and the dark energy scale emerge from the same super-Hubble quantum fluctuations of a non-minimally coupled ultra-light scalar field during primordial inflation. The current cosmic and solar-system observations constrain the non-minimal coupling to be small.

We formulate the two-point correlation function of primordial black holes (PBHs) at their formation time, based on the functional integration approach which has often been used in the context of halo clustering. We find that PBH clustering on super-Hubble scales could never be induced in the case where the initial primordial fluctuations are Gaussian, while it can be enhanced by the so-called local-type trispectrum (four-point correlation function) of the primordial curvature perturbations.

Within the framework that primordial black holes are formed by the direct gravitational collapse of large primordial density perturbations in the radiation dominated stage, we derive the threshold of the density contrast for the formation of rotating primordial black holes based on the simple Jeans criterion. It is found that the threshold value increases in proportion to the square of the angular momentum. We then apply the recently refined analysis on the formation threshold for non-rotating black holes to the case of rotating black holes, and contrast the derived threshold with the former. Caveats and effects ignored in our analysis are also presented, which suggests that the uncertainties of our result can be addressed only by means of numerical relativity.

We perform a fully relativistic analysis of even-parity linear perturbations around a static and spherically symmetric solution in the most general scalar-tensor theory with second-order field equations. This paper is a sequel to Kobayashi {\em et al.} (2012), in which the linear perturbation analysis for the odd-parity modes is presented. Expanding the Horndeski action to second order in perturbations and eliminating auxiliary variables, we derive the quadratic action for even-parity perturbations written solely in terms of two dynamical variables. The two perturbations can be interpreted as the gravitational and scalar waves. Correspondingly, we obtain two conditions to evade ghosts and two conditions for the absence of gradient instabilities. Only one in each pair of conditions yields a new stability criterion, as the conditions derived from the stability of the gravitational-wave degree of freedom coincide with those in the odd-parity sector. Similarly, the propagation speed of one of the two modes is the same as that for the odd-parity mode, while the other differs in general from them. Our result is applicable to all the theories of gravitation with an extra single scalar degree of freedom such as the Brans-Dicke theory, $f(R)$ models, and Galileon gravity.

We point out that the gravitational-wave event GW150914 observed by the LIGO detectors can be explained by the coalescence of primordial black holes (PBHs). It is found that the expected PBH merger rate would exceed the rate estimated by the LIGO scientific Collaboration and Virgo Collaboration if PBHs were the dominant component of dark matter, while it can be made compatible if PBHs constitute a fraction of dark matter. Intriguingly, the abundance of PBHs required to explain the suggested lower bound on the event rate, $> 2$ events ${\rm Gpc}^{-3} {\rm yr}^{-1}$, roughly coincides with the existing upper limit set by the nondetection of the cosmic microwave background spectral distortion. This implies that the proposed PBH scenario may be tested in the not-too-distant future.

As the first step to extend our understanding of higher-derivative theories, within the framework of analytic mechanics of point particles, we construct a ghost-free theory involving third-order time derivatives in Lagrangian. While eliminating linear momentum terms in the Hamiltonian is necessary and sufficient to kill the ghosts associated with higher derivatives for Lagrangian with at most second-order derivatives, we find that this is necessary but not sufficient for the Lagrangian with higher than second-order derivatives. We clarify a set of ghost-free conditions under which we show that the Hamiltonian is bounded, and that equations of motion are reducible into a second-order system.

This is a review article on the primordial black holes (PBHs), with particular focus on the massive ones ($\gtrsim 10^{15}{\rm g}$) which have not evaporated by the present epoch by the Hawking radiation. By the detections of gravitational waves by LIGO, we have gained a completely novel tool to observationally search for PBHs complementary to the electromagnetic waves. Based on the perspective that gravitational-wave astronomy will make a significant progress in the next decades, a purpose of this article is to give a comprehensive review covering a wide range of topics on PBHs. After discussing PBH formation as well as several inflation models leading to PBH production, we summarize various existing and future observational constraints. We then present topics on formation of PBH binaries, gravitational waves from PBH binaries, various observational tests of PBHs by using gravitational waves.

It has been proposed that primordial black holes (PBHs) form binaries in the radiation dominated era. Once formed, some fraction of them may merge within the age of the Universe by gravitational radiation reaction. We investigate the merger rate of the PBH binaries when the PBHs have a distribution of masses around ${\cal O}(10) M_\odot$, which is a generalization of the previous studies where the PBHs are assumed to have the same mass. After deriving a formula for the merger time probability distribution in the PBH mass plane, we evaluate it under two different approximations. We identify a quantity constructed from the mass-distribution of the merger rate density per unit cosmic time and comoving volume $\mathcal{R}(m_1,m_2)$, $\alpha = -{(m_1+m_2)}^2\partial^2 \ln\mathcal{R}/\partial m_1\partial m_2 $, which universally satisfies $0.97 \lesssim \alpha \lesssim 1.05$ for all binary masses independently of the PBH mass function. This result suggests that the measurement of this quantity is useful for testing the PBH scenario.

If cosmic strings are formed in the early universe, their associated loops emit gravitational waves during the whole cosmic history and contribute to the stochastic gravitational wave background at all frequencies. We provide a new estimate of the stochastic gravitational wave spectrum by considering a realistic cosmological loop distribution, in scaling, as it can be inferred from Nambu-Goto numerical simulations. Our result takes into account various effects neglected so far. We include both gravitational wave emission and backreaction effects on the loop distribution and show that they produce two distinct features in the spectrum. Concerning the string microstructure, in addition to the presence of cusps and kinks, we show that gravitational wave bursts created by the collision of kinks could dominate the signal for wiggly strings, a situation which may be favoured in the light of recent numerical simulations. In view of these new results, we propose four prototypical scenarios, within the margin of the remaining theoretical uncertainties, for which we derive the corresponding signal and estimate the constraints on the string tension put by both the LIGO and European Pulsar Timing Array (EPTA) observations. The less constrained of these scenarios is shown to have a string tension GU < 7.2 x 10^{-11}, at 95% of confidence. Smooth loops carrying two cusps per oscillation verify the two-sigma bound GU < 1.0 x 10^{-11} while the most constrained of all scenarios describes very kinky loops and satisfies GU < 6.7 x 10^{-14} at 95% of confidence.

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.

We perform a fully relativistic analysis of odd-type linear perturbations around a static and spherically symmetric solution in the most general scalar-tensor theory with second-order field equations in four-dimensional spacetime. It is shown that, as in the case of general relativity, the quadratic action for the perturbations reduces to the one having only a single dynamical variable, from which concise formulas for no-ghost and no-gradient instability conditions are derived. Our result is applicable to all the theories of gravity with an extra scalar degree of freedom. We demonstrate how the generic formulas can be applied to some particular examples such as the Brans-Dicke theory, $f(R)$ models, and Galileon gravity.

We investigate the internal structure and the mass-radius relation of neutron stars in a recently proposed scalar-tensor theory dubbed asymmetron in which a massive scalar field undergoes spontaneous scalarization inside neutron stars. We focus on the case where the Compton wavelength is shorter than 10 km, which has not been investigated in the literature. By solving the modified Einstein equations, either purely numerically or by partially using a semianalytic method, we find that not only the weakening of gravity by spontaneous scalarization but also the scalar force affect the internal structure significantly in the massive case. We also find that the maximum mass of neutron stars is larger for certain parameter sets than that in general relativity and reaches 2 solar mass even if the effect of strange hadrons is taken into account. There is even a range of parameters where the maximum mass of neutron stars largely exceeds the threshold that violates the causality bound in general relativity.

We present an impact of coupling between dark matter and a scalar field, which might be responsible for dark energy, on measurements of redshift-space distortions. We point out that, in the presence of conformal and/or disformal coupling, linearized continuity and Euler equations for total matter fluid significantly deviate from the standard ones even in the sub-horizon scales. In such a case, a peculiar velocity of total matter field is determined not only by a logarithmic time derivative of its density perturbation but also by density perturbations for both dark matter and baryon, leading to a large modification of the physical interpretation of observed data obtained by measurements of redshift-space distortions. We reformulate galaxy two-point correlation function in the redshift space based on the modified continuity and Euler equations. We conclude from the resultant formula that the true value of the linear growth rate of large-scale structure cannot be necessarily constrained by single-redshift measurements of the redshift-space distortions, unless one observes the actual time-evolution of structure.

An invertible field transformation is such that the old field variables correspond one-to-one to the new variables. As such, one may think that two systems that are related by an invertible transformation are physically equivalent. However, if the transformation depends on field derivatives, the equivalence between the two systems is nontrivial due to the appearance of higher derivative terms in the equations of motion. To address this problem, we prove the following theorem on the relation between an invertible transformation and Euler-Lagrange equations: If the field transformation is invertible, then any solution of the original set of Euler-Lagrange equations is mapped to a solution of the new set of Euler-Lagrange equations, and vice versa. We also present applications of the theorem to scalar-tensor theories.

Regardless of the long history of gauge theories, it is not well recognized under which condition gauge fixing at the action level is legitimate. We address this issue from the Lagrangian point of view, and prove the following theorem on the relation between gauge fixing and Euler-Lagrange equations: In any gauge theory, if a gauge fixing is complete, i.e., the gauge functions are determined uniquely by the gauge conditions, the Euler-Lagrange equations derived from the gauge-fixed action are equivalent to those derived from the original action supplemented with the gauge conditions. Otherwise, it is not appropriate to impose the gauge conditions before deriving Euler-Lagrange equations as it may in general lead to inconsistent results. The criterion to check whether a gauge fixing is complete or not is further investigated. We also provide applications of the theorem to scalar-tensor theories and make comments on recent relevant papers on theories of modified gravity, in which there are confusions on gauge fixing and counting physical degrees of freedom.

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.

In the context of classical mechanics, we study the conditions under which higher-order derivative theories can evade the so-called Ostrogradsky instability. More precisely, we consider general Lagrangians with second order time derivatives, of the form $L(\ddot\phi^a,\dot\phi^a,\phi^a;\dot q^i,q^i)$ with $a = 1,\cdots, n$ and $i = 1,\cdots, m$. For $n=1$, assuming that the $q^i$'s form a nondegenerate subsystem, we confirm that the degeneracy of the kinetic matrix eliminates the Ostrogradsky instability. The degeneracy implies, in the Hamiltonian formulation of the theory, the existence of a primary constraint, which generates a secondary constraint, thus eliminating the Ostrogradsky ghost. For $n>1$, we show that, in addition to the degeneracy of the kinetic matrix, one needs to impose extra conditions to ensure the presence of a sufficient number of secondary constraints that can eliminate all the Ostrogradsky ghosts. When these conditions that ensure the disappearance of the Ostrogradsky instability are satisfied, we show that the Euler-Lagrange equations, which involve a priori higher order derivatives, can be reduced to a second order system.

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 systematically show that in potential driven generalized G-inflation models, quantum corrections coming from new physics at the strong coupling scale can be avoided, while producing observable tensor modes. The effective action can be approximated by the tree level action, and as a result, these models are internally consistent, despite the fact that we introduced new mass scales below the energy scale of inflation. Although observable tensor modes are produced with sub-strong coupling scale field excursions, this is not an evasion of the Lyth bound, since the models include higher-derivative non-canonical kinetic terms, and effective rescaling of the field would result in super-Planckian field excursions. We argue that the enhanced kinetic term of the inflaton screens the interactions with other fields, keeping the system weakly coupled during inflation.

We propose a new scalar-tensor model which induces significant deviation from general relativity inside dense objects like neutron stars, while passing solar-system and terrestrial experiments, extending a model proposed by Damour and Esposito-Farese. Unlike their model, we employ a massive scalar field dubbed asymmetron so that it not only realizes proper cosmic evolution but also can account for the cold dark matter. In our model, asymmetron undergoes spontaneous scalarization inside dense objects, which results in reduction of the gravitational constant by a factor of order unity. This suggests that observational tests of constancy of the gravitational constant in high density phase are the effective ways to look into the asymmetron model.

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.

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.

Scalar field with generalized kinetic interactions metamorphoses depending on its field value, ranging from cosmological constant to stiff matter. We show that such a scalar field can give rise to temporal enhancement of the curvature perturbation in the primordial Universe, leading to efficient production of primordial black holes while the enhancement persists. If the inflation energy scale is high, those mini-black holes evaporate by the Hawking radiation much before Big Bang nucleosynthesis and the effective reheating of the Universe is achieved by the black hole evaporation. Dominance of PBHs and the reheating by their evaporation modify the expansion history of the primordial Universe. This results in a characteristic feature of the spectrum of primordial tensor modes in the DECIGO frequency band, opening an interesting possibility of testing PBH reheating scenario by measuring the primordial tensor modes. If the inflation energy scale is low, the PBH mass can be much larger than the solar mass. In this case, PBH is an interesting candidate for seeds for supermassive black holes residing in present galaxies.

Within the framework of the standard quantum electrodynamics, we compute contribution of vacuum polarization at one-loop order to the power spectrum of the magnetic field on inflationary (de Sitter) background. It is found that the one-loop term exhibits the infrared secular growth that is proportional to the number of $e$-folds. The use of the dynamical renormalization group method, which amounts to partial resummation of higher loop diagrams, shows that the resummed power spectrum is free from the secular growth and the loop effect only changes the power from $(k/a)^4$ at the tree level to ${(k/a)}^{4-\nu}$, where $\nu$ represents the contribution from vacuum polarization. The parameter $\nu$, being proportional to the square of the gauge coupling constant as well as the number of fermion species, is a simple function of a ratio of fermion mass to the Hubble parameter and is positive irrespective of the fermion mass. Thus, the loop effect always enhances the infrared magnetic field strength. We find that $\nu \simeq 5 \times 10^{-3}$ is the possible maximum contribution to $\nu$ from a single fermion. This estimate suggests that either large number of fermion species or large coupling constant is a necessary condition for the loop effect to be responsible for the seed of the cosmic magnetic fields.

We study the dynamics of the oscillating gauged scalar field in a thermal bath. A Langevin type equation of motion of the scalar field, which contains both dissipation and fluctuation terms, is derived by using the real-time finite temperature effective action approach. The existence of the quantum fluctuation-dissipation relation between the non-local dissipation term and the Gaussian stochastic noise terms is verified. We find the noise variables are anti-correlated at equal-time. The dissipation rate for the each mode is also studied, which turns out to depend on the wavenumber.

We revisit the method of using the photon-graviton conversion mechanism in the presence of the external magnetic field to probe small-scale primordial magnetic fields that may exist between the last scattering surface and present. Specifically, we investigate impacts on the conversion efficiency due to the presence of matter, including the plasma collective effect and the atomic polarizability. In general, these effects tend to reduce the conversion probability. Under this more realistic picture and based on the precision of COBE's measurement of CMB (cosmic microwave background) blackbody spectrum, we find an upper bound for the primordial magnetic field strength, B < 30G, at the time of recombination. Although at present the bound based on the photon-graviton conversion mechanism is not as tight as that obtained by the direct use of CMB temperature anisotropy, it nevertheless provides an important independent constraint on primordial magnetic fields and at epochs in addition to the recombination. The bound can be significantly improved if the CMB blackbody spectrum measurement becomes more precise in future experiments such as PIXIE.

We discuss implications of Planck results for models with local type non-Gaussianity. In light of the recent results of the Planck satellite, we constrain model parameters of several representative models and give the prediction of trispectrum, in particular, gNL. We also consider interesting possibilities that trispectrum appears as the first signature of the non-Gaussianities of the curvature perturbations, that is, fNL is small while gNL can be significantly large.

Effects invoked by the passage of the cosmic string through the Earth are investigated. The cosmic string induces global oscillations of the Earth whose amplitude and acceleration both linearly depend on the string line density. For the line density maximally allowed by the cosmic observations, the oscillations are perceivable even to human beings and may cause serious damages to the environment. Use of the sophisticated accelerograph enables us to detect the string of a line density down to ten orders of magnitude smaller than the cosmologically relevant value.

We provide a general formula for calculating correlators of arbitrary function of a Gaussian field. This work extends the standard leading-order approximation based on the delta N formalism to the case where truncation of the delta N at some low order does not yield the correct answer. As an application of this formula, we investigate 2, 3 and 4-point functions of the primordial curvature perturbation generated in the massless preheating model by approximating the mapping between the curvature perturbation and the Gaussian field as a sum of the many spiky normal distribution functions as suggested by lattice calculations. We also discuss observational consequences of this case and show that trispectrum would be a key observable to search signature of preheating in the CMB map. It is found the forms of the curvature correlation functions for any delta N, at the leading order in the correlator of the Gaussian field, coincide with the standard local type ones. Within this approximation, it is also found that the standard formula for the non-linearity parameters given by the product of the derivatives of the e-folding number still holds after we replace the bare e-folding number appearing in the original delta N expansion with the one smoothed in the field space with a Gaussian window function.

As additional perturbative degrees of freedom, it is known that magnetic fields of inflationary origin can source curvature perturbations on super-Hubble scales. By requiring the magnetic generated curvature to remain smaller than its inflationary adiabatic counterpart during inflation and reheating, we derive new constraints on the maximal field value today, the reheating energy scale and its equation of state parameter. These bounds end up being stronger by a few order of magnitude than those associated with a possible backreaction of the magnetic field onto the background. Our results are readily applicable to any slow-roll single field inflationary models and any magnetic field having its energy density scaling as a^gamma during inflation. As an illustrative example, massive inflation is found to remain compatible with a magnetic field today Bo = 5 x 10^(-15) G for some values of gamma only if a matter dominated reheating takes place at energies larger than 10^5 GeV. Conversely, assuming gamma=-1, massive inflation followed by a matter dominated reheating cannot explain large scale magnetic fields larger than 10^(-20) G today.

We give a concise formula for the non-Gaussianity of the primordial curvature perturbation generated on super-horizon scales in multi-scalar inflation model without assuming slow-roll conditions. This is an extension of our previous work. Using this formula, we study the generation of non-Gaussianity for the double inflation models in which the slow-roll conditions are temporarily violated after horizon exit, and we show that the non-linear parameter $f_{NL}$ for such models is suppressed by the slow-roll parameters evaluated at the time of horizon exit.

Ultra slow-roll inflation has recently been used to challenge the non-Gaussianity consistency relation. We show that this inflationary scenario belongs to a one parameter class of models and we study its properties and observational predictions. We demonstrate that the power spectrum remains scale-invariant and that the bi-spectrum is of the local type with fnl=5(3-ns)/4 which, indeed, represents a modification of the consistency relation. However, we also show that the system is unstable and suffers from many physical problems among which is the difficulty to correctly WMAP normalize the model. We conclude that ultra slow-roll inflation remains a very peculiar case, the physical relevance of which is probably not sufficient to call into question the validity of the consistency relation.

We provide a general formalism to calculate the infrared correlators of multiple interacting scalar fields in the de Sitter space by means of the stochastic approach. These scalar fields are treated as test fields and hence our result is applicable to the models such as the curvaton scenario where the fields that yield initially isocurvature modes do not contribute to the cosmic energy density during inflationary expansion. The stochastic formalism combined with the argument of conformal invariance of the correlators reflecting the de Sitter isometries allows us to fix the form and amplitude of the three-point functions completely and partially for the four-point functions in terms of calculable quantities. It turns out that naive scaling argument employed in the previous literature does not necessarily hold and we derive the necessary and sufficient condition for the correlator to obey the naive scaling. We also find that correlation functions can in principle exhibit more complicated structure than argued in the literature.

Within the framework of the recently proposed ghost-free massive gravity, a cosmological constant-type self-accelerating solution has been obtained for Minkowski and de Sitter reference metrics. We ease the assumption on the reference metric and find the self-accelerating solution for the reference metric respecting only isotropy, thus considerably extending the range of known solutions.

Chern-Simons gravitational theories are extensions of general relativity in which the parity is violated due to the Chern-Simons term. We study linear perturbations on the static and spherically symmetric background spacetime both for nondynamical and dynamical Chern-Simons theories. We do not make an assumption that the background Chern-Simons scalar field vanishes, which has been adopted in the literature. By eliminating nondynamical variables using their constraint equations, we derive the reduced second order action from which a set of closed evolution equations containing only dynamical variables are immediately obtained and therefore the number of propagating degrees of freedom as well. It is found that ghost is present both for the nondynamical case and for the dynamical case unless the background Chern-Simons scalar field vanishes. It is also found that if the background scalar field vanishes, ghost degrees of freedom are killed and all the modes propagate at the speed of light.

Geometrical cross sections of dust aggregates determine their coupling with disk gas, which governs their motions in protoplanetary disks. Collisional outcomes also depend on geometrical cross sections of initial aggregates. In the previous paper, we performed three-dimensional N-body simulations of sequential collisions of aggregates composed of a number of sub-micron-sized icy particles and examined radii of gyration (and bulk densities) of the obtained aggregates. We showed that collisional compression of aggregates is not efficient and that aggregates remain fluffy. In the present study, we examine geometrical cross sections of the aggregates. Their cross sections decreases due to the compression as well as their gyration radii. It is found that a relation between the cross section and the gyration radius proposed by Okuzumi et al. is valid for the compressed aggregates. We also refine the compression model proposed in our previous paper. The refined model enables us to calculate the evolution of both gyration radii and cross sections of growing aggregates and reproduces well our numerical results of sequential aggregate collisions. The refined model can describe non-equal-mass collisions as well as equal-mass case. Although we do not take into account oblique collisions in the present study, oblique collisions would further hinder compression of aggregates.

We consider the evolution of an initially FLRW universe after the formation of a long, straight, cosmic string with arbitrary tension and mass per unit length. The birth of the string sources scalar and tensor-type perturbations in the background metric and both density and velocity perturbations in the background fluid, which compensate for the string mass and maintain energy conservation. The former generate the deficit angle within the light cone of the string and a gravitational shock front at the cosmological horizon, whereas the latter are confined within the sound cone. We study the properties of the metric within each region of the resulting spacetime and give the explicit coordinate transformations which demonstrate non-violation of causality. This paper generalizes the work of previous studies for the Nambu-Goto string.

There is an observational indication of extragalactic magnetic fields. No known astrophysical process can explain the origin of such large scale magnetic fields, which motivates us to look for their origin in primordial inflation. By solving the linearized Einstein equations, we study metric perturbations sourced by magnetic fields that are produced during inflation. This leads to a simple but robust bound on the inflation models by requiring that the induced metric perturbation should not exceed the observed value 10^-5. In case of the standard single field inflation model, the bound can be converted into a lower bound on the Hubble parameter during inflation.

We explicitly show the fully non-linear equivalence of the $\delta$N and the covariant formalisms for the superhorizon curvature perturbations, which enables us to safely evaluate the non-Gaussian quantities of the curvature perturbation in either formalism. We also discuss isocurvature perturbations in the covariant formalism and clarify the relation between the fully non-linear evolution of the curvature covector and that of the curvature perturbation for multiple interacting fluids.

We study linear perturbations around the static and spherically symmetric spacetime for the gravitational theories whose Lagrangian depends on Ricci scalar and the parity violating Chern-Simons term. By an explicit construction, we show that Hamiltonian for the perturbation variables is not bounded from below in general, suggesting that such a background spacetime is unstable against perturbations. This gives a strong limit on a phenomenological gravitational model which violates parity. We also provide a necessary and sufficient condition for the theory to belong to a special class in which no such instability occurs. For such theories, the number of propagating modes for $\ell \ge 2$ is three, one from the odd and the other two from the even. Unlike in the case of $f(R)$ theories, those modes are coupled each other, which can be used as a distinctive feature to test the parity violating theories from observations. All the modes propagate at the speed of light. No-ghost condition and no-tachyon condition are the same as those in $f(R)$ theories. For the dipole perturbations, the odd and the even modes completely decouple. The odd mode gives a slowly-rotating BH solution whose metric is linearized in its angular momentum. We provide an integral expression of such a solution. On the other hand, the even mode propagates at the speed of light. For the monopole perturbation, in addition to a mode which just shifts the mass of the background BH, there is also one even mode that propagates at the speed of light.

We propose a new variant model of the modulated reheating. If particles have large scale fluctuations on their velocities, or equivalently their Lorentz factors, the decay rate also fluctuates and the curvature perturbation is induced via their decay processes in analogy with modulated reheating. For example, if they are produced nonthermally by the decay of another field with its mass fluctuating on large scales, such a situation is realized. We explicitly calculate the resulting curvature perturbation and non-linearity parameters and show that the effect of velocity-modulation is not negligible if the particles are semi-relativistic at the decay.

If more than one curvaton dominate the Universe at different epochs from each other, curvature perturbations can be temporarily enhanced to a value much larger than the observed one 10^{-5}. The traces of the enhancement may be left as higher order correlation functions, that is, as non-Gaussianity, the stochastic gravitational waves that are sourced by scalar-scalar mode couplings, as well as the primordial black holes that are formed by the gravitational collapse of the enhanced curvature perturbations. We first confirm that such a temporal enhancement indeed occurs by solving the linearized perturbation equations both numerically and analytically. We then derive an analytic expression of the full-order curvature perturbation which does not rely on the frequently used sudden decay approximation and is exact on super-horizon scales. By using this analytic formula, we provide expressions of the non-linearity parameters f_nl, tau_nl and g_nl. If both two curvatons contribute to the final curvature perturbations, the strong non-Gaussianity appears in the trispectrum rather than in the bispectrum. We also find a unique consistency relation between tau_nl and g_nl without f_nl. By using the second-order perturbation theory, we numerically show that the spectrum of the induced gravitational waves has a plateau corresponding to duration of the enhancement and such gravitational waves can be probed by ultimate-DECIGO and space-based atomic interferometers. We finally calculate the abundance of the primordial black holes and put a constraint on the amplitude of the enhanced curvature perturbations.

For the local-type primordial perturbation, it is known that there is an inequality between the bispectrum and the trispectrum. By using the diagrammatic method, we develop a general formalism to systematically construct the similar inequalities up to any order correlation function. As an application, we explicitly derive all the inequalities up to six and eight-point functions.

We derive the equation of matter density perturbations on sub-horizon scales around a flat Friedmann-Lema\^\i tre-Robertson-Walker background for the general Lagrangian density $f(R,\GB)$ that is a function of a Ricci scalar $R$ and a Gauss-Bonnet term $\GB$. We find that the effective gravitational constant generically scales as distance squared at small distances. The effect of this diminishing of the gravitational constant might be important in the gravitational dynamics of cosmic objects such as galaxies, which can be in principle tested by observations. We also provide the general expressions for the effective anisotropic stress, which is useful to constrain modified gravity models from observations of large-scale structure and weak lensing. We also find that there is a special class of theories which evade this unusual behaviour and that the condition to belong to this special class is exactly the same as the one for not having super-luminal modes with propagation speed proportional to their wavenumber.

We study linear metric perturbations around a spherically symmetric static spacetime for general f(R,G) theories, where R is the Ricci scalar and G is the Gauss-Bonnet term. We find that unless the determinant of the Hessian of f(R,G) is zero, even-type perturbations have a ghost for any multi-pole mode. In order for these theories to be plausible alternatives to General Relativity, the theory should satisfy the condition that the ghost is massive enough to effectively decouple from the other fields. We study the requirement on the form of f(R,G) which satisfies this condition. We also classify the number of propagating modes both for the odd-type and the even-type perturbations and derive the propagation speeds for each mode.

We classify models generating large local-type non-Gaussianity into some categories by using some "consistency relations" among the non-linearity parameters f_{NL}^{local}, \tau_{NL}^{local} and g_{NL}^{local}, which characterize the size of bispectrum for the former and trispectrum for the later two. Then we discuss how one can discriminate models of large local-type non-Gaussianity with such relations. We first classify the models by using the ratio of \tau_{NL}^{local}/(6f_{NL}^{local}/5)^2, which is unity for "single-source" models and deviates from unity for "multi-source" ones. We can make a further classification of models in each category by utilizing the relation between f_{NL}^{local} and g_{NL}^{local}. Our classification suggests that observations of trispectrum would be very helpful to distinguish models of large non-Gaussianity and may reveal the generation mechanism of primordial fluctuations.

We show that current cosmic acceleration can be explained by an almost massless scalar field experiencing quantum fluctuations during primordial inflation. Provided its mass does not exceed the Hubble parameter today, this field has been frozen during the cosmological ages to start dominating the universe only recently. By using supernovae data, completed with baryonic acoustic oscillations from galaxy surveys and cosmic microwave background anisotropies, we infer the energy scale of primordial inflation to be around a few TeV, which implies a negligible tensor-to-scalar ratio of the primordial fluctuations. Moreover, our model suggests that inflation lasted for an extremely long period. Dark energy could therefore be a natural consequence of cosmic inflation close to the electroweak energy scale.

In order to classify modified gravity models according to their physical properties, we analyze the cosmological linear perturbations for f(R,G) theories (R being the Ricci scalar and G, the Gauss-Bonnet term) with a minimally coupled perfect fluid. For the scalar type perturbations, we identify in general six degrees of freedom. We find that two of these physical modes obey the same dispersion relation as the one for a non-relativistic de Broglie wave. This means that spacetime is either highly unstable or its fluctuations undergo a scale-dependent super-luminal propagation. Two other modes correspond to the degrees of freedom of the perfect fluid, and propagate with the sound speed of such a fluid. The remaining two modes correspond to the entropy and temperature perturbations of the perfect fluid, and completely decouple from the other modes for a barotropic equation of state. We then provide a concise condition on f(R,G) theories, that both f(R) and R+f(G) do fulfill, to avoid the de Broglie type dispersion relation. For the vector type perturbation, we find that the perturbations decay in time. For the tensor type perturbation, the perturbations can be either super-luminal or sub-luminal, depending on the model. No-ghost conditions are also obtained for each type of perturbation.

We discuss generation of non-Gaussianity in density perturbation through the super-horizon evolution during inflation by using the so-called $\delta N$ formalism. We first provide a general formula for the non-linearity parameter generated during inflation. We find that it is proportional to the slow-roll parameters, multiplied by the model dependent factors that may enhance the non-gaussianity to the observable ranges. Then we discuss three typical examples to illustrate how difficult to generate sizable non-Gaussianity through the super-horizon evolution. First example is the double inflation model, which shows that temporal violation of slow roll conditions is not enough for the generation of non-Gaussianity. Second example is the ordinary hybrid inflation model, which illustrates the importance of taking into account perturbations on small scales. Finally, we discuss Kadota-Stewart model. This model gives an example in which we have to choose rather unnatural initial conditions even if large non-Gaussianity can be generated.

We describe the linear cosmological perturbations of a perfect fluid at the level of an action, providing thus an alternative to the standard approach based only on the equations of motion. This action is suited not only to perfect fluids with a barotropic equation of state, but also to those for which the pressure depends on two thermodynamical variables. By quantizing the system we find that (1) some perturbation fields exhibit a noncommutativity quite analogous to the one observed for a charged particle moving in a strong magnetic field, (2) local curvature and pressure perturbations cannot be measured simultaneously, (3) ghosts appear if the null energy condition is violated.

We provide an analytical expression for the trispectrum of the Cosmic Microwave Background (CMB) temperature anisotropies induced by cosmic strings. Our result is derived for the small angular scales under the assumption that the temperature anisotropy is induced by the Gott-Kaiser-Stebbins effect. The trispectrum is predicted to decay with a non-integer power-law exponent l^(-r) with 6<r<7, depending on the string microstructure, and thus on the string model. For Nambu-Goto strings, this exponent is related to the string mean square velocity and the loop distribution function. We then explore two classes of wavenumber configuration in Fourier space, the kite and trapezium quadrilaterals. The trispectrum can be of any sign and appears to be strongly enhanced for all squeezed quadrilaterals.

The Cosmic Microwave Background (CMB) bispectrum of the temperature anisotropies induced by a network of cosmic strings is derived for small angular scales, under the assumption that the principal cause of temperature fluctuations is the Gott-Kaiser-Stebbins (GKS) effect. We provide analytical expressions for all isosceles triangle configurations in Fourier space. Their overall amplitude is amplified as the inverse cube of the angle and diverges for flat triangles. The isosceles configurations generically lead to a negative bispectrum with a power law decay l^(-6) for large multipole l. However, collapsed triangles are found to be associated with a positive bispectrum whereas the squeezed triangles still exhibit negative values. We then compare our analytical estimates to a direct computation of the bispectrum from a set of 300 statistically independent temperature maps obtained from Nambu-Goto cosmic string simulations in a Friedmann-Lemaitre-Robertson-Walker (FLRW) universe. We find good agreement for the overall amplitude, the power law behaviour and angle dependency of the various triangle configurations. At l~500 the cosmic string GKS effect contributes approximately the same equilateral CMB bispectrum amplitude as an inflationary model with |fNL|~10^3, if the strings contribute about 10% of the temperature power spectrum at l=10. Current bounds on fNL are not derived using cosmic string bispectrum templates, and so our fNL estimate cannot be used to derive bounds on strings. However it does suggest that string bispectrum templates should be included in the search of CMB non-Gaussianities.

We study the propagation of the scalar modes around a Friedmann-Lemaitre-Robertson-Walker universe for general modifications of gravity in the presence of a real scalar field. In general, there will be two propagating scalar perturbation fields, which will have in total four degrees of freedom. Two of these degrees will have a superluminal propagation--with k-dependent speed of propagation--whereas the other two will travel with the speed of light. Therefore, the scalar degrees of freedom do not modify the general feature of modified gravity models: the appearance of modes whose frequency depends on the second power of the modulus of the wave vector. Constraints are given and special cases are discussed.

We have found for the general class of Modified Gravity Models f(R,G) a new instability which can arise in vacuum for the scalar modes of the cosmological perturbations if the background is not de Sitter. In particular, the short-wavelength modes, if stable, in general have a group velocity which depends linearly in k, the wave number. Therefore these modes will be in general superluminal. We have also discussed the condition for which in general these scalar modes will be ghost-like. There is a subclass of these models, defined out of properties of the function f(R,G) and to which the f(R) and f(G) models belong, which however does not have this feature.

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.

We present a new efficient method for computing the non-linearity parameters of the higher order correlation functions of local type curvature perturbations in inflation models having a $\cal N$-component scalar field, focusing on the non-Gaussianity generated during the evolution on super-horizon scales. In contrast to the naive expectation that the number of operations necessary to compute the $n$-point functions is proportional to ${\cal N}^n$, it grows only linearly in ${\cal N}$ in our formalism. Hence, our formalism is particularly powerful for the inflation models composed of a multi-component scalar field. Explicit formulas obtained by applying our method are provided for $n=2,3,4$ and 5, which correspond to power-, bi-, tri- and {\it quad}-spectra, respectively. We also discuss how many parameters we need to parameterize the amplitude and the shape of the higher order correlation functions of local type.

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.

We study non-Gaussianity, the spectral index of primordial scalar fluctuations and tensor modes in models where fluctuations from the inflaton and the curvaton can both contribute to the present cosmic density fluctuations. Even though simple single-field inflation models generate only tiny non-Gaussianity, if we consider such a mixed scenario, large non-Gaussianity can be produced. Furthermore, we study the inflationary parameters such as the spectral index and the tensor-to-scalar ratio in this kind of models and discuss in what cases models predict the spectral index and tensor modes allowed by the current data while generating large non-Gaussianity, which may have many implications for model-buildings of the inflationary universe.

We point out that the results by Brandenberger et al. in Phys.Rev.D77:083502(2008) that the geometry around the straight cosmic strings with stationary junctions is flat to linear order in the string tension can be immediately extended to any order.

We investigate non-Gaussianity in the modulated reheating scenario where fluctuations of the decay rate of the inflaton generate adiabatic perturbations, paying particular attention to the non-linearity parameters $f_{\rm NL}, \tau_{\rm NL}$ and $g_{\rm NL}$ as well as the scalar spectral index and tensor-to-scalar ratio which characterize the nature of the primordial power spectrum. We also take into account the pre-existing adiabatic perturbations produced from the inflaton fluctuations. It has been known that the non-linearity between the curvature perturbations and the fluctuations of the decay rate can yield non-Gaussianity at the level of $f_{\rm NL} \sim \mathcal{O}(1)$, but we find that the non-linearity between the decay rate and the modulus field which determines the decay rate can generate much greater non-Gaussianity. We also discuss a consistency relation among non-linearity parameters which holds in the scenario and find that the modulated reheating yields a different one from that of the curvaton model. In particular, they both can yield a large positive $f_{\rm NL}$ but with a different sign of $g_{\rm NL}$. This provides a possibility to discriminate these two competitive models by looking at the sign of $g_{\rm NL}$. Furthermore, we work on some concrete inflation models and investigate in what cases models predict the spectral index and the tensor-to-scalar ratio allowed by the current data while generating large non-Gaussianity, which may have many implications for model-buildings of the inflationary universe.

We investigate the non-Gaussianity of primordial curvature perturbation in the modulated reheating scenario where the primordial perturbation is generated due to the spacial fluctuation of the inflaton decay rate to radiation. We use the $\delta N$ formalism to evaluate the trispectrum of curvature perturbation as well as its bispectrum. We give expressions for three non-linear parameters $f_{NL}, \tau_{NL}$ and $g_{NL}$ in the modulated reheating scenario. If the intrinsic non-Gaussianity of scalar field fluctuations and third derivative of the decay rate with respect to scalar fields are negligibly small, $g_{NL}$ has at least the same order of magnitude as $f_{NL}$. We also give general inequality between $f_{NL}$ and $\tau_{NL}$ which is true for other inflationary scenarios as long as primordial non-Gaussianity comes from super-horizon evolution.

We analyze the non-Gaussianity for primordial curvature perturbations generated in multi-scalar slow-roll inflation model including the model with non-separable potential by making use of $\delta N$ formalism. Many authors have investigated the possibility of large non-Gaussianity for the models with separable potential, and they have found that the non-linear parameter, $f_{NL}$, is suppressed by the slow-roll parameters. We show that for the non-separable models $f_{NL}$ is given by the product of a factor which is suppressed by the slow-roll parameters and a possible enhancement factor which is given by exponentials of quantities of O(1).

We show that the primordial curvature perturbations may originate not from the quantum fluctuations of inflaton during inflation but from isocurvature perturbations which are amplified during preheating after inflation. We consider a simple preheating model, whose potential is given by $V(\phi, \chi)={1/4}\lambda \phi^4+{1/2}g^2 \phi^2 \chi^2$ with $g^2/\lambda=2$, as a possible realization of generating curvature perturbations during preheating. We make use of the $\delta N$ formalism which requires only knowledge of the homogeneous background solutions in order to evaluate the evolution of curvature perturbations on super-horizon scales. We solve the background equations numerically and find that the amplitude and the spectral index of curvature perturbations originating from preheating can be tuned to the observed values if the isocurvature perturbations at the end of inflation is not suppressed on super-horizon scales. We also point out that the tensor to scalar ratio in ${1/4}\lambda \phi^4$ inflation model can be significantly lowered, hence letting the ${1/4}\lambda \phi^4$ model, which is ruled out by the Wilkinson Microwave Anisotropy Probe (WMAP) data combined with SDSS data, get back into the observationally allowed region.

We present fully non-linear simulations of a self-interacting scalar field in the early universe undergoing tachyonic preheating. We find that density perturbations on sub-horizon scales which are amplified by tachyonic instability maintain long range correlations even during the succeeding parametric resonance, in contrast to the standard models of preheating dominated by parametric resonance. As a result the final spectrum exhibits memory and is not universal in shape. We find that throughout the subsequent era of parametric resonance the equation of state of the universe is almost dust-like, hence the Jeans wavelength is much smaller than the horizon scale. If our 2D simulations are accurate reflections of the situation in 3D, then there are wide regions of parameter space ruled out by over-production of black holes. It is likely however that realistic parameter values, consistent with COBE/WMAP normalisation, are safetly outside this black hole over-production region.

We present exact solutions of the massless Klein-Gordon equation in a spacetime in which an infinite straight cosmic string resides. The first solution represents a plane wave entering perpendicular to the string direction. We also present and analyze a solution with a static point-like source. In the short wavelength limit these solutions approach the results obtained by using the geometrical optics approximation: magnification occurs if the observer lies in front of the string within a strip of angular width $8\pi G\mu$, where $\mu$ is the string tension. We find that when the distance from the observer to the string is less than $ 10^{-3} {(G \mu)}^{-2}\lambda \sim 150 {\rm Mpc} (\lambda/{\rm AU}) (G\mu/10^{-8})^{-2}$, where $\lambda$ is the wave length, the magnification is significantly reduced compared with the estimate based on the geometrical optics due to the diffraction effect. For gravitational waves from neutron star(NS)-NS mergers the several lensing events per year may be detected by DECIGO/BBO.

We consider the scattering of the gravitational waves by the weak gravitational fields of lens objects. We obtain the scattered gravitational waveform by treating the gravitational potential of the lens to first order, i.e. using the Born approximation. We find that the effect of scattering on the waveform is roughly given by the Schwarzschild radius of the lens divided by the wavelength of gravitational wave for a compact lens object. If the lenses are smoothly distributed, the effect of scattering is of the order of the convergence field $\kappa$ along the line of sight to the source. In the short wavelength limit, the amplitude is magnified by $1+\kappa$, which is consistent with the result in weak gravitational lensing.

Wave effects can be important for the gravitational lensing of gravitational waves. In such a case, wave optics must be used in stead of geometric optics. We consider a plane wave entering a lens object and solve numerically the wave equation for three lens models: the uniform density sphere, the singular isothermal sphere, and the Hernquist model. By comparing our numerical solutions with the analytical solutions under the thin lens approximation, we evaluate the error of this approximation. The results show that the relative error of the thin lens approximation is small if the geometrical thickness of the lens is much smaller than the distance between the lens and the observer.

We provide a simple but robust argument that primordial black hole (PBH) production generically does {\em not} exceed astrophysical bounds during the resonant preheating phase after inflation. This conclusion is supported by fully nonlinear lattice simulations of various models in two and three dimensions which include rescattering but neglect metric perturbations. We examine the degree to which preheating amplifies density perturbations at the Hubble scale and show that at the end of the parametric resonance, power spectra are universal, with no memory of the power spectrum at the end of inflation. In addition we show how the probability distribution of density perturbations changes from exponential on very small scales to Gaussian when smoothed over the Hubble scale -- the crucial length for studies of primordial black hole formation -- hence justifying the standard assumption of Gaussianity.