Almost 30 years have passed since the successful detection of supernova neutrinos from SN 1987A. In the last decades, remarkable progress has been made in neutrino detection technique, through which it may be possible to detect neutrinos from a new source, pre-supernova (pre-SN) neutrinos. They are emitted from a massive star prior to core bounce. Because neutrinos escape from the core freely, they carry information about the stellar physics directly. Pre-SN neutrinos may play an important role in verifying our understanding of stellar evolution for massive stars. Observations of pre-SN neutrinos, moreover, may serve as an alarm regarding a supernova explosion a few days in advance if the progenitor is located in our vicinity, enabling us to observe the next galactic supernova. In this review, we summarize the current status of pre-SN neutrino studies from both of the theoretical and observational points of view.

We examine the sensitivity of neutrino emissions to stellar evolution models for a 15$M_\odot$ progenitor, paying particular attention to a phase prior to the collapse. We demonstrate that the number luminosities in both electron-type neutrinos ($\nu_e$) and their anti-partners ($\bar{\nu}_e$) differ by more than an order of magnitude by changing spatial resolutions and nuclear network sizes on stellar evolution models. We also develop a phenomenological model to capture the essential trend of the diversity, in which neutrino luminosities are expressed as a function of central density, temperature and electron fraction. In the analysis, we show that neutrino luminosity can be well characterized by these central quantities. This analysis also reveals that the most influential quantity to the time evolution of $\nu_e$ luminosity is matter density, while it is temperature for $\bar{\nu}_e$. These qualitative trends will be useful and applicable to constrain the physical state of progenitors at the final stages of stellar evolution from future neutrino observations, although more detailed systematic studies including various mass progenitors are required to assess the applicability.

The small energy exchange via nucleon recoils in neutrino-nucleon scattering is now supposed to be one of the important factors for successful explosion of core-collapse supernovae (CCSNe) as they can change neutrino spectra through accumulation of a large number of scatterings. In finite-difference methods employed for neutrino transport in CCSN simulations, we normally can not afford to deploy a large enough number of energy bins needed to resolve this small energy exchange and sub-grid techniques are employed one way or another. In this paper we study quantitatively with the Monte Carlo (MC) method how well such a treatment performs. We first investigate the effects of nucleon recoils on the neutrino spectra and confirm that the average energy is reduced by $\sim$15% for heavy-lepton neutrinos and by much smaller quantities for other types of neutrinos in a typical post-bounce situation. It is also observed that the nucleon scattering dominates the electron scattering in the thermalization of neutrino spectra in all flavors. We then study possible artifacts that the coarse energy grid may produce in the finite-difference methods. In order to mimic the latter calculation, we re-distribute MC particles in each energy bin after a certain interval in a couple of ways and study how the results are affected and depend on the energy-resolution. We also discuss possible implications of our results for the finite-difference methods.

We present a linear stability analysis of the fast-pairwise neutrino flavor conversion based on a result of our latest axisymmetric core-collapse supernova (CCSN) simulation with full Boltzmann neutrino transport. In the CCSN simulation, coherent asymmetric neutrino emissions of electron-type neutrinos ($\nu_{\rm e}$) and their anti-particles ($\bar{\nu}_{\rm e}$), in which the asymmetry of $\nu_{\rm e}$ and $\bar{\nu}_{\rm e}$ is anti-correlated with each other, occur at almost the same time as the onset of aspherical shock expansion. We find that the asymmetric neutrino emissions play a crucial role on occurrences of fast flavor conversions. The linear analysis shows that unstable modes appear in both pre- and post-shock flows; for the latter they appear only in the hemisphere of higher $\bar{\nu}_{\rm e}$ emissions (the same hemisphere with stronger shock expansion). We analyze in depth the characteristics of electron-lepton-number (ELN) crossing by closely inspecting the angular distributions of neutrinos in momentum space. The ELN crossing happens in various ways, and the property depends on the radius: in the vicinity of neutron star, $\bar{\nu}_{\rm e}$ ($\nu_{\rm e}$) dominates over $\nu_{\rm e}$ ($\bar{\nu}_{\rm e}$) in the forward (backward) direction: at the larger radius the ELN crossing occurs in the opposite way. We also find that the non-radial ELN crossing occurs at the boundary between no ELN crossing and the radial one, which is an effect of genuine multi-D transport. Our findings indicate that the collective neutrino oscillation may occur more commonly in CCSNe and suggest that the CCSN community needs to accommodate these oscillations self-consistently in the modelling of CCSNe.

We make a strong case that the fast neutrino-flavor conversion, one of the collective flavor oscillation modes, commonly occurs in core-collapse supernovae (CCSNe). It is confirmed in the numerical data obtained in realistic simulations of CCSNe but the argument is much more generic and applicable universally: the coherent neutrino-nucleus scattering makes the electron lepton number (ELN) change signs at some inward direction and trigger the flavor conversion in the outward direction in the pre-shock region. Although the ELN crossing is tiny and that is why it has eluded recognition so far, it is still large enough to induce the flavor conversion. Our findings will have an important observational consequences for CCSNe neutrinos.

