We develop a theory of nonlinear cosmological perturbations on superhorizon scales where a characteristic length scale of perturbations is longer than the Hubble radius, in general theoretical frameworks. Our formalism is based on the spatial gradient expansion approach by adopting the ADM decomposition. Nonlinear superhorizon perturbation including both scalar (curvature perturbation) and tensor (gravitational waves) modes can be dealt with valid up to a second-order in the expansion. First we will review the formalism for a standard general relativity (GR) gravity plus a general kinetic single scalar (k-inflation) with a general form of the potential in the context of inflationary cosmology. That is the basic overview of our procedure. Then it can be extended to more general framework, that is (1) beyond k-inflation (Galileon inflation), (2) a multi-component scalar field with a general kinetic term and a general form of the potential and also (3) beyond Einstein gravity (general scalar-tensor theory), which can lead to several kinds of modified gravity. These theories are motivated not only inflation, but also the topic of dark energy. We provide a formalism to obtain the solution and construct nonlinear curvature perturbation in such general theoretical situation and it can be applied to the calculation of the superhorizon evolution of a primordial non-Gaussianity beyond the so-called $\delta N$ formalism, showing fully nonlinear interaction of both scalar and tensor modes.

We explore the collision of two cylindrical bubbles in classical general relativity with a scalar field stress-energy tensor. Inside each bubble the field rests at a local minimum of the potential with non-negative energy density. Outside the field rests at zero potential, the global minimum. The calculation resolves the connection from the inner de-Sitter region to the asymptotically flat Minkowski spacetime. We choose initial conditions such that the two bubbles collide and study the full nonlinear evolution by means of a two-dimensional numerical simulation of Einstein's equations. The collision generates a strongly interacting region with spatially varying fields and potentials. These circumstances promote dynamical exploration of the potential's landscape. No horizon is present and the scalar curvature invariants eventually diverge. We speculate that Schwarzschild-like horizons will encompass only part of the complicated, interesting regions of spacetime in the analogous case of colliding spherical bubbles.

We argue a scenario motivated by the context of string landscape, where our universe is produced by a new vacuum bubble embedded in an old bubble and these bubble universes have not only different cosmological constants, but also their own different gravitational constants. We study these effects on the primordial curvature perturbations. In order to construct a model of varying gravitational constants, we use the Jordan-Brans-Dicke (JBD) theory where different expectation values of scalar fields produce difference of constants. In this system, we investigate the nucleation of bubble universe and dynamics of the wall separating two spacetimes. In particular, the primordial curvature perturbation on superhorizon scales can be affected by the wall trajectory as the boundary effect. We show the effect of gravitational constant in the exterior bubble universe can provide a peak like a bump feature at a large scale in a modulation of power spectrum.

We focus on the evolution of curvature perturbation on superhorizon scales by adopting the spatial gradient expansion and show that the nonlinear theory, called the beyond $\delta N$-formalism as the next-leading order in the expansion. As one application of our formalism for a single scalar field, we investigate the case of varying sound speed. In our formalism, we can deal with the time evolution in contrast to $\delta N$-formalism, where curvature perturbations remain just constant, and nonlinear curvature perturbation follows the simple master equation whose form is similar as one in linear theory. So the calculation of bispectrum can be done in the next-leading order in the expansion as similar as the case of deriving the power spectrum. We discuss localized features of both primordial power and bispectrum generated by the effect of varying sound speed with a finite duration time. We can see a local feature like a bump in the equilateral bispectrum.

We discuss the possibility that we could obtain some hints of the heavy physics during inflation by analyzing local features of the primordial bispectrum. A heavy scalar field could leave large signatures in the primordial spectra through the parametric resonance between its background oscillation and the fluctuations. Since the duration of the heavy-mode oscillations is finite, the effect of the resonance is localized in momentum space. In this paper, we show that the bispectrum is amplified when such a resonance occurs, and that the peak amplitude of the feature can be O(10^{1-2}), or as large as O(10^5) depending on the type of interactions. In particular, the resonance can give large contributions in finitely squeezed configurations, while the bispectrum cannot be large in the exact squeezed limit. We also find that there is a relation between the scales at which the features appear in the bispectrum and the power spectrum, and that the feature in the bispectrum can be much larger than that in the power spectrum. If correlated features are observed at characteristic scales in the primordial spectra, it will indicate the presence of heavy degrees of freedom. By analyzing these features, we may be able to obtain some information on the physics behind inflation.

