In our previous work, we found new types of the cosmic string solutions in the Abelian-Higgs model with an enhanced $U(1)$ global symmetry. We dubbed those solutions as the compensated/uncompensated strings. The compensated string is similar to the conventional cosmic string in the Abrikosov-Nielsen-Olesen (ANO) string, around which only the would-be Nambu-Goldstone (NG) boson winds. Around the uncompensated string, on the other hand, the physical NG boson also winds, where the physical NG boson is associated with the spontaneous breaking of the enhanced symmetry. Our previous simulation in the 2+1 dimensional spacetime confirmed that both the compensated/uncompensated strings are formed at the phase transition of the symmetry breaking. Non-trivial winding of the physical NG boson around the strings potentially causes the so-called axion domain-wall problem when the model is applied to the axion model. In this paper, we perform simulation in the 3+1 dimensional spacetime to discuss the fate of the uncompensated strings. We observe that the evolution of the string-network is highly complicated in the 3+1 dimensional simulation compared with that seen in the previous simulation. Despite such complications, we find that the number of the uncompensated strings which could cause can be highly suppressed at late times. Our observation suggests that the present setup can be applied to the axion model without suffering from the axion domain-wall problem.

We study the cosmic microwave background (CMB) in the framework of the degenerate higher-order scalar-tensor (DHOST) theory to test gravity theories. This theoretical framework includes the wide class of dark energy models such as the Horndeski theory and its extensions as certain limits, and the general relativity can be also recovered. In this study, to test gravity theories with CMB, we formulate the linear perturbations of gravity and matter in the theory and their effective description parameterised by time-dependent effective field theory (EFT) parameters, $\alpha_i$ $(i=B,K,T,M,H,L)$ and $\beta_i$ $(i=1,2,3)$. Based on the resultant DHOST framework, we develop a numerical code to solve Boltzmann equations consistently. We then show that the angular power spectra of the CMB temperature anisotropies, E-mode and lensing potential as a demonstration and find that the parameter characterising the DHOST theory, $\beta_1$, provides the larger modifications of the spectra, compared with other EFT parameters. We also show the results in the case of a specific model in which the cosmic expansion as well as the EFT parameters are consistently determined.

We present the possibility that the seesaw mechanism with thermal leptogenesis can be tested using the stochastic gravitational background. Achieving neutrino masses consistent with atmospheric and solar neutrino data, while avoiding non-perturbative couplings, requires right-neutrinos lighter than the typical scale of grand unification. This scale separation suggests a symmetry protecting the right handed neutrinos from getting a mass. Thermal leptogenesis would then require that such a symmetry be broken below the reheating temperature. We enumerate all such possible symmetries consistent with these minimal assumptions and their corresponding defects, finding that in many cases, gravitational waves from the network of cosmic strings should be detectable. Estimating the predicted gravitational wave background we find that future space-borne missions could probe the entire range relevant for thermal leptogenesis.

We explore the formation and the evolution of the string network in the Abelian Higgs model with two complex scalar fields. A special feature of the model is that it possesses a global $U(1)$ symmetry in addition to the $U(1)$ gauge symmetry. Both symmetries are spontaneously broken by the vacuum expectation values of the two complex scalar fields. As we will show the dynamics of the string network is quite rich compared with that in the ordinary Abelian Higgs model with a single complex scalar field. In particular, we find a new type of string solutions in addition to the conventional Abrikosov-Nielsen-Olesen (local) string solution. We call this the uncompensated string. An isolated uncompensated string has a logarithmic divergent string tension as in the case of the global strings, although it is accompanied by a non-trivial gauge field configuration. We also perform classical lattice simulations in the $2+1$ dimensional spacetime, which confirms the formation of the uncompensated strings at the phase transition. We also find that most of the uncompensated strings evolve into the local strings at later time when the gauge charge of the scalar field with a smaller vacuum expectation value is larger than that of the scalar field with a larger vacuum expectation value.

