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.

A combined interpretation of the CALET $e^+\, +\, e^-$ spectrum up to 3 TeV and the AMS-02 positron spectrum up to 500 GeV was performed and the results are discussed. To parametrize the background electron flux, we assume a smoothly broken power-law spectrum with an exponential cut-off for electrons and fit this parametrization to the measurements, with either a pulsar or 3-body decay of fermionic Dark Matter as the extra electron-positron pair source responsible for the positron excess. We found that depending on the parameters for the background, both Dark Matter decay and the pulsar model can explain the combined measurements. While the Dark Matter decay scenario is constrained by the Fermi-LAT $\gamma$-ray measurement, we show that 3-body decay of a 800 GeV Dark Matter can be compatible with the $\gamma$-ray flux measurement. We discuss the capability of CALET to discern decaying Dark Matter models from a generic pulsar source scenario, based on simulated data for five years of data-taking.

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.

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.

We present the detection of a blackbody component in GRB 160107A emission by using the combined spectral data of the CALET Gamma-ray Burst Monitor (CGBM) and the MAXI Gas Slit Camera (GSC). The MAXI/GSC detected the emission $\sim$45 s prior to the main burst episode observed by the CGBM. The MAXI/GSC and the CGBM spectrum of this prior emission period is well fit by a blackbody with the temperature of $1.0^{+0.3}_{-0.2}$ keV plus a power-law with the photon index of $-1.6 \pm 0.3$. We discuss the radius to the photospheric emission and the main burst emission based on the observational properties. We stress the importance of the coordinated observations via various instruments collecting the high quality data over a broad energy coverage in order to understand the GRB prompt emission mechanism.

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 ISS-based CALET (CALorimetric Electron Telescope) detector can play an important role in indirect search for Dark Matter (DM), measuring the electron+positron flux in the TeV region for the first time directly. With its fine energy resolution of approximately $2\%$ and good proton rejection ratio ($1:10^5$) it has the potential to search for fine structures in the Cosmic Ray (CR) electron spectrum. In this context we discuss the ability of CALET to discern between signals originating from astrophysical sources and DM decay or annihilation. We fit a parametrization of the local interstellar electron and positron spectra to current measurements, with either a pulsar or 3-body decay of fermionic DM as the extra source causing the positron excess. The expected CALET data for scenarios in which DM decay explains the excess are calculated and analyzed. The signal from this particular 3-body DM decay which can explain the recent measurements of from the AMS$-02$ experiment is shown to be distinguishable from a single pulsar source causing the positron excess by 5 years of observation with CALET, based on the shape of the spectrum. We also study the constraints from extra-galactic diffuse $\gamma$-ray data on this DM-only explanation of the positron excess and show that especially for the possibly remaining parameter space a clearly identifiable signature in the CR electron spectrum exists.

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.

CALET (Calorimetric Electron Telescope), installed on the ISS in August 2015, directly measures the electron+positron cosmic rays flux up to 20 TeV. With its proton rejection capability of 1:10^5 and an aperture of 1200 cm^2 sr, it will provide good statistics even well above one TeV, while also featuring an energy resolution of 2%, which allows it to detect fine structures in the spectrum. Such structures may originate from Dark Matter annihilation or decay, making indirect Dark Matter search one of CALET's main science objectives among others such as identification of signatures from nearby supernova remnants, study of the heavy nuclei spectra and gamma astronomy. The latest results from AMS-02 on positron fraction and total electron+positron flux can be fitted with a parametrization including a single pulsar as an extra power law source with exponential cut-off, which emits an equal amount of electrons and positrons. This single pulsar scenario for the positron excess is extrapolated into the TeV region and the expected CALET data for this case are simulated. Based on this prediction for CALET data, the sensitivity of CALET to Dark Matter annihilation in the galactic halo has been calculated. It is shown that CALET could significantly improve the limits compared to current data, especially for those Dark Matter candidates that feature a large fraction of annihilation directly into electron+positron, such as the LKP (Lightest Kaluza-Klein Particle).

