We study the effect of primordial scalar curvature perturbations on the propagation of gravitational waves over cosmic distances. We point out that such curvature perturbations deform the isotropic spectrum of any stochastic background of gravitational waves of primordial origin through the (integrated) Sachs-Wolfe effect. Computing the changes in the amplitude and frequency of the propagating gravitational wave induced at linear order by scalar curvature perturbations, we show that the resulting deformation of each frequency bin of the gravitational wave spectrum is described by a linearly biased Gaussian with the variance $\sigma^2 \simeq \int d\ln k \Delta_{\mathcal R}^2$, where $\Delta_{\mathcal R}^2(k)$ denotes the amplitude of the primordial curvature perturbations. The linear bias encodes the correlations between the changes induced in the frequency and amplitude of the gravitational waves. Taking into account the latest bounds on $\Delta_{\mathcal R}^2$ from primordial black hole and gravitational wave searches, we demonstrate that the resulting ${\mathcal O}(\sigma)$ deformation can be significant for extremely peaked gravitational wave spectra. We further provide an order of magnitude estimate for broad spectra, for which the net distortion is ${\mathcal O}(\sigma^2)$.

We point out that there is a high-frequency tail of the stochastic inflationary gravitational wave background that scales as $f^{-1/2}$ with frequency $f$. This contribution comes from the graviton vacuum fluctuation amplified by the inflaton coherent oscillation during the reheating stage. It contains information on inflaton properties such as the inflaton mass as well as the thermal history of the early Universe.

We investigate the sensitivity of Higgs(-like) inflation to higher dimensional operators in the nonminimal couplings and in the potential, both in the metric and Palatini formalisms. We find that, while inflationary predictions are relatively stable against the higher dimensional operators around the attractor point in the metric formalism, they are extremely sensitive in the Palatini one: for the latter, inflationary predictions are spoiled by $|\xi_4| \gtrsim 10^{-6}$ in the nonminimal couplings $(\xi_2 \phi^2 + \xi_4 \phi^4 + \cdots)R$, or by $|\lambda_6| \gtrsim 10^{-16}$ in the Jordan-frame potential $\lambda_4 \phi^4 + \lambda_6 \phi^6 + \cdots$ (both in Planck units). This extreme sensitivity results from the absence of attractor in the Palatini formalism. Our study underscores the challenge of realizing inflationary models with the nonminimal coupling in the Palatini formalism.

It has been argued that the mixed Higgs-$R^2$ model acts as the UV extension of the Higgs inflation, pushing up its cut-off scale in the vacuum close up to the Planck scale. In this letter, we study the inflaton oscillation stage after inflation, focusing on the effective mass of the phase direction of the Higgs field, which can cause a violent preheating process. We find that the "spikes" in the effective mass of the phase direction observed in the Higgs inflation still appear in the mixed Higgs-$R^2$ model. While the spikes appear above the cut-off scale in the Higgs-only case, they appear below the cut-off scale when the model is extended with $R^2$ term though reheating cannot be completed in the violent particle production regime since the spikes get milder.

Galactic-scale structure is of particular interest since it provides important clues to dark matter properties and its observation is improving. Weakly interacting massive particles (WIMPs) behave as cold dark matter on galactic scales, while beyond-WIMP candidates suppress galactic-scale structure formation. Suppression in the linear matter power spectrum has been conventionally characterized by a single parameter, the thermal warm dark matter mass. On the other hand, the shape of suppression depends on the underlying mechanism. It is necessary to introduce multiple parameters to cover a wide range of beyond-WIMP models. Once multiple parameters are introduced, it becomes harder to share results from one side to the other. In this work, we propose adopting neural network technique to facilitate the communication between the two sides. To demonstrate how to work out in a concrete manner, we consider a simplified model of light feebly interacting massive particles.

We study scalar bubble collisions in first-order phase transitions focusing on the relativistic limit. We propose 'trapping equation' which describes the wall behavior after collision, and test it with numerical simulations in several setups. We also examine the energy dynamics after collision and discuss its implications to gravitational wave production.