We investigate a possibility of constraining statistical anisotropies of the primordial tensor perturbations by using future observations for the Cosmic Microwave Background (CMB) B-mode polarization. By parameterizing a statistically-anisotropic tensor power spectrum as $P_h ({\boldsymbol{k}}) = P_h (k) \sum_n g_n \cos^n \theta_{\boldsymbol{k}}$, where $\theta_{\boldsymbol{k}}$ is an angle of the direction of $\hat{k}={\boldsymbol{k}}/k$ from a preferred direction, we find that it would be possible for future B-mode observations such as CMB-S4 to detect the tensor statistical anisotropy at the level of $g_n \sim {\mathcal O} (0.1)$.

Given observations of B-mode polarization power spectrum of the cosmic microwave background (CMB), we can reconstruct power spectra of primordial tensor modes from the early Universe without assuming their functional form such as a power-law spectrum. Shape of the reconstructed spectra can then be used to probe the origin of tensor modes in a model-independent manner. We use the Fisher matrix to calculate the covariance matrix of tensor power spectra reconstructed in bins. We find that the power spectra are best reconstructed at wavenumbers in the vicinity of $k\approx 6\times 10^{-4}$ and $5\times 10^{-3}~{\rm Mpc}^{-1}$, which correspond to the "reionization bump" at $\ell\lesssim 6$ and "recombination bump" at $\ell\approx 80$ of the CMB B-mode power spectrum, respectively. The error bar between these two wavenumbers is larger because of lack of the signal between the reionization and recombination bumps. The error bars increase sharply towards smaller (larger) wavenumbers because of the cosmic variance (CMB lensing and instrumental noise). To demonstrate utility of the reconstructed power spectra we investigate whether we can distinguish between various sources of tensor modes including those from the vacuum metric fluctuation and SU(2) gauge fields during single-field slow-roll inflation, open inflation and massive gravity inflation. The results depend on the model parameters, but we find that future CMB experiments are sensitive to differences in these models. We make our calculation tool available on-line.

We study relativistic stars in degenerate higher-order scalar-tensor theories that evade the constraint on the speed of gravitational waves imposed by GW170817. It is shown that the exterior metric is given by the usual Schwarzschild solution if the lower order Horndeski terms are ignored in the Lagrangian and a shift symmetry is assumed. However, this class of theories exhibits partial breaking of Vainshtein screening in the stellar interior and thus modifies the structure of a star. Employing a simple concrete model, we show that for high-density stars the mass-radius relation is altered significantly even if the parameters are chosen so that only a tiny correction is expected in the Newtonian regime. We also find that, depending on the parameters, there is a maximum central density above which solutions cease to exist.

Under certain conditions the collision and intercommutation of two cosmic strings can result in the formation of a third string, with the three strings then remaining connected at Y-junctions. The kinematics and dynamics of collisions of this type have been the subject of analytical and numerical analyses in the special case in which the strings are Nambu-Goto. Cosmic strings, however, may well carry currents, in which case their dynamics is not given by the Nambu-Goto action. Our aim is to extend the kinematic analysis to more general kinds of string model. We focus in particular on the collision of strings described by conservative elastic string models, characteristic of current carrying strings, and which are expected to form in a cosmological context. As opposed to Nambu-Goto strings collisions, we show that in this case the collision cannot lead to the formation of a third elastic string: if dynamically such a string forms then the joining string must be described by a more general equation of state. This process will be studied numerically in a forthcoming publication.

We consider a cosmological model in which the tensor mode becomes massive during inflation, and study the Cosmic Microwave Background (CMB) temperature and polarization bispectra arising from the mixing between the scalar mode and the massive tensor mode during inflation. The model assumes the existence of a preferred spatial frame during inflation. The local Lorentz invariance is already broken in cosmology due to the existence of a preferred rest frame. The existence of a preferred spatial frame further breaks the remaining local SO(3) invariance and in particular gives rise to a mass in the tensor mode. At linear perturbation level, we minimize our model so that the vector mode remains non-dynamical, while the scalar mode is the same as the one in single-field slow-roll inflation. At non-linear perturbation level, this inflationary massive graviton phase leads to a sizeable scalar-scalar-tensor coupling, much greater than the scalar-scalar-scalar one, as opposed to the conventional case. This scalar-scalar-tensor interaction imprints a scale dependent feature in the CMB temperature and polarization bispectra. Very intriguingly, we find a surprizing similarity between the predicted scale dependence and the scale-dependent non-Gaussianities at low multipoles hinted in the WMAP and Planck results.