The CALET Gamma-ray Burst Monitor (CGBM) is the secondary scientific instrument of the CALET mission on the International Space Station (ISS), which is scheduled for launch by H-IIB/HTV in 2014. The CGBM provides a broadband energy coverage from 7 keV to 20 MeV, and simultaneous observations with the primary instrument Calorimeter (CAL) in the GeV - TeV gamma-ray range and Advanced Star Camera (ASC) in the optical for gamma-ray bursts (GRBs) and other X-gamma-ray transients. The CGBM consists of two kinds of scintillators: two LaBr$_3$(Ce) (7 keV - 1 MeV) and one BGO (100 keV - 20 MeV) each read by a single photomultiplier. The LaBr$_3$(Ce) crystal, used in space for the first time here for celestial gamma-ray observations, enables GRB observations over a broad energy range from low energy X-ray emissions to gamma rays. The detector performance and structures have been verified using the bread-board model (BBM) via vibration and thermal vacuum tests. The CALET is currently in the development phase of the proto-flight model (PFM) and the pre-flight calibration of the CGBM is planned for August 2013. In this paper, we report on the current status and expected performance of CALET for GRB observations.

The neutron-star Low-Mass X-ray Binary Aquila X-1 was observed by Suzaku for seven times, from 2007 September 28 to October 30. The observations successfully traced an outburst decay in which the source luminosity decreased almost monotonically from $\sim 10^{37}$ erg s$^{-1}$ to $\sim 10^{34}$ erg s$^{-1}$, by $\sim 3$ orders of magnitude. To investigate luminosity-dependent changes in the accretion geometry, five of the seven data sets with a typical exposure of $\sim 18$ ks each were analyzed; the other two were utilized in a previous work \citep{Sakurai2012}. The source was detected up to 100 keV in the 2nd to the 4th observations, to 40 keV in the 5th, and to 10 keV on the last two occasions. All spectra were reproduced successfully by Comptonized blackbody model with relatively high ($\gtrsim 2.0$) optical depths, plus an additional softer optically-thick component. The faintest three spectra were reproduced alternatively by a single Comptonized blackbody model with a relatively low ($\lesssim 0.8$) optical depth. The estimated radius of the blackbody emission, including seed photons for the Comptonization, was $10 \pm 2$ km at a 0.8--100 keV luminosity of $2.4\times 10^{36}$ erg s$^{-1}$ (the 2nd to the 4th observations). In contrast, it decreased to $7 \pm 1$ km and further to $3 \pm 1$ km, at a luminosity of $(4.8-5.2)\times 10^{35}$ erg s$^{-1}$ (the 5th observation) and $\sim 2\times 10^{34}$ erg s$^{-1}$ (the 6th and 7th), respectively, regardless of the above model ambiguity. This can be taken as evidence for the emergence of a weak magnetosphere of the neutron star.

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.

Unified X-ray spectral and timing studies of Cygnus X-1 in the low/hard and hard intermediate state were conducted in a model-independent manner, using broadband Suzaku data acquired on 25 occasions from 2005 to 2009 with a total exposure of ~ 450 ks. The unabsorbed 0.1--500 keV source luminosity changed over 0.8--2.8% of the Eddington limit for 14.8 solar masses. Variations on short (1--2 seconds) and long (days to months) time scales require at least three separate components: a constant component localized below ~2 keV, a broad soft one dominating in the 2--10 keV range, and a hard one mostly seen in 10--300 keV range. In view of the truncated disk/hot inner flow picture, these are respectively interpreted as emission from the truncated cool disk, a soft Compton component, and a hard Compton component. Long-term spectral evolution can be produced by the constant disk increasing in temperature and luminosity as the truncation radius decreases. The soft Compton component likewise increases, but the hard Compton does not, so that the spectrum in the hard intermediate state is dominated by the soft Compton component; on the other hand, the hard Compton component dominates the spectrum in the dim low/hard state, probably associated with a variable soft emission providing seed photons for the Comptonization.

