Evidence for a large-scale supergalactic cosmic ray multiplet (arrival directions correlated with energy) structure is reported for ultra-high energy cosmic ray (UHECR) energies above 10$^{19}$ eV using seven years of data from the Telescope Array (TA) surface detector and updated to 10 years. Previous energy-position correlation studies have made assumptions regarding magnetic field shapes and strength, and UHECR composition. Here the assumption tested is that, since the supergalactic plane is a fit to the average matter density of the local Large Scale Structure (LSS), UHECR sources and intervening extragalactic magnetic fields are correlated with this plane. This supergalactic deflection hypothesis is tested by the entire field-of-view (FOV) behavior of the strength of intermediate-scale energy-angle correlations. These multiplets are measured in spherical cap section bins (wedges) of the FOV to account for coherent and random magnetic fields. The structure found is consistent with supergalactic deflection, the previously published energy spectrum anisotropy results of TA (the hotspot and coldspot), and toy-model simulations of a supergalactic magnetic sheet. The seven year data post-trial significance of this supergalactic structure of multiplets appearing by chance, on an isotropic sky, is found by Monte Carlo simulation to be 4.2$\sigma$. The ten years of data post-trial significance is 4.1$\sigma$. Furthermore, the starburst galaxy M82 is shown to be a possible source of the TA Hotspot, and an estimate of the supergalactic magnetic field using UHECR measurements is presented.

Motivated by the detection of a significant dipole structure in the arrival directions of ultrahigh-energy cosmic rays above 8 EeV reported by the Pierre Auger Observatory (Auger), we search for a large-scale anisotropy using data collected with the surface detector array of the Telescope Array Experiment (TA). With 11 years of TA data, a dipole structure in a projection of the right ascension is fitted with an amplitude of 3.3+- 1.9% and a phase of 131 +- 33 degrees. The corresponding 99% confidence-level upper limit on the amplitude is 7.3%. At the current level of statistics, the fitted result is compatible with both an isotropic distribution and the dipole structure reported by Auger.

Ultra high energy cosmic rays provide the highest known energy source in the universe to measure proton cross sections. Though conditions for collecting such data are less controlled than an accelerator environment, current generation cosmic ray observatories have large enough exposures to collect significant statistics for a reliable measurement for energies above what can be attained in the lab. Cosmic ray measurements of cross section use atmospheric calorimetry to measure depth of air shower maximum ($X_{\mathrm{max}}$), which is related to the primary particle's energy and mass. The tail of the $X_{\mathrm{max}}$ distribution is assumed to be dominated by showers generated by protons, allowing measurement of the inelastic proton-air cross section. In this work the proton-air inelastic cross section measurement, $\sigma^{\mathrm{inel}}_{\mathrm{p-air}}$, using data observed by Telescope Array's Black Rock Mesa and Long Ridge fluorescence detectors and surface detector array in hybrid mode is presented. $\sigma^{\mathrm{inel}}_{\mathrm{p-air}}$ is observed to be $520.1 \pm 35.8$[Stat.] $^{+25.0}_{-40}$[Sys.]~mb at $\sqrt{s} = 73$ TeV. The total proton-proton cross section is subsequently inferred from Glauber formalism and is found to be $\sigma^{\mathrm{tot}}_{\mathrm{pp}} = 139.4 ^{+23.4}_{-21.3}$ [Stat.]$ ^{+15.0}_{-24.0}$[Sys.]~mb.

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 performed Population III (Pop III) binary evolution by using population synthesis simulations for seven different models. We found that Pop III binaries tend to be binary black holes (BBHs) with chirp mass $M_{\rm chirp} \sim 30~M_{\odot}$ and they can merge at present day due to long merger time. The merger rate densities of Pop III BBHs at $z=0$ ranges 3.34--21.2 $\rm /yr/Gpc^3$ which is consistent with the aLIGO/aVIRGO result of 9.7--101 $\rm /yr/Gpc^3$. These Pop III binaries might contribute to some part of the massive BBH gravitational wave (GW) sources detected by aLIGO/aVIRGO. We also calculated the redshift dependence of Pop III BBH mergers. We found that Pop III low spin BBHs tend to merge at low redshift, while Pop III high spin BBHs do at high redshift, which can be confirmed by future GW detectors such as ET, CE, and DECIGO. These detectors can also check the redshift dependence of BBH merger rate and spin distribution. Our results show that except for one model, the mean effective spin $\left\langle \chi_{\rm eff} \right\rangle$ at $z=0$ ranges $0.02$--$0.3$ while at $z=10$ it does $0.16$--$0.64$. Therefore, massive stellar-mass BBH detection by GWs will be a key for the stellar evolution study in the early universe.

We present an upper limit on the flux of ultra-high-energy down-going neutrinos for $E > 10^{18}\ \mbox{eV}$ derived with the nine years of data collected by the Telescope Array surface detector (05-11-2008 -- 05-10-2017). The method is based on the multivariate analysis technique, so-called Boosted Decision Trees (BDT). Proton-neutrino classifier is built upon 16 observables related to both the properties of the shower front and the lateral distribution function.

A negative ion micro time projection chamber (NI$\mu$TPC) was developed and its performance studied. An NI$\mu$TPC is a novel technology that enables the measurement of absolute $z$ coordinates for self-triggering TPCs. This technology provides full-fiducialization analysis, which is not possible with conventional gaseous TPCs, and is useful for directional dark matter searches in terms of background rejection and the improvement of the angular resolution. The developed NI$\mu$TPC prototype had a detection volume of 12.8 $\times$ 25.6 $\times$ 144 mm$^{3}$. The absolute $z$ coordinate was determined with a location accuracy of 16 mm using minority carrieres of SF$_{5}^{-}$. Simultaneously, there was a successful reconstruction of the three-dimensional (3D) tracks with a spatial resolution of 130 $\mu\rm{m}$. This is the first demonstration of 3D tracking with the detection of absolute $z$ coordinates, and it is an important step in improving the sensitivity of directional dark matter searches.

The surface detector (SD) of the Telescope Array (TA) experiment allows one to indirectly detect photons with energies of order $10^{18}$ eV and higher and to separate photons from the cosmic-ray background. In this paper we present the results of a blind search for point sources of ultra-high energy (UHE) photons in the Northern sky using the TA SD data. The photon-induced extensive air showers (EAS) are separated from the hadron-induced EAS background by means of a multivariate classifier based upon 16 parameters that characterize the air shower events. No significant evidence for the photon point sources is found. The upper limits are set on the flux of photons from each particular direction in the sky within the TA field of view, according to the experiment's angular resolution for photons. Average 95% C.L. upper limits for the point-source flux of photons with energies greater than $10^{18}$, $10^{18.5}$, $10^{19}$, $10^{19.5}$ and $10^{20}$ eV are $0.094$, $0.029$, $0.010$, $0.0073$ and $0.0058$ km$^{-2}$yr$^{-1}$ respectively. For the energies higher than $10^{18.5}$ eV, the photon point-source limits are set for the first time. Numerical results for each given direction in each energy range are provided as a supplement to this paper.

We report a search for astronomical neutrinos in the energy region from several GeV to TeV in the direction of the blazar TXS0506+056 using the Super-Kamiokande detector following the detection of a 100 TeV neutrino from the same location by the IceCube collaboration. Using Super-Kamiokande neutrino data across several data samples observed from April 1996 to February 2018 we have searched for both a total excess above known backgrounds across the entire period as well as localized excesses on smaller time scales in that interval. No significant excess nor significant variation in the observed event rate are found in the blazar direction. Upper limits are placed on the electron and muon neutrino fluxes at 90\% confidence level as $6.03 \times 10^{-7}$ and $4.52 \times 10^{-7}$ to $9.26 \times 10^{-10}$ [${\rm erg}/{\rm cm}^2/{\rm s}$], respectively.

Magnetospheric Multiscale (MMS) observations of plasma turbulence generated by a Kelvin-Helmholtz (KH) event at the Earth's magnetopause are compared with a high-resolution two-dimensional (2D) hybrid direct numerical simulation (DNS) of decaying plasma turbulence driven by large-scale balanced Alfv\'enic fluctuations. The simulation, set up with four observation-driven physical parameters (ion and electron betas, turbulence strength, and injection scale) exhibits a quantitative agreement on the spectral, intermittency, and cascade-rate properties with in situ observations, despite the different driving mechanisms. Such agreement demonstrates a certain universality of the turbulent cascade from magnetohydrodynamic (MHD) to sub-ion scales, whose properties are mainly determined by the selected parameters, also indicating that the KH instability-driven turbulence has a quasi-2D nature. The validity of the Taylor hypothesis in the sub-ion spatial range suggests that the fluctuations at sub-ion scales have predominantly low frequencies, consistent with a kinetic Alfv\'en wave-like nature or with quasi-static structures. Finally, the third-order structure function analysis indicates that the cascade rate of the turbulence generated by a KH event in the magnetopause is an order of magnitude larger than in the ambient magnetosheath.