We investigate the imprint of reheating on the gravitational wave spectrum produced by self-ordering of multi-component scalar fields after a global phase transition. The equation of state of the Universe during reheating, which usually has different behaviour from that of a radiation-dominated Universe, affects the evolution of gravitational waves through the Hubble expansion term in the equations of motion. This gives rise to a different power-law behavior of frequency in the gravitational wave spectrum. The reheating history is therefore imprinted in the shape of the spectrum. We perform $512^3$ lattice simulations to investigate how the ordering scalar field reacts to the change of the Hubble expansion and how the reheating effect arises in the spectrum. We also compare the result with inflation-produced gravitational waves, which has a similar spectral shape, and discuss whether it is possible to distinguish the origin between inflation and global phase transition by detecting the shape with future direct detection gravitational wave experiments such as DECIGO.

The mechanism of thermal inflation, a relatively short period of accelerated expansion after primordial inflation, is a desirable ingredient for a certain class of particle physics models if they are not to be in contention with the cosmology of the early Universe. Though thermal inflation is most simply described in terms of a thermal effective potential, a thermal environment also gives rise to thermal fluctuations that must be taken into account. We numerically study the effects of these thermal fluctuations using lattice simulations. We conclude that though they do not ruin the thermal inflation scenario, the phase transition at the end of thermal inflation proceeds through phase mixing and is therefore not accompanied by the formations of bubbles nor appreciable amplitude of gravitational waves.

In this paper, we introduce a new approach to a treatment of the gravitational effects (redshift, time delay and lensing) on the observed cosmic microwave background (CMB) anisotropies based on the Boltzmann equation. From the Liouville's theorem in curved spacetime, the intensity of photons is conserved along a photon geodesic when non-gravitational scatterings are absent. Motivated by this fact, we derive a second-order line-of-sight formula by integrating the Boltzmann equation along a perturbed geodesic (curve) instead of a background geodesic (line). In this approach, the separation of the gravitational and intrinsic effects are manifest. This approach can be considered as a generalization of the remapping approach of CMB lensing, where all the gravitational effects can be treated on the same footing.

Redshift-space distortions (RSD) offer an attractive method to measure the growth of cosmic structure on large scales, and combining with the measurement of the cosmic expansion history, it can be used as cosmological tests of gravity. With the advent of future galaxy redshift surveys aiming at precisely measuring the RSD, an accurate modeling of RSD going beyond linear theory is a critical issue in order to detect or disprove small deviations from general relativity (GR). While several improved models of RSD have been recently proposed based on the perturbation theory (PT), the framework of these models heavily relies on GR. Here, we put forward a new PT prescription for RSD in general modified gravity models. As a specific application, we present theoretical predictions of the redshift-space power spectra in f(R) gravity model, and compare them with N-body simulations. Using the PT template that takes into account the effects of both modifications of gravity and RSD properly, we successfully recover the fiducial model parameter in N-body simulations in an unbiased way. On the other hand, we found it difficult to detect the scale dependence of the growth rate in a model-independent way based on GR templates.

Based on the multi-point propagator expansion, we present resummed perturbative calculations for cosmological power spectra and correlation functions in the context of modified gravity. In a wide class of modified gravity models that have a screening mechanism to recover general relativity (GR) on small scales, we apply the eikonal approximation to derive the governing equation for resummed propagator that partly includes the non-perturbative effect in the high-k limit. The resultant propagator in the high-k limit contains the new corrections arising from the screening mechanism as well as the standard exponential damping. We explicitly derive the expression for new high-k contributions in specific modified gravity models, and find that in the case of f(R) gravity for a currently constrained model parameter, the corrections are basically of the sub-leading order and can be neglected. Thus, in f(R) gravity, similarly to the GR case, we can analytically construct the regularized propagator that reproduces both the resummed high-k behavior and the low-k results computed with standard perturbation theory, consistently taking account of the nonlinear modification of gravity valid at large scales. With the regularized multi-point propagators, we give predictions for power spectrum and correlation function at one-loop order, and compare those with N-body simulations in f(R) gravity model. As an important application, we also discuss the redshift-space distortions and compute the anisotropic power spectra and correlation functions.

