We investigate gravitational waves with sub-nanoHz frequencies ($10^{-11}$ Hz $\lesssim f_{\rm GW} \lesssim 10^{-9}$ Hz) from the spatial distribution of the spin-down rates of milli-second pulsars. As we suggested in Yonemaru et al. 2018, gravitational waves from a single source induces the bias in the observed spin-down rates of pulsars depending on the relative direction between the source and pulsar. To improve the constraints on the time derivative of gravitational-wave amplitude obtained in our previous work (Kumamoto et al. 2019), we adopt a more sophisticated statistical method called the Mann-Whitney U test. Applying our method to the ATNF pulsar catalogue, we first found that the current data set is consistent with no GW signal from any direction in the sky. Then, we estimate the effective angular resolution of our method to be $(66~{\rm deg})^2$ by studying the probability distribution of the test statistic. Finally, we investigate gravitational-wave signal from the Galactic Centre and M87 and, comparing simulated mock data sets with the real pulsar data, we obtain the upper bounds on the time derivative as $\dot{h}_{\rm GC} < 8.9 \times 10^{-19} {\rm s}^{-1}$ for the Galactic Centre and $\dot{h}_{\rm M87} < 3.3 \times 10^{-19} {\rm s}^{-1}$ for M87, which are stronger than the ones obtained in Kumamoto et al. 2019 by factors of 7 and 25, respectively.

Pulsar search with time-domain observation is very computationally expensive and data volume will be enormous with the next generation telescopes such as the Square Kilometre Array. We apply artificial neural networks (ANNs), a machine learning method, for efficient selection of pulsar candidates from radio continuum surveys, which are much cheaper than time-domain observation. With observed quantities such as radio fluxes, sky position and compactness as inputs, our ANNs output the "score" that indicates the degree of likeliness of an object to be a pulsar. We demonstrate ANNs based on existing survey data by the TIFR GMRT Sky Survey (TGSS) and the NRAO VLA Sky Survey (NVSS) and test their performance. Precision, which is the ratio of the number of pulsars classified correctly as pulsars to that of any objects classified as pulsars, is about 96$\%$. Finally, we apply the trained ANNs to unidentified radio sources and our fiducial ANN with five inputs (the galactic longitude and latitude, the TGSS and NVSS fluxes and compactness) generates 2,436 pulsar candidates from 456,866 unidentified radio sources. These candidates need to be confirmed if they are truly pulsars by time-domain observations. More information such as polarization will narrow the candidates down further.

The globular cluster Omega Centauri is the most massive and luminous cluster in the Galaxy. The $\gamma$-ray source FL8Y J1326.7$-$4729 is coincident with the core of the cluster, leading to speculation that hitherto unknown radio pulsars or annihilating dark matter may be present in the cluster core. Here we report on the discovery of five millisecond pulsars in Omega Centauri following observations with the Parkes radio telescope. Four of these pulsars are isolated with spin periods of 4.1, 4.2, 4.6 and 6.8 ms. The fifth has a spin period of 4.8 ms and is in an eclipsing binary system with an orbital period of 2.1 hours. Deep radio continuum images of the cluster centre with the Australian Telescope Compact Array reveal a small population of compact radio sources making it likely that other pulsars await discovery. We consider it highly likely that the millisecond pulsars are the source of the $\gamma$-ray emission. The long-term timing of these pulsars opens up opportunities to explore the dynamics and interstellar medium of the cluster.

We probe ultra-low-frequency gravitational waves (GWs) with statistics of spin-down rates of milli-second pulsars (MSPs) by a method proposed in our prevous work (Yonemaru et al. 2016). The considered frequency range is $10^{-12}{\rm Hz} \lesssim f_{\rm GW} \lesssim 10^{-10}$Hz, which cannot be accessed by the conventional pulsar timing array. The effect of such low-frequency GWs appears as a bias to spin-down rates which has a quadrupole pattern in the sky. We use the skewness of the spin-down rate distribution and the number of MSPs with negative spin-down rates to search for the bias induced by GWs. Applying this method to 149 MSPs selected from the ATNF pulsar catalog, we derive upper bounds on the time derivative of the GW amplitudes of $\dot{h} < 6.2 \times 10^{-18}~{\rm sec}^{-1}$ and $\dot{h} < 8.1 \times 10^{-18}~{\rm sec}^{-1}$ in the directions of the Galactic Center and M87, respectively. Approximating the GW amplitude as $\dot{h} \sim 2 \pi f_{\rm GW} h$, the bounds translate into $h < 3 \times 10^{-9}$ and $h < 4 \times 10^{-9}$, respectively, for $f_{\rm GW} = 1/(100~{\rm yr})$. Finally, we give the implications to possible super-massive black hole binaries at these sites.

