We investigate the squeezing of primordial gravitational waves (PGWs) in terms of quantum discord. We construct a classical state of PGWs without quantum discord and compare it with the Bunch-Davies vacuum. Then it is shown that the oscillatory behavior of the angular-power spectrum of the cosmic microwave background (CMB) fluctuations induced by PGWs can be the signature of the quantum discord of PGWs. In addition, we discuss the effect of quantum decoherence on the entanglement and the quantum discord of PGWs for super-horizon modes. For the state of PGWs with decoherence effect, we examine the decoherence condition and the correlation condition introduced by C. Kiefer et al. (Class. Quantum Grav. 24 (2007) 1699). We show that the decoherence condition is not sufficient for the separability of PGWs and the correlation condition implies that the PGWs in the matter-dominated era have quantum discord.

The superradiant scattering of Alfv\'en waves (Alfv\'enic superradiance) in a forcefree magnetosphere is discussed to reveal the relationship between the Blandford-Znajek (BZ) process and superradiance. For simplicity, we consider a four-dimensional rotating black string spacetime of which each $z=\text{const}$ slice is a Ba\~nados-Teitelboim-Zanelli solution as an analogy of the equatorial plane of the Kerr spacetime. Then, it is confirmed that the condition for Alfv\'enic superradiance coincides with that for the BZ process, and the wave amplification can be very large due to a resonant scattering for some parameter sets of the wave frequency and the angular velocity of the magnetic field line. Moreover, by analysis of the Poynting flux, we first show that the BZ process can be interpreted as the long wavelength limit of Alfv\'enic superradiance.

The so-called black hole shadow is a dark region which is expected to appear in a fine image of optical observation of black holes. It is essentially an absorption cross section of black hole, and the boundary of shadow is determined by unstable circular orbits of photons (UCOP). If there exists a compact object possessing UCOP but no black hole horizon, it can provide us with the same shadow image with black holes, and a detection of shadow image cannot be a direct evidence of black hole existence. Then, this paper examine whether or not such compact objects can exist under some suitable conditions. We investigate thoroughly the static spherical polytropic ball of perfect fluid with single polytrope index, and then investigate a representative example of the piecewise polytropic ball. Our result is that the spherical polytropic ball which we have investigated cannot possess UCOP, if the sound speed at center is subluminal (slower-than-light). This means that, if the polytrope treated in this paper is a good model of stellar matter in compact objects, the detection of shadow image is regarded as a good evidence of black hole existence. As a by-product, we have found the upper bound of the mass-to-radius radio (M/R) of polytropic ball with single index, M/R < 0.281, under the subluminal-sound-speed condition.

The ratio of total mass $M$ to surface radius $R$ of spherical perfect fluid ball has an upper bound, $M/R < B$. Buchdahl obtained $B = 4/9$ under the assumptions; non-increasing mass density in outward direction, and barotropic equation of states. Barraco and Hamity decreased the Buchdahl's bound to a lower value $B = 3/8$ $(< 4/9)$ by adding the dominant energy condition to Buchdahl's assumptions. In this paper, we further decrease the Barraco-Hamity's bound to $B \simeq 0.3636403$ $(< 3/8)$ by adding the subluminal (slower-than-light) condition of sound speed. In our analysis, we solve numerically Tolman-Oppenheimer-Volkoff equations, and the mass-to-radius ratio is maximized by variation of mass, radius and pressure inside the fluid ball as functions of mass density.

We study scattering of waves by black holes. Solving a massless scalar field with a point source in the Schwarzschild spacetime, waves scattered by the black hole is obtained numerically. We then reconstruct images of the black hole from scattered wave data for specified scattering angles. For the forward and the backward directions, obtained wave optical images of black holes show rings that correspond to the black hole glories associated with existence of the unstable circular photon orbit in the Schwarzschild spacetime.

Void models provide a possible explanation of the "accelerated expansion" of the Universe without dark energy. To make the conventional void models more realistic, we allow the void, an underdense region around us, to be anisotropic and consider an average of the distance-redshift relations over the solid angle subtended at the observer. We first show that after taking the average of a form of the optical scalar equation (distance equation), the effective distance equation we obtain coincides with the one for the Lemaitre-Tolman-Bondi universe with a Dyer-Roeder-like extension. We then numerically solve the equation to compare with observational data of Type Ia supernovae. We find that anisotropy allows smaller size of void and larger Omega_m.

We investigate quantum entanglement of a scalar field in the inflationary universe. By introducing a bipartite system using a lattice model of scalar field, we apply the separability criterion based on the partial transpose operation and numerically calculate the bipartite entanglement between separate spatial regions. We find that the initial entangled state becomes separable or dis-entangled when the size of the spatial regions exceed the Hubble horizon. This is necessary condition for the appearance of classicality of the quantum fluctuation. We further investigate the condition for the appearance of the classical distribution function and find that the condition is given by the inequality for the symplectic eigenvalue of the covariance matrix of the scalar field.

Motivated by the inverse problem for the Lemaitre-Tolman-Bondi dust solution, in which problem the luminosity distance function $D_L(z)$ is taken as an input to select a specific model, we compute the function $D_L(z)$ of the LTB solution up to the third order of $z$. To perform the otherwise cumbersome computation, we introduce a new convenient form of the LTB solution, in which the solution is explicit and unified. With this form of the LTB solution we obtain the luminosity distance function with the full generality. We in particular find that the function exactly coincides with that of a homogeneous and isotropic dust solution up to second order, if we demand that the solution be regular at the center.