We revisit the production of gravitational waves from unstable domain walls analyzing their spectrum by the use of field theoretic lattice simulations with grid size $1024^3$, which is larger than the previous study. We have recognized that there exists an error in the code used in the previous study, and the correction of the error leads to the suppression of the spectrum of gravitational waves at high frequencies. The peak of the spectrum is located at the scale corresponding to the Hubble radius at the time of the decay of domain walls, and its amplitude is consistent with the naive estimation based on the quadrupole formula. Using the numerical results, the magnitude and the peak frequency of gravitational waves at the present time are estimated. It is shown that for some choices of parameters the signal of gravitational waves is strong enough to be probed in the future gravitational wave experiments.

The breaking of $U(1)_R$ symmetry plays a crucial role in modeling the breaking of supersymmetry (SUSY). In the models that possess both SUSY preserving and SUSY breaking vacua, tube-like cosmic strings called R-tubes, whose surfaces are constituted by domain walls interpolating a false and a true vacuum with some winding numbers, can exist. Their (in)stability can strongly constrain SUSY breaking models theirselves. In the present study, we investigate the dynamical (in)stability of two colliding metastable tube-like strings by field-theoretic simulations. From them, we find that the strings become unstable, depending on the relative collision angle and speed of two strings, and the false vacuum is eventually filled out by the true vacuum owing to rapid expansion of the strings or unstable bubbles created as remnants of the collision.

We study the network of Type-I cosmic strings using the field-theoretic numerical simulations in the Abelian-Higgs model. For Type-I strings, the gauge field plays an important role, and thus we find that the correlation length of the strings is strongly dependent upon the parameter \beta, the ratio between the self-coupling constant of the scalar field and the gauge coupling constant, namely, \beta=\lambda/2e^2. In particular, if we take the cosmic expansion into account, the network becomes densest in the comoving box for a specific value of \beta for \beta<1.

In massive gravity, galileon, and braneworld explanations of cosmic acceleration, force modifications are screened by nonlinear derivative self-interactions of the scalar field mediating that force. Interactions between the field of a central body ("A") and an orbiting body ("B") imply that body B does not move as a test body in the field of body A if the orbit is smaller than the Vainshtein radius of body B. We find through numerical solutions of the joint field at the position of B that the A-field Laplacian is nearly perfectly screened by the B self-field, whereas first derivative or net forces are reduced in a manner that scales with the mass ratio of the bodies as (M_B/M_A)^{3/5}. The latter causes mass-dependent reductions in the universal perihelion precession rate, with deviations for the Earth-Moon system at the ~4% level. In spite of universal coupling, which preserves the microscopic equivalence principle, the motion of macroscopic screened bodies depends on their mass providing in principle a means for testing the Vainshtein mechanism.

We investigate the cosmological constraints on axion models where the domain wall number is greater than one. In these models, multiple domain walls attached to strings are formed, and they survive for a long time. Their annihilation occurs due to the effects of explicit symmetry breaking term which might be raised by Planck-scale physics. We perform three-dimensional lattice simulations and compute the spectra of axions and gravitational waves produced by long-lived domain walls. Using the numerical results, we estimated relic density of axions and gravitational waves. We find that the existence of long-lived domain walls leads to the overproduction of cold dark matter axions, while the density of gravitational waves is too small to observe at the present time. Combining the results with other observational constraints, we find that the whole parameter region of models are excluded unless an unacceptable fine-tuning exists.

We analyze the spectrum of axions radiated from collapse of domain walls, which have received less attention in the literature. The evolution of topological defects related to the axion models is investigated by performing field-theoretic lattice simulations. We simulate the whole process of evolution of the defects, including the formation of global strings, the formation of domain walls and the annihilation of the defects due to the tension of walls. The spectrum of radiated axions has a peak at the low frequency, which implies that axions produced by the collapse of domain walls are not highly relativistic. We revisit the relic abundance of cold dark matter axions and find that the contribution from the decay of defects can be comparable with the contribution from strings. This result leads to a more severe upper bound on the axion decay constant.

