We study evolution of a hierarchical four-body (2+2) system composed by a pair of mass-transferring white dwarf binaries. Applying a simplified model around the synchronous state of two inner orbital periods, we newly find that the four body system could settle down to a limit cycle with a small period gap. The period gap generates an amplitude variation of emitted gravitational waves, as a beat effect. Depending on model parameters, the beat period could be 1-10 yr and a large amplitude variation might be observed by space gravitational-wave detectors.

We apply the independent component analysis (ICA) to the real data from a gravitational wave detector for the first time. Specifically we use the iKAGRA data taken in April 2016, and calculate the correlations between the gravitational wave strain channel and 35 physical environmental channels. Using a couple of seismic channels which are found to be strongly correlated with the strain, we perform ICA. Injecting a sinusoidal continuous signal in the strain channel, we find that ICA recovers correct parameters with enhanced signal-to-noise ratio, which demonstrates usefulness of this method. Among the two implementations of ICA used here, we find the correlation method yields the optimal result for the case environmental noises act on the strain channel linearly.

We discuss the prospects of high precision pointing of our transmitter to habitable planets around Galactic main sequence stars. For an efficient signal delivery, the future sky positions of the host stars should be appropriately extrapolated with accuracy better than the beam opening angle $\Theta$ of the transmitter. Using the latest data release (DR2) of Gaia, we estimate the accuracy of the extrapolations individually for $4.7\times 10^7$ FGK stars, and find that the total number of targets could be $\sim 10^7$ for the accuracy goal better than 1". Considering the pairwise nature of communication, our study would be instructive also for SETI (Search for Extraterrestrial Intelligence), not only for sending signals outward.

We discuss the prospects of LISA for detecting neutron star binaries (NSBs) in the local group galaxies such as LMC and M31. Using the recently estimated merger rate $ {\rm 1540 \, Gpc^{-3}yr^{-1}}$ and inversely applying the conventional arguments based on the $B$-band galaxy luminosities, we estimate the frequency distributions of NSBs in the local galaxies. We find that, after 10 year observation with its current design sensitivity, LISA might detect $\sim 5$ NSBs both in LMC and M31 with signal-to-noise ratios larger than 10. Some of the NSBs might be three-dimensionally localized well within LMC. These binaries will be useful for studying various topics including the origin of r-process elements.

We discuss a Galaxy-wide coordinated signaling scheme with which a SETI observer needs to examine a tiny fraction of the sky. The target sky direction is determined as a function of time, based on high-precision measurements of a progenitor of a conspicuous astronomical event such as a coalescence of a double neutron star binary. In various respects, such a coordinated scheme would be advantageous for both transmitters and receivers, and might be widely prevailing as a tacit adjustment. For this scheme, the planned space gravitational-wave detector LISA and its follow-on missions have a potential to narrow down the target sky area by a factor of $10^{3\textit{-}4}$, and could have a large impact on future SETI experiments.

We discuss the possibility of receiving a radio signal from extra-Galactic intelligence, around the time when we observe a binary neutron star merger in their galaxy. High-precision measurements of the binary parameters allow them to send the signal ~ $10^4$ years before they themselves observe the merger signal. Using the SKA, we might receive ~ $10^4$ bits of data, transmitted from 40 Mpc distance with the output power of ~ 1 TW. We also discuss related topics for GW170817 and mention potential roles of future gravitational wave detectors in relation to this transmission scheme.

KAGRA is a second-generation interferometric gravitational-wave detector with 3-km arms constructed at Kamioka, Gifu in Japan. It is now in its final installation phase, which we call bKAGRA (baseline KAGRA), with scientific observations expected to begin in late 2019. One of the advantages of KAGRA is its underground location of at least 200 m below the ground surface, which brings small seismic motion at low frequencies and high stability of the detector. Another advantage is that it cools down the sapphire test mass mirrors to cryogenic temperatures to reduce thermal noise. In April-May 2018, we have operated a 3-km Michelson interferometer with a cryogenic test mass for 10 days, which was the first time that km-scale interferometer was operated at cryogenic temperatures. In this article, we report the results of this "bKAGRA Phase 1" operation. We have demonstrated the feasibility of 3-km interferometer alignment and control with cryogenic mirrors.

We discuss a multimessenger strategy to detect radio pulses from Galactic binary neutron stars in a very tight orbit with the period shorter than 10 min. On one hand, all-sky surveys by radio instruments are inefficient for detecting faint pulsars in very tight binaries due partly to the rarity of targets and primarily to the need of correction for severe Doppler smearing. On the other hand, the Laser Interferometer Space Antenna (LISA) will detect these binaries with a very large signal-to-noise ratio and determine the orbital frequency, binary parameters, and sky location to high accuracy. The information provided by LISA will reduce the number of required pointings by two to six orders of magnitude and that of required trials for the corrections by about nine orders of magnitude, increasing the chance of discovering radio pulsars. For making full use of this strategy, it is desirable to operate high-sensitivity radio instruments such as Square Kilometer Array Phase 2 simultaneously with LISA.