We investigate the evolution of the non-linear long wavelength fluctuations during preheating after inflation. By using the separate universe approach, the temporal evolution of the power spectrum of the scalar fields and the curvature variable is obtained numerically. We found that the amplitude of the large scale fluctuations is suppressed after non-linear evolution during preheating.

We present a model of an inhomogeneous universe that leads to accelerated expansion after taking spatial averaging. The model universe is the Tolman-Bondi solution of the Einstein equation and contains both a region with positive spatial curvature and a region with negative spatial curvature. We find that after the region with positive spatial curvature begins to re-collapse, the deceleration parameter of the spatially averaged universe becomes negative and the averaged universe starts accelerated expansion. We also discuss the generality of the condition for accelerated expansion of the spatially averaged universe.

We investigate the backreaction of cosmological long wavelength perturbations on the evolution of the Universe. By applying the renormalization group method to a Friedmann-Robertson-Walker universe with long wavelength fluctuations, we demonstrate that the renormalized solution with the backreaction effect is equivalent to that of the separate universe. Then, using the effective Friedmann equation, we show that only non-adiabatic mode of long wavelength fluctuations affects the expansion law of the spatially averaged universe.

We investigate the back reaction of cosmological perturbations on an inflationary universe using the renormalization-group method. The second-order zero mode solution which appears by the nonlinearity of the Einstein equation is regarded as a secular term of a perturbative expansion, we renormalized a constant of integration contained in the background solution and absorbed the secular term to this constant in a gauge-invariant manner. The resultant renormalization-group equation describes the back reaction effect of inhomogeneity on the background universe. For scalar type classical perturbation, by solving the renormalization-group equation, we find that the back reaction of the long wavelength fluctuation works as a positive spatial curvature, and the short wavelength fluctuation works as a radiation fluid. For the long wavelength quantum fluctuation, the effect of back reaction is equivalent to a negative spatial curvature.

We investigated the back reaction of cosmological perturbations on the evolution of the universe using the second order perturbation of the Einstein's equation. To incorporate the back reaction effect due to the inhomogeneity into the framework of the cosmological perturbation, we used the renormalization group method. The second order zero mode solution which appears by the non-linearities of the Einstein's equation is regarded as a secular term of the perturbative expansion, we renormalized a constant of integration contained in the background solution and absorbed the secular term to this constant. For a dust dominated universe, using the second order gauge invariant quantity, we derived the renormalization group equation which determines the effective dynamics of the Friedman-Robertson-Walker universe with the back reaction effect in a gauge invariant manner. We obtained the solution of the renormalization group equation and found that perturbations of the scalar mode and the long wavelength tensor mode works as positive spatial curvature, and the short wavelength tensor mode as radiation fluid.

Using the renormalization group method, we improved the first order solution of the long-wavelength expansion of the Einstein equation. By assuming that the renormalization group transformation has the property of Lie group, we can regularize the secular divergence caused by the spatial gradient terms and absorb it to the background seed metric. The solution of the renormalization group equation shows that the renormalized metric describes the behavior of gravitational collapse in the expanding universe qualitatively well.

We investigate gauge invariant cosmological perturbations in a spatially flat Friedman-Robertson-Walker universe with scalar fields. It is well known that the evolution equation for the gauge invariant quantities has exact solutions in the long-wavelength limit. We find that these gauge invariant solutions can be obtained by differentiating the background solution with respect to parameters contained in the background system. This method is very useful when we analyze the long-wavelength behavior of cosmological perturbation with multiple scalar fields.

We investigate the cosmological perturbation of two-scalar field model during the reheating phase after inflation. Using the exact solution of the perturbation in long-wavelength limit, which is expressed in terms of the background quantities, we analyze the behavior of the metric perturbation. The oscillating inflaton field gives rise to the parametric resonance of the massless scalar field and this leads to the amplification of the iso-curvature mode of the metric perturbations.

The evolution of the cosmological perturbation during the oscillatory stage of the scalar field is investigated. For the power law potential of the inflaton field, the evolution equation of the Mukhanov's gauge invariant variable is reduced to the Mathieu equation and the density perturbation grows by the parametric resonance.

As an application of Gradient Expansion (long-wavelength) approximation, we studied the inhomogeneous universe including the gravitational wave(GW). For a plane-symmetric cosmological model, we could implement the 2nd order expansion and analyze the purely non-linear perturbation of GW. The result indicates that the non-linear effect plays an important role in anisotropic stage of the universe. We also confirm that the perturbation from the linear source term in gradient expansion is consistent with the result of linear perturbation.

Using the long wave perturbation scheme(gradient expansion), the effect of inhomogeneity on the inflationary phase is investigated. We solved the perturbation equation of which source term comes from inhomogeneity of a scalar field and a seed metric. The result indicates that sub-horizon scale inhomogeneity strongly affects the onset of inflation.

We investigate the evolution of the apparent horizons in a numerically gererated worm hole spacetime. The behavior of the apparent horizons is affected by the dynamics of the matter field. By using the local mass of the system, we interpret the evolution of the worm hole structure. Figures are available by mail to author.

We investigate the global structure of the space time with a spherically symmetric inhomogeneity using a metric junction, and classify all possible types. We found that a motion with a negative gravitational mass is possible although the energy condition of the matter is not violated. Using the result, formation of black hole and worm hole during the inflationary era is discussed.