We study the cosmological evolution of domain walls bounded by strings which arise naturally in axion models. If we introduce a bias in the potential, walls become metastable and finally disappear. We perform two dimensional lattice simulations of domain wall networks and estimate the decay rate of domain walls. By using the numerical results, we give a constraint for the bias parameter and the Peccei-Quinn scale. We also discuss the possibility to probe axion models by direct detection of gravitational waves produced by domain walls.

Cosmological evolution of axionic string network is analyzed in terms of field-theoretic simulations in a box of 512^3 grids, which are the largest ever, using a new and more efficient identification scheme of global strings. The scaling parameter is found to be \xi=0.87 +- 0.14 in agreement with previous results. The energy spectrum is calculated precisely using a pseudo power spectrum estimator which significantly reduces the error in the mean reciprocal comoving momentum. The resultant constraint on the axion decay constant leads to f_a <= 3*10^11 GeV. We also discuss implications for the early Universe.

We study Q-ball formation in the expanding universe on 1D, 2D and 3D lattice simulations. We obtain detailed Q-ball charge distributions, and find that the distribution is peaked at Q^{3D}_{peak} \simeq 1.9\times 10^{-2} (|\Phi_{in}|/m)^2, which is greater than the existing result by about 60%. Based on the numerical simulations, we discuss how the Q-ball formation proceeds. Also we make a comment on possible deviation of the charge distributions from what was conjectured in the past.

We study the production of gravitational waves from cosmic domain walls created during phase transition in the early universe. We investigate the process of formation and evolution of domain walls by running three dimensional lattice simulations. If we introduce an approximate discrete symmetry, walls become metastable and finally disappear. We calculate the spectrum of gravitational waves produced by collapsing metastable domain walls. Extrapolating the numerical results, we find the signal of gravitational waves produced by domain walls whose energy scale is around 10^10-10^12GeV will be observable in the next generation gravitational wave interferometers.

We study the non-linear evolution of baryon acoustic oscillations in the matter power spectrum and correlation function from the improved perturbation theory (PT). Based on the framework of renormalized PT, we apply the {\it closure approximation} that truncates the infinite series of loop contributions at one-loop order, and obtain a closed set of integral equations for power spectrum and non-linear propagator. The resultant integral expressions keep important non-perturbative properties which can dramatically improve the prediction of non-linear power spectrum. Employing the Born approximation, we then derive the analytic expressions for non-linear power spectrum and the predictions are made for non-linear evolution of baryon acoustic oscillations in power spectrum and correlation function. A detailed comparison between improved PT results and N-body simulations shows that a percent-level agreement is achieved in a certain range in power spectrum and in a rather wider range in correlation function. Combining a model of non-linear redshift-space distortion, we also evaluate the power spectrum and correlation function in correlation function. In contrast to the results in real space, the agreement between N-body simulations and improved PT predictions tends to be worse, and a more elaborate modeling for redshift-space distortion needs to be developed. Nevertheless, with currently existing model, we find that the prediction of correlation function has a sufficient accuracy compared with the cosmic-variance errors for future galaxy surveys with volume of a few (Gpc/h)^3 at z>=0.5.

We present a formalism to calculate the non-linear matter power spectrum in modified gravity models that explain the late-time acceleration of the Universe without dark energy. Any successful modified gravity models should contain a mechanism to recover General Relativity (GR) on small scales in order to avoid the stringent constrains on deviations from GR at solar system scales. Based on our formalism, the quasi non-linear power spectrum in the Dvali-Gabadadze-Porratti (DGP) braneworld models and $f(R)$ gravity models are derived by taking into account the mechanism to recover GR properly. We also extrapolate our predictions to fully non-linear scales using the Parametrized Post Friedmann (PPF) framework. In $f(R)$ gravity models, the predicted non-linear power spectrum is shown to reproduce N-body results. We find that the mechanism to recover GR suppresses the difference between the modified gravity models and dark energy models with the same expansion history, but the difference remains large at weakly non-linear regime in these models. Our formalism is applicable to a wide variety of modified gravity models and it is ready to use once consistent models for modified gravity are developed.