The General Antiparticle Spectrometer (GAPS) will carry out a sensitive dark matter search by measuring low-energy ($\mathrm{E} < 0.25 \mathrm{GeV/nucleon}$) cosmic ray antinuclei. The primary targets are low-energy antideuterons produced in the annihilation or decay of dark matter. At these energies antideuterons from secondary/tertiary interactions are expected to have very low fluxes, significantly below those predicted by well-motivated, beyond the standard model theories. GAPS will also conduct low-energy antiproton and antihelium searches. Combined, these observations will provide a powerful search for dark matter and provide the best observations to date on primordial black hole evaporation on Galactic length scales. The GAPS instrument detects antinuclei using the novel exotic atom technique. It consists of a central tracker with a surrounding time-of-flight (TOF) system. The tracker is a one cubic meter volume containing 10 cm-diameter lithium-drifted silicon (Si(Li)) detectors. The TOF is a plastic scintillator system that will both trigger the Si(Li) tracker and enable better reconstruction of particle tracks. After coming to rest in the tracker, antinuclei will form an excited exotic atom. This will then de-excite via characteristic X-ray transitions before producing a pion/proton star when the antiparticle annihilates with the nucleus. This unique event topology will give GAPS the nearly background-free detection capability required for a rare-event search. Here we present the scientific motivation for the GAPS experiment, its design and its current status as it prepares for flight in the austral summer of 2021-22.

We report on the highest energy photons from the Crab Nebula observed by the Tibet air shower array with the underground water-Cherenkov-type muon detector array. Based on the criterion of muon number measured in an air shower, we successfully suppress 99.92% of the cosmic-ray background events with energies $E>100$ TeV. As a result, we observed 24 photon-like events with $E>100$ TeV against 5.5 background events, which corresponds to 5.6$\sigma$ statistical significance. This is the first detection of photons with $E>100$ TeV from an astrophysical source.

In this paper, we present the analysis and results of a direct measurement of the cosmic-ray proton spectrum with the CALET instrument onboard the International Space Station, including the detailed assessment of systematic uncertainties. The observation period used in this analysis is from October 13, 2015 to August 31, 2018 (1054 days). We have achieved the very wide energy range necessary to carry out measurements of the spectrum from 50 GeV to 10 TeV covering, for the first time in space, with a single instrument the whole energy interval previously investigated in most cases in separate subranges by magnetic spectrometers (BESS-TeV, PAMELA, and AMS-02) and calorimetric instruments (ATIC, CREAM, and NUCLEON). The observed spectrum is consistent with AMS-02 but extends to nearly an order of magnitude higher energy, showing a very smooth transition of the power-law spectral index from -2.81 +- 0.03 (50--500 GeV) neglecting solar modulation effects (or -2.87 +- 0.06 including solar modulation effects in the lower energy region) to -2.56 +- 0.04 (1--10 TeV), thereby confirming the existence of spectral hardening and providing evidence of a deviation from a single power law by more than 3 sigma.

The CALorimetric Electron Telescope (CALET) is a high-energy astroparticle physics space experiment installed on the International Space Station (ISS), developed and operated by Japan in collaboration with Italy and the United States. The CALET mission goals include the investigation of possible nearby sources of high-energy electrons, of the details of galactic particle acceleration and propagation, and of potential signatures of dark matter. CALET measures the cosmic-ray electron + positron flux up to 20 TeV, gamma-rays up to 10 TeV, and nuclei with Z=1 to 40 up to 1,000 TeV for the more abundant elements during a long-term observation aboard the ISS. Starting science operation in mid-October 2015, CALET performed continuous observation without major interruption with close to 20 million triggered events over 10 GeV per month. Based on the data taken during the first two-years, we present an overview of CALET observations: uses w/o major interruption 1) Electron + positron energy spectrum, 2) Nuclei analysis, 3) Gamma-ray observation including a characterization of on-orbit performance. Results of the electromagnetic counterpart search for LIGO/Virgo gravitational wave events are discussed as well.

We have developed large-area lithium-drifted silicon (Si(Li)) detectors to meet the unique requirements of the General Antiparticle Spectrometer (GAPS) experiment. GAPS is an Antarctic balloon-borne mission scheduled for the first flight in late 2020. The GAPS experiment aims to survey low-energy cosmic-ray antinuclei, particularly antideuterons, which are recognized as essentially background-free signals from dark matter annihilation or decay. The GAPS Si(Li) detector design is a thickness of 2.5 mm, diameter of 10 cm and 8 readout strips. The energy resolution of <4 keV (FWHM) for 20 to 100 keV X-rays at temperature of -35 to -45 C, far above the liquid nitrogen temperatures frequently used to achieve fine energy resolution, is required. We developed a high-quality Si crystal and Li-evaporation, diffusion and drift methods to form a uniform Li-drifted layer. Guard ring structure and optimal etching of the surface are confirmed to suppress the leakage current, which is a main source of noise. We found a thin un-drifted layer retained on the p-side effectively suppresses the leakage current. By these developments, we succeeded in developing the GAPS Si(Li) detector. As the ultimate GAPS instrument will require >1000 10-cm diameter Si(Li) detectors to achieve high sensitivity to rare antideuteron events, high-yield production is also a key factor for the success of the GAPS mission.