The recent detections of gravitational waves (GWs) reported by LIGO/Virgo collaborations have made significant impact on physics and astronomy. A global network of GW detectors will play a key role to solve the unknown nature of the sources in coordinated observations with astronomical telescopes and detectors. Here we introduce KAGRA (former name LCGT; Large-scale Cryogenic Gravitational wave Telescope), a new GW detector with two 3-km baseline arms arranged in the shape of an "L", located inside the Mt. Ikenoyama, Kamioka, Gifu, Japan. KAGRA's design is similar to those of the second generations such as Advanced LIGO/Virgo, but it will be operating at the cryogenic temperature with sapphire mirrors. This low temperature feature is advantageous for improving the sensitivity around 100 Hz and is considered as an important feature for the third generation GW detector concept (e.g. Einstein Telescope of Europe or Cosmic Explorer of USA). Hence, KAGRA is often called as a 2.5 generation GW detector based on laser interferometry. The installation and commissioning of KAGRA is underway and its cryogenic systems have been successfully tested in May, 2018. KAGRA's first observation run is scheduled in late 2019, aiming to join the third observation run (O3) of the advanced LIGO/Virgo network. In this work, we describe a brief history of KAGRA and highlights of main feature. We also discuss the prospects of GW observation with KAGRA in the era of O3. When operating along with the existing GW detectors, KAGRA will be helpful to locate a GW source more accurately and to determine the source parameters with higher precision, providing information for follow-up observations of a GW trigger candidate.

We study the possibility of orbital synchronization capture for a hierarchical quadrupole stellar system composed by two binaries emitting gravitational waves. Based on a simple model including the mass transfer for white dwarf binaries, we find that the capture might be realized for inter-binary distances less than their gravitational wavelength. We also discuss related intriguing phenomena such as a parasitic relation between the coupled white dwarf binaries and significant reductions of gravitational and electromagnetic radiations.

Following the detection of the GW170817 signal and its associated electromagnetic emissions, we discuss the prospects of the local Hubble parameter measurement using double neutron stars (DNSs). The kilonova emissions of GW170817 are genuinely unique in terms of the rapid evolution of color and magnitude and we expect that, for a good fraction $\gtrsim 50\%$ of the DNS events within $\sim 200$Mpc, we could identify their host galaxies, using their kilonovae. At present, the estimated DNS merger rate $(1.5^{+3.2}_{-1.2})\times 10^{-6} {\rm Mpc^{-3} yr^{-1}}$ has a large uncertainty. But, if it is at the high end, we could measure the local Hubble parameter $H_L$ with the level of $\Delta H_L/H_L\sim 0.042$ ($1\sigma$ level), after the third observational run (O3). This accuracy is four times better than that obtained from GW170817 alone, and we will be able to examine the Hubble tension at $2.1\sigma$ level.

KAGRA is a 3-km interferometric gravitational wave telescope located in the Kamioka mine in Japan. It is the first km-class gravitational wave telescope constructed underground to reduce seismic noise, and the first km-class telescope to use cryogenic cooling of test masses to reduce thermal noise. The construction of the infrastructure to house the interferometer in the tunnel, and the initial phase operation of the interferometer with a simple 3-km Michelson configuration have been completed. The first cryogenic operation is expected in 2018, and the observing runs with a full interferometer are expected in 2020s. The basic interferometer configuration and the current status of KAGRA are described.

Based on the secular theory, we discuss the orbital evolution of stars in a nuclear star cluster to which a secondary massive black hole is infalling with vanishing eccentricity. We find that the eccentricities of the stars could show sharp transitions, depending strongly on their initial conditions. By examining the phase-space structure of an associated Hamiltonian, we show that these characteristic behaviors are partly due to a probabilistic bifurcation at a separatrix crossing, resulting from the retrograde apsidal precession by the cluster potential. We also show that separatrix crossings are closely related to realization of a large eccentricity and could be important for astrophysical phenomena such as tidal disruption events or gravitational wave emissions.

We propose that stellar-mass binary black holes like GW150914 will become a tool to explore the local Universe within ~100Mpc in the era of the Laser Interferometer Space Antenna (LISA). High calibration accuracy and annual motion of LISA could enable us to localize up to ~60 binaries more accurately than the error volume of ~100Mpc^3 without electromagnetic counterparts under moderately optimistic assumptions. This accuracy will give us a fair chance to determine the host object solely by gravitational waves. By combining the luminosity distance extracted from gravitational waves with the cosmological redshift determined from the host, the local value of the Hubble parameter will be determined up to a few % without relying on the empirically constructed distance ladder. Gravitational-wave cosmography would pave the way for resolution of the disputed Hubble tension, where the local and global measurements disagree in the value of the Hubble parameter at 3.4sigma level, which amounts to ~9%.

The primary target for the planned space-borne gravitational wave interferometers DECIGO/BBO is a primordial gravitational wave background (PGWB). However there exist astrophysical foregrounds and among them, gravitational waves from neutron star (NS) binaries are the solid and strong component that must be identified and subtracted. In this paper, we discuss the geometry of detector configurations preferable for identifying the NS/NS binary signals. As a first step, we analytically estimate the minimum signal-to-noise ratios (SNRs) of the binaries for several static detector configurations that are characterized by adjustable geometrical parameters, and determine the optimal values for these parameters. Next we perform numerical simulations to take into account the effect of detector motions, and find reasonable agreements with the analytical results. We show that, with the standard network formed by 4 units of triangle detectors, the proposed BBO sensitivity would be sufficient in receiving gravitational waves from all the NS/NS binaries at $z\le 5$ with SNRs higher than 25. We also discuss the minimum sensitivity of DECIGO required for the foreground identification.

We discuss the gravitational wave emission and the orbital evolution of a hierarchical triple system composed of an inner binary black hole (BBH) and an outer tertiary. Depending on the kick velocity at the merger, the merged BBH could tidally disrupt the tertiary. Even though the fraction of BBH mergers accompanied by such disruptions is expected to be much smaller than unity, the existence of a tertiary and its basic parameters (e.g. semimajor axis, projected mass) can be examined for more than 1000 BBHs with the space GW detector LISA and its follow-on missions. This allows us to efficiently prescreen the targets for the follow-up searches for the tidal disruption events (TDEs). The TDE probability would be significantly higher for triple systems with aligned orbital- and spin-angular momenta, compared with random configurations.

