XMASS-I is a single-phase liquid xenon detector whose purpose is direct detection of dark matter. To achieve the low background requirements necessary in the detector, a new model of photomultiplier tubes (PMTs), R10789, with a hexagonal window was developed based on the R8778 PMT used in the XMASS prototype detector. We screened the numerous component materials for their radioactivity. During development, the largest contributions to the reduction of radioactivity came from the stem and the dynode support. The glass stem was exchanged to the Kovar alloy one and the ceramic support were changed to the quartz one. R10789 is the first model of Hamamatsu Photonics K. K. that adopted these materials for low background purposes and provided a groundbreaking step for further reductions of radioactivity in PMTs. Measurements with germanium detectors showed 1.2$\pm$0.3 mBq/PMT of $^{226}$Ra, less than 0.78 mBq/PMT of $^{228}$Ra, 9.1$\pm$2.2 mBq/PMT of $^{40}$K, and 2.8$\pm$0.2 mBq/PMT of $^{60}$Co. In this paper, the radioactive details of the developed R10789 are described together with our screening methods and the components of the PMT.

A search for dark matter using an underground single-phase liquid xenon detector was conducted at the Kamioka Observatory in Japan, particularly for Weakly Interacting Massive Particles (WIMPs). We have used 705.9 live days of data in a fiducial volume containing 97 kg of liquid xenon at the center of the detector. The event rate in the fiducial volume after the data reduction was ${\rm (4.2 \pm 0.2) \times 10^{-3} \, day^{-1}kg^{-1} keV_{ee}^{-1}}$ at ${\rm 5 \, keV_{ee}}$, with a signal efficiency of ${\rm 20\%}$. All the remaining events are consistent with our background evaluation, mostly of the "mis-reconstructed events" originated from $^{210}$Pb in the copper plates lining the detector's inner surface. The obtained upper limit on a spin-independent WIMP-nucleon cross section was ${\rm 2.2 \times 10^{-44} \, cm^{2}}$ for a WIMP mass of ${\rm 60 \, GeV/c^{2}}$ at the $90\%$ confidence level, which was the most stringent limit among results from single-phase liquid xenon detectors.

An annual modulation signal due to the Earth orbiting around the Sun would be one of the strongest indications of the direct detection of dark matter. In 2016, we reported a search for dark matter by looking for this annual modulation with our single-phase liquid xenon XMASS-I detector. That analysis resulted in a slightly negative modulation amplitude at low energy. In this work, we included more than one year of additional data, which more than doubles the exposure to 800 live days with the same 832 kg target mass. When we assume weakly interacting massive particle (WIMP) dark matter elastically scattering on the xenon target, the exclusion upper limit for the WIMP-nucleon cross section was improved by a factor of 2 to 1.9$\times$10$^{-41}$cm$^2$ at 8 GeV/c$^2$ at 90\% confidence level with our newly implemented data selection through a likelihood method. For the model-independent case, without assuming any specific dark matter model, we obtained more consistency with the null hypothesis than before with a $p$-value of 0.11 in the 1$-$20 keV energy region. This search probed this region with an exposure that was larger than that of DAMA/LIBRA. We also did not find any significant amplitude in the data for periodicity with periods between 50 and 600 days in the energy region between 1 to 6 keV.

In theories with the large extra dimensions beyond the standard 4-dimensional spacetime, axions could propagate in such extra dimensions, and acquire Kaluza-Klein (KK) excitations. These KK axions are produced in the Sun and could solve unexplained heating of the solar corona. While most of the solar KK axions escape from the solar system, a small fraction is gravitationally trapped in orbits around the Sun. They would decay into two photons inside a terrestrial detector. The event rate is expected to modulate annually depending on the distance from the Sun. We have searched for the annual modulation signature using $832\times 359$ kg$\cdot$days of XMASS-I data. No significant event rate modulation is found, and hence we set the first experimental constraint on the KK axion-photon coupling of $4.8 \times 10^{-12}\, \mathrm{GeV}^{-1}$ at 90% confidence level for a KK axion number density of $\bar{n}_\mathrm{a} = 4.07 \times 10^{13}\, \mathrm{m}^{-3}$, the total number of extra dimensions $n = 2$, and the number of extra dimensions $\delta = 2$ that axions can propagate in.

We report the measurement of the emission time profile of scintillation from gamma-ray induced events in the XMASS-I 832 kg liquid xenon scintillation detector. Decay time constant was derived from a comparison of scintillation photon timing distributions between the observed data and simulated samples in order to take into account optical processes such as absorption and scattering in liquid xenon. Calibration data of radioactive sources, $^{55}$Fe, $^{241}$Am, and $^{57}$Co were used to obtain the decay time constant. Assuming two decay components, $\tau_1$ and $\tau_2$, the decay time constant $\tau_2$ increased from 27.9 ns to 37.0 ns as the gamma-ray energy increased from 5.9 keV to 122 keV. The accuracy of the measurement was better than 1.5 ns at all energy levels. A fast decay component with $\tau_1 \sim 2$ ns was necessary to reproduce data. Energy dependencies of $\tau_2$ and the fraction of the fast decay component were studied as a function of the kinetic energy of electrons induced by gamma-rays. The obtained data almost reproduced previously reported results and extended them to the lower energy region relevant to direct dark matter searches.