We present a new numerical scheme to treat the non-linear evolution of cosmological power spectra. Governing equations for matter power spectra have been previously derived by a non-perturbative technique with closure approximation. Solutions of the resultant closure equations just correspond to the resummation of an infinite class of perturbation corrections, and they consistently reproduce the one-loop results of standard perturbation theory. We develop a numerical algorithm to solve closure evolutions in both perturbative and non-perturbative regimes. The present numerical scheme is particularly suited for examining non-linear matter power spectrum in general cosmological models, including modified theory of gravity. As a demonstration, we study weakly non-linear evolution of power spectrum in a class of modified gravity models, as well as various dark energy models.

We apply a non-linear statistical method in turbulence to the cosmological perturbation theory and derive a closed set of evolution equations for matter power spectra. The resultant closure equations consistently recover the one-loop results of standard perturbation theory and beyond that, it is still capable of treating the non-linear evolution of matter power spectra. We find the exact integral expressions for the solutions of closure equations. These analytic expressions coincide with the renormalized one-loop results presented by Crocce & Scoccimarro (2006,2007). By constructing the non-linear propagator, we analytically evaluate the non-linear matter power spectra based on the first-order Born approximation of the integral expressions and compare it with those of the renormalized perturbation theory.

We study the behaviour of scalar perturbations in the radiation-dominated era of Randall-Sundrum braneworld cosmology by numerically solving the coupled bulk and brane master wave equations. We find that density perturbations with wavelengths less than a critical value (set by the bulk curvature length) are amplified during horizon re-entry. This means that the radiation era matter power spectrum will be at least an order of magnitude larger than the predictions of general relativity (GR) on small scales. Conversely, we explicitly confirm from simulations that the spectrum is identical to GR on large scales. Although this magnification is not relevant for the cosmic microwave background or measurements of large scale structure, it will have some bearing on the formation of primordial black holes in Randall-Sundrum models.

In this paper we quantise scalar perturbations in a Randall-Sundrum-type model of inflation where the inflaton field is confined to a single brane embedded in five-dimensional anti-de Sitter space-time. In the high energy regime, small-scale inflaton fluctuations are strongly coupled to metric perturbations in the bulk and gravitational back-reaction has a dramatic effect on the behaviour of inflaton perturbations on sub-horizon scales. This is in contrast to the standard four-dimensional result where gravitational back-reaction can be neglected on small scales. Nevertheless, this does not give rise to significant particle production, and the correction to the power spectrum of the curvature perturbations on super-horizon scales is shown to be suppressed by a slow-roll parameter. We calculate the complete first order slow-roll corrections to the spectrum of primordial curvature perturbations.

We discuss a robust data analysis method to detect a stochastic background of gravitational waves in the presence of non-Gaussian noise. In contrast to the standard cross-correlation (SCC) statistic frequently used in the stochastic background searches, we consider a {\it generalized cross-correlation} (GCC) statistic, which is nearly optimal even in the presence of non-Gaussian noise. The detection efficiency of the GCC statistic is investigated analytically, particularly focusing on the statistical relation between the false-alarm and the false-dismissal probabilities, and the minimum detectable amplitude of gravitational-wave signals. We derive simple analytic formulae for these statistical quantities. The robustness of the GCC statistic is clarified based on these formulae, and one finds that the detection efficiency of the GCC statistic roughly corresponds to the one of the SCC statistic neglecting the contribution of non-Gaussian tails. This remarkable property is checked by performing the Monte Carlo simulations and successful agreement between analytic and simulation results was found.

We discuss the cosmological evolution of the inflationary gravitational wave background (IGWB) in the Randall-Sundrum single-brane model. In the braneworld cosmology, in which three-dimensional space-like hypersurface that we live in is embedded in five-dimensional anti de Sitter (AdS_5) space-time, the evolution of gravitational wave (GW) modes is affected by the non-standard expansion of the universe and the excitation of the Kaluza-Klein modes (KK-modes), which are significant in the high-energy regime of the universe. We numerically evaluate these two effects by solving the evolution equation for GWs propagating through the AdS_5 space-time. Using a plausible initial condition from inflation, we find that the excitation of KK-modes can be characterized by a simple scaling relation above the critical frequency f_crit determined from the length scale of the fifth dimension \ell. The remarkable point is that this relation generally holds as long as the matter content of the universe is described by the perfect fluid with the EOS p=w\rho for 0\leq w\leq 1. The resultant scaling relation is translated into the energy spectrum of the IGWB as \Omega_GW\propto f^{(3w-1)/(3w+2)} for f>f_crit. This indicates that, in the radiation dominant case (w=1/3), the two high-energy effects accidentally compensate each other and the spectrum becomes almost the same as the one predicted in the four-dimensional theory, i.e., \Omega_GW\propto f^0.

