Search

The north-south (NS) anisotropy of galactic cosmic rays (GCRs) is dominated by a diamagnetic drift flow of GCRs in the interplanetary magnetic field (IMF), allowing us to derive key parameters of cosmic-ray propagation, such as the density gradient and diffusion coefficient. We propose a new method to analyze the rigidity spectrum of GCR anisotropy and reveal a solar cycle variation of the NS anisotropy's spectrum using ground-based muon detectors in Nagoya, Japan, and Hobart, Australia. The physics-based correction method for the atmospheric temperature effect on muons is used to combine the different-site detectors free from local atmospheric effects. NS channel pairs in the multi-directional muon detectors are formed to enhance sensitivity to the NS anisotropy, and in this process, general graph matching in graph theory is introduced to survey optimized pairs. Moreover, Bayesian estimation with the Gaussian process allows us to unfold the rigidity spectrum without supposing any analytical function for the spectral shape. Thanks to these novel approaches, it has been discovered that the rigidity spectrum of the NS anisotropy is dynamically varying with solar activity every year. It is attributed to a rigidity-dependent variation of the radial density gradient of GCRs based on the nature of the diamagnetic drift in the IMF. The diffusion coefficient and mean-free-path length of GCRs as functions of the rigidity are also derived from the diffusion-convection flow balance. This analysis expands the estimation limit of the mean-free-path length into $\le200$ GV rigidity region from $<10$ GV region achieved by solar energetic particle observations., Comment: 22 pages, 10 figures, submitted to The Astrophysical Journal

This paper presents the global analysis of two extended decreases of the galactic cosmic ray intensity observed by world-wide networks of ground-based detectors in 2012. This analysis is capable of separately deriving the cosmic ray density (or omnidirectional intensity) and anisotropy each as a function of time and rigidity. A simple diffusion model along the spiral field line between Earth and a cosmic-ray barrier indicates the long duration of these events resulting from about 190$^\circ$ eastern extension of a barrier such as an IP-shock followed by the sheath region and/or the corotating interaction region (CIR). It is suggested that the coronal mass ejection merging and compressing the preexisting CIR at its flank can produce such the extended barrier. The derived rigidity spectra of the density and anisotropy both vary in time during each event period. In particular we find that the temporal feature of the ``phantom Forbush decrease'' reported in an analyzed period is dependent on rigidity, looking quite different at different rigidities. From these rigidity spectra of the density and anisotropy, we derive the rigidity spectrum of the average parallel mean-free-path of pitch angle scattering along the spiral field line and infer the power spectrum of the magnetic fluctuation and its temporal variation. Possible physical cause of the strong rigidity dependence of the ``phantom Forbush decrease'' is also discussed. These results demonstrate the high-energy cosmic rays observed at Earth responding to remote space weather., Comment: Accepted for publication in ApJ

Compilation of papers presented by the GAPS Collaboration at the 38th International Cosmic Ray Conference (ICRC), held July 26 through August 3, 2023 in Nagoya, Japan.

Gamma rays from HESS J1849$-$000, a middle-aged TeV pulsar wind nebula (PWN), are observed by the Tibet air shower array and the muon detector array. The detection significance of gamma rays reaches $4.0\, \sigma$ and $4.4\, \sigma$ levels above 25 TeV and 100 TeV, respectively, in units of Gaussian standard deviation $\sigma$. The energy spectrum measured between $40\, {\rm TeV} < E < 320\, {\rm TeV}$ for the first time is described with a simple power-law function of ${\rm d}N/{\rm d}E = (2.86 \pm 1.44) \times 10^{-16}(E/40\, {\rm TeV})^{-2.24 \pm 0.41}\, {\rm TeV}^{-1}\, {\rm cm}^{-2}\, {\rm s}^{-1}$. The gamma-ray energy spectrum from the sub-TeV ($E < 1\, {\rm TeV}$) to sub-PeV ($100\, {\rm TeV} < E < 1\, {\rm PeV}$) ranges including the results of previous studies can be modeled with the leptonic scenario, inverse Compton scattering by high-energy electrons accelerated by the PWN of PSR J1849$-$0001. On the other hand, the gamma-ray energy spectrum can also be modeled with the hadronic scenario in which gamma rays are generated from the decay of neutral pions produced by collisions between accelerated cosmic-ray protons and the ambient molecular cloud found in the gamma-ray emitting region. The cutoff energy of cosmic-ray protons $E_{\rm p\, cut}$, cut is estimated at ${\rm log}_{10}(E_{\rm p,\, cut}/{\rm TeV}) = 3.73^{+2.98}_{-0.66}$, suggesting that protons are accelerated up to the PeV energy range. Our study thus proposes that HESS J1849$-$000 should be further investigated as a new candidate for a Galactic PeV cosmic-ray accelerator, PeVatron., Comment: 10 pages, 2 figures, Accepted for publication from the Astrophysical Journal