We discuss the scalar mode of gravitational waves emerging in the context of $F(R)$ gravity by taking into account the chameleon mechanism. Assuming a toy model with a specific matter distribution to reproduce the environment of detection experiment by a ground-based gravitational wave observatory, we find that chameleon mechanism remarkably suppresses the scalar wave in the atmosphere of Earth, compared with the tensor modes of the gravitational waves. We also discuss the possibility to detect and constrain scalar waves by the current gravitational observatories and advocate a necessity of the future space-based observations.

Supernova detection is a major objective of the Super-Kamiokande (SK) experiment. In the next stage of SK (SK-Gd), gadolinium (Gd) sulfate will be added to the detector, which will improve the ability of the detector to identify neutrons. A core-collapse supernova will be preceded by an increasing flux of neutrinos and anti-neutrinos, from thermal and weak nuclear processes in the star, over a timescale of hours; some of which may be detected at SK-Gd. This could provide an early warning of an imminent core-collapse supernova, hours earlier than the detection of the neutrinos from core collapse. Electron anti-neutrino detection will rely on inverse beta decay events below the usual analysis energy threshold of SK, so Gd loading is vital to reduce backgrounds while maximising detection efficiency. Assuming normal neutrino mass ordering, more than 200 events could be detected in the final 12 hours before core collapse for a 15-25 solar mass star at around 200 pc, which is representative of the nearest red supergiant to Earth, $\mathrm{\alpha}$Ori (Betelgeuse). At a statistical false alarm rate of 1 per century, detection could be up to 10 hours before core collapse, and a pre-supernova star could be detected by SK-Gd up to 600 pc away. A pre-supernova alert could be provided to the astrophysics community following gadolinium loading.

In this paper, we report the implementation and measured performance of our extreme-scale global simulation code on Sunway TaihuLight and two PEZY-SC2 systems: Shoubu System B and Gyoukou. The numerical algorithm is the parallel Barnes-Hut tree algorithm, which has been used in many large-scale astrophysical particle-based simulations. Our implementation is based on our FDPS framework. However, the extremely large numbers of cores of the systems used (10M on TaihuLight and 16M on Gyoukou) and their relatively poor memory and network bandwidth pose new challenges. We describe the new algorithms introduced to achieve high efficiency on machines with low memory bandwidth. The measured performance is 47.9, 10.6 PF, and 1.01PF on TaihuLight, Gyoukou and Shoubu System B (efficiency 40\%, 23.5\% and 35.5\%). The current code is developed for the simulation of planetary rings, but most of the new algorithms are useful for other simulations, and are now available in the FDPS framework.

In keeping with the Luxembourg government's initiative to support the future use of space resources, ASIME 2018 was held in Belval, Luxembourg on April 16-17, 2018. The goal of ASIME 2018: Asteroid Intersections with Mine Engineering, was to focus on asteroid composition for advancing the asteroid in-space resource utilisation domain. What do we know about asteroid composition from remote-sensing observations? What are the potential caveats in the interpretation of Earth-based spectral observations? What are the next steps to improve our knowledge on asteroid composition by means of ground-based and space-based observations and asteroid rendez-vous and sample return missions? How can asteroid mining companies use this knowledge? ASIME 2018 was a two-day workshop of almost 70 scientists and engineers in the context of the engineering needs of space missions with in-space asteroid utilisation. The 21 Questions from the asteroid mining companies were sorted into the four asteroid science themes: 1) Potential Targets, 2) Asteroid-Meteorite Links, 3) In-Situ Measurements and 4) Laboratory Measurements. The Answers to those Questions were provided by the scientists with their conference presentations and collected by A. Graps or edited directly into an open-access collaborative Google document or inserted by A. Graps using additional reference materials. During the ASIME 2018, first day and second day Wrap-Ups, the answers to the questions were discussed further. New readers to the asteroid mining topic may find the Conversation boxes and the Mission Design discussions especially interesting.

We present the results of the search for ultra-high-energy photons with nine years of data from the Telescope Array surface detector. A multivariate classifier is built upon 16 reconstructed parameters of the extensive air shower. These parameters are related to the curvature and the width of the shower front, the steepness of the lateral distribution function, and the timing parameters of the waveforms sensitive to the shower muon content. A total number of two photon candidates found in the search is fully compatible with the expected background. The $95\%\,$CL limits on the diffuse flux of the photons with energies greater than $10^{18.0}$, $10^{18.5}$, $10^{19.0}$, $10^{19.5}$ and $10^{20.0}$ eV are set at the level of $0.067$, $0.012$, $0.0036$, $0.0013$, $0.0013~\mbox{km}^{-2}\mbox{yr}^{-1}\mbox{sr}^{-1}$ correspondingly.

Gravitational waves from a merger of two neutron stars (NSs) were discovered for the first time in GW170817, together with diverse electromagnetic counterparts, providing a direct clue to the origin of short gamma-ray bursts (sGRBs). The associated sGRB 170817A was much fainter than typical, suggesting off-axis emission from a relativistic jet. However the observed prompt spectrum is inconsistent with the spectral (Amati) relation and causes the compactness problem in the simplest off-axis model. We focus on this spectral puzzle and suggest a possible key that the off-axis emission generally comes from the off-center jet, neither the jet core nor the line-of-sight jet, if the jet structure is exponentially faint outward as inferred from observations. The off-center jet could be loaded with baryon or cocoon. The off-axis model predicts that roughly $\sim 10\%$ events are brighter at smaller viewing angles than sGRB 170817A, although the exact event rate sensitively depends on uncertainties of the off-center structure. The model also predicts outliers to Amati relation, providing future tests to reveal the central engine activities.

Ultrarefractory (UR) phases in CAIs could have formed at higher T compared to common CAI minerals and thus they potentially provide constraints on very high-T processes in the solar nebula. We report a detailed characterization of an UR phase davisite bearing CAI from a reduced type CV chondrite. Absence of secondary iron- and/or alkali-rich phases and occurrence of LIME olivine indicate that primitive chemical and isotopic compositions are preserved in the CAI. Davisite occur only in one lithological unit that consists of three chemically and isotopically distinct parts: i) $^{16}$O-poor regions with reversely-zoned melilite and davisite; ii) $^{16}$O-rich regions consisting of unzoned, gehlenitic melilite, diopside and spinel; and iii) spinel framboids composed of $^{16}$O-rich spinel and $^{16}$O-poor melilite. Random distribution of chemical and isotopic heterogeneities with sharp boundaries in the CAI indicates its formation by an aggregation of mineral assemblages formed and processed separately at different time and/or space. Although isotope exchange prior to the final agglomeration of the CAI cannot be ruled out, we suggest that modification of chemical and isotopic composition of porous CAI precursors or aggregation of isotopically distinct mineral assemblages are alternative scenarios for the origin of O-isotopic heterogeneity in CAIs. In either case, coexistence of spatially and/or temporally distinct $^{16}$O-rich and -poor gaseous reservoirs at the earliest stage of the solar system formation is required. The grain-scale oxygen isotopic disequilibrium in the CAI indicate that post-formation heating of the CAI was short, which can be achieved by rapid outward transport of the CAI. High Ti$^{3+}$/Ti$^{\mathrm{tot}}$ ratios of pyroxene and presence of LIME olivine document that the entire CAI formation process took place under highly reducing conditions.

The results on ultra-high-energy cosmic rays (UHECR) mass composition obtained with the Telescope Array surface detector are presented. The analysis employs the boosted decision tree (BDT) multivariate analysis built upon 14 observables related to both the properties of the shower front and the lateral distribution function. The multivariate classifier is trained with Monte-Carlo sets of events induced by the primary protons and iron. An average atomic mass of UHECR is presented for energies $10^{18.0}-10^{20.0}\ \mbox{eV}$. The average atomic mass of primary particles shows no significant energy dependence and corresponds to $\langle \ln A \rangle = 2.0 \pm 0.1 (stat.) \pm 0.44 (syst.)$. The result is compared to the mass composition obtained by the Telescope Array with $\mbox{X}_{\mbox{max}}$ technique along with the results of other experiments. Possible systematic errors of the method are discussed.

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.