The General AntiParticle Spectrometer (GAPS) is a balloon-borne instrument designed to detect cosmic-ray antimatter using the novel exotic atom technique, obviating the strong magnetic fields required by experiments like AMS, PAMELA, or BESS. It will be sensitive to primary antideuterons with kinetic energies of $\approx0.05-0.2$ GeV/nucleon, providing some overlap with the previously mentioned experiments at the highest energies. For $3\times35$ day balloon flights, and standard classes of primary antideuteron propagation models, GAPS will be sensitive to $m_{\mathrm{DM}}\approx10-100$ GeV c$^{-2}$ WIMPs with a dark-matter flux to astrophysical flux ratio approaching 100. This clean primary channel is a key feature of GAPS and is crucial for a rare event search. Additionally, the antiproton spectrum will be extended with high statistics measurements to cover the $0.07 \leq E \leq 0.25 $ GeV domain. For $E>0.2$ GeV GAPS data will be complementary to existing experiments, while $E<0.2$ GeV explores a new regime. The first flight is scheduled for late 2020 in Antarctica. These proceedings will describe the astrophysical processes and backgrounds relevant to the dark matter search, a brief discussion of detector operation, and construction progress made to date.

We present results on searches for gamma-ray counterparts of the LIGO/Virgo gravitational-wave events using CALorimetric Electron Telescope ({\sl CALET}) observations. The main instrument of {\sl CALET}, CALorimeter (CAL), observes gamma-rays from $\sim1$ GeV up to 10 TeV with a field of view of nearly 2 sr. In addition, the {\sl CALET} gamma-ray burst monitor (CGBM) views $\sim$3 sr and $\sim2\pi$ sr of the sky in the 7 keV -- 1 MeV and the 40 keV -- 20 MeV bands, respectively, by using two different crystal scintillators. The {\sl CALET} observations on the International Space Station started in October 2015, and here we report analyses of events associated with the following gravitational wave events: GW151226, GW170104, GW170608, GW170814 and GW170817. Although only upper limits on gamma-ray emission are obtained, they correspond to a luminosity of $10^{49}\sim10^{53}$ erg s$^{-1}$ in the GeV energy band depending on the distance and the assumed time duration of each event, which is approximately the order of luminosity of typical short gamma-ray bursts. This implies there will be a favorable opportunity to detect high-energy gamma-ray emission in further observations if additional gravitational wave events with favorable geometry will occur within our field-of-view. We also show the sensitivity of {\sl CALET} for gamma-ray transient events which is the order of $10^{-7}$~erg\,cm$^{-2}$\,s$^{-1}$ for an observation of 100~s duration.

We analyze the short cosmic ray intensity increase ("cosmic ray burst": CRB) on June 22, 2015 utilizing a global network of muon detectors and derive the global anisotropy of cosmic ray intensity and the density (i.e. the omnidirectional intensity) with 10-minute time resolution. We find that the CRB was caused by a local density maximum and an enhanced anisotropy of cosmic rays both of which appeared in association with Earth's crossing of the heliospheric current sheet (HCS). This enhanced anisotropy was normal to the HCS and consistent with a diamagnetic drift arising from the spatial gradient of cosmic ray density, which indicates that cosmic rays were drifting along the HCS from the north of Earth. We also find a significant anisotropy along the HCS, lasting a few hours after the HCS crossing, indicating that cosmic rays penetrated into the inner heliosphere along the HCS. Based on the latest geomagnetic field model, we quantitatively evaluate the reduction of the geomagnetic cut-off rigidity and the variation of the asymptotic viewing direction of cosmic rays due to a major geomagnetic storm which occurred during the CRB and conclude that the CRB is not caused by the geomagnetic storm, but by a rapid change in the cosmic ray anisotropy and density outside the magnetosphere.

Extended results on the cosmic-ray electron + positron spectrum from 11 GeV to 4.8 TeV are presented based on observations with the Calorimetric Electron Telescope (CALET) on the International Space Station utilizing the data up to November 2017. The analysis uses the full detector acceptance at high energies, approximately doubling the statistics compared to the previous result. CALET is an all-calorimetric instrument with a total thickness of 30 $X_0$ at normal incidence and fine imaging capability, designed to achieve large proton rejection and excellent energy resolution well into the TeV energy region. The observed energy spectrum in the region below 1 TeV shows good agreement with Alpha Magnetic Spectrometer (AMS-02) data. In the energy region below $\sim$300 GeV, CALET's spectral index is found to be consistent with the AMS-02, Fermi Large Area Telescope (Fermi-LAT) and Dark Matter Particle Explorer (DAMPE), while from 300 to 600 GeV the spectrum is significantly softer than the spectra from the latter two experiments. The absolute flux of CALET is consistent with other experiments at around a few tens of GeV. However, it is lower than those of DAMPE and Fermi-LAT with the difference increasing up to several hundred GeV. The observed energy spectrum above $\sim$1 TeV suggests a flux suppression consistent within the errors with the results of DAMPE, while CALET does not observe any significant evidence for a narrow spectral feature in the energy region around 1.4 TeV. Our measured all-electron flux, including statistical errors and a detailed breakdown of the systematic errors, is tabulated in the Supplemental Material in order to allow more refined spectral analyses based on our data.

We examine the possible influence of Earth-directed coronal mass ejections (ECMEs) on the Sun's shadow in the 3~TeV cosmic-ray intensity observed by the Tibet-III air shower (AS) array. We confirm a clear solar-cycle variation of the intensity deficit in the Sun's shadow during ten years between 2000 and 2009. This solar-cycle variation is overall reproduced by our Monte Carlo (MC) simulations of the Sun's shadow based on the potential field model of the solar magnetic field averaged over each solar rotation period. We find, however, that the magnitude of the observed intensity deficit in the Sun's shadow is significantly less than that predicted by MC simulations, particularly during the period around solar maximum when a significant number of ECMEs is recorded. The $\chi^2$ tests of the agreement between the observations and the MC simulations show that the difference is larger during the periods when the ECMEs occur, and the difference is reduced if the periods of ECMEs are excluded from the analysis. This suggests the first experimental evidence of the ECMEs affecting the Sun's shadow observed in the 3~TeV cosmic-ray intensity.