We study prospects for detecting extragalactic binary black holes similar to GW150914 by evolved Laser Interferometer Space Antenna (eLISA). We find that the majority of detected binary black holes will not merge within reasonable observation periods of eLISA in any configuration. While long-arm detectors are highly desired for promoting multiband gravitational-wave astronomy by increasing the detections of merging binaries, the number of total detections can be increased also by improving the acceleration noise. A monochromatic approximation works well to derive semiquantitative features of observational prospects for non-merging binaries with clearly indicating the parameter dependence. Our estimate also suggests that the number of galaxies in the error volume is so small that the host galaxy may be determined uniquely with high confidence.

Pre-DECIGO consists of three spacecraft arranged in an equilateral triangle with 100km arm lengths orbiting 2000km above the surface of the earth. It is hoped that the launch date will be in the late 2020s. Pre-DECIGO has one clear target: binary black holes (BBHs) like GW150914 and GW151226. Pre-DECIGO can detect $\sim 30M_\odot-30M_\odot$ BBH mergers up to redshift $z\sim 30$. The cumulative event rate is $\sim 1.8\times 10^{5}\,{\rm events~yr^{-1}}$ in the Pop III origin model of BBHs like GW150914, and it saturates at $z\sim 10$, while in the primordial BBH (PBBH) model, the cumulative event rate is $ \sim 3\times 10^{4}\,{\rm events~ yr^{-1}}$ at $z=30$ even if only $0.1\%$ of the dark matter consists of PBHs, and it is still increasing at $z=30$. In the Pop I/II model of BBHs, the cumulative event rate is $(3-10)\times10^{5}\,{\rm events~ yr^{-1}}$ and it saturates at $z \sim 6$. We present the requirements on orbit accuracy, drag free techniques, laser power, frequency stability, and interferometer test mass. For BBHs like GW150914 at 1Gpc, SNR$\sim 90$ is achieved with the definition of Pre-DECIGO in the $0.01-100$Hz band. Pre-DECIGO can measure the mass spectrum and the $z$-dependence of the merger rate to distinguish various models of BBHs like GW150914. Pre-DECIGO can also predict the direction of BBHs at $z=0.1$ with an accuracy of $\sim 0.3\,{\rm deg}^2$ and a merging time accuracy of $\sim 1$s at about a day before the merger so that ground-based GW detectors further developed at that time as well as electromagnetic follow-up observations can prepare for the detection of merger in advance. For intermediate mass BBHs at a large redshift $z > 10$, the QNM frequency after the merger can be within the Pre-DECIGO band so that the ringing tail can also be detectable to confirm the Einstein theory of general relativity with SNR$\sim 35$. [abridged]

Based on a simple geometrical approach, we analyze the evolution of the Kozai-Lidov mechanism for stars around shrinking massive black hole binaries on circular orbits. We find that, due to a peculiar bifurcation pattern induced by the Newtonian potential of stellar clusters, the orbit of stars could become highly eccentric. This transition occurs abruptly for stars with small initial eccentricities. Our approach would be also useful for studying the Kozai-Lidov mechanism in various astrophysical contexts.

We discuss the prospects of eLISA for detecting gravitational waves (GWs) from Galactic binary black holes (BBHs) similar to GW150914. For a comoving merger rate that is consistent with current observation, eLISA is likely to identify at least one BBH with a sufficient signal-to-noise ratio. In addition, eLISA has a potential to measure the eccentricity of the BBH as small as $e\sim 0.02$, corresponding to the residual value $e\sim 10^{-6}$ at 10Hz. Therefore, eLISA could provide us with a crucial information to understand the formation processes of relatively massive BBHs like GW150914. We also derive a simple scaling relation for the expected number of detectable Galactic BBHs.

A period of inflation in the early universe produces a nearly scale-invariant spectrum of gravitational waves over a huge range in wavelength. If the amplitude of this gravitational wave background is large enough to be detectable with microwave background polarization measurements, it will also be detectable directly with a space-based laser interferometer. Using a Monte Carlo sampling of inflation models, we demonstrate that the combination of these two measurements will strongly constrain the expansion history during inflation and the physical mechanism driving it.

In the light of the history of researches on electromagnetic wave spectrum, a sharp emission line of gravitational-wave background (GWB) would be an interesting observational target. Here we study an efficient method to detect a line GWB by correlating data of multiple ground-based detectors. We find that the width of frequency bin for coarse graining is a critical parameter, and the commonly-used value 0.25 Hz is far from optimal, decreasing the signal-to-noise ratio by up to a factor of seven. By reanalyzing the existing data with a smaller bin width, we might detect a precious line signal from the early universe.

We analytically discuss probability distribution function (PDF) for inclinations of merging compact binaries whose gravitational waves are coherently detected by a network of ground based interferometers. The PDF would be useful for studying prospects of (1) simultaneously detecting electromagnetic signals (such as gamma-ray-bursts) associated with binary mergers and (2) statistically constraining the related theoretical models from the actual observational data of multi-messenger astronomy. Our approach is similar to Schutz (2011), but we explicitly include the dependence of the polarization angles of the binaries, based on the concise formulation given in Cutler and Flanagan (1994). We find that the overall profiles of the PDFs are similar for any networks composed by the second generation detectors (Advanced-LIGO, Advanced-Virgo, KAGRA, LIGO-India). For example, 5.1% of detected binaries would have inclination angle less than 10 degree with at most 0.1% differences between the potential networks. A perturbative expression is also provided for generating the PDFs with a small number of parameters given by directional averages of the quantity $\epsilon$ that characterises the asymmetry of network sensitivities to incoming two orthogonal polarization modes.