The Pierre Auger Collaboration (Auger) recently reported a correlation between the arrival directions of cosmic rays with energies above 39 EeV and the flux pattern of 23 nearby starburst galaxies (SBGs). In this Letter, we tested the same hypothesis using cosmic rays detected by the Telescope Array experiment (TA) in the 9-year period from May 2008 to May 2017. Unlike the Auger analysis, we did not optimize the parameter values but kept them fixed to the best-fit values found by Auger, namely 9.7% for the anisotropic fraction of cosmic rays assumed to originate from the SBGs in the list and 12.9{\deg} for the angular scale of the correlations. The energy threshold we adopted is 43 EeV, corresponding to 39 EeV in Auger when taking into account the energy-scale difference between two experiments. We find that the TA data is compatible with isotropy to within 1.1{\sigma} and with the Auger result to within 1.4{\sigma}, meaning that it is not capable to discriminate between these two hypotheses.

We investigate the validity of screening mechanism of the fifth force for chameleon field in highly inhomogeneous density profile. For simplicity, we consider a static and spherically symmetric system which is composed of concentric infinitely thin shells. We calculate the fifth force profile by using a numerical method for a relatively large Compton wavelength of the chameleon field. An approximate solution is also derived for the small Compton wavelength limit. Our results show that, if the thin-shell condition for the corresponding smoothed density profile is satisfied, the fifth force is safely screened outside the system irrespective of the configuration of the shells inside the system. In contrast to the screening outside the system, we find that the fifth force can be comparable to the Newtonian gravitational force inside the system. If the system is highly inhomogeneous, the chameleon field cannot trace the potential minimum varying with the density and repeats being kicked, climbing up and rolling down the potential even when the effective mass of the chameleon field is sufficiently large in the system on average. One should not feel complacent about the wellbehavedness of the fifth force field with an averaged density distribution when we consider inhomogeneous objects.

In the faint short gamma-ray burst sGRB 170817A followed by the gravitational waves (GWs) from a merger of two neutron stars (NSs) GW170817, the spectral peak energy is too high to explain only by canonical off-axis emission. We investigate off-axis appearance of an sGRB prompt emission scattered by a cocoon, which is produced through the jet-merger-ejecta interaction, with either sub-relativistic or mildly-relativistic velocities. We show that the observed properties of sGRB 170817A, in particular the high peak energy, can be consistently explained by the Thomson-scattered emission with a typical sGRB jet, together by its canonical off-axis emission, supporting that an NS-NS merger is the origin of sGRBs. The scattering occurs at $\lesssim 10^{10}$--$10^{12}\,{\rm cm}$ not far from the central engine, implying the photospheric or internal shock origin of the sGRB prompt emission. The boundary between the jet and cocoon is sharp, which could be probed by future observations of off-axis afterglows. The scattering model predicts a distribution of the spectral peak energy that is similar to the observed one but with a cutoff around $\sim$ MeV energy, and its correlations with the luminosity, duration, and time lag from GWs, providing a way to distinguish it from alternative models.

An evolving Japanese gravitational-wave (GW) mission in the deci-Hz band: B-DECIGO (DECihertz laser Interferometer Gravitational wave Observatory) will enable us to detect GW150914-like binary black holes, GW170817-like binary neutron stars, and intermediate-mass binary black holes out to cosmological distances. The B-DECIGO band slots in between the aLIGO-Virgo-KAGRA-IndIGO (hecto-Hz) and LISA (milli-Hz) bands for broader bandwidth; the sources described emit GWs for weeks to years across the multiband to accumulate high signal-to-noise ratios. This suggests the possibility that joint detection would greatly improve the parameter estimation of the binaries. We examine B-DECIGO's ability to measure binary parameters and assess to what extent multiband analysis could improve such measurement. Using non-precessing post-Newtonian waveforms with the Fisher matrix approach, we find for systems like GW150914 and GW170817 that B-DECIGO can measure the mass ratio to within $< 0.1\%$, the individual black-hole spins to within $< 10\%$, and the coalescence time to within $< 5\,$s about a week before alerting aLIGO and electromagnetic facilities. Prior information from B-DECIGO for aLIGO can further reduce the uncertainty in the measurement of, e.g., certain neutron star tidally-induced deformations by factor of $\sim 6$, and potentially determine the spin-induced neutron star quadrupole moment. Joint LISA and B-DECIGO measurement will also be able to recover the masses and spins of intermediate-mass binary black holes at percent-level precision. However, there will be a large systematic bias in these results due to post-Newtonian approximation of exact GW signals.

We demonstrate the dragging of the magnetic field by the super-Alfvenic shear flows out of the reconnection plane can strongly localize the reconnection x-line in collisionless plasmas, reversing the current direction at the x-line. Reconnection with this new morphology, which is impossible in resistive-magnetohydrodynamic (MHD), is enabled by electron inertia. Surprisingly, the quasi-steady reconnection rate remains of order 0.1 even though the aspect ratio of the local x-line geometry is larger than unity. We explain this by examining the transport of the reconnected magnetic flux and the opening angle made by the upstream magnetic field, concluding that the reconnection rate is still limited by the constraint imposed at the inflow region. This study further suggests the nearly universal fast rate value of order 0.1 cannot be explained by the physics of tearing modes, nor can it be explained by a universal localization mechanism.

Bursts of gamma ray showers have been observed in coincidence with downward propagating negative leaders in lightning flashes by the Telescope Array Surface Detector (TASD). The TASD is a 700~square kilometer cosmic ray observatory located in southwestern Utah, U.S.A. In data collected between 2014 and 2016, correlated observations showing the structure and temporal development of three shower-producing flashes were obtained with a 3D lightning mapping array, and electric field change measurements were obtained for an additional seven flashes, in both cases co-located with the TASD. National Lightning Detection Network (NLDN) information was also used throughout. The showers arrived in a sequence of 2--5 short-duration ($\le$10~$\mu$s) bursts over time intervals of several hundred microseconds, and originated at an altitude of $\simeq$3--5 kilometers above ground level during the first 1--2 ms of downward negative leader breakdown at the beginning of cloud-to-ground lightning flashes. The shower footprints, associated waveforms and the effect of atmospheric propagation indicate that the showers consist primarily of downward-beamed gamma radiation. This has been supported by GEANT simulation studies, which indicate primary source fluxes of $\simeq$$10^{12}$--$10^{14}$ photons for $16^{\circ}$ half-angle beams. We conclude that the showers are terrestrial gamma-ray flashes (TGFs), similar to those observed by satellites, but that the ground-based observations are more representative of the temporal source activity and are also more sensitive than satellite observations, which detect only the most powerful TGFs.

The Telescope Array observatory utilizes fluorescence detectors and surface detectors to observe air showers produced by ultra high energy cosmic rays in the Earth's atmosphere. Cosmic ray events observed in this way are termed hybrid data. The depth of air shower maximum is related to the mass of the primary particle that generates the shower. This paper reports on shower maxima data collected over 8.5 years using the Black Rock Mesa and Long Ridge fluorescence detectors in conjunction with the array of surface detectors. We compare the means and standard deviations of the observed $X_{\mathrm{max}}$ distributions with Monte Carlo $X_{\mathrm{max}}$ distributions of unmixed protons, helium, nitrogen, and iron, all generated using the QGSJet~II-04 hadronic model. We also perform an unbinned maximum likelihood test of the observed data, which is subjected to variable systematic shifting of the data $X_{\mathrm{max}}$ distributions to allow us to test the full distributions, and compare them to the Monte Carlo to see which elements are not compatible with the observed data. For all energy bins, QGSJet~II-04 protons are found to be compatible with Telescope Array hybrid data at the 95% confidence level after some systematic $X_{\mathrm{max}}$ shifting of the data. Three other QGSJet~II-04 elements are found to be compatible using the same test procedure in an energy range limited to the highest energies where data statistics are sparse.

One of the uncertainties in interpretation of ultra-high energy cosmic ray (UHECR) data comes from the hadronic interaction models used for air shower Monte Carlo (MC) simulations. The number of muons observed at the ground from UHECR-induced air showers is expected to depend upon the hadronic interaction model. One may therefore test the hadronic interaction models by comparing the measured number of muons with the MC prediction. In this paper, we present the results of studies of muon densities in UHE extensive air showers obtained by analyzing the signal of surface detector stations which should have high $\it{muon \, purity}$. The muon purity of a station will depend on both the inclination of the shower and the relative position of the station. In 7 years' data from the Telescope Array experiment, we find that the number of particles observed for signals with an expected muon purity of $\sim$65% at a lateral distance of 2000 m from the shower core is $1.72 \pm 0.10{\rm (stat.)} \pm 0.37 {\rm (syst.)}$ times larger than the MC prediction value using the QGSJET II-03 model for proton-induced showers. A similar effect is also seen in comparisons with other hadronic models such as QGSJET II-04, which shows a $1.67 \pm 0.10 \pm 0.36$ excess. We also studied the dependence of these excesses on lateral distances and found a slower decrease of the lateral distribution of muons in the data as compared to the MC, causing larger discrepancy at larger lateral distances.