Rapid spectral changes in the hard X-ray on a time scale down to ~0.1 s are studied by applying "shot analysis" technique to the Suzaku observations of the black hole binary Cygnus X-1, performed on 2008 April 18 during the low/hard state. We successfully obtained the shot profiles covering 10--200 keV with the Suzaku HXD-PIN and HXD-GSO detector. It is notable that the 100-200 keV shot profile is acquired for the first time owing to the HXD-GSO detector. The intensity changes in a time-symmetric way, though the hardness does in a time-asymmetric way. When the shot-phase-resolved spectra are quantified with the Compton model, the Compton y-parameter and the electron temperature are found to decrease gradually through the rising phase of the shot, while the optical depth appears to increase. All the parameters return to their time-averaged values immediately within 0.1 s past the shot peak. We have not only confirmed this feature previously found in energies below ~60 keV, but also found that the spectral change is more prominent in energies above ~100 keV, implying the existence of some instant mechanism for direct entropy production. We discuss possible interpretations on the rapid spectral changes in the hard X-ray band.

We present observations of a transient He-like Fe K alpha absorption line in Suzaku observations of the black hole binary Cygnus X-1 on 2011 October 5 near superior conjunction during the high/soft state, which enable us to map the full evolution from the start and the end of the episodic accretion phenomena or dips for the first time. We model the X-ray spectra during the event and trace their evolution. The absorption line is rather weak in the first half of the observation, but instantly deepens for ~10 ks, and weakens thereafter. The overall change in equivalent width is a factor of ~3, peaking at an orbital phase of ~0.08. This is evidence that the companion stellar wind feeding the black hole is clumpy. By analyzing the line with a Voigt profile, it is found to be consistent with a slightly redshifted Fe XXV transition, or possibly a mixture of several species less ionized than Fe XXV. The data may be explained by a clump located at a distance of ~10^(10-12) cm with a density of ~10^((-13)-(-11)) g cm^-3, which accretes onto and/or transits the line-of-sight to the black hole, causing an instant decrease in the observed degree of the ionization and/or an increase in density of the accreting matter. Continued monitoring for individual events with future X-ray calorimeter missions such as ASTRO-H and AXSIO will allow us to map out the accretion environment in detail and how it changes between the various accretion states.

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 report on the geometry of accretion disk and high energy coronae in the strong Comptonization state (the very high/steep power law/hard intermediate state) based on a Suzaku observation of the famous Galactic black hole GX 339-4. These data were taken just before the peak of the 2006-2007 outburst, and the average X-ray luminosity in the 0.7-200 keV band is estimated to be 2.9E38 erg/s for a distance of 8 kpc. We fit the spectrum with both simple (independent disk and corona) and sophisticated (energetically coupled disk and corona) models, but all fits imply that the underlying optically thick disk is truncated significantly before the innermost stable circular orbit around the black hole. We show this directly by a comparison with similarly broadband data from a disk dominated spectrum at almost the same luminosity observed by XMM-Newton and RXTE 3 days after the Suzaku observation. During the Suzaku observation, the QPO frequency changes from 4.3 Hz to 5.5 Hz, while the spectrum softens. The energetically coupled model gives a corresponding 5+/- 8 % decrease in derived inner radius of the disk. While this is not significant, it is consistent with the predicted change in QPO frequency from Lense-Thirring precession of the hot flow interior to the disk and/or a deformation mode of this flow, as a higher QPO frequency implies a smaller size scale for the corona. This is consistent with the truncated disk extending further inwards towards the black hole.

