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, without assuming asymptotic flatness, a gravitational lens for an observer and source that are within a finite distance from a lens object. The proposed lens equation is consistent with the deflection angle of light that is defined for nonasymptotic observer and source by Takizawa et al. [Phys. Rev. D 101, 104032 (2020)] based on the Gauss-Bonnet theorem with using the optical metric. This lens equation, though it is shown to be equivalent to the lens equation by Bozza [Phys. Rev. D 78, 103005 (2008)], is linear in the deflection angle. Therefore, the proposed equation is more convenient for the purpose of doing an iterative analysis. As an explicit example of an asymptotically nonflat spacetime, we consider a static and spherically symmetric solution in Weyl conformal gravity, especially a case that $\gamma$ parameter in the Weyl gravity model is of the order of the inverse of the present Hubble radius. For this case, we examine iterative solutions for the finite-distance lens equation up to the third order. The effect of the Weyl gravity on the lensed image position begins at the third order and it is linear in the impact parameter of light. The deviation of the lensed image position from the general relativistic one is $\sim 10^{-2}$ microarcsecond for the lens and source with a separation angle of $\sim 1$ arcminute, where we consider a cluster of galaxies with $10^{14} M_{\odot}$ at $\sim 1$ Gpc for instance. The deviation becomes $\sim 10^{-1}$ microarcseconds, even if the separation angle is $\sim 10$ arcminutes. Therefore, effects of the Weyl gravity model are negligible in current and near-future observations of gravitational lensing. On the other hand, the general relativistic corrections at the third order $\sim 0.1$ milliarcseconds can be relevant with VLBI observations.

The gravitational deflection angle of light for an observer and source at finite distance from a lens object has been studied by Ishihara et al. [Phys. Rev. D, 94, 084015 (2016)], based on the Gauss-Bonnet theorem with using the optical metric. Their approach to finite-distance cases is limited within an asymptotically flat spacetime. By making several assumptions, we give an interpretation of their definition from the observer's viewpoint: The observer assumes the direction of a hypothetical light emission at the observer position and makes a comparison between the fiducial emission direction and the direction along the real light ray. The angle between the two directions at the observer location can be interpreted as the deflection angle by Ishihara et al. The present interpretation does not require the asymptotic flatness. Motivated by this, we avoid such asymptotic regions to discuss another integral form of the deflection angle of light. This form makes it clear that the proposed deflection angle can be used not only for asymptotically flat spacetimes but also for asymptotically nonflat ones. We examine the proposed deflection angle in two models for the latter case; Kottler (Schwarzschild-de Sitter) solution in general relativity and a spherical solution in Weyl conformal gravity. Effects of finite distance on the light deflection in Weyl conformal gravity result in an extra term in the deflection angle, which may be marginally observable in a certain parameter region. On the other hand, those in Kottler spacetime are beyond reach of the current technology.

General metric theories in a four-dimensional spacetime allow at most six polarization states (two spin-0, two spin-1 and two spin-2) of gravitational waves (GWs). If a sky location of a GW source with the electromagnetic counterpart satisfies a single equation that we propose in this paper, both the spin-1 modes and spin-2 ones can be eliminated from a certain combination of strain outputs at four ground-based GW interferometers (e.g. a network of aLIGO-Hanford, aLIGO-Livingston, Virgo and KAGRA), where this equation describes curves on the celestial sphere. This means that, if a GW source is found in the curve (or its neighborhood practically), a direct test of scalar (spin-0) modes separately from the other (vector and tensor) modes become possible in principle. The possibility of such a direct test is thus higher than an earlier expectation (Hagihara et al. PRD, 100, 064010, 2019), in which they argued that the vector modes could not be completely eliminated. We discuss also that adding the planned LIGO-India detector as a fifth detector will increase the feasibility of scalar polarization tests.

General metric theories of gravity in four-dimensional spacetimes can contain at most six polarization modes (two spin-0, two spin-1 and two spin-2) of gravitational waves (GWs). It has been recently shown that, with using four GW non-coaligned detectors, a direct test of the spin-1 modes can be done in principle separately from the spin-0 and spin-2 modes for a GW source in particular sky positions [Hagihara et al., Phys. Rev. D 98, 064035 (2018)]. They have found particular sky positions that satisfy a condition of killing completely the spin-0 modes in a so-called null stream which is a linear combination of the signal outputs to kill the spin-2 modes. The present paper expands the method to discuss a possibility that the spin-0 modes are not completely but effectively suppressed in the null streams to test the spin-1 modes separately from the other modes, especially with an expected network of Advanced LIGO, Advanced Virgo and KAGRA. We study also a possibility that the spin-1 modes are substantially suppressed in the null streams to test the spin-0 modes separately from the other modes, though the spin-1 modes for any sky position cannot be completely killed in the null streams. Moreover, we find that the coefficient of the spin-0 modes in the null stream is significantly small for the GW170817 event, so that an upper bound can be placed on the amplitude of the spin-1 modes as $< 6 \times 10^{-23}$.

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.

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.

By using a method improved with a generalized optical metric, the deflection of light for an observer and source at finite distance from a lens object in a stationary, axisymmetric and asymptotically flat spacetime has been recently discussed [Ono, Ishihara, Asada, Phys. Rev. D 96, 104037 (2017)]. By using this method, in the weak field approximation, we study the deflection angle of light for an observer and source at finite distance from a rotating Teo wormhole, especially by taking account of the contribution from the geodesic curvature of the light ray in a space associated with the generalized optical metric. Our result of the deflection angle of light is compared with a recent work on the same wormhole but limited within the asymptotic source and observer [Jusufi, Ovgun, Phys. Rev. D 97, 024042, (2018)], in which they employ another approach proposed by Werner with using the Nazim's osculating Riemannian construction method via the Randers-Finsler metric. We show that the two different methods give the same result in the asymptotic limit. We obtain also the corrections to the deflection angle due to the finite distance from the rotating wormhole.

Assuming that, for a given source of gravitational waves (GWs), we know its sky position, as is the case of GW events with an electromagnetic counterpart such as GW170817, we discuss a null stream method to probe GW polarizations including spin-0 (scalar) GW modes and spin-1 (vector) modes, especially with an expected network of Advanced LIGO, Advanced Virgo and KAGRA. For two independent null streams for four non-co-aligned GW detectors, we study a location on the sky, exactly at which the spin-0 modes of GWs vanish in any null stream for the GW detector network, though the strain output at a detector may contain the spin-0 modes. Our numerical calculations show that there exist seventy sky positions that satisfy this condition of killing the spin-0 modes in the null streams. If a GW source with an electromagnetic counterpart is found in one of the seventy sky positions, the spin-1 modes will be testable separately from the spin-0 modes by the null stream method. In addition, we study a superposition of the two null streams to show that any one of the three modes (one combined spin-0 and two spin-1 modes) can be eliminated by suitably adjusting a weighted superposition of the null streams and thereby a set of the remaining polarization modes can be experimentally tested.