The CALorimetric Electron Telescope (CALET), launched for installation on the International Space Station (ISS) in August, 2015, has been accumulating scientific data since October, 2015. CALET is intended to perform long-duration observations of high-energy cosmic rays onboard the ISS. CALET directly measures the cosmic-ray electron spectrum in the energy range of 1 GeV to 20 TeV with a 2% energy resolution above 30 GeV. In addition, the instrument can measure the spectrum of gamma rays well into the TeV range, and the spectra of protons and nuclei up to a PeV. In order to operate the CALET onboard ISS, JAXA Ground Support Equipment (JAXA-GSE) and the Waseda CALET Operations Center (WCOC) have been established. Scientific operations using CALET are planned at WCOC, taking into account orbital variations of geomagnetic rigidity cutoff. Scheduled command sequences are used to control the CALET observation modes on orbit. Calibration data acquisition by, for example, recording pedestal and penetrating particle events, a low-energy electron trigger mode operating at high geomagnetic latitude, a low-energy gamma-ray trigger mode operating at low geomagnetic latitude, and an ultra heavy trigger mode, are scheduled around the ISS orbit while maintaining maximum exposure to high-energy electrons and other high-energy shower events by always having the high-energy trigger mode active. The WCOC also prepares and distributes CALET flight data to collaborators in Italy and the United States. As of August 31, 2017, the total observation time is 689 days with a live time fraction of the total time of approximately 84%. Nearly 450 million events are collected with a high-energy (E>10 GeV) trigger. By combining all operation modes with the excellent-quality on-orbit data collected thus far, it is expected that a five-year observation period will provide a wealth of new and interesting results.

We analyze the Sun's shadow observed with the Tibet-III air shower array and find that the shadow's center deviates northward (southward) from the optical solar disc center in the "Away" ("Toward") IMF sector. By comparing with numerical simulations based on the solar magnetic field model, we find that the average IMF strength in the "Away" ("Toward") sector is $1.54 \pm 0.21_{\rm stat} \pm 0.20_{\rm syst}$ ($1.62 \pm 0.15_{\rm stat} \pm 0.22_{\rm syst}$) times larger than the model prediction. These demonstrate that the observed Sun's shadow is a useful tool for the quantitative evaluation of the average solar magnetic field.

In August 2015, the CALorimetric Electron Telescope (CALET), designed for long exposure observations of high energy cosmic rays, docked with the International Space Station (ISS) and shortly thereafter began tocollect data. CALET will measure the cosmic ray electron spectrum over the energy range of 1 GeV to 20 TeV with a very high resolution of 2% above 100 GeV, based on a dedicated instrument incorporating an exceptionally thick 30 radiation-length calorimeter with both total absorption and imaging (TASC and IMC) units. Each TASC readout channel must be carefully calibrated over the extremely wide dynamic range of CALET that spans six orders of magnitude in order to obtain a degree of calibration accuracy matching the resolution of energy measurements. These calibrations consist of calculating the conversion factors between ADC units and energy deposits, ensuring linearity over each gain range, and providing a seamless transition between neighboring gain ranges. This paper describes these calibration methods in detail, along with the resulting data and associated accuracies. The results presented in this paper show that a sufficient accuracy was achieved for the calibrations of each channel in order to obtain a suitable resolution over the entire dynamic range of the electron spectrum measurement.

First results of a cosmic-ray electron + positron spectrum, from 10 GeV to 3 TeV, is presented based upon observations with the CALET instrument on the ISS starting in October, 2015. Nearly a half million electron + positron events are included in the analysis. CALET is an all-calorimetric instrument with total vertical thickness of 30 $X_0$ and a fine imaging capability designed to achieve a large proton rejection and excellent energy resolution well into the TeV energy region. The observed energy spectrum over 30 GeV can be fit with a single power law with a spectral index of -3.152 $\pm$ 0.016 (stat.+ syst.). Possible structure observed above 100 GeV requires further investigation with increased statistics and refined data analysis.

The GAPS experiment is designed to carry out a sensitive dark matter search by measuring low-energy cosmic ray antideuterons and antiprotons. GAPS will provide a new avenue to access a wide range of dark matter models and masses that is complementary to direct detection techniques, collider experiments and other indirect detection techniques. Well-motivated theories beyond the Standard Model contain viable dark matter candidates which could lead to a detectable signal of antideuterons resulting from the annihilation or decay of dark matter particles. The dark matter contribution to the antideuteron flux is believed to be especially large at low energies (E < 1 GeV), where the predicted flux from conventional astrophysical sources (i.e. from secondary interactions of cosmic rays) is very low. The GAPS low-energy antiproton search will provide stringent constraints on less than 10 GeV dark matter, will provide the best limits on primordial black hole evaporation on Galactic length scales, and will explore new discovery space in cosmic ray physics. Unlike other antimatter search experiments such as BESS and AMS that use magnetic spectrometers, GAPS detects antideuterons and antiprotons using an exotic atom technique. This technique, and its unique event topology, will give GAPS a nearly background-free detection capability that is critical in a rare-event search. GAPS is designed to carry out its science program using long-duration balloon flights in Antarctica. A prototype instrument was successfully flown from Taiki, Japan in 2012. GAPS has now been approved by NASA to proceed towards the full science instrument, with the possibility of a first long-duration balloon flight in late 2020. Here we motivate low-energy cosmic ray antimatter searches and discuss the current status of the GAPS experiment and the design of the payload.

