Parity-violating (PV) gravity has recently attracted interest in several aspects. One of them is the axion-graviton coupling to test axion-dark matter model. Moreover, by extending Chern-Simons (CS) gravity to include derivatives of a scalar field up to second order, a more general class of PV gravity theory, which we call the CNCL model, has been proposed~[M. Crisostomi et al., Phys. Rev. D, 97, 044034 (2018)]. The model can be further extended by including even higher derivatives of the scalar field and/or higher curvature terms. In this paper, we discuss the effect of parity violation in gravitational sector on the propagation of gravitational waves from binary coalescence by introducing a model-independent parametrization of modification. Our parametrization includes the CNCL model as well as CS gravity. The effect of parity violation on the gravitational waveform is maximum when the source binary orientation to our line of sight is edge-on, while the modified waveform reduces to the parity-symmetric one when the source is face-on. We perform a search for the signature of such modification by using LIGO/Virgo O1/O2 catalog. We find that the catalog data is consistent with general relativity and obtain constraints on parity violation in gravity for various post-Newtonian order modifications for the first time. The obtained constraint on CS gravity is consistent with the results in previous works. On the other hand, the constraint on the CNCL model that we obtain is tighter than the previous results by roughly 7 orders of magnitude.

The direct detection of gravitational waves now provides a new channel of testing gravity theories. Despite that the parametrized post-Einsteinian framework is a powerful tool to quantitatively investigate effects of modification of gravity theory, the gravitational waveform in this framework is still extendable. One of such extensions is to take into account the gradual activation of dipole radiation due to massive fields, which are still only very weakly constrained if their mass $m$ is greater than $10^{-16}$ eV from pulsar observations. Ground-based gravitational-wave detectors, LIGO, Virgo, and KAGRA, are sensitive to this activation in the mass range, $10^{-14}$ eV $\lesssim m \lesssim 10^{-13}$ eV. Hence, we discuss a dedicated test for dipole radiation due to a massive field using the LIGO-Virgo collaboration's open data. In addition, assuming Einstein-dilaton-Gauss-Bonnet (EdGB) type coupling, we combine the results of the analysis of the binary black hole events to obtain the 90\% confidence level constraints on the coupling parameter $\alpha_{\rm EdGB}$ as $\sqrt{\alpha_{\rm EdGB}} \lesssim 2.47$ km for any mass less than $6 \times 10^{-14}$ eV for the first time, including $\sqrt{\alpha_{\rm EdGB}} \lesssim 1.85$ km in the massless limit.

We promote the precise gamma-ray observation project Gamma-Ray Astro-Imager with Nuclear Emulsion (GRAINE), which uses balloon-borne emulsion gamma-ray telescopes. The emulsion telescope realizes observations with high angular resolution, polarization sensitivity, and large aperture area in the 0.01--100 GeV energy region. Herein, we report the data analysis of emulsion tracks and the first demonstration of gamma-ray imaging via an emulsion telescope by using the flight data from the balloon experiment performed in 2015 (GRAINE 2015). The emulsion films were scanned by the latest read-out system for a total area of 41 m$^2$ in three months, and then the gamma-ray event selection was automatically processed. Millions of electron-pair events are accumulated in the balloon-borne emulsion telescope. The emulsion telescope detected signals from a calibration source (gamma rays from the interaction of cosmic rays with an aluminum plate) with a high significance during the balloon observation and created a gamma-ray image consistent with the source size and the expected angular resolution in the energy range of 100--300 MeV. The flight performance obtained in the GRAINE 2015 experiment proves that balloon-borne emulsion telescope experiments with larger area are feasible while maintaining expected imaging performance.

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.

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.

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.

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.

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.

Spherical solar dynamo simulations are performed. Self-consistent, fully compressible magnetohydrodynamic system with a stably stratified layer below the convective envelope is numerically solved with a newly developed simulation code based on the Yin-Yang grid. The effects of penetrative convection are studied by comparing two models with and without the stable layer. The differential rotation profile in both models is reasonably solar-like with equatorial acceleration. When considering the penetrative convection, a tachocline-like shear layer is developed and maintained beneath the convection zone without assuming any forcing. While turbulent magnetic field becomes predominant in the region where the convective motion is vigorous, mean-field component is preferentially organized in the region where the convective motion is less vigorous. Especially in the stable layer, the strong large-scale field with a dipole symmetry is spontaneously built up. The polarity reversal of the mean-field component takes place globally and synchronously throughout the system regardless the presence of the stable layer. Our results suggest that the stably stratified layer is a key component for organizing the large-scale strong magnetic field, but is not essential for the polarity reversal.

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.

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 present imaging simulations of the Sunyaev-Zel'dovich effect of galaxy clusters for the Atacama Large Millimeter/submillimeter Array (ALMA) including the Atacama Compact Array (ACA). In its most compact configuration at 90GHz, ALMA will resolve the intracluster medium with an effective angular resolution of 5 arcsec. It will provide a unique probe of shock fronts and relativistic electrons produced during cluster mergers at high redshifts, that are hard to spatially resolve by current and near-future X-ray detectors. Quality of image reconstruction is poor with the 12m array alone but improved significantly by adding ACA; expected sensitivity of the 12m array based on the thermal noise is not valid for the Sunyaev-Zel'dovich effect mapping unless accompanied by an ACA observation of at least equal duration. The observations above 100 GHz will become excessively time-consuming owing to the narrower beam size and the higher system temperature. On the other hand, significant improvement of the observing efficiency is expected once Band 1 is implemented in the future.