Continuing work initiated in an earlier publication [H. Asada, Phys. Rev. D {\bf 80}, 064021 (2009)], the gravitational radiation reaction to Lagrange's equilateral triangular solution of the three-body problem is investigated in an analytic method. The previous work is based on the energy balance argument, which is sufficient for a two-body system because the number of degrees of freedom (the semimajor axis and the eccentricity in quasi-Keplerian cases, for instance) equals that of the constants of motion such as the total energy and the orbital angular momentum. In a system with three (or more) bodies, however, the number of degrees of freedom is more than that of the constants of motion. Therefore, the present paper discusses the evolution of the triangular system by directly treating the gravitational radiation reaction force to each body. The perturbed equations of motion are solved by using the Laplace transform technique. It is found that the triangular configuration is adiabatically shrinking and is kept in equilibrium by increasing the orbital frequency due to the radiation reaction if the mass ratios satisfy the Newtonian stability condition. Long-term stability involving the first post-Newtonian corrections is also discussed.

Major construction and initial-phase operation of a second-generation gravitational-wave detector KAGRA has been completed. The entire 3-km detector is installed underground in a mine in order to be isolated from background seismic vibrations on the surface. This allows us to achieve a good sensitivity at low frequencies and high stability of the detector. Bare-bones equipment for the interferometer operation has been installed and the first test run was accomplished in March and April of 2016 with a rather simple configuration. The initial configuration of KAGRA is named {\it iKAGRA}. In this paper, we summarize the construction of KAGRA, including the study of the advantages and challenges of building an underground detector and the operation of the iKAGRA interferometer together with the geophysics interferometer that has been constructed in the same tunnel.

By using the Gauss-Bonnet theorem, the bending angle of light in a static, spherically symmetric and asymptotically flat spacetime has been recently discussed, especially by taking account of the finite distance from a lens object to a light source and a receiver [Ishihara, Suzuki, Ono, Asada, Phys. Rev. D 95, 044017 (2017)]. We discuss a possible extension of the method of calculating the bending angle of light to stationary, axisymmetric and asymptotically flat spacetimes. For this purpose, we consider the light rays on the equatorial plane in the axisymmetric spacetime. We introduce a spatial metric to define the bending angle of light in the finite-distance situation. We show that the proposed bending angle of light is coordinate-invariant by using the Gauss-Bonnet theorem. The non-vanishing geodesic curvature of the photon orbit with the spatial metric is caused in gravitomagnetism, even though the light ray in the four-dimensional spacetime follows the null geodesic. Finally, we consider Kerr spacetime as an example in order to examine how the bending angle of light is computed by the present method. The finite-distance correction to the gravitomagnetic deflection angle due to the Sun's spin is around a pico-arcsecond level. The finite-distance corrections for Sgr A$^{\ast}$ also are estimated to be very small. Therefore, the gravitomagnetic finite-distance corrections for these objects are unlikely to be observed with present technology.

This paper reviews a phenomenological approach to the gravitational lensing by exotic objects such as the Ellis wormhole lens, where exotic lens objects may follow a non-standard form of the equation of state or may obey a modified gravity theory. A gravitational lens model is proposed in the inverse powers of the distance, such that the Schwarzschild lens and exotic lenses can be described in a unified manner as a one parameter family. As observational implications, the magnification, shear, photo-centroid motion and time delay in this lens model are discussed.

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.

Continuing work initiated in an earlier publication [Ishihara, Suzuki, Ono, Kitamura, Asada, Phys. Rev. D {\bf 94}, 084015 (2016) ], we discuss a method of calculating the bending angle of light in a static, spherically symmetric and asymptotically flat spacetime, especially by taking account of the finite distance from a lens object to a light source and a receiver. For this purpose, we use the Gauss-Bonnet theorem to define the bending angle of light, such that the definition can be valid also in the strong deflection limit. Finally, this method is applied to Schwarzschild spacetime in order to discuss also possible observational implications. The proposed corrections for Sgr A$^{\ast}$ for instance are able to amount to $\sim 10^{-5}$ arcseconds for some parameter range, which may be within the capability of near-future astronomy, while also the correction for the Sun in the weak field limit is $\sim 10^{-5}$ arcseconds.

We discuss a possible extension of calculations of the bending angle of light in a static, spherically symmetric and asymptotically flat spacetime to a non-asymptotically flat case. We examine a relation between the bending angle of light and the Gauss-Bonnet theorem by using the optical metric. A correspondence between the deflection angle of light and the surface integral of the Gaussian curvature may allow us to take account of the finite distance from a lens object to a light source and a receiver. Using this relation, we propose a method for calculating the bending angle of light for such cases. Finally, this method is applied to two examples of the non-asymptotically flat spacetimes to suggest finite-distance corrections: Kottler (Schwarzschild-de Sitter) solution to the Einstein equation and an exact solution in Weyl conformal gravity.

We study linear nonradial perturbations and stability of a marginal stable circular orbit (MSCO) such as the innermost stable circular orbit (ISCO) of a test particle in stationary axisymmetric spacetimes which possess a reflection symmetry with respect to the equatorial plane. The proposed approach is applied to Kerr solution and Majumdar-Papapetrou solution to Einstein equation. Finally, we reexamine MSCOs for a modified metric of a rapidly spinning black hole that has been recently proposed by Johannsen and Psaltis [PRD, 83, 124015 (2011)]. We show that, for the Johannsen and Psaltis's model, circular orbits that are stable against radial perturbations for some parameter region become unstable against vertical perturbations. This suggests that the last circular orbit for this model may be larger than the ISCO.

In terms of Sturm's theorem, we reexamine a marginal stable circular orbit (MSCO) such as the innermost stable circular orbit (ISCO) of a timelike geodesic in any spherically symmetric and static spacetime. MSCOs for some of exact solutions to the Einstein's equation are discussed. Strum's theorem is explicitly applied to the Kottler (often called Schwarzschild-de Sitter) spacetime. Moreover, we analyze MSCOs for a spherically symmetric, static and vacuum solution in Weyl conformal gravity.