We study the gravitational wave (GW) production in extremely strong first order phase transitions where the latent heat density dominates the plasma energy density, $\alpha \gtrsim 1$. In such transitions, bubbles develop extremely thin and relativistic fluid configurations, resulting in strong shock waves after collisions. We first propose a strategy to understand the GW production in such a system by separating the problem into the propagation part and the collision part. Focusing on the former, we next develop an effective theory for the propagation of the relativistic fluid shells. Using this effective theory, we finally calculate the expected duration of the relativistic fluid configurations and discuss its implications to the GW production.

We study gravitational-wave production from bubble dynamics (bubble collisions and sound waves) during a cosmic first-order phase transition with an analytic approach. We first propose modeling the system with the thin-wall approximation but without the envelope approximation often adopted in the literature, in order to take bubble propagation after collisions into account. The bubble walls in our setup are considered as modeling the scalar field configuration and/or the bulk motion of the fluid. We next write down analytic expressions for the gravitational-wave spectrum, and evaluate them with numerical methods. It is found that, in the long-lasting limit of the collided bubble walls, the spectrum grows from $\propto f^3$ to $\propto f^1$ in low frequencies, showing a significant enhancement compared to the one with the envelope approximation. It is also found that the spectrum saturates in the same limit, indicating a decrease in the correlation of the energy-momentum tensor at late times. We also discuss the implications of our results to gravitational-wave production both from bubble collisions (scalar dynamics) and sound waves (fluid dynamics).

We investigate the sensitivity of future space-based interferometers such as LISA and DECIGO to the parameters of new particle physics models which drive a first-order phase transition in the early Universe. We first perform a Fisher matrix analysis on the quantities characterizing the gravitational wave spectrum resulting from the phase transition, such as the peak frequency and amplitude. We next perform a Fisher analysis for the quantities which determine the properties of the phase transition, such as the latent heat and the time dependence of the bubble nucleation rate. Since these quantities are determined by the model parameters of the new physics, we can estimate the expected sensitivities to such parameters. We illustrate this point by taking three new physics models for example: (1) models with additional isospin singlet scalars (2) a model with an extra real Higgs singlet, and (3) a classically conformal $B-L$ model. We find that future gravitational wave observations play complementary roles to future collider experiments in pinning down the parameters of new physics models driving a first-order phase transition.

We propose a realization of cosmic inflation with the Higgs field when the Higgs potential has degenerate vacua by employing the recently proposed idea of hillclimbing inflation. The resultant inflationary predictions exhibit a sizable deviation from those of the ordinary Higgs inflation.

We study gravitational-wave production from bubble collisions in a cosmic first-order phase transition, focusing on the possibility of model separation by the bubble nucleation rate dependence of the resulting gravitational-wave spectrum. By using the method of relating the spectrum with the two-point correlator of the energy-momentum tensor $\left< T(x)T(y) \right>$, we first write down analytic expressions for the spectrum with a Gaussian correction to the commonly used nucleation rate, $\Gamma \propto e^{\beta t}\rightarrow e^{\beta t-\gamma^2t^2}$, under the thin-wall and envelope approximations. Then we quantitatively investigate how the spectrum changes with the size of the Gaussian correction. It is found that the spectral shape shows ${\mathcal O}(10)\%$ deviation from $\Gamma \propto e^{\beta t}$ case for some physically motivated scenarios. We also briefly discuss detector sensitivities required to distinguish different spectral shapes.

We study particle production at the preheating era in inflation models with nonminimal coupling $\xi \phi^2R$ and quartic potential $\lambda \phi^4/4$ for several cases: real scalar inflaton, complex scalar inflaton and Abelian Higgs inflaton. We point out that the preheating proceeds much more violently than previously thought. If the inflaton is a complex scalar, the phase degree of freedom is violently produced at the first stage of preheating. If the inflaton is a Higgs field, the longitudinal gauge boson production is similarly violent. This is caused by a spike-like feature in the time dependence of the inflaton field, which may be understood as a consequence of the short time scale during which the effective potential or kinetic term changes suddenly. The produced particles typically have very high momenta $k \lesssim \sqrt{\lambda}M_\text{P}$. The production might be so strong that almost all the energy of the inflaton is carried away within one oscillation for $\xi^2\lambda \gtrsim {\mathcal O}(100)$. This may partly change the conventional understandings of the (p)reheating after inflation with the nonminimal coupling to gravity such as Higgs inflation. We also discuss the possibility of unitarity violation at the preheating stage.

