We present a search for low energy antineutrino events coincident with the gravitational wave events GW150914 and GW151226, and the candidate event LVT151012 using KamLAND, a kiloton-scale antineutrino detector. We find no inverse beta-decay neutrino events within $\pm 500$ seconds of either gravitational wave signal. This non-detection is used to constrain the electron antineutrino fluence and the luminosity of the astrophysical sources.

In the late stages of nuclear burning for massive stars ($M>8~M_{\sun}$), the production of neutrino-antineutrino pairs through various processes becomes the dominant stellar cooling mechanism. As the star evolves, the energy of these neutrinos increases and in the days preceding the supernova a significant fraction of emitted electron anti-neutrinos exceeds the energy threshold for inverse beta decay on free hydrogen. This is the golden channel for liquid scintillator detectors because the coincidence signature allows for significant reductions in background signals. We find that the kiloton-scale liquid scintillator detector KamLAND can detect these pre-supernova neutrinos from a star with a mass of $25~M_{\sun}$ at a distance less than 690~pc with 3$\sigma$ significance before the supernova. This limit is dependent on the neutrino mass ordering and background levels. KamLAND takes data continuously and can provide a supernova alert to the community.

We report a measurement of the neutrino-electron elastic scattering rate of 862 keV 7Be solar neutrinos based on a 165.4 kton-day exposure of KamLAND. The observed rate is 582 +/- 90 (kton-day)^-1, which corresponds to a 862 keV 7Be solar neutrino flux of (3.26 +/- 0.50) x 10^9 cm^-2s^-1, assuming a pure electron flavor flux. Comparing this flux with the standard solar model prediction and further assuming three flavor mixing, a nu_e survival probability of 0.66 +/- 0.14 is determined from the KamLAND data. Utilizing a global three flavor oscillation analysis, we obtain a total 7Be solar neutrino flux of (5.82 +/- 0.98) x 10^9 cm^-2s^-1, which is consistent with the standard solar model predictions.

We search for electron anti-neutrinos ($\overline{\nu}_e$) from long and short-duration gamma-ray bursts~(GRBs) using data taken by the KamLAND detector from August 2002 to June 2013. No statistically significant excess over the background level is found. We place the tightest upper limits on $\overline{\nu}_e$ fluence from GRBs below 7 MeV and place first constraints on the relation between $\overline{\nu}_e$ luminosity and effective temperature.

We report a measurement of the neutrino-electron elastic scattering rate from 8B solar neutrinos based on a 123 kton-day exposure of KamLAND. The background-subtracted electron recoil rate, above a 5.5 MeV analysis threshold is 1.49+/-0.14(stat)+/-0.17(syst) events per kton-day. Interpreted as due to a pure electron flavor flux with a 8B neutrino spectrum, this corresponds to a spectrum integrated flux of 2.77+/-0.26(stat)+/-0.32(syst) x 10^6 cm^-2s^-1. The analysis threshold is driven by 208Tl present in the liquid scintillator, and the main source of systematic uncertainty is due to background from cosmogenic 11Be. The measured rate is consistent with existing measurements and with Standard Solar Model predictions which include matter enhanced neutrino oscillation.

The vertical absolute fluxes of atmospheric muons and muon charge ratio have been measured precisely at different geomagnetic locations by using the BESS spectrometer. The observations had been performed at sea level (30 m above sea level) in Tsukuba, Japan, and at 360 m above sea level in Lynn Lake, Canada. The vertical cutoff rigidities in Tsukuba (36.2 N, 140.1 E) and in Lynn Lake (56.5 N, 101.0 W) are 11.4 GV and 0.4 GV, respectively. We have obtained vertical fluxes of positive and negative muons in a momentum range from 0.6 to 20 GeV/c with systematic errors less than 3 % in both measurements. By comparing the data collected at two different geomagnetic latitudes, we have seen an effect of cutoff rigidity. The dependence on the atmospheric pressure and temperature, and the solar modulation effect have been also clearly observed. We also clearly observed the decrease of charge ratio of muons at low momentum side with at higher cutoff rigidity region.

The energy spectra of cosmic-ray low-energy antiprotons and protons have been measured by BESS in 1999 and 2000, during a period covering the solar magnetic field reversal. Based on these measurements, a sudden increase of the antiproton to proton flux ratio following the solar magnetic field reversal was observed, and it generally agrees with a drift model of the solar modulation.

We report cosmic-ray proton and helium spectra in energy ranges of 1 to 120 GeV and 1 to 54 GeV/nucleon, respectively, measured by a balloon flight of the BESS spectrometer in 1998. The magnetic-rigidity of the cosmic-rays was reliably determined by highly precise measurement of the circular track in a uniform solenoidal magnetic field of 1 Tesla. Those spectra were determined within overall uncertainties of +-5 % for protons and +- 10 % for helium nuclei including statistical and systematic errors.

The energy spectrum of cosmic-ray antiprotons has been measured in the range 0.18 to 3.56 GeV, based on 458 antiprotons collected by BESS in recent solar-minimum period. We have detected for the first time a distinctive peak at 2 GeV of antiprotons originating from cosmic-ray interactions with the interstellar gas. The peak spectrum is reproduced by theoretical calculations, implying that the propagation models are basically correct and that different cosmic-ray species undergo a universal propagation. Future BESS flights toward the solar maximum will help us to study the solar modulation and the propagation in detail and to search for primary antiproton components.

We calculate the solar-modulated energy spectra of cosmic-ray antiprotons ($\bar{p}$'s) from two candidate primary sources, i.e., evaporating primordial black holes and the annihilation of neutralino dark matter, as well as for the secondary $\bar{p}$'s produced by cosmic-ray interactions with interstellar gas. A large enhancement toward the solar minimum phase emerges in the low-energy flux of $\bar{p}$'s from the primary sources, whereas the flux of the secondary $\bar{p}$'s, falling steeply below 2 GeV, does not significantly vary. This enables us to conduct a very sensitive search for primary $\bar{p}$ components by precisely measuring the $\bar{p}$ spectrum, especially at low energies, throughout the forthcoming solar minimum phase.

We investigate low-energy cosmic-ray antiprotons ($\bar{p}$'s) arising from the fragmentation of quarks and gluons emitted from evaporating primordial black holes (PBHs). To calculate the local interstellar flux of these $\bar{p}$'s, their propagation in the Galaxy is described by a 3-D Monte Carlo simulation based on the diffusion model. This flux is used with recent observations to derive new upper limits on (i) the local PBH explosion rate ${\cal R}<1.7\times10^{-2}$ pc$^{-3}$yr$^{-1}$, (ii) the fraction of the Universe's mass going into PBHs with particular mass, and (iii) the average density of PBHs in the Universe.