We study a marginally stable circular orbit (MSCO) such as the innermost stable circular orbit (ISCO) of a timelike geodesic in any spherically symmetric and static spacetime. It turns out that the metric components are separable from the constants of motion along geodesics. We show also that a metric component $g_{rr}$ with a radial coordinate $r$ does not affect MSCOs. This suggests that, as a test of gravity, any ISCO measurement may be put into the same category as gravitational redshift experiments. MSCOs for exact solutions to the Einstein's equation are also mentioned.

Continuing work initiated in earlier publications [Ichita, Yamada and Asada, Phys. Rev. D {\bf 83}, 084026 (2011); Yamada and Asada, Phys. Rev. D {\bf 86}, 124029 (2012)], we examine the post-Newtonian (PN) effects on the stability of the triangular solution in the relativistic three-body problem for general masses. For three finite masses, a condition for stability of the triangular solution is obtained at the first post-Newtonian (1PN) order, and it recovers previous results for the PN restricted three-body problem when one mass goes to zero. The stability regions still exist even at the 1PN order, though the PN triangular configuration for general masses is less stable than the PN restricted three-body case as well as the Newtonian one.

The gravitational lensing effects in the weak gravitational field by exotic lenses have been investigated intensively to find nonluminous exotic objects. Gravitational lensing based on 1/r^n fall-off metric, as a one-parameter model that can treat by hand both the Schwarzschild lens (n=1) and the Ellis wormhole (n=2) in the weak field, has been recently studied. Only for n=1 case, however, it has been explicitly shown that effects of relativistic lens images by the strong field on the light curve can be neglected. We discuss whether relativistic images by the strong field can be neglected for n>1 in the Tangherlini spacetime which is one of the simplest models for our purpose. We calculate the divergent part of the deflection angle for arbitrary n and the regular part for n=1, 2 and 4 in the strong field limit, the deflection angle for arbitrary n under the weak gravitational approximation. We also compare the radius of the Einstein ring with the radii of the relativistic Einstein rings for arbitrary n. We conclude that the images in the strong gravitational field have little effect on the total light curve and that the time-symmetric demagnification parts in the light curve will appear even after taking account of the images in the strong gravitational field for n>1.

Toward a test of parity violation in a gravity theory, possible effects of Chern-Simons (CS) gravity on an interferometer have been recently discussed. Continuing work initiated in an earlier publication [Okawara, Yamada and Asada, Phys. Rev. Lett. 109, 231101 (2012)], we study possible altitudinal and directional dependence of relativistic Sagnac effect in CS modified gravity. We compare the CS effects on Sagnac interferometers with the general relativistic Lense-Thirring (LT) effects. Numerical calculations show that the eastbound Sagnac interferometer might be preferred for testing CS separately, because LT effects on this interferometer cancel out. The size of the phase shift induced in the CS model might have an oscillatory dependence also on the altitude of the interferometer through the CS mass parameter $m_{CS}$. Therefore, the international space station site as well as a ground-based experiment is also discussed.

A moment approach for orbit determinations of astrometric binaries from astrometric observations alone has been recently studied for a low signal-to-noise ratio (Iwama et al. 2013, PASJ, 65, 2). With avoiding a direct use of the time-consuming Kepler equation, temporal information is taken into account to increase the accuracy of statistical moments. As numerical tests, 100 realizations are done and the mean and the standard deviation are also evaluated. For a semi-major axis, the difference between the mean of the recovered values and the true value decreases to less than a tenth in the case of $10000$ observed points. Therefore, the present moment approach works better than the previous one for the orbit determinations when one has a number of the observed points. The present approach is thus applicable to Cyg X-1.

Gravitational lens models, some of which might act as if a concave lens, have been recently investigated by using a static and spherically symmetric modified spacetime metric that depends on the inverse distance to the $n$-th power [Kitamura, Nakajima and Asada, PRD 87, 027501 (2013)]. We reexamine the time delay of light in a gravitational concave lens as well as a gravitational convex one. The frequency shift due to the time delay is also investigated. We show that the sign of the time delay in the lens models is the same as that of the deflection angle of light. The size of the time delay decreases with increase in the parameter $n$. We discuss also possible parameter ranges that are relevant to pulsar timing measurements in our galaxy.

Gravitational lens models with negative convergence inspired by modified gravity theories, exotic matter and energy have been recently examined, in such a way that a static and spherically symmetric modified spacetime metric depends on the inverse distance to the $n$-th power ($n=1$ for Schwarzschild metric, $n=2$ for Ellis wormhole, and $n \neq 1$ for an extended spherical distribution of matter such as an isothermal sphere) in the weak-field approximation [Kitamura, Nakajima and Asada, PRD 87, 027501 (2013), Izumi et al. PRD 88 024049 (2013)]. Some of the models act as if a convex lens, whereas the others are repulsive on light rays like a concave lens. The present paper considers microlensed image centroid motions by the exotic lens models. Numerical calculations show that, for large $n$ cases in the convex-type models, the centroid shift from the source position might move on a multiply-connected curve like a bow tie, while it is known to move on an ellipse for $n=1$ case and to move on an oval curve for $n=2$. The distinctive feature of the microlensed image centroid may be used for searching (or constraining) localized exotic matter or energy with astrometric observations. It is shown also that the centroid shift trajectory for concave-type repulsive models might be elongated vertically to the source motion direction like a prolate spheroid, whereas that for convex-type models such as the Schwarzschild one is tangentially elongated like an oblate spheroid.

Gravitational lens models with negative convergence(surface mass density projected onto the lens plane) inspired by modified gravity theories, exotic matter and energy have been recently discussed in such a way that a static and spherically-symmetric modified spacetime metric depends on the inverse distance to the power of positive $n$(n=1 for Schwarzschild metric, n=2 for Ellis wormhole) in the weak-field approximation [Kitamura, Nakajima and Asada, PRD 87, 027501 (2013)], and it has been shown that demagnification of images could occur for $n>1$ lens models associated with exotic matter (and energy), though they cause the gravitational pull on light rays. The present paper considers gravitational lensing shear by the demagnifying lens models and other models such as negative-mass compact objects causing the gravitational repulsion on light rays like a concave lens. It is shown that images by the lens models for the gravitational pull are tangentially elongated, whereas those by the repulsive ones are radially distorted. This feature of lensed image shapes may be used for searching(or constraining) localized exotic matter or energy with gravitational lensing surveys. It is suggested also that an underdense region such as a cosmic void might produce radially elongated images of background galaxies rather than tangential ones.