We consider gravitational wave production by bubble collisions during a cosmological first-order phase transition. In the literature, such spectra have been estimated by simulating the bubble dynamics, under so-called thin-wall and envelope approximations in a flat background metric. However, we show that, within these assumptions, the gravitational wave spectrum can be estimated in an analytic way. Our estimation is based on the observation that the two-point correlator of the energy-momentum tensor $\langle T(x)T(y)\rangle$ can be expressed analytically under these assumptions. Though the final expressions for the spectrum contain a few integrations that cannot be calculated explicitly, we can easily estimate it numerically. As a result, it is found that the most of the contributions to the spectrum come from single-bubble contribution to the correlator, and in addition the fall-off of the spectrum at high frequencies is found to be proportional to $f^{-1}$. We also provide fitting formulae for the spectrum.

We study the cosmological history of the classical conformal $B-L$ gauge extension of the standard model, in which the physical scales are generated via the Coleman-Weinberg-type symmetry breaking. Especially, we consider the thermal phase transition of the U$(1)_{B-L}$ symmetry in the early universe and resulting gravitational-wave production. Due to the classical conformal invariance, the phase transition tends to be a first-order one with ultra-supercooling, which enhances the strength of the produced gravitational waves. We show that, requiring (1) U$(1)_{B-L}$ is broken after the reheating, (2) the $B-L$ gauge coupling does not blow up below the Planck scale, (3) the thermal phase transition completes in almost all the patches in the universe, the gravitational wave spectrum can be as large as $\Omega_{\rm GW} \sim 10^{-8}$ at the frequency $f \sim 0.01$-$1$Hz for some model parameters, and a vast parameter region can be tested by future interferometer experiments such as eLISA, LISA, BBO and DECIGO.

We study production of light particles due to oscillation of the Hubble parameter or the scale factor. Any coherently oscillating scalar field, irrespective of its energy fraction in the universe, imprints such an oscillating feature on them. Not only the Einstein gravity but extended gravity models, such as models with non-minimal (derivative) coupling to gravity and $f(R)$ gravity, lead to oscillation of the scale factor. We present a convenient way to estimate the gravitational particle production rate in these circumstances. Cosmological implications of gravitational particle production, such as dark matter/radiation and moduli problem, are discussed. For example, if the theory is described solely by the standard model plus the Peccei-Quinn sector, the Starobinsky $R^2$ inflation may lead to observable amount of axion dark radiation.

In a wide class of new physics models, there exist scalar fields that obtain vacuum expectation values of high energy scales. We study the possibility that the standard model Higgs field has experienced first order phase transition at the high energy scale due to the couplings with these scalar fields. We estimate the amount of gravitational waves produced by the phase transition, and discuss observational consequences.

We study cosmological evolution after inflation in models with non-minimal derivative coupling to gravity. The background dynamics is solved and particle production associated with rapidly oscillating Hubble parameter is studied in detail. In addition, production of gravitons through the non-minimal derivative coupling with the inflaton is studied. We also find that the sound speed squared of the scalar perturbation oscillates between positive and negative values when the non-minimal derivative coupling dominates over the minimal kinetic term. This may lead to an instability of this model. We point out that the particle production rates are the same as those in the Einstein gravity with the minimal kinetic term, if we require the sound speed squared is positive definite.

We point out that the inflaton inevitably couples to all non-conformally coupled matters gravitationally through an oscillation in the Hubble parameter or the cosmic scale factor. It leads to particle production during the inflaton oscillation regime, which is most efficient just after inflation. Moreover, the analysis is extended to the model with non-minimal inflaton couplings to gravity, in which the Hubble parameter oscillates more violently. We apply our results to the graviton production by the inflaton: gravitons are also produced just after inflation, but the non-minimal coupling does not induce inflaton decay into the graviton pair.