The neutron-star Low-Mass X-ray Binary Aquila X-1 was observed seven times in total with the Suzaku X-ray observatory from September 28 to October 30 in 2007, in the decaying phase of an outburst. In order to constrain the flux-dependent accretion geometry of this source over wider energy bands than employed in most of previous works, the present study utilized two out of the seven data sets. The 0.8-31 keV spectrum on September 28, taken with the XIS and HXD-PIN for an exposure of 13.8 ks, shows an absorbed 0.8-31 keV flux of $3.6\times 10^{-9}$ erg s$^{-1}$ cm$^{-2}$, together with typical characteristics of the soft state of this type of objects. The spectrum was successfully explained by an optically-thick disk emission plus a Comptonized blackbody component. Although these results are in general agreement with previous studies, the significance of a hard tail recently reported using the same data was inconclusive in our analysis. The spectrum acquired on October 9 for an exposure of 19.7 ks was detected over a 0.8-100 keV band with the XIS, HXD-PIN, and HXD-GSO, at an absorbed flux of $8.5\times 10^{-10}$ erg s$^{-1}$ cm$^{-2}$ (in 0.8-100 keV). It shows characteristics of the hard state, and was successfully explained by the same two continuum components but with rather different parameters including much stronger thermal Comptonization, of which the seed photon source was identified with blackbody emission from the neutron-star surface. As a result, the accretion flow in the hard state is inferred to take a form of an optically-thick and geometrically-thin disk down to a radius of $21\pm 4$ km from the neutron star, and then turn into an optically-thin nearly-spherical hot flow.

We present the result of a systematic study of pileup phenomena seen in the X-ray Imaging Spectrometer, an X-ray CCD instrument, onboard the Suzaku observatory. Using a data set of observed sources in a wide range of brightness and spectral hardness, we characterized the pileup fraction, spectral hardening, and grade migration as a function of observed count rate in a frame per pixel. Using the pileup fraction as a measure of the degree of pileup, we found that the relative spectral hardening (the hardness ratio normalized to the intrinsic spectral hardness), branching ratio of split events, and that of detached events increase monotonically as the pileup fraction increases, despite the variety of brightness and hardness of the sample sources. We derived the pileup fraction as a function of radius used for event extraction. Upon practical considerations, we found that events outside of the radius with a pileup fraction of 1% or 3% are useful for spectral analysis. We present relevant figures, tables, and software for the convenience of users who wish to apply our method for their data reduction of piled-up sources.

The Type I Seyfert galaxy Markarian 509 was observed with Suzaku in 2010 November, for a gross time span of 2.2 days. Timing and spectral properties of the 0.5-45 keV X-rays, detected with the XIS and HXD, consistently revealed the presence of a soft spectral component that remained constant while the total X- ray intensity varied by \pm10%. This stable soft component, found in the 0.5-3.0 keV range, was interpreted as a result of thermal Comptonization in a corona with a temperature of ~ 0.5 keV and an optical depth of ~ 18. The time-avearged 0.5-45 keV Suzaku spectrum was reproduced successfully, as a combination of this thermal Comptonization component, a harder power-law of photon index ~ 1.8, moderate reflection, and an iron K-emission line. By analyzing four archival Suzaku datasets of the same object obtained in 2006, the thermal Comptonization component, which was stable during the 2.2 day pointing in 2010, was found to vary on time scales of a few weeks, independently of the power-law component. Implications of these results are discussed in terms of the "multi-zone Comptonization" view, obtained with Suzaku from the black hole binary Cygnus X-1.

Improvements of in-orbit calibration of GSO scintillators in the Hard X-ray Detector on board Suzaku are reported. To resolve an apparent change of the energy scale of GSO which appeared across the launch for unknown reasons, consistent and thorough re-analyses of both pre-launch and in-orbit data have been performed. With laboratory experiments using spare hardware, the pulse height offset, corresponding to zero energy input, was found to change by ~0.5 of the full analog voltage scale, depending on the power supply. Furthermore, by carefully calculating all the light outputs of secondaries from activation lines used in the in-orbit gain determination, their energy deposits in GSO were found to be effectively lower, by several percent, than their nominal energies. Taking both these effects into account, the in-orbit data agrees with the on-ground measurements within ~5%, without employing the artificial correction introduced in the previous work (Kokubun et al. 2007). With this knowledge, we updated the data processing, the response, and the auxiliary files of GSO, and reproduced the HXD-PIN and HXD-GSO spectra of the Crab Nebula over 12-300 keV by a broken powerlaw with a break energy of ~110 keV.