We develop a theory of nonlinear cosmological perturbations on superhorizon scales for generic single-field inflation. Our inflaton is described by the Lagrangian of the form $W(X,\phi)-G(X,\phi)\Box\phi$ with $X=-\partial^{\mu}\phi\partial_{\mu}\phi/2$, which is no longer equivalent to a perfect fluid. This model is more general than k-inflation, and is called G-inflation. A general nonlinear solution for the metric and the scalar field is obtained at second order in gradient expansion. We derive a simple master equation governing the large-scale evolution of the nonlinear curvature perturbation. It turns out that the nonlinear evolution equation is deduced as a straightforward extension of the corresponding linear equation for the curvature perturbation on uniform $\phi$ hypersurfaces.

We develop a theory of nonlinear cosmological perturbations on superhorizon scales for a multi-component scalar field with a general kinetic term and a general form of the potential in the context of inflationary cosmology. We employ the ADM formalism and the spatial gradient expansion approach, characterised by O(\epsilon^2), where \epsilon=1/(HL) is a small parameter representing the ratio of the Hubble radius to the characteristic length scale L of perturbations. We provide a formalism to obtain the solution in the multi-field case. This formalism can be applied to the superhorizon evolution of a primordial non-Gaussianity beyond the so-called \delta N formalism which is equivalent to O(\epsilon^0) of the gradient expansion. In doing so, we also derive fully nonlinear gauge transformation rules valid through O(\epsilon^2). These fully nonlinear gauge transformation rules can be used to derive the solution in a desired gauge from the one in a gauge where computations are much simpler. As a demonstration, we consider an analytically solvable model and construct the solution explicitly.

We investigate the possibility that a heavy scalar field, whose mass exceeds the Hubble scale during inflation, could leave non-negligible signatures in the Cosmic Microwave Background (CMB) temperature anisotropy power spectrum through the parametric resonance between its background oscillations and the inflaton fluctuations. By assuming the heavy scalar field couples with the inflaton derivatively, we show that the resonance can be efficient without spoiling the slow-roll inflation. The primordial power spectrum modulated by the resonance has a sharp peak at a specific scale and could be an origin of the anomalies observed in the angular power spectrum of the CMB. In some values of parameters, the modulated spectrum can fit the observed data better than the simple power-law power spectrum, though the resultant improvement of the fit is not large enough and hence other observations such as non-Gaussianity are necessary to confirm that the CMB anomalies are originated from the resonant effect of the heavy scalar field. The resonant signatures can provide an opportunity to observe heavy degrees of freedom during inflation and improve our understanding of physics behind inflation.

We study the effects of sudden change in the sound velocity on primordial curvature perturbation spectrum in inflationary cosmology, assuming that the background evolution satisfies the slow-roll condition throughout. It is found that the power spectrum acquires oscillating features which are determined by the ratio of the sound speed before and after the transition and the wavenumeber which crosses the sound horizon at the transition, and their analytic expression is given. In some values of those parameters, the oscillating primordial power spectrum can better fit the observed Cosmic Microwave Background temperature anisotropy power spectrum than the simple power-law power spectrum, although introduction of such a new degree of freedom is not justified in the context of Akaike's Information Criterion.

We clarify the behavior of curvature perturbations in a nonlinear theory in case the inflaton temporarily stops during inflation. We focus on the evolution of curvature perturbation on superhorizon scales by adopting the spatial gradient expansion and show that the nonlinear theory, called the {\it beyond} $\delta N$-formalism for a general single scalar field as the next-leading order in the expansion. Both the leading-order in the expansion ($\delta N$-formalism) and our nonlinear theory include the solutions of full-nonlinear orders in the standard perturbative expansion. Additionally, in our formalism, we can deal with the time evolution in contrast to $\delta N$-formalism, where curvature perturbations remain just constant, and show decaying modes do not couple with growing modes as similar to the case with linear theory. We can conclude that although the decaying mode diverges when $\dot{\phi}$ vanishes, there appears no trouble for both the linear and nonlinear theory since these modes will vanish at late times.