We consider background dynamics of generalized Galileon theories in the context of inflation, where gravity and inflaton are non-minimally coupled to each other. In the inflaton oscillation regime, the Hubble parameter and energy density oscillate violently in many cases, in contrast to the Einstein gravity with minimally coupled inflaton. However, we find that there is an adiabatic invariant in the inflaton oscillation regime in any generalized Galileon theory. This adiabatic invariant is useful in estimating the expansion law of the universe and also the particle production rate due to the oscillation of the Hubble parameter.

Motivated by recent progress in our understanding of the $B$-mode polarization of cosmic microwave background (CMB), which provides important information about the inflationary gravitational waves (IGWs), we study the possibility to acquire information about the early universe using future space-based gravitational wave (GW) detectors. We perform a detailed statistical analysis to estimate how well we can determine the reheating temperature after inflation as well as the amplitude, the tensor spectral index, and the running of the inflationary gravitational waves. We discuss how the accuracies depend on noise parameters of the detector and the minimum frequency available in the analysis. Implication of such a study on the test of inflation models is also discussed.

We study cosmological scenarios in which high-energy neutrinos are emitted from the decay of long-lived massive particles at the cosmic time later than a redshift of 10^6. The high-energy neutrino events recently observed by the IceCube experiment suggest a new source of high-energy cosmic-ray neutrinos; decay of a heavy particle can be one of the possibilities. We calculate the spectrum of the high-energy neutrinos emitted from the decay of long-lived particles, taking account of the neutrino scattering processes with background neutrinos. Then, we derive bounds on the scenario using the observation of high-energy cosmic-ray neutrino flux. We also study constraints from the spectral distortions of the cosmic microwave background and the big-bang nucleosynthesis. In addition, we show that the PeV neutrinos observed by the IceCube experiment can originate from the decay of a massive particle with its mass as large as O(10^10 GeV).

Motivated by the recent IceCube result, we study high energy cosmic-ray neutrino flux from the decay of a long-lived particle. Because neutrinos are so transparent, high energy neutrinos produced in the past may also contribute to the present neutrino flux. We point out that the PeV neutrino events observed by IceCube may originate in the decay of a particle much heavier than PeV if its lifetime is shorter than the present cosmic time. It is shown that the mass of the particle responsible for the IceCube event can be as large as $\sim 10^{10}\ {\rm GeV}$. We also discuss several possibilities to acquire information about the lifetime of the long-lived particle.

We study the expansion law of the universe dominated by the oscillating scalar field with non-minimal derivative coupling to gravity as G^{\mu \nu} \partial_{\mu} \phi \partial_{\nu} \phi. In this system the Hubble parameter oscillates with a frequency of the effective mass of the scalar field, which formerly caused a difficulty in analyzing how the universe expands. We find an analytical solution for power law potentials and interpret the solution in an intuitive way by using a new invariant of the system. As a result, we find marginally accelerated expansion for the quadratic potential and no accelerated expansion for the potential with higher power.

We study the effects of various phenomena which may have happened in the early universe on the spectrum of inflationary gravitational waves. The phenomena include phase transitions, entropy productions from non-relativistic matter, the production of dark radiation, and decoupling of dark matter/radiation from thermal bath. These events can create several characteristic signatures in the inflationary gravitational wave spectrum, which may be direct probes of the history of the early universe and the nature of high-energy physics.

Recent cosmological observations indicate the existence of extra light species, i.e., dark radiation. In this paper we show that signatures of the dark radiation are imprinted in the spectrum of inflationary gravitational waves. If the dark radiation is produced by the decay of a massive particle, high frequency mode of the gravitational waves are suppressed. In addition, due to the effect of the anisotropic stress caused by the dark radiation, a dip in the gravitational wave spectrum may show up at the frequency which enters the horizon at the time of the dark radiation production. Once the gravitational wave spectrum is experimentally studied in detail, we can infer the information on how and when the dark radiation was produced in the Universe.

We discuss the effects of cosmic phase transition on the spectrum of primordial gravitational waves generated during inflation. The energy density of the scalar condensation responsible for the phase transition may become sizable at the epoch of phase transition, which significantly affects the evolution of the universe. As a result, the amplitudes of the gravitational waves at high frequency modes are suppressed. Thus the gravitational wave spectrum can be a probe of phase transition in the early universe.