A new detection method for gravitational waves (GWs) with ultra-low frequencies ($f_{\rm GW} \lesssim 10^{-10}~{\rm Hz}$), which is much lower than the range of pulsar timing arrays (PTAs), was proposed in Yonemaru et al. (2016). This method utilizes the statistical properties of spin-down rates of milli-second pulsars (MSPs) and the sensitivity was evaluated in Yonemaru et al. (2018). There, some simplifying assumptions, such as neglect of the "pulsar term" and spatially uniform distribution of MSPs, were adopted and the sensitivity on the time derivative of GW amplitude was estimated to be $10^{-19}~{\rm s}^{-1}$ independent of the direction, polarization and frequency of GWs. In this paper, extending the previous analysis, realistic simulations are performed to evaluate the sensitivity more reasonably. We adopt a model of 3-dimensional pulsar distribution in our Galaxy and take the pulsar term into account. As a result, we obtain expected sensitivity as a function of the direction, polarization and frequency of GWs. The dependence on GW frequency is particularly significant and the sensitivity becomes worse by a few orders for $< 10^{-12}~{\rm Hz}$ compared to the previous estimates.

A new detection method for ultra-low frequency gravitational waves (GWs) with a frequency much lower than the observational range of pulsar timing arrays (PTAs) was suggested in Yonemaru et al. (2016). In the PTA analysis, ultra-low frequency GWs ($\lesssim 10^{-10}$ Hz) which evolve just linearly during the observation time span are absorbed by the pulsar spin-down rates since both have the same effect on the pulse arrival time. Therefore, such GWs cannot be detected by the conventional method of PTAs. However, the bias on the observed spin-down rates depends on relative direction of a pulsar and GW source and shows a quadrupole pattern in the sky. Thus, if we divide the pulsars according to the position in the sky and see the difference in the statistics of the spin-down rates, ultra-low frequency GWs from a single source can be detected. In this paper, we evaluate the potential of this method by Monte-Carlo simulations and estimate the sensitivity, considering only the "Earth term" while the "pulsar term" acts like random noise for GW frequencies $10^{-13}-10^{-10}$ Hz. We find that with 3,000 milli-second pulsars, which are expected to be discovered by a future survey with the Square Kilometre Array, GWs with the derivative of amplitude of about $3 \times 10^{-19}~\rm{s}^{-1}$ can in principle be detected. Implications for possible supermassive binary black holes in Sgr${}^*$ and M87 are also given.

In this paper, we study gravitational-wave (GW) emission from a hypothetical supermassive black-hole (SMBH) binary at the center of M87. The existence of a SMBH other than that usually identified with the central AGN is a possible explanation for the observed displacement ($\sim O(1)~{\rm pc}$) between the AGN and the galactic centroid, and it is reasonable to assume consid- ering the evolution of SMBHs through galaxy mergers. Because the period of the binary and the resulting GWs is much longer than the observational time span, we calculate the variation of the GW amplitude, rather than the amplitude itself. We investigate the dependence on the orbital elements and the second BH mass taking the observational constraints into account. The frequency of the GWs is too low to be detected with the conventional pulsar timing array and we propose a new method to detect such low-frequency GWs with the distribution func- tion of pulsar spin-down rates. Although the GWs from a SMBH binary which explains the observed displacement is extremely hard to be detected even with the new method, GWs are still a useful way to probe the M87 center.

Pulsar timing arrays are one of the powerful tools to test the existence of cosmic strings through searching for the gravitational wave background. The amplitude of the background connects to information on cosmic strings such as the tension and string network properties. In addition, one may be able to extract more information on properties of cosmic strings by measuring anisotropies in the gravitational wave (GW) background. In this paper, we provide estimates of the level of anisotropy expected in the GW background generated by cusps on cosmic strings. We find that the anisotropy level strongly depends on the initial loop size $\alpha$, and thus we may be able to put constraint on $\alpha$ by measuring the anisotropy of the GW background. We also find that certain regions of the parameter space can be probed by shifting the observation frequency of GWs.

The Square Kilometre Array will revolutionize pulsar studies with its wide field-of-view, wide-band observation and high sensitivity, increasing the number of observable pulsars by more than an order of magnitude. Pulsars are of interest not only for the study of neutron stars themselves but for their usage as tools for probing fundamental physics such as general relativity, gravitational waves and nuclear interaction. In this article, we summarize the activity and interests of SKA-Japan Pulsar Science Working Group, focusing on an investigation of modified gravity theory with the supermassive black hole in the Galactic Centre, gravitational-wave detection from cosmic strings and binary supermassive black holes, a study of the physical state of plasma close to pulsars using giant radio pulses and determination of magnetic field structure of Galaxy with pulsar pairs.