This paper is a sequel to our previous one (Kato et al.2015), which calculated the luminosities and spectra of electron-type anti-neutrinos ($\bar{\nu}_e$'s) from the progenitors of core-collapse supernovae. Expecting that a capability to detect electron-type neutrinos ($\nu_e$'s) will increase dramatically with the emergence of liquid-argon detectors such as DUNE, we broaden the scope in this study to include all-flavors of neutrinos emitted from the pre-bounce phase. We pick up three progenitor models of an electron capture supernova (ECSN) and iron-core collapse supernovae (FeCCSNe). We find that the number luminosities reach $\sim10^{57} \mathrm{s^{-1}}$ and $\sim10^{53} \mathrm{s^{-1}}$ at maximum for $\nu_e$ and $\bar{\nu}_e$, respectively. We also estimate the numbers of detection events at terrestrial neutrino detectors including DUNE, taking flavor oscillations into account and assuming the distance to the progenitors to be 200 pc. It is demonstrated that $\bar{\nu}_e$'s from the ECSN-progenitor will be undetected at almost all detectors, whereas we will be able to observe $\gtrsim$15900 $\nu_e$'s at DUNE for the inverted mass hierarchy. From the FeCCSN-progenitors, the number of $\bar{\nu}_e$ events will be largest for JUNO, 200-900 $\bar{\nu}_e$'s, depending on the mass hierarchy whereas the number of $\nu_e$ events at DUNE is $\gtrsim$2100 for the inverted mass hierarchy. These results imply that the detection of $\bar{\nu}_e$'s is useful to distinguish FeCCSN- from ECSN-progenitors, while $\nu_e$'s will provide us with detailed information on the collapse phase regardless of the type and mass of progenitor.

The distribution and isotopic composition of volatile elements in planetary materials holds a key to the characterization of the early solar system and the Moon's formation. The Moon and Earth are chemically and isotopically very similar. However, the Moon is highly depleted in volatile elements and the origin of this depletion is still debated. We present gallium isotopic and elemental measurements in a large set of lunar samples to constrain the origin of this volatile depletion. We show that while Ga has a geochemical behavior different from zinc, both elements show a systematic enrichment in the heavier isotopes in lunar mare basalts and Mg-suite rocks compared to the silicate Earth, pointing to a global-scale depletion event. On the other hand, the ferroan anorthosites are isotopically heterogeneous, suggesting a secondary distribution of Ga at the surface of the Moon by volatilization and condensation. The isotopic difference of Ga between Earth and the Moon and the isotopic heterogeneity of the crustal ferroan anorthosites suggest that the volatile depletion occurred following the giant impact and during the lunar magma ocean phase. These results point toward a Moon that has lost its volatile elements during a whole-scale evaporation event and that is now relatively dry compared to Earth.

X-ray emission in young stellar objects (YSOs) is orders of magnitude more intense than in main sequence stars1,2, suggestive of cosmic ray irradiation of surrounding accretion disks. Protoplanetary disk irradiation has been detected around YSOs by HERSCHEL3. In our solar system, short-lived 10Be (half-life = 1.39 My4), which cannot be produced by stellar nucleosynthesis, was discovered in the oldest solar system solids, the calcium-aluminium-rich inclusions (CAIs)5. The high 10Be abundance, as well as detection of other irradiation tracers6,7, suggest 10Be likely originates from cosmic ray irradiation caused by solar flares8. Nevertheless, the nature of these flares (gradual or impulsive), the target (gas or dust), and the duration and location of irradiation remain unknown. Here we use the vanadium isotopic composition, together with initial 10Be abundance to quantify irradiation conditions in the early Solar System9. For the initial 10Be abundances recorded in CAIs, 50V excesses of a few per mil relative to chondrites have been predicted10,11. We report 50V excesses in CAIs up to 4.4 per mil that co-vary with 10Be abundance. Their co-variation dictates that excess 50V and 10Be were synthesised through irradiation of refractory dust. Modelling of the production rate of 50V and 10Be demonstrates that the dust was exposed to solar cosmic rays produced by gradual flares for less than 300 years at about 0.1 au from the protoSun.

We investigate the influences of the nuclear composition on the weak interaction rates of heavy nuclei during the core collapse of massive stars. The nuclear abundances in nuclear statistical equilibrium (NSE) are calculated by some equation of state (EOS) models including in-medium effects on nuclear masses. We systematically examine the sensitivities of electron capture and neutrino-nucleus scattering on heavy nuclei to the nuclear shell effects and the single nucleus approximation. We find that the washout of shell effects at high temperatures brings significant change to weak rates by smoothing the nuclear abundance distribution: the electron capture rate decreases by $\sim$20$\%$ in the early phase and increases by $\sim$40$\%$ in the late phase at most, while the cross section for neutrino-nucleus scattering is reduced by $\sim$15$\%$. This is because the open-shell nuclei become abundant instead of those with closed neutron shells as the shell effects disappear. We also find that the single-nucleus description based on the average values leads to underestimations of weak rates. Electron captures and neutrino coherent scattering on heavy nuclei are reduced by $\sim$80$\%$ in the early phase and by $\sim$5$\%$ in the late phase, respectively. These results indicate that NSE like EOS accounting for shell washout is indispensable for the reliable estimation of weak interaction rates in simulations of core-collapse supernovae.