The latest result in the Sloan Digital Sky Survey Quasar Lens Search (SQLS) has set the first cosmological constraints on negative-mass compact objects and Ellis wormholes. There are no multiple images lensed by the above two exotic objects for $\sim 50000$ distant quasars in the SQLS data. Therefore, an upper bound is put on the cosmic abundances of these lenses. The number density of negative mass compact objects is $n<10^{-8} (10^{-4}) h^3 {\rm Mpc}^{-3}$ at the mass scale $|M| > 10^{15} (10^{12}) M_\odot$, which corresponds to the cosmological density parameter $|\Omega| < 10^{-4}$ at the galaxy-scale mass range $|M|=10^{12-15}M_\odot$. The number density of the Ellis wormhole is $n<10^{-4} h^3 {\rm Mpc}^{-3}$ for a range of the throat radius $a = 10^{1-4}$pc, which is much smaller than the Einstein ring radius.

It has been recently suggested that possible effects of Chern-Simons gravity on a quantum interferometer are dependent on the latitude and direction of the interferometer on Earth in orbital motion around Sun. Continuing work initiated in the earlier publication [Okawara, Yamada and Asada, Phys. Rev. Lett. 109, 231101 (2012)], we perform numerical calculations of time variation in the induced phase shifts for nonequatorial cases. We show that the maximum phase shift at any latitude might occur at 6, 0 (and 12), and 18 hours (in local time) of each day, when the normal vector to the interferometer is vertical, eastbound and northbound, respectively. If two identical interferometers were located at different latitudes, the difference between two phase shifts that are measured at the same local time would be $O(\sin \delta\varphi)$ for a small latitude difference $\delta\varphi$. It might thus become maximally $\sim 20$ percents for $\delta\varphi \sim 10$ degrees, for instance.

We examine a gravitational lens model inspired by modified gravity theories and exotic matter and energy. We study an asymptotically flat, static, and spherically symmetric spacetime that is modified in such a way that the spacetime metric depends on the inverse distance to the power of positive $n$ in the weak-field approximation. It is shown analytically and numerically that there is a lower limit on the source angular displacement from the lens object to get demagnification. Demagnifying gravitational lenses could appear, provided the source position $\beta$ and the power $n$ satisfy $\beta > 2/(n+1)$ in the units of the Einstein ring radius under a large-$n$ approximation. Unusually, the total amplification of the lensed images, though they are caused by the gravitational pull, could be less than unity. Therefore, time-symmetric demagnification parts in numerical light curves by gravitational microlensing (F.Abe, Astrophys. J. 725, 787, 2010) may be evidence of an Ellis wormhole (being an example of traversable wormholes), but they do not always prove it. Such a gravitational demagnification of the light might be used for hunting a clue of exotic matter and energy that are described by an equation of state more general than the Ellis wormhole case. Numerical calculations for the $n=3$ and 10 cases show maximally $\sim 10$ and $\sim 60$ percent depletion of the light, when the source position is $\beta \sim 1.1$ and $\beta \sim 0.7$, respectively.

Continuing work initiated in an earlier publication [Ichita, Yamada and Asada, Phys. Rev. D 83, 084026 (2011)], we reexamine the post-Newtonian effects on Lagrange's equilateral triangular solution for the three-body problem. For three finite masses, it is found that a triangular configuration satisfies the post-Newtonian equation of motion in general relativity, if and only if it has the relativistic corrections to each side length. This post-Newtonian configuration for three finite masses is not always equilateral and it recovers previous results for the restricted three-body problem when one mass goes to zero. For the same masses and angular velocity, the post-Newtonian triangular configuration is always smaller than the Newtonian one.

A moment approach for orbit determinations of an astrometric binary with low signal-to-noise ratio from astrometric observations alone is proposed, especially aiming at a close binary system with a short orbital period such as Cyg-X1 and also at a star wobbled by planets. As an exact solution to the nonlinearly coupled equation system, the orbital elements are written in terms of the second and third moments of projected positions that are measured by astrometry. This may give a possible estimation of the true orbit.

We consider a coplanar system comprised of a massive central body (a star), a less massive secondary (a planet) on a circular orbit, and a test particle on a bound orbit exterior to that of the secondary. The gravitational pull exerted on the test particle by the secondary acts as a small perturbation, wherefore the trajectory of the particle can be described as an ellipse of a precessing perihelion. While the apsidal motion is defined overwhelmingly by the Newtonian portion of the secondary's gravity, the post-Newtonian portion, too, brings its tiny input. We explore whether this input may be of any astrophysical relevance in the next few decades. We demonstrate that the overall post-Newtonian input of the secondary's gravity can be split into two parts. One can be expressed via the orbital angular momentum of the secondary, another via its orbital radius. Despite some moderately large numerical factors showing up in the expressions for these two parts, the resulting post-Newtonian contributions from the secondary's gravity into the apsidal motion of the test particle turn out to be small enough to be neglected in the near-future measurements.

We reexamine the light deflection by an Ellis wormhole. The bending angle as a function of the ratio between the impact parameter and the throat radius of the wormhole is obtained in terms of a complete elliptic integral of the first kind. This result immediately yields asymptotic expressions in the weak field approximation. It is shown that an expression for the deflection angle derived (and used) in recent papers is valid at the leading order but it breaks down at the next order because of the nontrivial spacetime topology.

The relativistic perihelion precession due to the three-body interaction is derived. We consider a hierarchical coplanar three-body system, such as the Sun, Jupiter and Saturn, in which both the secondary object as the largest planet corresponding to Jupiter (mass $m_2$) and the third one corresponding to Saturn (mass $m_3$) orbit around the primary object corresponding to Sun (mass $m_1 \gg m_2 \gg m_3$), where the mean orbital radius of the third body is larger than that of the secondary one (denoted as $\ell$). We investigate the post-Newtonian effects on the motion of the third body (semimajor axis a, eccentricity e for the Keplerian orbital elements). Under some assumptions with a certain averaging, the relativistic perihelion precession of the third mass by the post-Newtonian three-body interaction is expressed as $6 G m_2 \ell^2 c^{-2} a^{-3} n (1+9e^2/16) (1-e^2)^{-3} $, where G and c denote the gravitational constant and the speed of light, respectively, and the mean motion for the third body is denoted as $n = 2\pi a^{3/2} G^{-1/2} (m_1+m_2)^{-1/2}$. For the Sun-Jupiter-Saturn system, it is $7.8 \times 10^{-6}$ arcsec/cy. This is larger than the Lense-Thirring effect by Sun but it cannot yet explain the recently reported value for the anomalous perihelion precession of Saturn as $-0.006 \pm 0.002$ arcsec/cy by Iorio (2009) based on the analyses by Pitjeva with the EPM2008 ephemerides.