Gravitational waves from a merger of two neutron stars (NSs) were discovered for the first time in GW170817, together with diverse electromagnetic (EM) counterparts. To make constraints on a relativistic jet from the NS merger, we calculate the EM signals in (1) the short gamma-ray burst sGRB 170817A from an off-axis jet, (2) the optical-infrared macronova (or kilonova), especially the blue macronova, from a jet-powered cocoon, and (3) the X-ray and radio afterglows from the interaction between the jet and interstellar medium. We find that a typical sGRB jet is consistent with these observations, and there is a parameter space to explain all the observations in a unified fashion with an isotropic energy $\sim 10^{51}$-$10^{52}$ erg, opening angle $\sim 20^{\circ}$, and viewing angle $\sim 30^{\circ}$. The off-axis emission is less de-beamed than the point-source case because the viewing angle is comparable to the opening angle. We also analytically show that the jet energy accelerates a fair fraction of the merger ejecta to a sub-relativistic velocity $\sim 0.3$-$0.4$c as a cocoon in a wide parameter range. The ambient density might be low $\sim 10^{-3}$-$10^{-6}$ cm$^{-3}$, which can be tested by future observations of radio flares and X-ray remnants.

We report the results of a neutrino search in Super-Kamiokande for coincident signals with the first detected gravitational wave produced by a binary neutron star merger, GW170817, which was followed by a short gamma-ray burst, GRB170817A, and a kilonova/macronova. We searched for coincident neutrino events in the range from 3.5 MeV to $\sim$100 PeV, in a time window $\pm$500 seconds around the gravitational wave detection time, as well as during a 14-day period after the detection. No significant neutrino signal was observed for either time window. We calculated 90% confidence level upper limits on the neutrino fluence for GW170817. From the upward-going-muon events in the energy region above 1.6 GeV, the neutrino fluence limit is $16.0^{+0.7}_{-0.6}$ ($21.3^{+1.1}_{-0.8}$) cm$^{-2}$ for muon neutrinos (muon antineutrinos), with an error range of $\pm5^{\circ}$ around the zenith angle of NGC4993, and the energy spectrum is under the assumption of an index of $-2$. The fluence limit for neutrino energies less than 100 MeV, for which the emission mechanism would be different than for higher-energy neutrinos, is also calculated. It is $6.6 \times 10^7$ cm$^{-2}$ for anti-electron neutrinos under the assumption of a Fermi-Dirac spectrum with average energy of 20 MeV.

We report on a measurement of the cosmic ray energy spectrum by the Telescope Array Low-Energy Extension (TALE) air fluorescence detector. The TALE air fluorescence detector is also sensitive to the Cherenkov light produced by shower particles. Low energy cosmic rays, in the PeV energy range, are detectable by TALE as "Cherenkov Events". Using these events, we measure the energy spectrum from a low energy of $\sim 2$ PeV to an energy greater than 100 PeV. Above 100 PeV TALE can detect cosmic rays using air fluorescence. This allows for the extension of the measurement to energies greater than a few EeV. In this paper, we will describe the detector, explain the technique, and present results from a measurement of the spectrum using $\sim 1000$ hours of observation. The observed spectrum shows a clear steepening near $10^{17.1}$ eV, along with an ankle-like structure at $10^{16.2}$ eV. These features present important constraints on galactic cosmic rays origin and propagation models. The feature at $10^{17.1}$ eV may also mark the end of the galactic cosmic rays flux and the start of the transition to extra-galactic sources.

The counter jet of a short gamma-ray burst (sGRB) has not yet been observed, while recent discoveries of gravitational waves (GWs) from a binary neutron star (NS) merger GW170817 and the associated sGRB 170817A have demonstrated that off-axis sGRB jets are detectable. We calculate the prompt emission from the counter jet of an sGRB and show that it is typically 23-26 mag in the optical-infrared band 10-10^3 sec after the GWs for an sGRB 170817A-like event, which is brighter than the early macronova (or kilonova) emission and detectable by LSST in the near future. We also propose a new method to constrain the unknown jet properties, such as the Lorentz factor, opening angle, emission radii, and jet launch time, by observing both the forward and counter jets. To scrutinize the counter jets, space GW detectors like DECIGO are powerful in forecasting the merger time (<~ 1 sec) and position (<~ 1 arcmin) (~ a week) before the merger.

We investigate the physical conditions of the CO gas near the young star cluster, Trumpler 14 of the Carina Nebula. The observations presented in this work are taken with the Fourier Transform Spectrometer (FTS) of the Spectral and Photometric Imaging REceiver (SPIRE) onboard the Herschel Space Observatory. Our field of view covers the edge of a cavity carved by Trumpler 14 about $1\,\mathrm{Myr}$ ago and marks the transition from HII regions to photo-dissociation regions. With the state-of-the-art Meudon PDR code, we successfully derive the physical conditions, which include the thermal pressure ($P$) and the scaling factor of radiation fields ($G_{\mathrm{UV}}$), from the observed CO spectral line energy distributions~(SLEDs) in the observed region. The derived $G_{\mathrm{UV}}$ values generally show an excellent agreement with the UV radiation fields created by nearby OB-stars and thus confirm that the main excitation source of the observed CO emission are the UV-photons provided by the massive stars. The derived thermal pressure is between $0.5-3\,\times\,10^{8}\,\mathrm{K\,cm^{-3}}$ with the highest values found along the ionization front in Car I-E region facing Trumpler 14, hinting that the cloud structure is similar to the recent observations of the Orion Bar. Comparing the derived thermal pressure with the radiation fields, we report the first observationally-derived and spatially-resolved $P \sim 2\times10^4\,G_{\mathrm{UV}}$ relationship. As direct comparisons of the modeling results to the observed $^{13}\mathrm{CO}$, [OI] $63\,\mathrm{\mu m}$, and [CII] $158\,\mathrm{\mu m}$ intensities are not straightforward, we urge the readers to be cautious when constraining the physical conditions of PDRs with combinations of $^{12}\mathrm{CO}$, $^{13}\mathrm{CO}$, [CI], [OI] $63\,\mathrm{\mu m}$, and [CII] $158\,\mathrm{\mu m}$ observations.

An intermediate-scale energy spectrum anisotropy has been found in the arrival directions of ultra-high energy cosmic rays of energies above $10^{19.2}$ eV in the northern hemisphere, using 7 years of data from the Telescope Array surface detector. A relative energy distribution test is done comparing events inside oversampled spherical caps of equal exposure, to those outside, using the Poisson likelihood ratio. The center of maximum significance is at $9^h$$16^m$, $45^{\circ}$. and has a deficit of events with energies $10^{19.2}$$\leq$$E$$<$$10^{19.75}$ eV and an excess for $E$$\geq$$10^{19.75}$ eV. The post-trial probability of this energy anisotropy, appearing by chance anywhere on an isotropic sky, is found by Monte Carlo simulation to be $9$$\times$$10^{-5}$ ($3.74$$\sigma_{global}$).

The energy of the ultrahigh energy spectral cutoff was measured, integrating over the northern hemisphere sky, by the Telescope Array (TA) collaboration, to be $10^{19.78 \pm 0.06}$ eV, in agreement with the High Resolution Fly's Eye (HiRes) experiment, whereas the Pierre Auger experiment, integrating over the southern hemisphere sky, measured the cutoff to be at 10$^{19.62 \pm 0.02}$ eV. An 11% energy scale difference between the TA and Auger does not account for this difference. However, in comparing the spectra of the Telescope Array and Pierre Auger experiments in the band of declination common to both experiments ($-15.7^{\circ} < \delta < 24.8^{\circ}$) we have found agreement in the energy of the spectral cutoff. While the Auger result is essentially unchanged, the TA cutoff energy has changed to $10^{19.59 \pm 0.06}$ eV. In this paper we argue that this is an astrophysical effect.

We present the results of a search in the Super-Kamiokande (SK) detector for excesses of neutrinos with energies above a few GeV that are in the direction of the track events reported in IceCube. Data from all SK phases (SK-I through SK-IV) were used, spanning a period from April 1996 to April 2016 and corresponding to an exposure of 225 kilotonne-years . We considered the 14 IceCube track events from a data set with 1347 livetime days taken from 2010 to 2014. We use Poisson counting to determine if there is an excess of neutrinos detected in SK in a 10 degree search cone (5 degrees for the highest energy data set) around the reconstructed direction of the IceCube event. No significant excess was found in any of the search directions we examined. We also looked for coincidences with a recently reported IceCube multiplet event. No events were detected within a $\pm$ 500 s time window around the first detected event, and no significant excess was seen from that direction over the lifetime of SK.