We study the torque acting on a planet embedded in an optically thick accretion disc, using global two-dimensional hydrodynamic simulations. The temperature of an optically thick accretion disc is determined by the energy balance between the viscous heating and the radiative cooling. The radiative cooling rate depends on the opacity of the disc. The opacity is expressed as a function of the temperature. We find the disc is divided into three regions that have different temperature distributions. The slope of the entropy distribution becomes steep in the inner region of the disc with the high temperature and the outer region of the disc with the low temperature, while it becomes shallow in the middle region with the intermediate temperature. Planets in the inner and outer regions move outward owing to the large positive corotation torque exerted on the planet by an adiabatic disc, on the other hand, a planet in the middle region moves inward toward the central star. Planets are expected to accumulate at the boundary between the inner and middle regions of the adiabatic disc. The positive corotation torque decreases with an increase in the viscosity of the disc. We find that the positive corotation torque acting on the planet in the inner region becomes too small to cancel the negative Lindblad torque when we include the large viscosity, which destroys the enhancement of the density in the horseshoe orbit of the planet. This leads to the inward migration of the planet in the inner region of the disc. A planet with 5 Earth masses in the inner region can move outward in a disc with the surface density of 100 g/cm^2 at 1 AU when the accretion rate of a disc is smaller than 2x10^{-8} solar mass/yr.

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.

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 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.

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.

We study Type I migration of a planet in a radiatively efficient disk using global two dimensional hydrodynamic simulations. The large positive corotation torque is exerted on a planet by an adiabatic disk at early times when the disk has the steep negative entropy gradient. The gas on the horseshoe orbit of the planet is compressed adiabatically during the change of the orbit from the slow orbit to the fast orbit, increasing its density and exerting the positive torque on the planet. The planet would migrate outward in the adiabatic disk before saturation sets in. We further study the effect of energy dissipation by radiation on Type I migration of the planet. The corotation torque decreases when the energy dissipates effectively because the density of the gas on the horseshoe orbit does not increase by the compression compared with the gas of the adiabatic disk. The total torque is mainly determined by the negative Lindblad torque and becomes negative. The planet migrates inward toward the central star in the radiatively efficient disk. The migration velocity is dependent on the radiative efficiency and greatly reduced if the radiative cooling works inefficiently.

We have undertaken a search for the infrared emission from the intracluster dust in the Coma cluster of galaxies by the Multiband Imaging Photometer for Spitzer. Our observations yield the deepest mid and far-infrared images of a galaxy cluster ever achieved. In each of the three bands, we have not detected a signature of the central excess component in contrast to the previous report on the detection by Infrared Space Observatory (ISO). We still find that the brightness ratio between 70 and 160 microns shows a marginal sign of the central excess, in qualitative agreement with the ISO result. Our analysis suggests that the excess ratio is more likely due to faint infrared sources lying on fluctuating cirrus foreground. Our observations yield the 2 sigma upper limits on the excess emission within 100 kpc of the cluster center as 5 x 10^-3 MJy/sr, 6 x 10^-2 MJy/sr, and 7 x 10^-2 MJy/sr, at 24, 70, and 160 microns, respectively. These values are in agreement with those found in other galaxy clusters and suggest that dust is deficient near the cluster center by more than 3 orders of magnitude compared to the interstellar medium.

For 117 clusters of galaxies, we explore the detectability of intracluster dust grains by current and future infrared facilities, taking into account both collisional heating and sputtering of grains by ambient plasma. If the dust grains are injected into the intergalactic space with the amount and size comparable to the Galactic values, the dust-to-gas ratio is typically $10^{-6}$ and the mean dust temperature is $\sim 30$ K near the cluster center. The predicted infrared intensities lie marginally above the detection thresholds for Spitzer Space Telescope, ASTRO-F, Herschel and SPICA missions. For some nearby clusters such as Perseus, A3571, A2319, A3112 and A2204, good detections of intracluster dust signal are expected in the 70 $\mu$m band. Given rather tight constraints on the dust temperature from observed electron density and temperature, the dust mass can be inferred directly from the infrared observations. Further constraints on the size distribution will be obtained once multi-band data are available. They will definitely provide a powerful probe of the dust injection processes and dust-gas interactions in the intergalactic space.