Continuing work initiated in an earlier publication (Asada, MNRAS. 394 (2009) 818), we make a systematic attempt to determine, as a function of lens and source parameters, the positions of images by multi-plane gravitational lenses. By extending the previous single-plane work, we present a method of Taylor-series expansion to solve the multi-plane lens equation in terms of mass ratios except for the neighborhood of the caustics. The advantage of this method is that it allows a systematic iterative analysis and clarifies the dependence on lens and source parameters. In concordance with the multi-plane lensed-image counting theorem that the lower bound on the image number is $2^N$ for N planes with a single point mass on each plane, our iterative results show how $2^N$ images are realized. Numerical tests are done to investigate if the Taylor expansion method is robust. The method with a small mass ratio works well for changing a plane separation, whereas it breaks down in the inner domain near the caustics.

Continuing work initiated in an earlier publication (Abe, ApJ, 725 (2010) 787), we study the gravitational microlensing effects of the Ellis wormhole in the weak-field limit. First, we find a suitable coordinate transformation, such that the lens equation and analytic expressions of the lensed image positions can become much simpler than the previous ones. Second, we prove that two images always appear for the weak-field lens by the Ellis wormhole. By using these analytic results, we discuss astrometric image centroid displacements due to gravitational microlensing by the Ellis wormhole. The astrometric image centroid trajectory by the Ellis wormhole is different from the standard one by a spherical lensing object that is expressed by the Schwarzschild metric. The anomalous shift of the image centroid by the Ellis wormhole lens is smaller than that by the Schwarzschild lens, provided that the impact parameter and the Einstein ring radius are the same. Therefore, the lensed image centroid by the Ellis wormhole moves slower. Such a difference, though it is very small, will be in principle applicable for detecting or constraining the Ellis wormhole by using future high-precision astrometry observations. In particular, the image centroid position gives us an additional information, so that the parameter degeneracy existing in photometric microlensing can be partially broken. The anomalous shift reaches the order of a few micro arcsec. if our galaxy hosts a wormhole with throat radius larger than $10^5$ km. When the source moves tangentially to the Einstein ring for instance, the maximum position shift of the image centroid by the Ellis wormhole is 0.18 normalized by the Einstein ring radius. For the same source trajectory, the maximum difference between the centroid displacement by the Ellis wormhole lens and that by the Schwarzschild one is -0.16 in the units of the Einstein radius.

Continuing work initiated in earlier publications [Yamada, Asada, Phys. Rev. D 82, 104019 (2010), 83, 024040 (2011)], we investigate the post-Newtonian effects on Lagrange's equilateral triangular solution for the three-body problem. For three finite masses, it is found that the equilateral triangular configuration satisfies the post-Newtonian equation of motion in general relativity, if and only if all three masses are equal. When a test mass is included, the equilateral configuration is possible for two cases: (1) one mass is finite and the other two are zero, or (2) two of the masses are finite and equal, and the third one is zero, namely a symmetric binary with a test mass. The angular velocity of the post-Newtonian equilateral triangular configuration is always smaller than the Newtonian one, provided that the masses and the side length are the same.

Continuing work initiated in an earlier publication [Yamada, Asada, Phys. Rev. D 82, 104019 (2010)], we investigate collinear solutions to the general relativistic three-body problem. We prove the uniqueness of the configuration for given system parameters (the masses and the end-to-end length). First, we show that the equation determining the distance ratio among the three masses, which has been obtained as a seventh-order polynomial in the previous paper, has at most three positive roots, which apparently provide three cases of the distance ratio. It is found, however, that, even for such cases, there exists one physically reasonable root and only one, because the remaining two positive roots do not satisfy the slow motion assumption in the post-Newtonian approximation and are thus discarded. This means that, especially for the restricted three-body problem, exactly three positions of a third body are true even at the post-Newtonian order. They are relativistic counterparts of the Newtonian Lagrange points L1, L2 and L3. We show also that, for the same masses and full length, the angular velocity of the post-Newtonian collinear configuration is smaller than that for the Newtonian case. Provided that the masses and angular rate are fixed, the relativistic end-to-end length is shorter than the Newtonian one.

Microlensing studies are usually based on the lens equation that is valid only to the first order in the gravitational constant G and lens mass M. We consider corrections to the conventional lens equation in terms of differentiable functions, so that they can express not only the second-order effects of GM in general relativity but also modified gravity theories. As a generalization of Ebina et al. (Prog. Theor. Phys. 104 (2000) 1317), we show that, provided that the spacetime is static, spherically symmetric and asymptotically flat, the total amplification by microlensing remains unchanged at the linear order of the correction to the deflection angle, if and only if the correction takes a particular form as the inverse square of the impact parameter, whereas the magnification factor for each image is corrected. It is concluded that the light curve shape by microlensing is inevitably changed and will thus allow us to probe modified gravity, unless a modification to the deflection angle takes the particular form. No systematic deviation in microlensing observations has been reported. For instance, therefore, the Yukawa-type correction is constrained as the characteristic length > $10^{14}$ m.

The three-body problem is reexamined in the framework of general relativity. The Newtonian three-body problem admits Euler's collinear solution, where three bodies move around the common center of mass with the same orbital period and always line up. The solution is unstable. Hence it is unlikely that such a simple configuration would exist owing to general relativistic forces dependent not only on the masses but also on the velocity of each body. However, we show that the collinear solution remains true with a correction to the spatial separation between masses. Relativistic corrections to the Sun-Jupiter Lagrange points L1, L2 and L3 are also evaluated.