The origin of Phobos and Deimos in a giant impact generated disk is gaining larger attention. Although this scenario has been the subject of many studies, an evaluation of the chemical composition of the Mars' moons in this framework is missing. The chemical composition of Phobos and Deimos is unconstrained. The large uncertainty about the origin of the mid-infrared features, the lack of absorption bands in the visible and near-infrared spectra, and the effects of secondary processes on the moons' surface make the determination of their composition very difficult from remote sensing data. Simulations suggest a formation of a disk made of gas and melt with their composition linked to the nature of the impactor and Mars. Using thermodynamic equilibrium we investigate the composition of dust (condensates from gas) and solids (from a cooling melt) that result from different types of Mars impactors (Mars-, CI-, CV-, EH-, comet-like). Our calculations show a wide range of possible chemical compositions and noticeable differences between dust and solids depending on the considered impactors. Assuming Phobos and Deimos as result of the accretion and mixing of dust and solids, we find that the derived assemblage (dust rich in metallic-iron, sulphides and/or carbon, and quenched solids rich in silicates) can be compatible with the observations. The JAXA's MMX (Martian Moons eXploration) mission will investigate the physical and chemical properties of the Maroons, especially sampling from Phobos, before returning to Earth. Our results could be then used to disentangle the origin and chemical composition of the pristine body that hit Mars and suggest guidelines for helping in the analysis of the returned samples.

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.

We present the results of near- to mid-infrared (NIR to MIR) imaging and NIR spectroscopic observations of two galaxy mergers, NGC 2782 (Arp 215) and NGC 7727 (Arp 222), with the Infrared Camera on board AKARI. NGC 2782 shows extended MIR emission in the eastern side of the galaxy, which corresponds to the eastern tidal tail seen in the HI 21 cm map, while NGC 7727 shows extended MIR emission in the north of the galaxy, which is similar to the plumes seen in the residual image at the K-band after subtracting a galaxy model. Both extended structures are thought to have formed associated with their merger events. They show excess emission at 7--15 micron, which can be attributed to emission from polycyclic aromatic hydrocarbons (PAHs), while the observed spectral energy distributions decline longward of 24 micron, suggesting that very small grains (VSGs) are deficient. These characteristics of the observed MIR spectral energy distribution may be explained if PAHs are formed by fragmentation of VSGs during merger events. The star formation rate is estimated from the MIR PAH emission in the eastern tail region of NGC 2782 and it is in fair agreement with those estimated from Halpha and [CII] 158 micron. MIR observations are efficient for the study of dust processing and structures formed during merger events.

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.

We propose a method for confirmation of the existence of Population III (Pop III) stars with massive black hole binaries as GW150914 in gravitational wave (GW) observation. When we get enough number of events, we want to determine which model is closer to reality, with and without Pop III stars. We need to prepare various "Pop I/II models" and various "Pop I/II/III models" and investigate which model is consistent with the events. To demonstrate our analysis, we simulate detections of GW events for some examples of population synthesis models with and without Pop III stars. We calculate the likelihood ratio with the realistic number of events and evaluate the probability of identifying the existence of Pop III stars. In typical cases, our analysis can distinguish between Pop I/II model and Pop I/II/III model with 90% probability by 22 GW signals from black hole-black hole binary mergers.

The Telescope Array (TA) experiment is located in the western desert of Utah, USA, and observes ultra high energy cosmic rays (UHECRs) in the Northern hemisphere. At the highest energies, $E>10$~EeV, the shape of cosmic ray energy spectrum may carry an imprint of the source density distribution along the line of sight different in different directions of the sky. In this study, we search for such directional variations in the shape of the energy spectrum using events observed with the Telescope Array's surface detector. We divide the TA field of view into two nearly equal-exposure regions: the "on-source" region which we define as $\pm 30^\circ$ of the supergalactic plane containing mostly nearby structures, and the complementary "off-source" region where the sources are further away on average. We compare the UHECR spectra in these regions by fitting them to the broken power law and comparing the resulting parameters. We find that the off-source spectrum has an earlier break at highest energies. The chance probability to obtain such or larger difference in statistically equivalent distributions is estimated as $6.2\pm1.1\times10^{-4}$ ($3.2\sigma$) by a Monte-Carlo simulation. The observed difference in spectra is in a reasonable quantitative agreement with a simplified model that assumes that the UHECR sources trace the galaxy distribution from the 2MRS catalogue, primary particles are protons and the magnetic deflections can be neglected.

The Sasaki-Nakamura transformation gives a short-ranged potential and a convergent source term for the master equation of perturbations in the Kerr space-time. In this paper, we study the asymptotic behavior of the transformation, and present a new relaxed necessary and sufficient condition of the transformation to obtain the short-ranged potential in the assumption that the transformation converges in the far distance. Also, we discuss quasinormal mode frequencies which are determined by the information around the peak of the potential in the WKB analysis. Finally, in the extreme Kerr limit, $a/M \to 1$, where $M$ and $a$ denote the mass and spin parameter of a Kerr black hole, respectively, we find the peak location of the potential, $r_p/M \lesssim 1 + 1.8 \,(1-a/M)^{1/2}$ by using the new transformation. The uncertainty of the location is as large as that expected from the equivalence principle.

Apollo-type NEA (3200) Phaethon, classified at the B/F-type taxonomy, probably the main mass of the Phaethon-Geminid stream complex (PGC), can be the most metamorphic C-complex asteroid in our solar system, since it is heated up to ~1000 K by the solar radiation around its perihelion passages. Hence, its surface material may be easily decomposed in near-sun environment. Phaethon's spectrum exhibits extremely blue-slope in the VIS-NIR region (so-called Phaethon Blue). Another candidate large member of the PGC, (155140) 2005 UD, shows a B/F-type color, however with a C-type-like red color over its ~1/4 rotational part, which implies an exposition of less metamorphosed primordial internal structure of the PGC precursor by a splitting or breakup event long ago. If so, some rotational part of Phaethon should show the C-type color as well as 2005 UD. Hence, we carried out the time-series VIS-spectroscopic observations of Phaethon using 1-m telescope in order to detect such a signature. Also, R-band photometries were simultaneously performed in order to complement our spectroscopy. Consequently, we obtained a total of 68 VIS-spectrophotometric data, 78% of which show the B-type blue-color, as against the rest of 22% showing the C-type red-color. We successfully acquired rotationally time-resolved spectroscopic data, of which particular rotational phase shows a red-spectral slope as the C-type color, as 2005 UD does, suggesting longitudinal inhomogeneity on Phaethon's surface. We constrained this C-type red-colored area in the mid-latitude in Phaethon's southern hemisphere based on the rotationally time-resolved spectroscopy along with Phaethon's axial rotation state, of which size suggests the impact-induced origin of the PGC. We also surveyed the meteoritic analog of Phaethon's surface blue-color, and found thermally metamorphosed CI/CM chondrites as likely candidates.

Ordinary chondrites (OCs) are by far the most abundant meteorites (80% of all falls). Their origin has long been the matter of a heated debate. About thirty years ago (e.g., Pellas 1988), it was proposed that OCs should originate from S-type bodies (the most abundant asteroid spectral types in the inner part of the asteroid belt), but the apparent discrepancy between S-type asteroid and OC reflectance spectra generated what was known as the S-type--OC conundrum. This paradox has gradually been resolved over the years. It is now understood that space weathering processes are responsible for the spectral mismatch between S-type bodies and OCs. Furthermore, both telescopic observations and the first asteroid sample return mission (Hayabusa) indicate that most S-type bodies have mineralogies similar to those of OCs. Importantly, the S-type/OC link, which has remained sterile for more than 30 years, has been delivering fundamental constraints on the formation and evolution of planetesimals over the recent years.

Cosmic rays in the energy range $10^{18.0}$ - $10^{18.5}$ eV are thought to have a light, probably protonic, composition. To study their origin one can search for anisotropy in their arrival directions. Extragalactic cosmic rays should be isotropic, but galactic cosmic rays of this type should be seen mostly along the galactic plane, and there should be a shortage of events coming from directions near the galactic anticenter. This is due to the fact that, under the influence of the galactic magnetic field, the transition from ballistic to diffusive behavior is well advanced, and this qualitative picture persists over the whole energy range. Guided by models of the galactic magnetic field that indicate that the enhancement along the galactic plane should have a standard deviation of about 20$^\circ$ in galactic latitude, and the deficit in the galactic anticenter direction should have a standard deviation of about 50$^\circ$ in galactic longitude, we use the data of the Telescope Array surface detector in $10^{18.0}$ to $10^{18.5}$ eV energy range to search for these effects. The data are isotropic. Neither an enhancement along the galactic plane nor a deficit in the galactic anticenter direction is found. Using these data we place an upper limit on the fraction of EeV cosmic rays of galactic origin at 1.3% at 95% confidence level.