In a recent paper, Schutz proposed an analytical approximation for simplifying treatment of polarization angle and conveniently evaluating relative detection rates of compact binary inspirals for various networks of ground-based interferometers. We derived relative event rates by strictly handling polarization angle and quantitatively examine validity of Schutz's approximation. The associated error of the approximation is rigorously shown to be less than 1.02\%, irrespective of details of the detector networks.

We study the possibility of pre-merger localization of eccentric compact binary coalescences by second-generation gravitational-wave detector networks. Gravitational waves from eccentric binaries can be regarded as a sequence of pulses, which are composed of various higher harmonic modes than ones with twice the orbital frequency. The higher harmonic modes from a very early inspiral phase will not only contribute to the signal-to-noise ratio, but also allow us to localize the gravitational-wave source before the merger sets in. This is due to the fact that high-frequency gravitational waves are essential for the source localization via triangulation by ground-based detector networks. We found that the single-detector signal-to-noise ratio exceeds 5 at 10 min before the merger for a 1.4--1.4 $M_\odot$ eccentric binary neutron stars at 100 Mpc in optimal cases, and it can be localized up to 10 deg$^2$ at half a minute before the merger by a four-detector network. We will even be able to achieve $\sim$ 10 deg$^2$ at 10 min before the merger for a face-on eccentric compact binary by a five-detector network.

The Kozai mechanism for a hierarchical triple system could reduce the merger time of inner eccentric binary emitting gravitational waves (GWs), and has been qualitatively explained with the secular theory that is derived by averaging short-term orbital revolutions. However, with the secular theory, the minimum value of the inner pericenter distance could be excessively limited by the averaging operation. Compared with traditional predictions, the actual evolution of an eccentric inner binary could be accompanied by (i) a higher characteristic frequency of the pulse-like GWs around its pericenter passages, and (ii) a larger residual eccentricity at its final inspiral phase. These findings would be important for GW astronomy with the forthcoming advanced detectors.

We study dynamical evolution of a resonant triple system formed by an inner EMRI and an additional outer MBH. The relevant resonant state (\lambda_2-\varpi_1 ~ const) is supported by the relativistic apsidal precession of the inner EMRI, and, unlike standard mean motion resonances, the triple system can have a hierarchical orbital configuration (but different from the Kozai process). In order to analyze this unusual resonant system, we extend the so-called Hamiltonian approach, and derive a mapping from the EMRI-MBH triple system to a simple one-dimensional Hamiltonian. With the derived mapping, we make analytical predictions for characteristic quantities of the resonance, such as the capture probability, and find that they reasonably agree with numerical simulations up to moderate eccentricities.

We investigate the capability of pulsar timing arrays (PTAs) as a probe of primordial black holes (PBHs), which might constitute the Galactic dark matter. A PBH passing nearby the Earth or a pulsar gives an impulse acceleration and induces residuals on otherwise orderly pulsar timing data. We show that the timing residuals induced at pulsars are optimal for searching heavier PBHs than those at the Earth, and the two probes are highly complemental. Future facilities like SKA could detect PBHs with masses around 10^(22-28) g even if only a small fraction (< 1%) of the Galactic dark matter consists of these PBHs.

The maximum likelihood method is often used for parameter estimation in gravitational wave astronomy. Recently, an interesting approach was proposed by Vallisneri to evaluate the distributions of parameter estimation errors expected for the method. This approach is to statistically analyze the local peaks of the likelihood surface, and works efficiently even for signals with low signal-to-noise ratios. Focusing special attention to geometric structure of the likelihood surface, we follow the proposed approach and derive formulae for a simplified model of data analysis where the target signal has only one intrinsic parameter, along with its overall amplitude. Then we apply our formulae to correlation analysis of stochastic gravitational wave background with a power-law spectrum. We report qualitative trends of the formulae using numerical results specifically obtained for correlation analysis with two Advanced-LIGO detectors.

One component of a massive black hole binary (MBHB) might capture a small third body, and then a hierarchical, inclined triple system would be formed. With the post-Newtonian approximation including radiation reaction, we analyzed the evolution of the triple initially with small eccentricities. We found that an essentially new resonant relation could arise in the triple system. Here relativistic effects are crucial. Relativistic resonances, including the new one, stably work even for an outer MBHB of comparable masses, and significantly change the orbit of the inner small body.

A massive black hole binary might resonantly trap a star (e.g. a white dwarf) and the gas released by its tidal disruption might emit electromagnetic wave signals around the coalescence of the binary. With post-Newtonian equations of motion including gravitational radiation reaction, we numerically studied resonant trappings by black hole binaries with mass ratio 1/100. It is found that 2:1 (and simultaneously 4:2) mean motion resonances of the binaries would be strong and could, in principle, draw small third objects deep into relativistic regimes (e.g. ~10 Schwarzschild radii). The inclinations of the trapped objects could increase significantly and, in some cases, retrograde orbits could be realized eventually.

We show that a massive black hole binary might resonantly trap a small third body (e.g. a neutron star) down to a stage when the binary becomes relativistic due to its orbital decay by gravitational radiation. The final fate of the third body would be quite interesting for relativistic astrophysics. For example, the parent binary could expel the third body with a velocity more than 10% of the speed of light. We also discuss the implications of this three-body system for direct gravitational wave observation.