HESS J1843$-$033 is a very-high-energy gamma-ray source whose origin remains unidentified. This work presents, for the first time, the energy spectrum of gamma rays beyond $100\, {\rm TeV}$ from the HESS J1843$-$033 region using the data recorded by the Tibet air shower array and its underground muon detector array. A gamma-ray source with an extension of $0.34^{\circ} \pm 0.12^{\circ}$ is successfully detected above $25\, {\rm TeV}$ at $(\alpha,\, \delta) = (281.09^{\circ}\pm 0.10^{\circ},\, -3.76^{\circ}\pm 0.09^{\circ})$ near HESS J1843$-$033 with a statistical significance of $6.2\, \sigma$, and the source is named TASG J1844$-$038. The position of TASG J1844$-$038 is consistent with those of HESS J1843$-$033, eHWC J1842$-$035, and LHAASO J1843$-$0338. The measured gamma-ray energy spectrum in $25\, {\rm TeV} < E < 130\, {\rm TeV}$ is described with ${\rm d}N/{\rm d}E = (9.70\pm 1.89)\times 10^{-16} (E/40\, {\rm TeV})^{-3.26\pm 0.30}\, {\rm TeV}^{-1} {\rm cm}^{-2} {\rm s}^{-1}$, and the spectral fit to the combined spectra of HESS J1843$-$033, LHAASO J1843$-$0338, and TASG J1844$-$038 implies the existence of a cutoff at $49.5\pm 9.0\, {\rm TeV}$. Associations of TASG J1844-038 with SNR G28.6$-$0.1 and PSR J1844-0346 are also discussed in detail for the first time., Comment: 11 pages, 4 figures, 1 table

We analyze the cosmic-ray variations during a significant Forbush decrease observed with world-wide networks of ground-based neutron monitors and muon detectors during November 3-5, 2021. Utilizing the difference between primary cosmic-ray rigidities monitored by neutron monitors and muon detectors, we deduce the rigidity spectra of the cosmic-ray density (or omnidirectional intensity) and the first- and second-order anisotropies separately, for each hour of data. A clear two-step decrease is seen in the cosmic-ray density with the first $\sim2\%$ decrease after the interplanetary shock arrival followed by the second $\sim5\%$ decrease inside the magnetic flux rope (MFR) at 15 GV. Most strikingly, a large bidirectional streaming along the magnetic field is observed in the MFR with a peak amplitude of $\sim5\%$ at 15 GV which is comparable to the total density decrease inside the MFR. The bidirectional streaming could be explained by adiabatic deceleration and/or focusing in the expanding MFR, which have stronger effects for pitch angles near 90$^\circ$, or by selective entry of GCRs along a leg of the MFR. The peak anisotropy and density depression in the flux rope both decrease with increasing rigidity. The spectra vary dynamically indicating that the temporal variations of density and anisotropy appear different in neutron monitor and muon detector data., Comment: Accepted for publication in the Astrophysical Journal