Supermassive black holes observed at high redshift $z\gtrsim6$ could grow from direct collapse black holes (DCBHs) with mass $\sim10^5\,M_{\odot}$, which result from the collapse of supermassive stars (SMSs). If a relativistic jet is launched from a DCBH, it can break out of the collapsing SMS and produce a gamma-ray burst (GRB). Although most of the GRB jets are off-axis from our line of sight, we show that the energy injected from the jet into a cocoon is huge $\sim10^{55-56}\,{\rm{erg}}$, so that the cocoon fireball is observed as ultra-luminous supernovae of $\sim10^{45-46}\rm{\,erg\,s^{-1}}$ for $\sim5000 [(1+z)/16] \rm{\,days}$. They are detectable by the future telescopes with near infrared bands, such as, $Euclid$, $WFIRST$, $WISH$, and $JWST$ up to $z\sim20$ and $\sim 100$ events per year, providing a direct evidence of the DCBH scenario.

In the population synthesis simulations of Pop III stars, many BH (Black Hole)-BH binaries with merger time less than the age of the Universe $(\tau_{\rm H})$ are formed, while NS (Neutron Star)-BH binaries are not. The reason is that Pop III stars have no metal so that no mass loss is expected. Then, in the final supernova explosion to NS, much mass is lost so that the semi major axis becomes too large for Pop III NS-BH binaries to merge within $\tau_{\rm H}$. However it is almost established that the kick velocity of the order of $200-500{\rm~ km~s^{-1}}$ exists for NS from the observation of the proper motion of the pulsar. Therefore, the semi major axis of the half of NS-BH binaries can be smaller than that of the previous argument for Pop III NS-BH binaries to decrease the merging time. We perform population synthesis Monte Carlo simulations of Pop III NS-BH binaries including the kick of NS and find that the event rate of Pop III NS-BH merger rate is $\sim 1 {\rm Gpc^{-3} yr^{-1}}$. This suggests that there is a good chance of the detection of Pop III NS-BH mergers in O2 of Advanced LIGO and Advanced Virgo from this autumn.

We report the results from a search in Super-Kamiokande for neutrino signals coincident with the first detected gravitational wave events, GW150914 and GW151226, using a neutrino energy range from 3.5 MeV to 100 PeV. We searched for coincident neutrino events within a time window of $\pm$500 seconds around the gravitational wave detection time. Four neutrino candidates are found for GW150914 and no candidates are found for GW151226. The remaining neutrino candidates are consistent with the expected background events. We calculated the 90\% confidence level upper limits on the combined neutrino fluence for both gravitational wave events, which depends on event energy and topologies. Considering the upward going muon data set (1.6 GeV - 100 PeV) the neutrino fluence limit for each gravitational wave event is 14 - 37 (19 - 50) cm$^{-2}$ for muon neutrinos (muon antineutrinos), depending on the zenith angle of the event. In the other data sets, the combined fluence limits for both gravitational wave events range from 2.4$\times 10^{4}$ to 7.0$\times 10^{9}$ cm$^{-2}$.

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

Focusing on the remnant black holes after merging binary black holes, we show that ringdown gravitational waves of Population III binary black holes mergers can be detected with the rate of $5.9-500~{\rm events~yr^{-1}}~({\rm SFR_p}/ (10^{-2.5}~M_\odot~{\rm yr^{-1}~Mpc^{-3}})) \cdot ({\rm [f_b/(1+f_b)]/0.33})$ for various parameters and functions. This rate is estimated for the events with SNR$>8$ for the second generation gravitational wave detectors such as KAGRA. Here, ${\rm SFR_p}$ and ${\rm f_b}$ are the peak value of the Population III star formation rate and the fraction of binaries, respectively. When we consider only the events with SNR$>35$, the event rate becomes $0.046-4.21~{\rm events~yr^{-1}}~({\rm SFR_p}/ (10^{-2.5}~M_\odot~{\rm yr^{-1}~Mpc^{-3}})) \cdot ({\rm [f_b/(1+f_b)]/0.33})$. This suggest that for remnant black hole's spin $q_f>0.95$ we have the event rate with SNR$>35$ less than $0.037~{\rm events~yr^{-1}}~({\rm SFR_p}/ (10^{-2.5}~M_\odot~{\rm yr^{-1}~Mpc^{-3}})) \cdot ({\rm [f_b/(1+f_b)]/0.33})$, while it is $3-30~{\rm events~yr^{-1}}~({\rm SFR_p}/ (10^{-2.5}~M_\odot~{\rm yr^{-1}~Mpc^{-3}})) \cdot ({\rm [f_b/(1+f_b)]/0.33})$ for the third generation detectors such as Einstein Telescope. If we detect many Population III binary black holes merger, it may be possible to constrain the Population III binary evolution paths not only by the mass distribution but also by the spin distribution.

TARA (Telescope Array Radar) is a cosmic ray radar detection experiment colocated with Telescope Array, the conventional surface scintillation detector (SD) and fluorescence telescope detector (FD) near Delta, Utah, U.S.A. The TARA detector combines a 40 kW, 54.1 MHz VHF transmitter and high-gain transmitting antenna which broadcasts the radar carrier over the SD array and within the FD field of view, towards a 250 MS/s DAQ receiver. TARA has been collecting data since 2013 with the primary goal of observing the radar signatures of extensive air showers (EAS). Simulations indicate that echoes are expected to be short in duration (~10 microseconds) and exhibit rapidly changing frequency, with rates on the order of 1 MHz/microsecond. The EAS radar cross-section (RCS) is currently unknown although it is the subject of over 70 years of speculation. A novel signal search technique is described in which the expected radar echo of a particular air shower is used as a matched filter template and compared to waveforms obtained by triggering the radar DAQ using the Telescope Array fluorescence detector. No evidence for the scattering of radio frequency radiation by EAS is obtained to date. We report the first quantitative RCS upper limits using EAS that triggered the Telescope Array Fluorescence Detector.

Recent numerical relativity simulations of mergers of binary black holes suggest that the maximum final value of $a/M$ is $\sim 0.95$ for the coalescence of two equal mass black holes with aligned spins of the same magnitude $a/M=0.994$ which is close to the upper limit $a/M=0.998$ of accretion spin-up shown by Thorne. Using the WKB method, we suggest that if quasinormal modes with $a/M \sim 0.95$ are detected by the second generation gravitational wave detectors, we could confirm the strong gravity space-time based on Einstein's general relativity up to $1.33M$ which is only $\sim 1.014$ times the event horizon radius and within the ergoregion. One more message about black hole geometry is expected here. If the quasinormal mode is different from that of general relativity, we need to find the true theory of gravity which deviates from general relativity only near the black hole horizon.

Using the WKB method, we show that the peak location ($r_{\rm peak}$) of the potential, which determines the quasinormal mode frequency of the Kerr black hole, obeys an accurate empirical relation as a function of the specific angular momentum $a$ and the gravitational mass $M$. If the quasinormal mode with $a/M \sim 1$ is observed by gravitational wave detectors, we can confirm the black-hole space-time around the event horizon, $r_{\rm peak}=r_+ +O(\sqrt{1-q})$ where $r_+$ is the event horizon radius. While if the quasinormal mode is different from that of general relativity, we are forced to seek the true theory of gravity and/or face to the existence of the naked singularity.

The existence of the ergoregion of the Kerr space-time has not been confirmed observationally yet. We show that the confirmation would be possible by observing the quasinormal mode in gravitational waves. As an example, using the recent population synthesis results of Pop III binary black holes, we find that the peak of the final merger mass ($M_f$) is about $50~\rm M_{\odot}$, while the fraction of the final spin $q_f = a_f/M_f > 0.7$ needed for the confirmation of a part of ergoregion is $\sim 77\%$. To confirm the frequency of the quasinormal mode, ${\rm SNR} > 35$ is needed. The standard model of Pop III population synthesis tells us that the event rate for the confirmation of more than $50\%$ of the ergoregion by the second generation gravitational wave detectors is $\sim 2.3$ ${\rm events\ yr^{-1}\ (SFR_p/(10^{-2.5}\ M_\odot yr^{-1}\ Mpc^{-3}))} \cdot (\rm [f_b/(1+f_b)]/0.33)$ where ${\rm SFR_p}$ and ${\rm f_b}$ are the peak value of the Pop III star formation rate and the fraction of binaries, respectively.