The merger of a super-massive binary black hole (SBBH) is one of the most extreme events in the universe with a huge amount of energy released by gravitational radiation. Although the characteristic gravitational wave (GW) frequency around the merger event is far higher than the nHz regime optimal for pulsar timing arrays (PTAs), nonlinear GW memory might be a critical smoking gun of the merger event detectable with PTAs. In this paper, basic aspects of this interesting observation are discussed for SBBHs, and the detection numbers of their memory and inspiral GWs are estimated for ongoing and planned PTAs. We find that the expected detection number would be smaller than unity for the two-types of signals even with the Square Kilometer Array. We also provide various scaling relations that would be useful to study detection probabilities of GWs from individual SBBHs with PTAs.

We study an observational method to analyze non-Gaussianity of a gravitational wave (GW) background made by superposition of weak burst signals. The proposed method is based on fourth-order correlations of data from four detectors, and might be useful to discriminate the origin of a GW background. With a formulation newly developed to discuss geometrical aspects of the correlations, it is found that the method provides us with linear combinations of two interesting parameters, I_2 and V_2 defined by the Stokes parameters of individual GW burst signals. We also evaluate sensitivities of specific detector networks to these parameters.

We propose a non-Gaussianity test for gravitational wave backgrounds by combining data streams of multiple detectors. This simple method allows us to check whether a detected background is "smooth" enough to be consistent with an inflation-type background, or is contaminated by individually undetectable weak burst signals. The proposed test would be quite useful for the Big Bang Observer or DECIGO whose primary target is a background from inflation at 0.1-1Hz where gravitational wave bursts from supernovae of population III stars might become a troublesome foreground.

I propose a method to detect planets around compact binaries that are strong sources of gravitational radiation. This approach is to measure gravitational-wave phase modulations induced by the planets, and its prospect is studied with a Fisher matrix analysis. I find that, using the Laser Interferometer Space Antenna (LISA), planets can be searched for around ~3000 Galactic double white dwarfs with detection limit ~4 M_J (M_J~ 2x10^{30}g: the Jupiter mass). With its follow-on missions, planets with mass ~1 M_J might be detected around double neutron stars even at cosmological distances z~1. In this manner, gravitational wave observation has potential to make interesting contributions to extra-solar planetary science.

The Stokes V parameter characterizes asymmetry of amplitudes between right- and left-handed waves, and non-vanishing value of the V parameter yields a circularly polarized signal. Cosmologically, V parameter may be a direct probe for parity violation in the universe. In this paper, we theoretically investigate a measurement of this parameter, particularly focusing on the gravitational-wave backgrounds observed via ground-based interferometers. In contrast to the traditional analysis that only considers the total amplitude (or equivalently $\Omega_{GW}$), the signal analysis including a circular-polarized mode has a rich structure due to the multi-dimensionality of target parameters. We show that, by using the network of next-generation detectors, separation between polarized and unpolarized modes can be performed with small statistical loss induced by their correlation.

We show that pairs of widely separated interferometers are advantageous for measuring the Stokes parameter V of a stochastic background of gravitational waves. This parameter characterizes asymmetry of amplitudes of right- and left-handed waves and generation of the asymmetry is closely related to parity violation in the early universe. The advantageous pairs include LIGO(Livingston)-LCGT and AIGO-Virgo that are relatively insensitive to Omega_GW (the simple intensity of the background). Using at least three detectors, information of the intensity Omega_GW and the degree of asymmetry V can be separately measured.

We discuss the possibility of detecting the presence of primordial black holes (PBHs), such as those that might account for galactic dark matter, using modification of pulsar timing residuals when PBHs pass within ~1000 AU and impart impulse accelerations to the Earth. With this technique, PBHs with masses around 10^{25} g (~0.1 lunar mass) can be detected. Currently, the constraints on the abundance of such dark matter candidates are weak. A 30 year-long monitoring campaign with the proposed Square Kilometer Array (SKA) can rule out a PBH fraction more than ~1/10 in the solar neighborhood in the form of dark matter with mass ~10^{25} g.

We show that isotropic component of circular polarization of stochastic gravitational wave background can be explored by breaking two dimensional configuration of multiple laser interferometers for correlation analysis. By appropriately selecting orbital parameters for the proposed BBO mission, the circular polarization degree Pi can be measured down to Pi ~ 0.08 (10^{-15}/Omega_{GW})(SNR/5) with slightly (~10%) sacrificing the detection limit for the total intensity Omega_{GW} compared to the standard plane symmetric configuration. This might allow us to detect signature of parity violation in the very early universe.

We discussed prospects for directly detecting circular polarization signal of gravitational wave background. We found it is generally difficult to probe the monopole mode of the signal due to broad directivity of gravitational wave detectors. But the dipole (l=1) and octupole (l=3) modes of the signal can be measured in a simple manner by combining outputs of two unaligned detectors, and we can dig them deeply under confusion and detector noises. Around f~0.1mHz LISA will provide ideal data streams to detect these anisotropic components whose magnitudes are as small as ~1 percent of the detector noise level in terms of the non-dimensional energy density \Omega_{GW}(f).

We study prospects of a method to constrain the inclination of a coalescing compact binary by detecting its gravitational waves associated with a three-dimensionally localized (direction and distance) short-hard gamma-ray burst. We take advantage of a synergy of these two observations, and our method can be applied even with a single interferometer. For a nearly face-on binary the inclination angle $I$ can be constrained in the range 1-1/SNR < cosI \le 1 (SNR: the signal to noise ratio of gravitational wave detection), provided that the error of the distance estimation is negligible. This method would help us to study properties of the short-hard bursts, including potentially collimated jet-like structures as indicated by recent observation.