Continuing work initiated in an earlier publication [Sato and Asada, PASJ, 61, L29 (2009)], we consider light curves influenced by the orbital inclination and eccentricity of a companion in orbit around a transiting extrasolar planet (in a planet-satellite system or a hypothetical true binary). We show that the semimajor axis, eccentricity and inclination angle of a `moon' orbit around the host planet can be determined by transit method alone. For this purpose, we present a formulation for the parameter determinations in a small-eccentricity approximation as well as in the exact form. As a result, the semimajor axis is expressed in terms of observables such as brightness changes, transit durations and intervals in light curves. We discuss also a narrow region of parameters that produce a mutual transit by an extrasolar satellite.

Continuing work initiated in an earlier publication [Torigoe et al. Phys. Rev. Lett. {\bf 102}, 251101 (2009)], gravitational wave forms for a three-body system in Lagrange's orbit are considered especially in an analytic method. First, we derive an expression of the three-body wave forms at the mass quadrupole, octupole and current quadrupole orders. By using the expressions, we solve a gravitational-wave {\it inverse} problem of determining the source parameters to this particular configuration (three masses, a distance of the source to an observer, and the orbital inclination angle to the line of sight) through observations of the gravitational wave forms alone. For this purpose, the chirp mass to a three-body system in the particular configuration is expressed in terms of only the mass ratios by deleting initial angle positions. We discuss also whether and how a binary source can be distinguished from a three-body system in Lagrange's orbit or others.

We consider the effects of mutual transits by extrasolar planet-companion systems (in a true binary or a planet-satellite system) on light curves. We show that induced changes in light curves depend strongly on a ratio between a planet-companion's orbital velocity around their host star and a planet-companion's spin speed around their common center of mass. In both the slow and fast spin cases (corresponding to long and short distances between them, respectively), a certain asymmetry appears in light curves. We show that, especially in the case of short distances, occultation of one faint object by the other, while the transit of the planet-companion system occurs in front of its parent star, causes an apparent increase in light curves and characteristic fluctuations appear as important evidence of mutual transits. We show also that extrasolar mutual transits provide a complementary method of measuring the radii of two transiting objects, their separation and mass, and consequently identifying them as a true binary, planet-satellite system or others. Monitoring $10^5$ stars for three years with Kepler may lead to a discovery of a second Earth-Moon-like system if the fraction of such systems for an averaged star is larger than 0.05, or it may put upper limits on the fraction as f < 0.05.

Different numbers of self-gravitating particles (in different types of periodic motion) are most likely to generate very different shapes of gravitational waves, some of which, however, can be accidentally almost the same. One such example is a binary and a three-body system for Lagrange's solution. To track the evolution of these similar waveforms, we define a chirp mass to the triple system. Thereby, we show that the quadrupole waveforms cannot distinguish the sources. It is suggested that waveforms with higher $\ell$-th multipoles will be important for classification of them (with a conjecture of $\ell \leq N$ for N particles).

As an extension of four point mass lenses at the vertices of a rhombus, we present five point mass lenses at the center and vertices of a diamond, which constitute, for a source behind the center, a soluble model giving expressions of all the image positions (with the maximum number of images as twenty for five lenses). For a source near the center, all the image positions are obtained in the linear approximation.

This paper makes the first systematic attempt to determine using perturbation theory the positions of images by gravitational lensing due to arbitrary number of coplanar masses without any symmetry on a plane, as a function of lens and source parameters. We present a method of Taylor-series expansion to solve the lens equation under a small mass-ratio approximation. First, we investigate perturbative structures of a single-complex-variable polynomial, which has been commonly used. Perturbative roots are found. Some roots represent positions of lensed images, while the others are unphysical because they do not satisfy the lens equation. This is consistent with a fact that the degree of the polynomial, namely the number of zeros, exceeds the maximum number of lensed images if N=3 (or more). The theorem never tells which roots are physical (or unphysical). In this paper, unphysical ones are identified. Secondly, to avoid unphysical roots, we re-examine the lens equation. The advantage of our method is that it allows a systematic iterative analysis. We determine image positions for binary lens systems up to the third order in mass ratios and for arbitrary N point masses up to the second order. This clarifies the dependence on parameters. Thirdly, the number of the images that admit a small mass-ratio limit is less than the maximum number. It is suggested that positions of extra images could not be expressed as Maclaurin series in mass ratios. Magnifications are finally discussed.

Simplified solutions to determine binary elements by astrometry were obtained in terms of elementary functions (Asada et al. 2004), and therefore require neither iterative nor numerical methods. In the framework of the simplified solution, this paper discusses the remaining two parameters of the time of periastron passage and the longitude of ascending node in order to complete the solution. We thus clarify a difference between the simplified solution and other analytical methods.

We reexamined the gravitational time delay of light, allowing for various models of modified gravity. We clarify the dependence of the time delay (and induced frequency shift) on modified gravity models and investigate how to distinguish those models, when light propagates in static spherically symmetric spacetimes. Thus experiments by radio signal from spacecrafts at very different distances from Sun and future space-borne laser interferometric detectors could be a probe of modified gravity in the solar system.

We revisit the three-body problem in the framework of general relativity. The Newtonian N-body problem admits choreographic solutions, where a solution is called choreographic if every massive particles move periodically in a single closed orbit. One is a stable figure-eight orbit for a three-body system, which was found first by Moore (1993) and re-discovered with its existence proof by Chenciner and Montgomery (2000). In general relativity, however, the periastron shift prohibits a binary system from orbiting in a single closed curve. Therefore, it is unclear whether general relativistic effects admit a choreographic solution such as the figure eight. We carefully examine general relativistic corrections to initial conditions so that an orbit for a three-body system can be closed and a figure eight. This solution is still choreographic. This illustration suggests that the general relativistic N-body problem also may admit a certain class of choreographic solutions.

None of N-body gravitating systems have been considered to emit periodic gravitational waves because of their chaotic orbits when N=3 (or more). We employ a figure-eight orbit as a specific model for a 3-body system in order to illustrate that some of triple stars are capable of generating periodic waves. This illustration would imply that a certain class of N-body gravitating systems may be relevant to the gravitational waves generation. We show also that the total angular momentum of this 3-body system is not carried away by gravitational waves. A waveform generated by this system is volcano-shaped and thus different from that of a binary system. Finally, by evaluating the radiation reaction time scale, we give an order-of-magnitude estimate of merging event rates. The estimate suggests that figure-eight sources, which require carefully prepared initial states, may be too rare to detect.

We present a simplified solution to orbit determination of a binary system from astrometric observations. An exact solution was found by Asada, Akasaka and Kasai by assuming no observational errors. We extend the solution considering observational data. The generalized solution is expressed in terms of elementary functions, and therefore requires neither iterative nor numerical methods.

We present an exact solution of the equations for orbit determination of a two body system in a hyperbolic or parabolic motion. In solving this problem, we extend the method employed by Asada, Akasaka and Kasai (AAK) for a binary system in an elliptic orbit. The solutions applicable to each of elliptic, hyperbolic and parabolic orbits are obtained by the new approach, and they are all expressed in an explicit form, remarkably, only in terms of elementary functions. We show also that the solutions for an open orbit are recovered by making a suitable transformation of the AAK solution for an elliptic case.

The backreaction of nonlinear inhomogeneities to the cosmic expansion is re-analyzed in the framework of general relativity. Apparent discrepancies regarding the effect of the nonlinear backreaction, which exist among the results of previous works in different gauges, are resolved. By defining the spatially averaged matter energy density as a conserved quantity in the large comoving volume, it is shown that the nonlinear backreaction neither accelerates nor decelerates the cosmic expansion in a matter-dominated universe. The present result in the Newtonian gauge is consistent with the previous results obtained in the comoving synchronous gauge. Although our work does not give a complete proof, it strongly suggests the following no-go theorem: No cosmic acceleration occurs as a result of the nonlinear backreaction via averaging.

We study analytically a gravitational lens due to a deformed star, which is modeled by using a monopole and a quadrupole moment. Positions of the images are discussed for a source on the principal axis. We present explicit expressions for the lens equation for this gravitational lens as a single real tenth-order algebraic equation. Furthermore, we compute an expression for the caustics as a discriminant for the polynomial. Another simple parametric representation of the caustics is also presented in a more tractable form. A simple expression for the critical curves is obtained to clarify a topological feature of the critical curves; the curves are simply connected if and only if the distortion is sufficiently large.

It is believed that some numerical technique must be employed for the determination of the system parameters of a visual binary or a star with a planet because the relevant equations are not only highly nonlinear but also transcendental owing to the Kepler's equation. Such a common sense, however, is not true; we discover an analytic inversion formula, in which the original orbital parameters are expressed as elementary functions of the observable quantities such as the location of four observed points and the time interval between these points. The key thing is that we use the time interval but not the time of each observation in order to avoid treating the Kepler's equation. The present formula can be applied even in cases where the observations cover a short arc of the orbit during less than one period. Thus the formula will be useful in the future astrometric missions such as SIM, GAIA and JASMINE.

We study the effect of the core radius of an isothermal ellipsoidal gravitational lens on the caustics and the critical curves. We derive an analytic expression of the caustics for an isothermal ellipsoidal gravitational lens via a sixth-order algebraic equation. Since the expression is too long, by using another method we obtain a parametric representation of the critical curves in order to show analytically that there exist three cases: There are two curves for a small core radius, one for a quite large one, and no curves appear for an extremely large one, though the latter two cases are not realistic. The caustics are represented also by the same parameter.

The Euclidean algorithm in algebra is applied to a class of gravitational lenses for which the lens equation consists of any set of coupled polynomial equations in the image position. In general, this algorithm allows us to reduce an apparently coupled system to a single polynomial in one variable (say $x$ in Cartesian coordinates) without the other component (say $y$), which is expressed as a function of the first component. This reduction enables us to investigate the lensing properties in an algebraic manner: For instance, we can obtain an analytic expression of the caustics by computing the discriminant of the polynomial equation. To illustrate this Euclidean algorithm, we re-examine a binary gravitational lens and show that the lens equation is reduced to a single real fifth-order equation, in agreement with previous works. We apply this algorithm also to the linearlized Kerr lens and find that the lens equation is reduced to a single real fifth-order one.

Einstein gravity with extra dimensions or alternative gravity theories might suggest that the gravity propagation speed can be different from the light speed. Such a difference may play a vital role in the primordial universe. In recent, Kopeikin and Fomalont claimed the first measurement of the gravity speed by VLBI. However, the measurement has no relevance with the speed of gravity as I had shown before the observation was done. It seems that our conclusion has been established well by re-examining recent papers with great care.

We find a simple expression for critical curves of a binary gravitational lens. On the basis of this, we present a parametric representation of such curves. The caustics can also be expressed with the same parameterization. The present result is helpful for efficiently constructing many templates of light curves due to binary systems, particularly extrasolar planets, which cause spikes in the light curves when a source crosses the caustics.

We derive the deflection angle up to $O(m^2a)$ due to a Kerr gravitational lens with mass $m$ and specific angular momentum $a$. It is known that at the linear order in $m$ and $a$ the Kerr lens is observationally equivalent to the Schwarzschild one because of the invariance under the global translation of the center of the lens mass. We show, however, nonlinear couplings break the degeneracy so that the rotational effect becomes in principle separable for multiple images of a single source. Furthermore, it is distinguishable also for each image of an extended source and/or a point source in orbital motion. In practice, the correction at $O(m^2a)$ becomes $O(10^{-10})$ for the supermassive black hole in our galactic center. Hence, these nonlinear gravitational lensing effects are too small to detect by near-future observations.

It has been recently shown that the lens equation for a binary gravitational lens, which is apparently a coupled system, can be reduced to a real fifth-order (quintic) algebraic equation. Some algebraic properties of the real quintic equation are revealed. We find that the number of images on each side of the separation axis is independent of the mass ratio and separation unless the source crosses the caustics. Furthermore, the discriminant of the quintic equation enables us to study changes in the number of solutions, namely in the number of images. It is shown that this discriminant can be factorized into two parts: One represents the condition that the lens equation can be reduced to a single quintic equation, while the other corresponds to the caustics.

We present explicit expressions for the lens equation for a cored isothermal ellipsoidal gravitational lens as a single real sixth-order algebraic equation in two approaches; 2-dimensional Cartesian coordinates and 3-dimensional polar ones. We find a condition for physical solutions which correspond to at most five images. For a singular isothermal ellipsoid, the sixth-order equation is reduced to fourth-order one for which analytic solutions are well-known. Furthermore, we derive analytic criteria for determining the number of images for the singular lens, which give us simple expressions for the caustics and critical curves. The present formulation offers a useful way for studying galaxy lenses frequently modeled as isothermal ellipsoids.

It is shown that the lens equation for a binary gravitational lens being a set of two coupled real fifth-order algebraic equations (equivalent to a single complex equation of the same order) can be reduced to a single real fifth-order algebraic equation, which provides a much simpler way to study lensing by binary objects.

We investigate the light-cone effect on the Shapiro time delay. The extra time delay caused by Jupiter on the 8th of September 2002 can be measured by advanced VLBI (very long baseline interferometry). Our expression for the delay is in complete agreement with that of Kopeikin (2001), in which he argued that the excess time delay was due to the propagation of gravity. The present letter, however, shows that the excess comes from nothing but the propagation of light, namely the light-cone effect. To make a robust confirmation of general relativity by the coming Jupiter event, it is important to take account of the light-cone effect on the Shapiro time delay.

We have developed a perturbative approach to microlensing due to an extrasolar planetary lens. In particular, we have found analytic formulae for triple images. We have used the formulae to investigate the astrometric microlensing due to the extrasolar planetary lens, in expectation of dramatic improvements in the precision of the future astrometric measurements. For a weak lensing case, we have shown how the maximum angular size and the typical time scale of the anomalous shift of the lightcentroid are dependent on the mass ratio and angular separation between the star and the planet.

Studies of gravitational microlensing are usually based on the lens equation, which is valid only to first order in the gravitational constant $G$. However, the amplification factor of microlensing is a second-order quantity with respect to $G$. Despite this fact, conventional studies are still based on the lowest-order lens equation. Then the question naturally arises: Why are these conventional studies justified? We carefully study the relativistic correction to the amplification factor at $O(G^2)$. We show that the amplification factor for each image is corrected. However, the total amplification remains unchanged at this order.

The Lagrangian description of fluid flow in relativistic cosmology is extended to the case of flow accelerated by pressure. In the description, the entropy and the vorticity are obtained exactly for the barotropic equation of state. In order to determine the metric, the Einstein equation is solved perturbatively, when metric fluctuations are small but entropy inhomogeneities are large. Thus, the present formalism is applicable to the case when the inhomogeneities are small in the large scale but locally nonlinear.

It is known that the rotation of a gravitational lens affects properties of images. We consider an inverse problem: If the lens is dark, can we infer its rotation from the observed images? We find that, up to the first order in the gravitational constant G, a rotating lens is not distinguishable from a non-rotating one. The quadrupole contribution to the deflection angle is also calculated for a general, extended lens object.

We present a new set of 3.5 Post-Newtonian equations in which Newtonian hydrodynamics is coupled to the nonconservative effects of gravitational radiation emission. Our formalism differs in two significant ways from a similar 3.5 Post-Newtonian approach proposed by Blanchet (1993, 1997). Firstly we concentrate only on the radiation-reaction effects produced by a time-varying mass-current quadrupole $S_{ij}$. Secondly, we adopt a gauge in which the radiation-reaction force densities depend on the fourth time derivative of $S_{ij}$, rather than on the fifth, as in Blanchet's approach. This difference makes our formalism particularly well-suited to numerical implementation and could prove useful in performing fully numerical simulations of the recently discovered $r$-mode instability for rotating neutron stars subject to axial perturbations.

Distances play important roles in cosmological observations, especially in gravitational lens systems, but there is a problem in determining distances because they are defined in terms of light propagation, which is influenced gravitationally by the inhomogeneities in the universe. In this paper we first give the basic optical relations and the definitions of different distances in inhomogeneous universes. Next we show how the observational relations depend quantitatively on the distances. Finally, we give results for the frequency distribution of different distances and the shear effect on distances obtained using various methods of numerical simulation.

We develop the Lagrangian perturbation theory in the general relativistic cosmology, which enables us to take into account the vortical effect of the dust matter. Under the Lagrangian representation of the fluid flow, the propagation equation for the vorticity as well as the density is exactly solved. Based on this, the coupling between the density and vorticity is clarified in a non-perturbative way. The relativistic correspondence to the Lagrangian perturbation theory in the Newtonian cosmology is also emphasized.

We discuss how inhomogeneities of the universe affect observations of the gravitational lensing; (1) the bending angle, (2) the lensing statistics and (3) the time delay. In order to take account of the inhomogeneities, the Dyer-Roeder distance is used, which includes a parameter representing the clumpiness of the matter along the line of sight. It is shown analytically that all three combinations of distances appearing in the above observations (1)-(3) are monotonic with respect to the clumpiness in general, for any given set of the density parameter, cosmological constant and redshifts of the lens and the source. Some implications of this result for the observation are presented; the clumpiness decreases both the bending angle and the lensing event rate, while it increases the time delay. We also discuss cosmological tests using the gravitational lensing in the clumpy universe.

We investigate the cosmological test recently proposed by B. Fort, Y. Mellier and M. Dantel-Fort (FMD), where the observed location of the critical line in gravitational lensing is used to determine the cosmological parameters, $\Omega$ and $\lambda$. Applying this method to the cluster of galaxies Cl0024+1654, FMD obtained a constraint on the cosmological constant, $\lambda > 0.6$, assuming the spatially flat universe. It plays a crucial role in this method that the angular diameter distance-redshift relation depends on the cosmological models through the cosmological parameters. First, using the angular diameter distance in the Friedmann-Lemaitre- Robertson-Walker universe, we show that one can hardly determine $\Omega$ by this method without the assumption of the spatially flat universe. We also investigate the effect of inhomogeneities of the universe by using the Dyer-Roeder angular diameter distance. It is shown that the effect of inhomogeneities can become too large to be ignored, particularly for a high density universe. As a result, this method cannot be taken as a clear cosmological test to determine $\Omega$ and $\lambda$, though it may provide a bound on $\Omega$ and $\lambda$. Moreover, we mention the uncertainty of the determination of the velocity dispersion, which is regarded as one of the most serious problems in this test.

We present a formalism to obtain equilibrium configurations of uniformly rotating fluid in the second post-Newtonian approximation of general relativity. In our formalism, we need to solve 29 Poisson equations, but their source terms decrease rapidly enough at the external region of the matter(i.e., at worst $O(r^{-4})$). Hence these Poisson equations can be solved accurately as the boundary value problem using standard numerical methods.This formalism will be useful to obtain nonaxisymmetric uniformly rotating equilibrium configurations such as synchronized binary neutron stars just before merging and the Jacobi ellipsoid.

Using the (3+1) formalism in general relativity, we perform the post-Newtonian(PN) approximation to clarify what sort of gauge condition is suitable for numerical analysis of coalescing compact binary neutron stars and gravitational waves from them. We adopt a kind of transverse gauge condition to determine the shift vector. On the other hand, for determination of the time slice, we adopt three slice conditions(conformal slice, maximal slice and harmonic slice) and discuss their properties. Using these conditions, the PN hydrodynamic equations are obtained up through the 2.5PN order including the quadrupole gravitational radiation reaction. In particular, we describe methods to solve the 2PN tensor potential which arises from the spatial 3-metric. It is found that the conformal slice seems appropriate for analysis of gravitational waves in the wave zone and the maximal slice will be useful for describing the equilibrium configurations. The PN approximation in the (3+1) formalism will be also useful to perform numerical simulations using various slice conditions and, as a result, to provide an initial data for the final merging phase of coalescing binary neutron stars which can be treated only by fully general relativistic simulations.