The Karin family is a very young asteroid family created by an asteroid breakup 5.8 Myr ago. Since the members of this family probably have not experienced significant orbital or collisional evolution yet, it is possible that they still preserve properties of the original family-forming event in terms of their spin state. As we carried out a series of photometric observations of the Karin family asteroids, here we report an analysis result of lightcurves including the rotation period of eleven members. The mean rotation rate of the Karin family members turned out to be much lower than those of NEAs or smaller MBAs (diameter D<12 km), and even lower than that of larger MBAs (D>130 km). We investigated a correlation between the peak-to-peak variation magnitude reduced to zero solar phase angle and the rotation period of the eleven Karin family asteroids, and found a possible trend that elongated members have lower spin rate, and less elongated members have higher spin rate. However, this has to be confirmed by another series of future observations.

Recent population synthesis simulations of Pop III stars suggest that the event rate of coalescence of $\sim 30M_\odot$--$30M_\odot$ binary black holes can be high enough for the detection by the second generation gravitational wave detectors. The frequencies of chirp signal as well as quasinormal modes are near the best sensitivity of these detectors so that it would be possible to confirm Einstein's general relativity. Using the WKB method, we suggest that for the typical value of spin parameter $a/M\sim 0.7$ from numerical relativity results of the coalescence of binary black holes, the strong gravity of the black hole space-time at around the radius $2M$, which is just $\sim 1.17$ times the event horizon radius, would be confirmed as predicted by general relativity. The expected event rate with the signal-to-noise ratio $> 35$ needed for the determination of the quasinormal mode frequency with the meaningful accuracy is $0.17$--$7.2$~${\rm events~yr^{-1}~(SFR_p/(10^{-2.5}~M_\odot~yr^{-1}~Mpc^{-3}))} \cdot (\rm [f_b/(1+f_b)]/0.33)$ where ${\rm SFR_p}$ and ${\rm f_b}$ are the peak value of the Pop III star formation rate and the fraction of binaries, respectively. As for the possible optical counter part, if the merged black hole of mass $M\sim 60M_\odot$ is in the interstellar matter with $n\sim 100~{\rm cm^{-3}}$ and the proper motion of black hole is $\sim 1~{\rm km~s^{-1}}$, the luminosity is $\sim 10^{40}~{\rm erg~ s^{-1}}$ which can be detected up to $\sim 300~{\rm Mpc}$, for example, by Subaru-HSC and LSST with the limiting magnitude 26.

This paper presents the results of different searches for correlations between very high-energy neutrino candidates detected by IceCube and the highest-energy cosmic rays measured by the Pierre Auger Observatory and the Telescope Array. We first consider samples of cascade neutrino events and of high-energy neutrino-induced muon tracks, which provided evidence for a neutrino flux of astrophysical origin, and study their cross-correlation with the ultrahigh-energy cosmic ray (UHECR) samples as a function of angular separation. We also study their possible directional correlations using a likelihood method stacking the neutrino arrival directions and adopting different assumptions on the size of the UHECR magnetic deflections. Finally, we perform another likelihood analysis stacking the UHECR directions and using a sample of through-going muon tracks optimized for neutrino point-source searches with sub-degree angular resolution. No indications of correlations at discovery level are obtained for any of the searches performed. The smallest of the p-values comes from the search for correlation between UHECRs with IceCube high-energy cascades, a result that should continue to be monitored.

Extended emission is a mystery in short gamma-ray bursts (SGRBs). By making time resolved spectral analyses of brightest nine events observed by ${\it Swift}$ XRT, we obviously classify the early X-ray emission of SGRBs into two types. One is the extended emission with exponentially rapid decay, which shows significant spectral softening during hundreds seconds since the SGRB trigger and is also detected by ${\it Swift}$-BAT. The other is a dim afterglow only showing power-law decay over $10^4$ s. The correlations between the temporal decay and spectral indices of the extended emissions are inconsistent with the $\alpha$-$\beta$ correlation expected for the high-latitude curvature emission from a uniform jet. The observed too-rapid decay suggests the emission from a photosphere or a patchy surface, and manifests the stopping central engine via such as magnetic reconnection at the black hole.

Using our population synthesis code, we found that the typical chirp mass of Pop III BH-BHs is $\sim30~\msun$ with the total mass of $\sim60~\msun$ so that the inspiral chirp signal as well as quasi normal mode (QNM) of the merging BH are interesting targets of KAGRA. The detection rate of the coalescing Pop III BH-BHs is $\sim$180 $\rm events~yr^{-1}$$(\rm SFR_p/(10^{-2.5}~\msun \rm~yr^{-1}~Mpc^{-3}))\cdot([f_b/(1+f_b)]/0.33)\cdot Err_{sys}$ in our standard model where $\rm SFR_{p},~f_b$ and $\rm Err_{sys}$ are the peak value of the Pop III star formation rate, the binary fraction and the systematic error with $\rm Err_{sys}=1$ for our standard model, respectively. To evaluate the robustness of chirp mass distribution and the range of $\rm Err_{sys}$, we examine the dependence of the results on the unknown parameters and the distribution functions. We found that the chirp mass has a peak at $\sim 30 ~\msun$ in most of parameters and distribution functions as well as $\rm Err_{sys}$ ranges from 0.046 to 4. The minimum rate corresponds to the worst model which we think unlikely so that unless $ {\rm ~(SFR_p/(10^{-2.5}~\msun~yr^{-1}~Mpc^{-3}))\cdot([f_b/(1+f_b)]/0.33) \ll 0.1}$, we expect the Pop III BH-BHs merger rate of at least one event per year by KAGRA. Nakano, Tanaka \& Nakamura (2015) show that if S/N of QNM is larger than 35, we can confirm or refute the General Relativity (GR) more than 5 sigma level. In our standard model, the detection rate of Pop III BH-BHs whose S/N is larger than 35 is $3.2~\rm events~yr^{-1}$$(\rm SFR_p/(10^{-2.5}~\msun \rm~yr^{-1}~Mpc^{-3}))\cdot([f_b/(1+f_b)]/0.33)\cdot Err_{sys}$. Thus, there is a good chance to check whether GR is correct or not in the strong gravity region.

Joint contributions of the Pierre Auger Collaboration and the Telescope Array Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.

We have conducted three searches for correlations between ultra-high energy cosmic rays detected by the Telescope Array and the Pierre Auger Observatory, and high-energy neutrino candidate events from IceCube. Two cross-correlation analyses with UHECRs are done: one with 39 cascades from the IceCube `high-energy starting events' sample and the other with 16 high-energy `track events'. The angular separation between the arrival directions of neutrinos and UHECRs is scanned over. The same events are also used in a separate search using a maximum likelihood approach, after the neutrino arrival directions are stacked. To estimate the significance we assume UHECR magnetic deflections to be inversely proportional to their energy, with values $3^\circ$, $6^\circ$ and $9^\circ$ at 100 EeV to allow for the uncertainties on the magnetic field strength and UHECR charge. A similar analysis is performed on stacked UHECR arrival directions and the IceCube sample of through-going muon track events which were optimized for neutrino point-source searches.

The existence of black holes (BHs) of mass ~ 10^{9} M_sun at z > 6 is a big puzzle in astrophysics because even optimistic estimates of the accretion time are insufficient for stellar mass BHs of ~ 10 M_sun to grow into such supermassive BHs. A resolution of this puzzle might be the direct collapse of supermassive stars with mass M ~ 10^{5} M_sun into massive seed BHs. We find that if a jet is launched from the accretion disk around the central BH, the jet can break out the star because of the structure of the radiation pressure-dominated envelope. Such ultra-long gamma-ray bursts with duration of ~ 10^{4} - 10^{6} s and flux of 10^{-11} - 10^{-8} erg s^{-1} cm^{-2} could be detectable by Swift. We estimate an event rate of < 1yr^{-1}. The total explosion energy is > 10^{55} - 10^{56} erg. The resulting negative feedback delays the growth of the remnant BH by about 70 Myr or evacuates the host galaxy completely.

There is a forbidden region in the parameter space of quasinormal modes of black holes in general relativity. Using both inspiral and ringdown phases of gravitational waves from binary black holes, we propose two methods to test general relativity. We also evaluate how our methods will work when we apply them to Pop III black-hole binaries with typical masses. Adopting the simple mean of the estimated range of the event rate, we have the expected rate of 500 ${\rm yr^{-1}}$. Then, the rates of events with signal-to-noise ratios greater than 20 and greater than 50 are 32 ${\rm yr^{-1}}$ and 2 ${\rm yr^{-1}}$, respectively. Therefore, there is a good chance to confirm (or refute) the Einstein theory in the strong gravity region by observing the expected quasinormal modes.