We discuss prospects for direct measurement of stochastic gravitational wave background around 0.1-1Hz with future space missions. It is assumed to use correlation analysis technique with the optimal TDI variables for two sets of LISA-type interferometers. The signal to noise for detection of the background and the estimation errors for its basic parameters (amplitude, spectral index) are evaluated for proposed missions.

The divergence of the momentum density field of the large scale structure generates a secondary anisotropy contribution to the Cosmic Microwave Background (CMB). While the effect is best described as a non-linear extension to the well-known integrated Sachs-Wolfe effect, due to mathematical coincidences, the anisotropy contribution is also described as the lensing of the dipole seen in the rest-frame of a moving mass. Given the closeness, there is a remote possibility that local concentrations of mass in the form of the Great Attractor and the Shapley concentration generate large angular scale fluctuations in CMB and could potentially be responsible, at least partly, for some of the low-multipole anomalies in WMAP data. While the local anisotropy contribution peaks at low multipoles, for reasonable models of the mass and velocity distributions associated with local super structures we find that the amplitude of temperature anisotropies is at most at a level of 10$^{-2}$ $\mu$K and is substantially smaller than primordial fluctuations. It is extremely unlikely that the momentum density of local mass concentrations is responsible for any of the large angular scale anomalies in WMAP data.

We introduce a new formalism to describe the polarization signal of galaxy clusters on the whole sky. We show that a sparsely sampled, half-sky map of the cluster polarization signal at $z\sim 1$ would allow to better characterize the very large scale density fluctuations. While the horizon length is smaller in the past, two other competing effects significantly remove the contribution of the small scale fluctuations from the quadrupole polarization pattern at $z\sim 1$. For the standard Lambda-CDM universe with vanishing tensor mode, the quadrupole moment of the temperature anisotropy probed by WMAP is expected to have a ~32% contribution from fluctuations on scales below 6.3h^{-1}Gpc. This percentage would be reduced to ~2% level for the quadrupole moment of polarization pattern at $z\sim 1$. A cluster polarization map at $z \sim 1$ would shed light on the potentially anomalous features of the largest scale structure in the observable universe.

We discuss an astrometric timing effect on data analysis of continuous gravitational waves from rapidly rotating isolated neutron stars. Special attention is directed to the possibility of determining their distances by measuring the curvature of the wave fronts. We predict that if continuous gravitational waves from an unknown neutron star with a stable rotation are detected around 1kHz within 1/3yr by initial LIGO detectors and the ellipticity parameter epsilon is smaller than 10^{-6}, the distance r to the source can be estimated with relative error \Delta r/r of \sim 10% by using the broad band configuration of advanced LIGO detectors over 3 years. By combining the observed amplitude of the waves with the estimated distance, information on the parameter $\epsilon$ can be obtained purely through gravitational wave measurements.

The Laser Interferometer Space Antenna (LISA) is expected to detect N ~ 22 [x10^(+/- 1)] close white dwarf binaries in the Large Magellanic Cloud (LMC) through their gravitational radiation with signal-to-noise ratios greater than 10 in observational durations of 3 years or more. In addition to chirp mass, location on the sky, and other binary parameters, the distance to each binary is an independent parameter that can be extracted from an analysis of gravitational waves from these binaries. Using a sample of binaries, one can establish the mean distance to the LMC as well as the variance of this distance. Assuming no confusion noise at frequencies above 2 mHz, LISA might determine the LMC distance to ~ 4.5 (N/22)^(1/2)% and the line of sight extent of LMC to ~ 15 (N/22)^(1/4)$%, relative to its distance, at the one-sigma confidence. These estimates are competitive to some of the proposed direct geometric techniques to measure LMC distance in future with missions such as SIM and GAIA.

While there are no strong observational constraints on the gravitational wave background across six or more orders of magnitude between 10^{-16} Hz and 10^{-10} Hz and it is difficult to get a constraint below 10^{-12}Hz using objects in our Galaxy, we suggest that the anisotropy pattern of time variation of the redshift related to a sample of high redshift objects can be used to constrain the gravitational wave background around 10^{-12} Hz. Useful observations for the monitoring of an anisotropy signal in a global redshift change include spectroscopic observations of the Ly-alpha forest in absorption towards a sample of quasars, redshifted 21 cm line observations either in absorption or emission towards a sample of neutral HI regions before or during reionization, and high frequency (0.1 Hz to 1 Hz) gravitational wave analysis of a sample of neutron star--neutron star binaries detected with gravitational wave instruments such as the Decihertz Interferometer Gravitational Wave Observatory (DECIGO). The low frequency background can also be constrained by arcsecond-scale anisotropy observations of the CMB. For reasonable observations in the future involving extragalactic sources, we find best limits at the level of \Omega_{GW} < 10^{-5} at a frequency around 10^{-12} Hz while the eventual ultimate limit one cannot beat is \Omega_{GW} < 10^{-11}.

The high sensitivity of upcoming space-based gravitational wave detectors suggests the possibility that if halo dark matter were composed of primordial black holes (PBHs) with mass between $10^{16}$ g and 10$^{20}$ g, the gravitational interaction with detector test masses will lead to a detectable pulse-like signal during the fly-by. For an improved version of the Laser Interferometer Space Antenna with a reduced acceleration noise at the low-end of its frequency spectrum, we find an event rate, with signal-to-noise ratios greater than 5, of $\sim$ a few per decade involving black holes of mass $\sim$ 10$^{17}$ g. The detection rate improves significantly for second generation space based interferometers that are currently envisioned, though these events must be distinguished from those involving perturbations due to near-Earth asteroids. While the presence of primordial black holes below a mass of $\sim$ 10$^{16}$ g is now constrained based on the radiation released during their evaporation, the gravitational wave detectors will extend the study of PBHs to a several orders of magnitude higher masses.