We present upper limits in the hard X-ray and gamma-ray bands at the time of the LIGO gravitational-wave event GW 151226 derived from the CALorimetric Electron Telescope (CALET) observation. The main instrument of CALET, CALorimeter (CAL), observes gamma-rays from ~1 GeV up to 10 TeV with a field of view of ~2 sr. The CALET gamma-ray burst monitor (CGBM) views ~3 sr and ~2pi sr of the sky in the 7 keV - 1 MeV and the 40 keV - 20 MeV bands, respectively, by using two different scintillator-based instruments. The CGBM covered 32.5% and 49.1% of the GW 151226 sky localization probability in the 7 keV - 1 MeV and 40 keV - 20 MeV bands respectively. We place a 90% upper limit of 2 x 10^{-7} erg cm-2 s-1 in the 1 - 100 GeV band where CAL reaches 15% of the integrated LIGO probability (~1.1 sr). The CGBM 7 sigma upper limits are 1.0 x 10^{-6} erg cm-2 s-1 (7-500 keV) and 1.8 x 10^{-6} erg cm-2 s-1 (50-1000 keV) for one second exposure. Those upper limits correspond to the luminosity of 3-5 x 10^{49} erg s-1 which is significantly lower than typical short GRBs.

We analyze the galactic cosmic ray (GCR) density and its spatial gradient in Forbush Decreases (FDs) observed with the Global Muon Detector Network (GMDN) and neutron monitors (NMs). By superposing the GCR density and density gradient observed in FDs following 45 interplanetary shocks (IP-shocks), each associated with an identified eruption on the sun, we infer the average spatial distribution of GCRs behind IP-shocks. We find two distinct modulations of GCR density in FDs, one in the magnetic sheath and the other in the coronal mass ejection (CME) behind the sheath. The density modulation in the sheath is dominant in the western flank of the shock, while the modulation in the CME ejecta stands out in the eastern flank. This east-west asymmetry is more prominent in GMDN data responding to $\sim$ 60 GV GCRs than in NM data responding to $\sim$ 10 GV GCRs, because of softer rigidity spectrum of the modulation in the CME ejecta than in the sheath. The GSE-y component of the density gradient, $G_y$ shows a negative (positive) enhancement in FDs caused by the eastern (western) eruptions, while $G_z$ shows a negative (positive) enhancement in FDs by the northern (southern) eruptions. This implies the GCR density minimum being located behind the central flank of IP-shock and propagating radially outward from location of the solar eruption. We also confirmed the average $G_z$ changing its sign above and below the heliospheric current sheet, in accord with the prediction of the drift model for the large-scale GCR transport in the heliosphere.

Aiming to distinguish two types of progenitors of core collapse supernovae, i.e., one with a core composed mainly of oxygen and neon (abbreviated as ONe core) and the other with an iron core (or Fe core), we calculated the luminosities and spectra of neutrinos emitted from these cores prior to gravitational collapse, taking neutrino oscillation into account. We found that the total energies emitted as $\bar{\nu}_e$ from the ONe core are $\lesssim 10^{46}\ {\rm erg}$, which is much smaller than $\sim 10^{47}\ {\rm erg}$ for Fe cores. The average energy, on the other hand, is twice as large for the ONe core as those for the Fe cores. The neutrinos produced by the plasmon decays in the ONe core are more numerous than those from the electron-positron annihilation in both cores but they have much lower average energies $\lesssim 1\ {\rm MeV}$. Although it is difficult to detect the pre-supernova neutrinos from the ONe core even if it is located within 200$\ $pc from the earth, we expect $ \sim 9 - 43$ and $\sim 7 - 61$ events for Fe cores at KamLAND and Super-Kamiokande, respectively, depending on the progenitor mass and neutrino-mass hierarchy. These numbers might be increased by an order of magnitude if we envisage next-generation detectors such as JUNO. We will hence be able to distinguish the two types of progenitors by the detection or non-detection of the pre-supernova neutrinos if they are close enough ($\lesssim1\ {\rm kpc}$).

We analyze the three dimensional anisotropy of the galactic cosmic ray (GCR) intensities observed independently with a muon detector at Nagoya in Japan and neutron monitors over four solar activity cycles. We clearly see the phase of the free-space diurnal anisotropy shifting toward earlier hours around solar activity minima in A>0 epochs, due to the reduced anisotropy component parallel to the mean magnetic field. The average parallel component is consistent with a rigidity independent spectrum, while the perpendicular component increases with GCR rigidity. We suggest that this harder spectrum of the perpendicular component is due to contribution from the drift streaming. We find that the bidirectional latitudinal density gradient is positive in A>0 epoch, while it is negative in A<0 epoch, in accord with the drift model prediction. The radial density gradient, on the other hand, varies with ~11-year cycle with maxima (minima) in solar maximum (minimum) periods, but there is no significant difference seen between average radial gradients in A>0 and A<0 epochs. The average parallel mean free path is larger in A<0 than in A>0. We also find, however, that parallel mean free path (radial gradient) appears to persistently increase (decreasing) in the last three cycles of weakening solar activity. We suggest that simple differences between these parameters in A>0 and A<0 epochs are seriously biased by these long-term trends.