In this work we are reporting on the measurement of the proton-air inelastic cross section $\sigma^{\rm inel}_{\rm p-air}$ using the Telescope Array (TA) detector. Based on the measurement of the $\sigma^{\rm inel}_{\rm p-air}$ the proton-proton cross section $\sigma_{\rm p-p}$ value is also determined at $\sqrt{s} = 95_{-8}^{+5}$ TeV. Detecting cosmic ray events at ultra high energies with Telescope Array enables us to study this fundamental parameter that we are otherwise unable to access with particle accelerators. The data used in this report is the hybrid events observed by the Middle Drum fluorescence detector together with the surface array detector collected over five years. The value of the $\sigma^{\rm inel}_{\rm p-air}$ is found to be equal to $567.0 \pm 70.5 [{\rm Stat.}] ^{+29}_{-25} [{\rm Sys.}]$ mb. The total proton-proton cross section is subsequently inferred from Glauber Formalism and Block, Halzen and Stanev QCD inspired fit and is found to be equal to $170_{-44}^{+48} [{\rm Stat.}] _{-17}^{+19} [{\rm Sys.}] $mb.

Neutron star binary mergers are strong sources of gravitational waves (GWs). Promising electromagnetic counterparts are short gamma-ray bursts (GRBs) but the emission is highly collimated. We propose that the scattering of the long-lasting plateau emission in short GRBs by the merger ejecta produces nearly isotropic emission for $\sim 10^4$ s with flux $10^{-13}-10^{-10}$ erg cm$^{-2}$ s$^{-1}$ at 100 Mpc in X-ray. This is detectable by Swift XRT and wide field X-ray detectors such as ISS-Lobster, Einstein Probe, eROSITA and WF-MAXI, which are desired by the infrared and optical follow-ups to localize and measure the distance to the host galaxy. The scattered X-rays obtain linear polarization, which correlates with the jet direction, X-ray luminosity and GW polarizations. The activity of plateau emission is also a natural energy source of a macronova (or kilonova) detected in short GRB 130603B without the $r$-process radioactivity.

We examine the bright radio synchrotron counterparts of low-luminosity gamma-ray bursts (llGRBs) and relativistic supernovae (SNe) and find that they can be powered by spherical hypernova (HN) explosions. Our results imply that radio-bright HNe are driven by relativistic jets that are choked deep inside the progenitor stars or quasi-spherical magnetized winds from fast-rotating magnetars. We also consider the optical synchrotron counterparts of radio-bright HNe and show that they can be observed as precursors several days before the SN peak with an r-band absolute magnitude of M_r ~ -14 mag. While previous studies suggested that additional trans-relativistic components are required to power the bright radio emission, we find that they overestimated the energy budget of the trans-relativistic component by overlooking some factors related to the minimum energy of non-thermal electrons. If an additional trans-relativistic component exists, then a much brighter optical precursor with M_r ~ -20 mag can be expected. Thus, the scenarios of radio-bright HNe can be distinguished by using optical precursors, which can be detectable from < 100 Mpc by current SN surveys like the Kiso SN Survey, Palomar Transient Factory, and Panoramic Survey Telescope & Rapid Response System.

The gravitational waveforms in the ghost-free bigravity theory exhibit deviations from those in general relativity. The main difference is caused by graviton oscillations in the bigravity theory. We investigate the prospects for the detection of the corrections to gravitational waveforms from coalescing compact binaries due to graviton oscillations and for constraining bigravity parameters with the gravitational wave observations. We consider the bigravity model discussed by the De Felice-Nakamura-Tanaka subset of the bigravity model, and the phenomenological model in which the bigravity parameters are treated as independent variables. In both models, the bigravity waveform shows strong amplitude modulation, and there can be a characteristic frequency of the largest peak of the amplitude, which depends on the bigravity parameters. We show that there is a detectable region of the bigravity parameters for the advanced ground-based laser interferometers, such as Advanced LIGO, Advanced Virgo, and KAGRA. This region corresponds to the effective graviton mass of $\mu \geq 10^{-17}~{\rm cm}^{-1}$ for $\tilde{c}-1 \geq 10^{-19}$ in the phenomenological model, while $\mu \geq 10^{-16.5}~{\rm cm}^{-1}$ for $\kappa\xi_c^2 \geq 10^{0.5}$ in the De Felice-Nakamura-Tanaka subset of the bigravity model, respectively, where $\tilde{c}$ is the propagation speed of the massive graviton and $\kappa\xi_c^2$ corresponds to the corrections to the gravitational constant in general relativity. These regions are not excluded by existing solar system tests. We also show that, in the case of $1.4-1.4M_{\rm sun}$ binaries at the distance of $200~{\rm Mpc}$, $\log\mu^2$ is determined with an accuracy of ${\cal O}$(0.1)\% at the 1$\sigma$ level for a fiducial model with $\mu^2=10^{-33}~{\rm cm}^{-2}$ in the case of the phenomenological model.

We report on the search for steady point-like sources of neutral particles around 10$^{18}$ eV between 2008 May and 2013 May with the scintillator surface detector of the Telescope Array experiment. We found overall no significant point-like excess above 0.5 EeV in the northern sky. Subsequently, we also searched for coincidence with the Fermi bright Galactic sources. No significant coincidence was found within the statistical uncertainty. Hence, we set an upper limit on the neutron flux that corresponds to an averaged flux of 0.07 km$^{-2}$ yr$^{-1}$ for $E>1$ EeV in the northern sky at the 95% confidence level. This is the most stringent flux upper limit in a northern sky survey assuming point-like sources. The upper limit at the 95% confidence level on the neutron flux from Cygnus X-3 is also set to 0.2 km$^{-2}$ yr$^{-1}$ for $E>0.5$ EeV. This is an order of magnitude lower than previous flux measurements.

The main accretion phase of star formation is investigated in clouds with different metallicities in the range of 0 \le Z \le Z_\odot, resolving the protostellar radius. Starting from a near-equilibrium prestellar cloud, we calculate the cloud evolution up to \sim100 yr after the first protostar formation. The star formation process considerably differs between clouds with lower (Z \le 10^-4 Z_\odot) and higher (Z > 10^-4 Z_\odot) metallicities. Fragmentation frequently occurs and many protostars appear without forming a stable circumstellar disc in lower-metallicity clouds. In these clouds, although protostars mutually interact and some are ejected from the cloud centre, many remain as a small stellar cluster. In contrast, higher-metallicity clouds produce a single protostar surrounded by a nearly stable rotation-supported disc. In these clouds, although fragmentation occasionally occurs in the disc, the fragments migrate inwards and finally fall onto the central protostar. The difference in cloud evolution is due to different thermal evolutions and mass accretion rates. The thermal evolution of the cloud determines the emergence and lifetime of the first core. The first core develops prior to the protostar formation in higher-metallicity clouds, whereas no (obvious) first core appears in lower-metallicity clouds. The first core evolves into a circumstellar disc with a spiral pattern, which effectively transfers the angular momentum outwards and suppresses frequent fragmentation. In lower-metallicity clouds, the higher mass accretion rate increases the disc surface density within a very short time, rendering the disc unstable to self-gravity and inducing vigorous fragmentation.

We have developed a near-infrared camera called ANIR (Atacama Near-InfraRed camera) for the University of Tokyo Atacama Observatory 1.0m telescope (miniTAO) installed at the summit of Cerro Chajnantor (5640 m above sea level) in northern Chile. The camera provides a field of view of 5'.1 $\times$ 5'.1 with a spatial resolution of 0".298 /pixel in the wavelength range of 0.95 to 2.4 $\mu$m. Taking advantage of the dry site, the camera is capable of hydrogen Paschen-$\alpha$ (Pa$\alpha$, $\lambda=$1.8751 $\mu$m in air) narrow-band imaging observations, at which wavelength ground-based observations have been quite difficult due to deep atmospheric absorption mainly from water vapor. We have been successfully obtaining Pa$\alpha$ images of Galactic objects and nearby galaxies since the first-light observation in 2009 with ANIR. The throughputs at the narrow-band filters ($N1875$, $N191$) including the atmospheric absorption show larger dispersion (~10%) than those at broad-band filters (a few %), indicating that they are affected by temporal fluctuations in Precipitable Water Vapor (PWV) above the site. We evaluate the PWV content via the atmospheric transmittance at the narrow-band filters, and derive that the median and the dispersion of the distribution of the PWV are 0.40+/-0.30 mm for $N1875$ and 0.37+/-0.21 mm for $N191$, which are remarkably smaller (49+/-38% for $N1875$ and 59+/-26% for $N191$) than radiometry measurements at the base of Cerro Chajnantor (5100 m alt.). The decrease in PWV can be explained by the altitude of the site when we assume that the vertical distribution of the water vapor is approximated at an exponential profile with scale heights within 0.3-1.9 km (previously observed values at night). We thus conclude that miniTAO/ANIR at the summit of Cerro Chajnantor indeed provides us an excellent capability for a "ground-based" Pa$\alpha$ observation.