We discuss spatial fluctuations in the gravitational wave background arising from unresolved Galactic binary sources, such as close white dwarf binaries, due to the fact the galactic binary source distribution is anisotropic. We introduce a correlation analysis of the two data streams of the Laser Interferometer Space Antenna (LISA) to extract spherical harmonic coefficients, in an independent manner, of the hexadecapole moment ($l=4$) related to the projected two-dimensional density distribution of the binary source population. The proposed technique complements and improves over previous suggestions in the literature to measure the gravitational wave background anisotropy based on the time modulation of data as LISA orbits around the Sun. Such techniques, however, are restricted only to certain combinations of spherical harmonic coefficients of the galaxy with no ability to separate them individually. With LISA, $m=2,3$ and 4 coefficients of the hexadecapole ($l=4$) can be measured with signal-to-noise ratios at the level of 10 and above in a certain coordinate system. In addition to the hexadecapole coefficients, when combined with the time modulation analysis, the correlation study can also be used, in principle, to measure quadrupole coefficients of the binary distribution.

We study the anisotropies of the Galactic confusion noise background and its effects on LISA data analysis. LISA has two data streams of the gravitational waves signals relevant for low frequency regime. Due to the anisotropies of the background, the matrix for their confusion noises has off-diagonal components and depends strongly on the orientation of the detector plane. We find that the sky-averaged confusion noise level $\sqrt {S(f)}$ could change by a factor of 2 in three months, and would be minimum when the orbital position of LISA is either around the spring or autumn equinox.

The arrival times of gravitational waves and optical light from orbiting binaries provide a mechanism to understand the propagation speed of gravity when compared to that of light or electromagnetic radiation. This is achieved with a measurement of any offset between optically derived orbital phase related to that derived from gravitational wave data, at a specified location of one binary component with respect to the other. Using a sample of close white dwarf binaries (CWDBs) detectable with the Laser Interferometer Space Antenna (LISA) and optical light curve data related to binary eclipses from meter-class telescopes for the same sample, we determine the accuracy to which orbital phase differences can be extracted. We consider an application of these measurements involving a variation to the speed of gravity, when compared to the speed of light, due to a massive graviton. For a subsample of $\sim$ 400 CWDBs with high signal-to-noise gravitational wave and optical data with magnitudes brighter than 25, the combined upper limit on the graviton mass is at the level of $\sim 6 \times 10^{-24}$ eV. This limit is two orders of magnitude better than the present limit derived by Yukawa-correction arguments related to the Newtonian potential and applied to the Solar-system.

We study how the angular resolution of LISA for merging massive black-hole binaries would be improved if we observe multiple gravitational wave ``images'' due to strong gravitational lensing. The correlation between fitting parameters is reduced by the additional information of the second image which significantly reduces the error box on the sky. This improvement would be very helpful for identifying the host galaxy of a binary. The angular resolution expected with multiple detectors is also discussed.

The Laser Interferometer Space Antenna (LISA) is expected to detect close white dwarf binaries (CWDBs) through their gravitational radiation. Around 3000 binaries will be spectrally resolved at frequencies > 3 mHz, and their positions on the sky will be determined to an accuracy ranging from a few tens of arcminutes to a degree or more. Due to the small binary separation, the optical light curves of >~ 30% of these CWDBs are expected to show eclipses, giving a unique signature for identification in follow-up studies of the LISA error boxes. While the precise optical location improves binary parameter determination with LISA data, the optical light curve captures additional physics of the binary, including the individual sizes of the stars in terms of the orbital separation. To optically identify a substantial fraction of CWDBs and thus localize them very accurately, a rapid monitoring campaign is required, capable of imaging a square degree or more in a reasonable time, at intervals of 10--100 seconds, to magnitudes between 20 and 25. While the detectable fraction can be up to many tens of percent of the total resolved LISA CWDBs, the exact fraction is uncertain due to unknowns related to the white dwarf spatial distribution, and potentially interesting physics, such as induced tidal heating of the WDs due to their small orbital separation.

Response of an interferometer becomes complicated for gravitational wave shorter than the arm-length of the detector, as nature of wave appears strongly. We have studied how parameter estimation for merging massive black hole binaries are affected by this complicated effect in the case of LISA. It is shown that three dimensional positions of some binaries might be determined much better than the past estimations that use the long wave approximation. For equal mass binaries this improvement is most prominent at $\sim 10^5\sol$.

We calculate how accurately parameters of the short-period binaries $(10^{-4} {Hz} \lsim f\lsim 10^{-2} {Hz})$ will be determined from the gravitational waves by LISA. In our analysis the chirp signal ${\dot f}$ is newly included as a fitting parameters and dependence on observational period or wave frequency is studied in detail. Implications for gravitational wave astronomy are also discussed quantitatively.

The binary confusion noise spectrum at LISA band depends strongly on observational period and abundance of Galactic close white dwarf binaries (CWDBs). We have investigated how the number of the resolved Galactic CWDBs varies with operation period of LISA, and found that the resolved number would typically grows by a factor of 5 when the operation period increases from 1yr to 10yr. We have also made a similar estimation for number of CWDBs whose chirp signal can be measured in matched filtering analysis.