The lithium-drifted silicon (Si(Li)) detector developed for the General Antiparticle Spectrometer (GAPS) experiment features a thick (~2.2 mm) sensitive layer, large (10 cm) diameter, and excellent energy resolution (~4 keV for 20-100 keV X-rays) at a relatively high operating temperature (approximately -40C). Mass production of GAPS Si(Li) detectors has been performed to construct a large-volume silicon tracker for GAPS. We achieved the first success of the mass production of large-area Si(Li) detectors with a high (~90%) yield rate. Valuable datasets related to detector fabrication, such as detector performance and manufacturing parameters, were recorded and collected during the mass production. This study analyzes the datasets using statistical methods with the aim of comprehensively examining the mass production and to gain valuable insight into the fabrication method. Sufficient uniformities of the performance parameters (leakage current and capacitance) between detectors and strips are found, demonstrating high-quality and stable mass production. We also search for correlations between detector performance and manufacturing parameters by using data-mining techniques. Conventional multivariate analysis (multiple regression analysis) and machine-learning techniques (regression tree analysis) are complementarily used, and it is found that the Li-drift process makes a significant contribution to the performance parameters of the finished detectors. Detailed investigation of the drift process is performed using environmental data, and physical interpretations are presented. Our results provide valuable insight into the fabrication methods for this kind of large-area Si(Li) detector, and encourages future projects that require large-volume silicon trackers., Comment: 41 pages, 18 figures, accepted for publication in Nuclear Instruments and Methods in Physics Research Section A

Cosmic rays (protons and other atomic nuclei) are believed to gain energies of petaelectronvolts (PeV) and beyond at astrophysical particle accelerators called 'PeVatrons' inside our Galaxy. Although a characteristic feature of a PeVatron is expected to be a hard gamma-ray energy spectrum that extends beyond 100 teraelectronvolts (TeV) without a cutoff, none of the currently known sources exhibits such a spectrum due to the low maximum energy of accelerated cosmic rays or insufficient detector sensitivity around 100 TeV. Here we report the observation of gamma-ray emission from the supernova remnant G106.3+2.7 above 10 TeV. This work provides flux data points up to and above 100 TeV and indicates that the very-high-energy gamma-ray emission above 10 TeV is well correlated with a molecular cloud rather than the pulsar PSR J2229+6114. Regarding the gamma-ray emission mechanism of G106.3+2.7, this morphological feature appears to favor a hadronic origin via the {\pi}0 decay caused by accelerated relativistic protons over a leptonic one via the inverse-Compton scattering by relativistic electrons. Furthermore, we point out that an X-ray flux upper limit on the synchrotron spectrum would provide important information to firmly establish the hadronic scenario as the mechanism of particle acceleration at the source., Comment: Published in Nature Astronomy

We report observations of gamma-ray emissions with energies in the 100 TeV energy region from the Cygnus region in our Galaxy. Two sources are significantly detected in the directions of the Cygnus OB1 and OB2 associations. Based on their positional coincidences, we associate one with a pulsar PSR J2032+4127 and the other mainly with a pulsar wind nebula PWN G75.2+0.1 with the pulsar moving away from its original birthplace situated around the centroid of the observed gamma-ray emission. This work would stimulate further studies of particle acceleration mechanisms at these gamma-ray sources., Comment: Accepted for publication in the Physical Review Letters

We report, for the first time, the long-awaited detection of diffuse gamma rays with energies between 100 TeV and 1 PeV in the Galactic disk. Particularly, all gamma rays above 398 TeV are observed apart from known TeV gamma-ray sources and compatible with expectations from the hadronic emission scenario in which gamma rays originate from the decay of $\pi^0$'s produced through the interaction of protons with the interstellar medium in the Galaxy. This is strong evidence that cosmic rays are accelerated beyond PeV energies in our Galaxy and spread over the Galactic disk., Comment: Accepted for publication in the Physical Review Letters

We demonstrate that global observations of high-energy cosmic rays contribute to understanding unique characteristics of a large-scale magnetic flux rope causing a magnetic storm in August 2018. Following a weak interplanetary shock on 25 August 2018, a magnetic flux rope caused an unexpectedly large geomagnetic storm. It is likely that this event became geoeffective because the flux rope was accompanied by a corotating interaction region and compressed by high-speed solar wind following the flux rope. In fact, a Forbush decrease was observed in cosmic-ray data inside the flux rope as expected, and a significant cosmic-ray density increase exceeding the unmodulated level before the shock was also observed near the trailing edge of the flux rope. The cosmic-ray density increase can be interpreted in terms of the adiabatic heating of cosmic rays near the trailing edge of the flux rope, as the corotating interaction region prevents free expansion of the flux rope and results in the compression near the trailing edge. A northeast-directed spatial gradient in the cosmic-ray density was also derived during the cosmic-ray density increase, suggesting that the center of the heating near the trailing edge is located northeast of Earth. This is one of the best examples demonstrating that the observation of high-energy cosmic rays provides us with information that can only be derived from the cosmic ray measurements to observationally constrain the three-dimensional macroscopic picture of the interaction between coronal mass ejections and the ambient solar wind, which is essential for prediction of large magnetic storms., Comment: 19 pages, 3 figures, accepted for publication in the Space Weather