We report on a clear solar-cycle variation of the Sun's shadow in the 10 TeV cosmic-ray flux observed by the Tibet air shower array during a full solar cycle from 1996 to 2009. In order to clarify the physical implications of the observed solar cycle variation, we develop numerical simulations of the Sun's shadow, using the Potential Field Source Surface (PFSS) model and the Current Sheet Source Surface (CSSS) model for the coronal magnetic field. We find that the intensity deficit in the simulated Sun's shadow is very sensitive to the coronal magnetic field structure, and the observed variation of the Sun's shadow is better reproduced by the CSSS model. This is the first successful attempt to evaluate the coronal magnetic field models by using the Sun's shadow observed in the TeV cosmic-ray flux.

A new hybrid experiment has been started by AS{\gamma} experiment at Tibet, China, since August 2011, which consists of a low threshold burst-detector-grid (YAC-II, Yangbajing Air shower Core array), the Tibet air-shower array (Tibet-III) and a large underground water Cherenkov muon detector (MD). In this paper, the capability of the measurement of the chemical components (proton, helium and iron) with use of the (Tibet-III+YAC-II) is investigated by means of an extensive Monte Carlo simulation in which the secondary particles are propagated through the (Tibet-III+YAC-II) array and an artificial neural network (ANN) method is applied for the primary mass separation. Our simulation shows that the new installation is powerful to study the chemical compositions, in particular, to obtain the primary energy spectrum of the major component at the knee.

We analyze the precursor of a Forbush Decrease (FD) observed with the Global Muon Detector Network on December 14, 2006. An intense geomagnetic storm is also recorded during this FD with the peak Kp index of 8+. By using the "two-dimensional map" of the cosmic ray intensity produced after removing the contribution from the diurnal anisotropy, we succeed in extracting clear signatures of the precursor. A striking feature of this event is that a weak loss-cone signature is first recorded more than a day prior to the Storm Sudden Commencement (SSC) onset. This suggests that the loss-cone precursor appeared only 7 hours after the Coronal Mass Ejection (CME) eruption from the Sun, when the Interplanetary (IP) shock driven by the Interplanetary Coronal Mass Ejection (ICME) located at 0.4 AU from the Sun. We find the precursor being successively observed with multiple detectors in the network according to the Earth's spin and confirmed that the precursor continuously exists in space. The long lead time (15.6 hours) of this precursor which is almost twice the typical value indicates that the IMF was more quiet in this event than a typical power spectrum assumed for the IMF turbulence. The amplitude (-6.45 %) of the loss-cone anisotropy at the SSC onset is more than twice the FD size, indicating that the maximum intensity depression behind the IP shock is much larger than the FD size recorded at the Earth in this event. We also find the excess intensity from the sunward IMF direction clearly observed during ~10 hours preceding the SSC onset. It is shown that this excess intensity is consistent with the measurement of the particles accelerated by the head-on collisions with the approaching shock. This is the first detailed observation of the precursor due to the shock reflected particles with muon detectors.

We analyze the temporal variation of the diurnal anisotropy of sub-TeV cosmic ray intensity observed with the Matsushiro (Japan) underground muon detector over two full solar activity cycles in 1985-2008. The average sidereal amplitude over the entire period is 0.034+-0.003 %, which is roughly one third of the amplitude reported from AS and deep-underground muon experiments monitoring multi-TeV GCR intensity suggesting a significant attenuation of the anisotropy due to the solar modulation. We find, on the other hand, only weak correlations either with the solar activity- or magnetic-cycles. We examine the temporal variation of the "single-band valley depth" (SBVD) quoted by the Milagro experiment and, by contrast with recent Milagro reports, we find no steady increase in the Matsushiro observations in a 7-year period between 2000 and 2007. We suggest, therefore, that the steady increase of the SBVD reported by the Milagro experiment is not caused by the decreasing solar modulation in the declining phase of the 23rd solar activity cycle.

We analyze the large-scale two-dimensional sidereal anisotropy of multi-TeV cosmic rays by Tibet Air Shower Array, with the data taken from 1999 November to 2008 December. To explore temporal variations of the anisotropy, the data set is divided into nine intervals, each in a time span of about one year. The sidereal anisotropy of magnitude about 0.1% appears fairly stable from year to year over the entire observation period of nine years. This indicates that the anisotropy of TeV Galactic cosmic rays remains insensitive to solar activities since the observation period covers more than a half of the 23rd solar cycle.

We deduce on hourly basis the spatial gradient of the cosmic ray density in three dimensions from the directional anisotropy of high-energy (~50 GeV) galactic cosmic ray (GCR) intensity observed with a global network of muon detectors on the Earth's surface. By analyzing the average features of the gradient in the corotational interaction regions (CIRs) recorded in successive two solar activity minimum periods, we find that the observed latitudinal gradient (Gz) changes its sign from negative to positive on the Earth's heliospheric current sheet (HCS) crossing from the northern to the southern hemisphere in A<0 epoch, while it changes from positive to negative in A>0 epoch. This is in accordance with the drift prediction. We also find a negative enhancement in Gx after the HCS crossing in both A<0 and A>0 epochs, but not in Gy. This asymmetrical feature of Gx and Gy indicates significant contributions from the parallel and perpendicular diffusions to the the gradient in CIRs in addition to the contribution from the drift effect.

