In the full nonlinear cosmological perturbation theory in the leading order of the gradient expansion, all the types of the gauge invariant perturbation variables are defined. The metric junction conditions across the spacelike transition hypersurface are formulated in a manifestly gauge invariant manner. It is manifestly shown that all the physical laws such as the evolution equations, the constraint equations, and the junction conditions can be written using the gauge invariant variables which we defined only. Based on the existence of the universal adiabatic growing mode in the nonlinear perturbation theory and the $\rho$ philosophy where the physical evolution are described using the energy density $\rho$ as the evolution parameter, we give the definitions of the adiabatic perturbation variable and the entropic perturbation variables in the full nonlinear perturbation theory. In order to give the analytic order estimate of the nonlinear parameter $f_{NL}$, we present the exponent evaluation method. As the models where $f_{NL}$ changes continuously and becomes large, using the $\rho$ philosophy, we investigate the non-Gaussianity induced by the entropic perturbation of the component which does not govern the cosmic energy density, and we show that in order to obtain the significant non-Gaussianity it is necessary that the scalar field which supports the entropic perturbation is extremely small compared with the scalar field which supports the adiabatic perturbation.

In the general matter composition where the multiple scalar fields and the multiple perfect fluids coexist, in the leading order of the gradient expansion, we construct all of the solutions of the nonlinear evolutions of the locally homogeneous universe. From the momentum constraint, we derive the constraints which the solution constants of the locally homogeneous universe must satisfy. We construct the gauge invariant perturbation variables in the arbitrarily higher order nonlinear cosmological perturbation theory around the spatially flat Friedmann-Robertson-Walker universe. We construct the nonlinear long wavelength limit formula representing the long wavelength limit of the evolution of the nonlinear gauge invariant perturbation variables in terms of perturbations of the evolutions of the locally homogeneous universe. By using the long wavelength limit formula, we investigate the evolution of nonlinear cosmological perturbations in the universe dominated by the multiple slow rolling scalar fields with an arbitrary potential. The tau function and the N potential introduced in this paper make it possible to write the evolution of the multiple slow rolling scalar fields with an arbitrary interaction potential and the arbitrarily higher order nonlinear Bardeen parameter at the end of the slow rolling phase analytically. It is shown that the nonlinear parameters such as fNL and gNL are suppressed by the slow rolling expansion parameters.

We present the LWL formula which represents the long wavelengh limit of the solutions of evolution equations of cosmological perturbations in terms of the exactly homogeneous solutions in the most general case where multiple scalar fields and multiple perfect fluids coexist. We find the conserved quantity which has origin in the adiabatic decaying mode, and by regarding this quantity as the source term we determine the correction term which corrects the discrepancy between the exactly homogeneous perturbations and the $k \to 0$ limit of the evolutions of cosmological perturbations. This LWL formula is useful for investigating the evolutions of cosmological perturbations in the early stage of our universe such as reheating after inflation and the curvaton decay in the curvaton scenario. When we extract the long wavelength limits of evolutions of cosmological perturbations from the exactly homogeneos perturbations by the LWL formula, it is more convenient to describe the corresponding exactly homogeneous system with not the cosmological time but the scale factor as the evolution parameter. By applying the LWL formula to the reheating model and the curvaton model with multiple scalar fields and multiple radiation fluids, we obtain the S formula representing the final amplitude of the Bardeen parameter in terms of the initial adiabatic and isocurvature perturbations Keywords:cosmological perturbations,long wavelength limit,reheating,curvaton PACS number(s):98.80.Cq