The General Antiparticle Spectrometer (GAPS) is an Antarctic balloon experiment designed for low-energy (0.1$-$0.3 GeV/$n$) cosmic antinuclei as signatures of dark matter annihilation or decay. GAPS is optimized to detect low-energy antideuterons, as well as to provide unprecedented sensitivity to low-energy antiprotons and antihelium nuclei. The novel GAPS antiparticle detection technique, based on the formation, decay, and annihilation of exotic atoms, provides greater identification power for these low-energy antinuclei than previous magnetic spectrometer experiments. This work reports the sensitivity of GAPS to detect antihelium-3 nuclei, based on full instrument simulation, event reconstruction, and realistic atmospheric influence simulations. The report of antihelium nuclei candidate events by AMS-02 has generated considerable interest in antihelium nuclei as probes of dark matter and other beyond the Standard Model theories. GAPS is in a unique position to detect or set upper limits on the cosmic antihelium nuclei flux in an energy range that is essentially free of astrophysical background. In three 35-day long-duration balloon flights, GAPS will be sensitive to an antihelium flux on the level of $1.3^{+4.5}_{-1.2}\cdot 10^{-6}\mathrm{m^{-2}sr^{-1}s^{-1}}(\mathrm{GeV}/n)^{-1}$ (95% confidence level) in the energy range of 0.11$-$0.3 GeV/$n$, opening a new window on rare cosmic physics., Comment: Accepted for publication at Astroparticle Physics, 13 pages, 5 figures

The precise measurement of cosmic-ray antinuclei serves as an important means for identifying the nature of dark matter and other new astrophysical phenomena, and could be used with other cosmic-ray species to understand cosmic-ray production and propagation in the Galaxy. For instance, low-energy antideuterons would provide a "smoking gun" signature of dark matter annihilation or decay, essentially free of astrophysical background. Studies in recent years have emphasized that models for cosmic-ray antideuterons must be considered together with the abundant cosmic antiprotons and any potential observation of antihelium. Therefore, a second dedicated Antideuteron Workshop was organized at UCLA in March 2019, bringing together a community of theorists and experimentalists to review the status of current observations of cosmic-ray antinuclei, the theoretical work towards understanding these signatures, and the potential of upcoming measurements to illuminate ongoing controversies. This review aims to synthesize this recent work and present implications for the upcoming decade of antinuclei observations and searches. This includes discussion of a possible dark matter signature in the AMS-02 antiproton spectrum, the most recent limits from BESS Polar-II on the cosmic antideuteron flux, and reports of candidate antihelium events by AMS-02; recent collider and cosmic-ray measurements relevant for antinuclei production models; the state of cosmic-ray transport models in light of AMS-02 and Voyager data; and the prospects for upcoming experiments, such as GAPS. This provides a roadmap for progress on cosmic antinuclei signatures of dark matter in the coming years., Comment: 45 pages, 14 figures

Cosmic rays are charged particles whose flux observed at Earth shows temporal variations related to space weather phenomena and may be an important tool to study them. The cosmic ray intensity recorded with ground-based detectors also shows temporal variations arising from atmospheric variations. In the case of muon detectors, the main atmospheric effects are related to pressure and temperature changes. In this work, we analyze both effects using data recorded by the Global Muon Detector Network (GMDN), consisting of four multidirectional muon detectors at different locations, in the period between 2007 and 2016. For each GMDN directional channel, we obtain coefficients that describe the pressure and temperature effects. We then analyze how these coefficients can be related to the geomagnetic cutoff rigidity and zenith angle associated with cosmic-ray particles observed by each channel. In the pressure effect analysis, we found that the observed barometric coefficients show a very clear logarithmic correlation with the cutoff rigidity divided by the zenith angle cosine. On the other hand, the temperature coefficients show a good logarithmic correlation with the product of the cutoff and zenith angle cosine after adding a term proportional to the sine of geographical latitude of the observation site. This additional term implies that the temperature effect measured in the northern hemisphere detectors is stronger than that observed in the southern hemisphere. The physical origin of this term and of the good correlations found in this analysis should be studied in detail in future works.

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., Comment: 8 pages, 7 figures, Proc. 36th International Cosmic Ray Conference (ICRC 2019), Madison, USA

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., Comment: April 4, 2019; Submitted to the Physical Review Letters