From 2005 to 2009, 25 observations of Cyg X-1 were performed with Suzaku, achieving a total exposure of 446 ks. In all observations, the source was found in the low/hard state, while the 1.5--12.0 keV count rate of the All-Sky Monitor onboard RXTE varied by a factor of ~ 3. In each observation, the 10--60 keV HXD-PIN spectrum and the 60--400 keV HXD-GSO spectrum were fitted successfully by a thermal Comptonization model plus reflection by a thick neutral material. As the soft X-ray intensity increased, the Compton y-parameter was found to decrease from 1.0 to 0.6, while the solid angle of reflection to increase by ~ 30%. Also conducted was timing analysis over a frequency range of 0.001--10 Hz. As the source became brighter in soft X-rays, the characteristic frequency of hard X-ray variation increased from 0.03 to 0.3 Hz, while the fractional hard X-ray variation integrated over 0.001--0.01 Hz decreased by a factor of ~ 5. The signals in the 60--200 keV band were generally found to vary on shorter time scales than those in the 10--60 keV band. These spectral and timing results can be consistently interpreted by presuming that increases in the mass accretion rate cause the Comptonizing hot corona to shrink, while the optically-thick disk to intrude deeper therein.

Based on tests of a tentative detector for observing geo-synchrotron hard X-rays generated by primary electrons, we study the feasibility of probing cosmic electrons above a few TeV to over 10 TeV. Such high energy electrons are expected to give proof of sources near the Earth (e.g. supernova remnants such as Vela: age < ~10^5 years located within <~1kpc). The idea itself is rather old; a high energy electron emits synchrotron X-rays successively in the geomagnetic field and thus gives several X-rays aligned on a meter scale. This feature is a clue to overcome the background problem encountered in other traditional observation methods. We critically examine the feasibility of this approach assuming a satellite altitudes observation, and find that it is difficult to derive a precise energy spectrum of electrons but is possible to get a clear signal of the existence of several TeV electrons if the flux is comparable to the level predicted by a class of plausible models. For such observations, an exposure >1 m^2year would be needed. It would be attractive to incorporate the present scheme in the gamma-ray burst observations.

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.

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.

During thunderstorms on 2008 September 20, a simultaneous detection of gamma rays and electrons was made at a mountain observatory in Japan located 2770 m above sea level. Both emissions, lasting 90 seconds, were associated with thunderclouds rather than lightning. The photon spectrum, extending to 10 MeV, can be interpreted as consisting of bremsstrahlung gamma rays arriving from a source which is 60 - 130 m in distance at 90% confidence level. The observed electrons are likely to be dominated by a primary population escaping from an acceleration region in the clouds.

We argue that we may be able to sort out dark matter models in which electrons are generated through the annihilation and/or decay of dark matter, by using a fact that the initial energy spectrum is reflected in the cosmic-ray electron flux observed at the Earth even after propagation through the galactic magnetic field. To illustrate our idea we focus on three representative initial spectra: (i)monochromatic (ii)flat and (iii)double-peak ones. We find that those three cases result in significantly different energy spectra, which may be probed by the Fermi satellite in operation or an up-coming cosmic-ray detector such as CALET.

We have observed cosmic-ray electrons from 10 GeV to 800 GeV by a long duration balloon flight using Polar Patrol Balloon (PPB) in Antarctica. The observation was carried out for 13 days at an average altitude of 35 km in January 2004. The detector is an imaging calorimeter composed of scintillating-fiber belts and plastic scintillators inserted between lead plates with 9 radiation lengths. The performance of the detector has been confirmed by the CERN-SPS beam test and also investigated by Monte-Carlo simulations. New telemetry system using a commercial satellite of Iridium, power supply by solar batteries, and automatic level control using CPU have successfully been developed and operated during the flight. From the long duration balloon observations, we derived the energy spectrum of cosmic-ray electrons in the energy range from 100 GeV to 800 GeV. In addition, for the first time we derived the electron arrival directions above 100 GeV, which is consistent with the isotropic distribution.

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.