We propose a new method for determining total masses of low frequency eccentric binaries (such as, neutron star binaries with orbital frequency $f\gsim 10^{-3}$Hz) from their gravitational waves. In this method we use the frequency shift caused by periastron advance, and it works even at low frequency band where chirp signal due to radiation reaction is difficult to be measured. It is shown that the total masses of several Galactic neutron star binaries might be measured accurately (within a few percent error) by LISA with operation period of $\sim 10$ years.

It may be possible to construct a laser interferometer gravitational wave antenna in space with $h_{rms}\sim 10^{-27}$ at $ f\sim 0.1{\rm Hz}$ in this century. We show possible specification of this antenna which we call DECIGO. Using this antenna we show that 1) typically $ 10^5$ ($10^4\sim 10^6$) chirp signals of coalescing binary neutron stars per year may be detected with S/N $\sim 10^4$. 2) We can directly measure the acceleration of the universe by ten years observation of binary neutron stars. 3) The stochastic gravitational waves of $\Omega_{GW}\gsim 10^{-20}$ predicted by the inflation may be detected by correlation analysis for which effects of the recent cosmic acceleration would become highly important. Our formula for phase shift due to accelerating motion might be also applied for binary sources of LISA.

We study evolution of various statistical quantities of smoothed cosmic density and velocity fields using N-body simulations. The parameter $C\equiv <V^2\delta>/(<V^2> <\delta^2>)$ characterizes nonlinear coupling of these two fields and determines behavior of bulk velocity dispersion as a function of local density contrast. It is found that this parameter depends strongly on the smoothing scale even in quasi-linear regimes where the skewness parameter $S_3$ is nearly constant and close to the predicted value by the second-order perturbation theory. We also analyze weakly nonlinear effects caused by an adaptive smoothing known as the gather approach.

Evolution of power spectrum is studied for non-Gaussian models of structure formation. We generalize the dark-matter-approach to these models and find that the evolved spectrum at weakly nonlinear regime is mainly determined by a simple integral of the dark-matter-halo mass function in this approach. We also study the change of the nonlinear spectrum within the current observational constraint of the primordial non-Gaussianity.

We examine the polarization signal of the cosmic microwave background radiation associated with distant clusters. The polarization is induced by the Thomson scattering of microwave photons with ionized gas of clusters and contains information of quadrupole temperature anisotropies observed at the clusters. The three-dimensional map of the signal are expressed in terms of the spin-weighted harmonics for its angular dependence. Its radial dependence is expanded perturbatively with respect to the distances (equivalently redshifts) to the clusters. The independent information that we can extract out from the map is clarified explicitly.

Cosmological linear perturbation theory predicts that the peculiar velocity $V(x)$ and the matter overdensity $\delta(x)$ at a same point $x$ are statistically independent quantities, as log as the initial density fluctuations are random Gaussian distributed. However nonlinear gravitational effects might change the situation. Using framework of second-order perturbation theory and the Edgeworth expansion method, we study local density dependence of bulk velocity dispersion that is coarse-grained at a weakly nonlinear scale. For a typical CDM model, the first nonlinear correction of this constrained bulk velocity dispersion amounts to $\sim 0.3\delta$ (Gaussian smoothing) at a weakly nonlinear scale with a very weak dependence on cosmological parameters. We also compare our analytical prediction with published numerical results given at nonlinear regimes.

Smoothing operation to make continuous density field from observed point-like distribution of galaxies is crucially important for topological or morphological analysis of the large-scale structure, such as, the genus statistics or the area statistics (equivalently the level crossing statistics). It has been pointed out that the adaptive smoothing filters are more efficient tools to resolve cosmic structures than the traditional spatially fixed filters. We study weakly nonlinear effects caused by two representative adaptive methods often used in smoothed hydrodynamical particle (SPH) simulations. Using framework of second-order perturbation theory, we calculate the generalized skewness parameters for the adaptive methods in the case of initially power-law fluctuations. Then we apply the multidimensional Edgeworth expansion method and investigate weakly nonlinear evolution of the genus statistics and the area statistics. Isodensity contour surfaces are often parameterized by the volume fraction of the regions above a given density threshold. We also discuss this parameterization method in perturbative manner.

If primordial fluctuation is Gaussian distributed, higher-order cumulants of the cosmic fields reflect nonlinear mode coupling and provide useful information of gravitational instability picture of structure formation. We show that their expected deviation (sample variance) from the universal values is nonvanishing even in linear theory in the case where observed volume is finite. As a result, we find that the relative sample variance of the skewness of the smoothed velocity divergence field remains as large as $\sim 30%$ even if the survey depth is as deep as $\sim 150h^{1}Mpc$.

Since the cosmic peculiar velocity field depends on small wave-number modes strongly, we cannot probe its universal properties unless we observe a sufficiently large region. We calculate the expected deviation (sample variance) of the peculiar velocity dispersion from its universal value in the case observed volume is finite. Using linear theory we show that the sample variance remains as large as $\sim 10%$, even if the observed region is as deep as $100h^{-1}Mpc$ and that it seriously affects the estimation of cosmological parameters from the peculiar velocity field.

The statistical distribution of the radial pairwise peculiar velocity of galaxies is known to have an exponential form as implied by observations and explicitly shown in N-body simulations. Here we calculate its statistical distribution function using the Zel'dovich approximation assuming that the primordial density fluctuations are Gaussian distributed. We show that the exponential distribution is realized as a transient phenomena on megaparsec scales in the standard cold-dark-matter model.

Evolution of genus density is calculated from Gaussian initial conditions using Zel'dovich approximation. A new approach is introduced which formulates the desired quantity in a rotationally invariant manner. It is shown that normalized genus density does not depend on the initial spectral shape but is a function of the fluctuation amplitude only.