This study presents a fabrication process for lithium-drifted silicon (Si(Li)) detectors that, compared to previous methods, allows for mass production at a higher yield, while providing a large sensitive area and low leakage currents at relatively high temperatures. This design, developed for the unique requirements of the General Antiparticle Spectrometer (GAPS) experiment, has an overall diameter of 10 cm, with ~9 cm of active area segmented into 8 readout strips, and an overall thickness of 2.5 mm, with $\gtrsim$2.2 mm ($\gtrsim$90%) sensitive thickness. An energy resolution $\lesssim$4 keV full-width at half-maximum (FWHM) for 20-100 keV X-rays is required at the operating temperature ~-40C, which is far above the liquid nitrogen temperatures conventionally used to achieve fine energy resolution. High-yield production is also required for GAPS, which consists of $\gtrsim$1000 detectors. Our specially-developed Si crystal and custom methods of Li evaporation, diffusion and drifting allow for a thick, large-area and uniform sensitive layer. We find that retaining a thin undrifted layer on the p-side of the detector drastically reduces the leakage current, which is a dominant component of the energy resolution at these temperatures. A guard-ring structure and optimal etching of the detector surface are also confirmed to suppress the leakage current. We report on the mass production of these detectors that is ongoing now, and demonstrate it is capable of delivering a high yield of ~90%., Comment: 28 pages, 11 figures, submitted in 6/13/2019; re-submitted in 8/23/2019

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., Comment: submitted to IEEE NSS/MIC 2018 Proceedings

The General Antiparticle Spectrometer (GAPS) is designed to carry out indirect dark matter search by measuring low-energy cosmic-ray antiparticles. Below a few GeVs the flux of antiparticles produced by cosmic-ray collisions with the interstellar medium is expected to be very low and several well-motivated beyond-standard models predict a sizable contribution to the antideuteron flux. GAPS is planned to fly on a long-duration balloon over Antarctica in the austral summer of 2020. The primary detector is a 1m3 central volume containing planes of Si(Li) detectors. This volume is surrounded by a time-of-flight system to both trigger the Si(Li) detector and reconstruct the particle tracks. The detection principle of the experiment relies on the identification of the antiparticle annihilation pattern. Low energy antiparticles slow down in the apparatus and they are captured in the medium to form exotic excited atoms, which de-excite by emitting characteristic X-rays. Afterwards they undergo nuclear annihilation, resulting in a star of pions and protons. The simultaneous measurement of the stopping depth and the dE/dx loss of the primary antiparticle, of the X-ray energies and of the star particle-multiplicity provides very high rejection power, that is critical in rare-event search. GAPS will be able to perform a precise measurement of the cosmic antiproton flux below 250 MeV, as well as a sensitive search for antideuterons., Comment: 7 pages, 3 figures, Low Energy Antiproton Physics Conference (LEAP) 2018, Paris (France)

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., Comment: Talk presented CIPANP2018. 13 pages, LaTeX, 7 figures. Higher resolution figs at http://gaps1.astro.ucla.edu/gaps/papers2018/cipanp18

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., Comment: accepted for the publication in the Astrophysical Journal

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., Comment: 11 pages, 3 figures, accepted for publication in the ApJ

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., Comment: 7 pages, 4 figures

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., Comment: 8 pags, 3 figures, Proc. 35th International Cosmic Ray Conference (ICRC 2017), Busan, Korea

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., Comment: 54 pages, 11 figures, 1 table, accepted for publication in the Astrophysical Journal

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., Comment: accepted for the publication in the Astrophysical Journal