Using the Tibet-III air shower array, we search for TeV gamma-rays from 27 potential Galactic sources in the early list of bright sources obtained by the Fermi Large Area Telescope at energies above 100 MeV. Among them, we observe 7 sources instead of the expected 0.61 sources at a significance of 2 sigma or more excess. The chance probability from Poisson statistics would be estimated to be 3.8 x 10^-6. If the excess distribution observed by the Tibet-III array has a density gradient toward the Galactic plane, the expected number of sources may be enhanced in chance association. Then, the chance probability rises slightly, to 1.2 x 10^-5, based on a simple Monte Carlo simulation. These low chance probabilities clearly show that the Fermi bright Galactic sources have statistically significant correlations with TeV gamma-ray excesses. We also find that all 7 sources are associated with pulsars, and 6 of them are coincident with sources detected by the Milagro experiment at a significance of 3 sigma or more at the representative energy of 35 TeV. The significance maps observed by the Tibet-III air shower array around the Fermi sources, which are coincident with the Milagro >=3sigma sources, are consistent with the Milagro observations. This is the first result of the northern sky survey of the Fermi bright Galactic sources in the TeV region.

We develop a model anisotropy best-fitting to the two-dimensional sky-map of multi-TeV galactic cosmic ray (GCR) intensity observed with the Tibet III air shower (AS) array. By incorporating a pair of intensity excesses in the hydrogen deflection plane (HDP) suggested by Gurnett et al., together with the uni-directional and bi-directional flows for reproducing the observed global feature, this model successfully reproduces the observed sky-map including the "skewed" feature of the excess intensity from the heliotail direction, whose physical origin has long remained unknown. These additional excesses are modeled by a pair of the northern and southern Gaussian distributions, each placed ~50 degree away from the heliotail direction. The amplitude of the southern excess is as large as ~0.2 %, more than twice the amplitude of the northern excess. This implies that the Tibet AS experiment discovered for the first time a clear evidence of the significant modulation of GCR intensity in the heliotail and the asymmetric heliosphere.

We analyze the sidereal daily variations observed between 1985 and 2006 at Matsushiro, Japan (MAT) and between 1993 and 2005 at Liapootah, Tasmania (LPT). These stations comprise the two hemisphere network (THN) of underground muon detectors in Japan and Australia. Yearly mean harmonic vectors at MAT and LPT are more or less stable without any significant change in phase and amplitude in correlation with either the solar activity- or magnetic-cycles. In this paper, therefore, we analyze the average anisotropy over the entire observation periods, i.e. 1985-2006 for MAT and 1993-2005 for LPT. We apply to the THN data a best-fitting analysis based on a model anisotropy in space identical to that adopted by Amenomori et al. (2007) for Tibet III data. The median energies of primary cosmic rays recorded are ~0.5 TeV for THN and ~5 TeV for the Tibet III experiment. It is shown that the intensity distribution of the best-fit anisotropy is quite similar to that derived from Tibet III data, regardless of the order of magnitude difference in energies of primary particles. This, together with the THN observations, confirms that the analysis by Amenomori et al. (2007) based on the Tibet III experiment in the northern hemisphere is not seriously biased. The best-fit amplitudes of the anisotropy, on the other hand, are only one third or less of those reported by the Tibet III experiment, indicating attenuation due to solar modulation. The rigidity dependence of the anisotropy amplitude in the sub-TeV region is consistent with the spectrum reported by Hall et al. (1999), smoothly extending to the Tibet III result in the multi-TeV region. The amplitude at higher energies appears almost constant or gradually decreasing with increasing rigidity.

We present for the first time hourly variations of the spatial density gradient of 50 GeV cosmic rays within a sample solar rotation period in 2006. By inversely solving the transport equation, including diffusion, we deduce the gradient from the anisotropy that is derived from the observation made by the Global Muon Detector Network (GMDN). The anisotropy obtained by applying a new analysis method to the GMDN data is precise and free from atmospheric temperature effects on the muon count rate recorded by ground based detectors. We find the derived north-south gradient perpendicular to the ecliptic plane is oriented toward the Helioshperic Current Sheet (HCS) (i.e. southward in the toward sector of the Interplanetary Magnetic Field (IMF) and northward in the away sector). The orientation of the gradient component parallel to the ecliptic plane remains similar in both sectors with an enhancement of its magnitude seen after the Earth crosses the HCS. These temporal features are interpreted in terms of a local maximum of the cosmic ray density at the HCS. This is consistent with the prediction of the drift model for the $A<0$ epoch. By comparing the observed gradient with the numerical prediction of a simple drift model, we conclude that particle drifts in the large-scale magnetic field play an important role in organizing the density gradient, at least in the present $A<0$ epoch. We also found that corotating interaction regions did not have such a notable effect. Observations with the GMDN provide us with a new tool for investigating cosmic ray transport in the IMF.

The amplitude of the Compton-Getting (CG) anisotropy contains the power-law index of the cosmic-ray energy spectrum. Based on this relation and using the Tibet air-shower array data, we measure the cosmic-ray spectral index to be $-3.03 \pm 0.55_{stat} \pm < 0.62_{syst}$ between 6 TeV and 40 TeV, consistent with $-$2.7 from direct energy spectrum measurements. Potentially, this CG anisotropy analysis can be utilized to confirm the astrophysical origin of the ``knee'' against models for non-standard hadronic interactions in the atmosphere.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We are planning to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV (10 - 1000 TeV) energy region by discriminating them from cosmic-ray hadrons. We report on the possibility of detection of gamma rays in the 100 TeV energy region in our field of view, based on the improved sensitivity of our air shower array deduced from the full Monte Carlo simulation.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We have a plan to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV energy region by discriminating them from cosmic-ray hadrons. We have deduced the attainable sensitivity of the muon detector array using our Monte Carlo simulation. We report here on the detailed procedure of our Monte Carlo simulation.