A search for dark matter was conducted by looking for an annual modulation signal due to the Earth's rotation around the Sun using XMASS, a single phase liquid xenon detector. The data used for this analysis was 359.2 live days times 832 kg of exposure accumulated between November 2013 and March 2015. When we assume Weakly Interacting Massive Particle (WIMP) dark matter elastically scattering on the target nuclei, the exclusion upper limit of the WIMP-nucleon cross section 4.3$\times$10$^{-41}$cm$^2$ at 8 GeV/c$^2$ was obtained and we exclude almost all the DAMA/LIBRA allowed region in the 6 to 16 GeV/c$^2$ range at $\sim$10$^{-40}$cm$^2$. The result of a simple modulation analysis, without assuming any specific dark matter model but including electron/$\gamma$ events, showed a slight negative amplitude. The $p$-values obtained with two independent analyses are 0.014 and 0.068 for null hypothesis, respectively. we obtained 90\% C.L. upper bounds that can be used to test various models. This is the first extensive annual modulation search probing this region with an exposure comparable to DAMA/LIBRA.

One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime.

We calculate the secular changes of the orbital parameters of a point particle orbiting a Kerr black hole, due to the gravitational radiation reaction. For this purpose, we use the post-Newtonian (PN) approximation in the first order black hole perturbation theory, with the expansion with respect to the orbital eccentricity. In this work, the calculation is done up to the fourth post-Newtonian (4PN) order and to the sixth order of the eccentricity, including the effect of the absorption of gravitational waves by the black hole. We confirm that, in the Kerr case, the effect of the absorption appears at the 2.5PN order beyond the leading order in the secular change of the particle's energy and may induce a superradiance, as known previously for circular orbits. In addition, we find that the superradiance may be suppressed when the orbital plane inclines with respect to the equatorial plane of the central black hole. We also investigate the accuracy of the 4PN formulae by comparing to numerical results. If we require that the relative errors in the 4PN formulae are less than $10^{-5}$, the parameter region to satisfy the condition will be $p\gtrsim 50$ for $e=0.1$, $p\gtrsim 80$ for $e=0.4$, and $p\gtrsim 120$ for $e=0.7$ almost irrespective of the inclination angle nor the spin of the black hole, where $p$ and $e$ are the semi-latus rectum and the eccentricity of the orbit. The region can further be extended using an exponential resummation method to $p\gtrsim 40$ for $e=0.1$, $p\gtrsim 60$ for $e=0.4$, and $p\gtrsim 100$ for $e=0.7$. Although we still need the higher order calculations of the PN approximation and the expansion with respect to the orbital eccentricity to apply for data analysis of gravitational waves, the results in this paper would be an important improvement from the previous work at the 2.5PN order, especially for large $p$ region.

We discuss the inspiral of a small body around a Kerr black hole. When the time scale of the radiation reaction is sufficiently longer than its orbital period, the leading order orbital evolution is described only by the knowledge of the averaged evolution of the constants of motion, i.e., the energy, azimuthal angular momentum and the Carter constant. Although there is no conserved current composed of the perturbation field corresponding to the Carter constant, it has been shown that the averaged rate of change of the Carter constant can be given by a simple formula, when there exists a simultaneous turning point of the radial and polar oscillations. However, an inspiralling orbit may cross a "resonance" point, where the frequencies of the radial and polar orbital oscillations are in a rational ratio. At the resonance point, one cannot find a simultaneous turning point in general. Hence, even for the averaged rate of change of the Carter constant, a direct computation of the 'self-force", which is still challenging especially in the case of the Kerr background, seems to be necessary. In this paper, we present a method of computing the averaged rate of change of the Carter constant in a relatively simple manner at the "resonance" point.

We derive the analytical solutions of the bound timelike geodesic orbits in Kerr spacetime. The analytical solutions are expressed in terms of the elliptic integrals using Mino time $\lambda$ as the independent variable. Mino time decouples the radial and polar motion of a particle and hence leads to forms more useful to estimate three fundamental frequencies, radial, polar and azimuthal motion, for the bound timelike geodesics in Kerr spacetime. This paper gives the first derivation of the analytical expressions of the fundamental frequencies. This paper also gives the first derivation of the analytical expressions of all coordinates for the bound timelike geodesics using Mino time. These analytical expressions should be useful not only to investigate physical properties of Kerr geodesics but more importantly to applications related to the estimation of gravitational waves from the extreme mass ratio inspirals.