We develop a theory of nonlinear cosmological perturbations on superhorizon scales for a single scalar field with a general kinetic term and a general form of the potential. We employ the ADM formalism and the spatial gradient expansion approach, characterised by $O(\epsilon^m)$, where $\epsilon=1/(HL)$ is a small parameter representing the ratio of the Hubble radius to the characteristic length scale $L$ of perturbations. We obtain the general solution for a full nonlinear version of the curvature perturbation valid up through second-order in $\epsilon$ ($m=2$). We find the solution satisfies a nonlinear second-order differential equation as an extension of the equation for the linear curvature perturbation on the comoving hypersurface. Then we formulate a general method to match a perturbative solution accurate to $n$-th-order in perturbation inside the horizon to our nonlinear solution accurate to second-order ($m=2$) in the gradient expansion on scales slightly greater than the Hubble radius. The formalism developed in this paper allows us to calculate the superhorizon evolution of a primordial non-Gaussianity beyond the so-called $\delta N$ formalism or separate universe approach which is equivalent to leading order ($m=0$) in the gradient expansion. In particular, it can deal with the case when there is a temporary violation of slow-roll conditions. As an application of our formalism, we consider Starobinsky's model, which is a single field model having a temporary non-slow-roll stage due to a sharp change in the potential slope. We find that a large non-Gaussianity can be generated even on superhorizon scales due to this temporary suspension of slow-roll inflation.

We develop a theory of non-linear cosmological perturbations at superhorizon scales for a scalar field with a Lagrangian of the form $P(X,\phi)$, where $X=-\partial^{\mu}\phi\partial_{\mu}\phi$ and $\phi$ is the scalar field. We employ the ADM formalism and the spatial gradient expansion approach to obtain general solutions valid up to the second order in the gradient expansion. This formulation can be applied to, for example, DBI inflation models to investigate superhorizon evolution of non-Gaussianities. With slight modification, we also obtain general solutions valid up to the same order for a perfect fluid with a general equation of state $P=P(\rho)$.

We study the dynamics of colliding domain walls including self-gravity. The initial data is set up by applying a BPS domain wall in five-dimensional supergravity, and we evolve the system determining the final outcome of collisions. After a collision, a spacelike curvature singularity covered by a horizon is formed in the bulk, resulting in a black brane with trapped domain walls. This is a generic consequence of collisions, except for non-relativistic weak field cases, in which the walls pass through one another or multiple bounces take place without singularity formation. These results show that incorporating the self-gravity drastically changes a naive picture of colliding branes.

We study the behaviour of five-dimensional fermions localized on branes, which we describe by domain walls, when two parallel branes collide in a five-dimensional Minkowski background spacetime. We find that most fermions are localized on both branes as a whole even after collision. However, how much fermions are localized on which brane depends sensitively on the incident velocity and the coupling constants unless the fermions exist on both branes.

The thermal equilibrium of string gas is necessary to activate the Brandenberger-Vafa mechanism, which makes our observed 4-dimensional universe enlarge. Nevertheless, the thermal equilibrium is not realized in the original setup, a problem that remains as a critical defect. We study thermal equilibrium in the Hagedorn universe, and explore possibilities for avoiding the issue aforementioned flaw. We employ a minimal modification of the original setup, introducing a dilaton potential. Two types of potential are investigated: exponential and double-well potentials. For the first type, the basic evolutions of universe and dilaton are such that both the radius of the universe and the dilaton asymptotically grow in over a short time, or that the radius converges to a constant value while the dilaton rolls down toward the weak coupling limit. For the second type, in addition to the above solutions, there is another solution in which the dilaton is stabilized at a minimum of potential and the radius grows in proportion to $t$. Thermal equilibrium is realized for both cases during the initial phase. These simple setups provide possible resolutions of the difficulty.

We study collision of two domain walls in 5-dimensional asymptotically Anti de Sitter spacetime. This may provide the reheating mechanism of an ekpyrotic (or cyclic) brane universe, in which two BPS branes collide and evolve into a hot big bang universe. We evaluate a change of scalar field making the domain wall and can investigate the effect of a negative cosmological term in the bulk to the collision process and the evolution of our universe.

We study a particle production at the collision of two domain walls in 5-dimensional Minkowski spacetime. This may provide the reheating mechanism of an ekpyrotic (or cyclic) brane universe, in which two BPS branes collide and evolve into a hot big bang universe. We evaluate a production rate of particles confined to the domain wall. The energy density of created particles is given as $\rho \approx 20 \bar{g}^4 N_b ~m_\eta^4 $ where $\bar{g}$ is a coupling constant of particles to a domain-wall scalar field, $N_b $ is the number of bounces at the collision and $m_\eta$ is a fundamental mass scale of the domain wall. It does not depend on the width $d$ of the domain wall, although the typical energy scale of created particles is given by $\omega\sim 1/d$. The reheating temperature is evaluated as $T_{\rm R}\approx 0.88 ~ \bar{g} ~ N_b^{1/4}$. In order to have the baryogenesis at the electro-weak energy scale, the fundamental mass scale is constrained as $m_\eta \gsim 1.1\times 10^7$ GeV for $\bar{g}\sim 10^{-5}$.