Future missions of gravitational-wave astronomy will be operated by space-based interferometers, covering very wide range of frequency. Search for stochastic gravitational-wave backgrounds (GWBs) is one of the main targets for such missions, and we here discuss the prospects for direct measurement of isotropic and anisotropic components of (primordial) GWBs around the frequency 0.1-10 Hz. After extending the theoretical basis for correlation analysis, we evaluate the sensitivity and the signal-to-noise ratio for the proposed future space interferometer missions, like Big-Bang Observer (BBO), Deci-Hertz Interferometer Gravitational-wave Observer (DECIGO) and recently proposed Fabry-Perot type DECIGO. The astrophysical foregrounds which are expected at low frequency may be a big obstacle and significantly reduce the signal-to-noise ratio of GWBs. As a result, minimum detectable amplitude may reach h^2 \ogw = 10^{-15} \sim 10^{-16}, as long as foreground point sources are properly subtracted. Based on correlation analysis, we also discuss measurement of anisotropies of GWBs. As an example, the sensitivity level required for detecting the dipole moment of GWB induced by the proper motion of our local system is closely examined.

We discuss the gravitational wave background (GWB) from a cosmological population of gamma-ray bursts (GRBs). Among various emission mechanisms for the gravitational waves (GWs), we pay a particular attention to the vast anisotropic neutrino emissions from the accretion disk around the black hole formed after the so-called failed supernova explosions. The produced GWs by such mechanism are known as burst with memory, which could dominate over the low-frequency regime below \sim 10Hz. To estimate their amplitudes, we derive general analytic formulae for gravitational waveform from the axisymmetric jets. Based on the formulae, we first quantify the spectrum of GWs from a single GRB. Then, summing up its cosmological population, we find that the resultant value of the density parameter becomes roughly \Omega_{GW} \approx 10^{-20} over the wide-band of the low-frequency region, f\sim 10^{-4}-10^1Hz. The amplitude of GWB is sufficiently smaller than the primordial GWBs originated from an inflationary epoch and far below the detection limit.

We discuss the cosmological evolution of gravitational waves (GWs) after inflation in a brane-world cosmology embedded in five-dimensional anti-de Sitter (AdS_5) bulk spacetime. In a brane-world scenario, the evolution of GWs is affected by the non-standard cosmological expansion and the excitation of the Kaluza-Klein modes (KK-modes), which are significant in the high-energy regime of the universe. We numerically solve the wave equation of GWs in the Poincare coordinates of the AdS_5 spacetime. Using a plausible initial condition from inflation, we find that, while the behavior of GWs in the bulk is sensitive to the transition time from inflation to the radiation dominated epoch, the amplitude of GWs on the brane is insensitive to this time if the transition occurs early enough before horizon re-entry. As a result, the amplitude of GWs is suppressed by the excitation of KK-modes which escape from the brane into the bulk, and the effect may compensate the enhancement of the GWs by the non-standard cosmological expansion. Based on this, the influence of the high-energy effects on the stochastic background of GWs is discussed.

We study the evolution of gravitational waves(GWs) after inflation in a brane-world cosmology embedded in five-dimensional anti-de Sitter spacetime. Contrary to the standard four-dimensional results, the GWs at the high-energy regime in brane-world model suffer from the effects of the non-standard cosmological expansion and the excitation of the Kaluza-Klein modes(KK-modes), which can affect the amplitude of stochastic gravitational wave background significantly. To investigate these two high-energy effects quantitatively, we numerically solve the wave equation of the GWs in the radiation dominated epoch at relatively low-energy scales. We show that the resultant GWs is suppressed by the excitation of the KK modes. The created KK modes are rather soft and escape away from the brane to the bulk gravitational field. The results are also compared to the semi-analytic prediction from the low-energy approximation and the evolved amplitude of GWs on the brane reasonably matches the numerical simulations.