We study the equivalence principle, regarded as the cornerstone of general relativity, by analyzing the deformation observable of black hole shadows. Such deformation can arise from new physics and may be expressed as a phenomenological violation of the equivalence principle. Specifically, we assume that there is an additional background vector field that couples to the photons around the supermassive black hole. This type of coupling yields impact on the way the system depends on initial conditions, and affects the black hole shadow at different wavelengths by a different amount, and therefore observations of the shadow in different wavelengths could constrain such couplings. This can be tested by future multi-band observations. Adopting a specific form of the vector field, we obtain constraints on model parameters from Event Horizon Telescope observations and measurements of gas/stellar orbits.

We run three-dimensional numerical simulations for the accretion flow around the white dwarf (WD) in the progenitor system of Tycho's supernova (SN). The mass of the WD, mass of the companion star, and the orbital period are set to be 1$M_{\odot}$, 1.6$M_{\odot}$, and 0.794 day, respectively, based on theoretical and observational researches of Tycho's SN remnant (SNR). We find that when the magnetic field in the accreted material is negligible, outflowing wind is concentrated near the equatorial plane. When the magnetic field has energy equipartition with internal energy, polar wind is comparable with the equatorial wind. A carefully chosen magnetic field between the above two cases ($B=5.44\times10^3 \rm{G}$) can roughly reproduce the latitude-dependent wind required to form the peculiar periphery of Tycho's SNR.

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.

We compile a sample of 93 long gamma-ray bursts (GRBs) from Fermi satellite and 131 from Konus-Wind, which have measured redshifts and well determined spectra, and estimate their pseudo Lorentz factors (${\Gamma_0}$) using the tight ${L_{{\rm{iso}}}}$-${E_{\rm p}}$-${\Gamma_0}$ correlation. The statistical properties and pair correlations of temporal and spectral parameters are studied in the observer frame, rest frame and comoving frame, respectively. We find that the distributions of the duration, peak energy, isotropic energy and luminosity in the different frames are basically lognormal, and their distributions in the comoving frame are narrow, clustering around $T'_{\rm 90}\sim 4000$ s, $E'_{\rm p,c}\sim 0.7$ keV, $E'_{\rm iso,c} \sim 8\times10^{49}$ erg and $L'_{\rm iso,c}\sim 2.5\times10^{46}$ erg s$^{-1}$, where the redshift evolution effect has been taken into account. We also find that the values of ${\Gamma_0}$ are broadly distributed between few tens and several hundreds with median values $\sim 270$. We further analyze the pair correlations of all the quantities, and well confirm ${E_{{\rm{iso}}}}$-${E_{\rm p}}$, ${L_{{\rm{iso}}}}$-${E_{\rm p}}$, ${L_{{\rm{iso}}}}$-${\Gamma_0}$ and $E_{\rm{iso}}$-${\Gamma_0}$ relations, and find that the corresponding relations in the comoving frame do still exist, but have large dispersions. This suggests not only the well-known spectrum-energy relations are intrinsic correlations, but also the observed correlations are governed by the Doppler effect. In addition, the peak energies of long GRBs are independent of durations both in the rest frame and in the comoving frame. And there is a weak anticorrelation between the peak energy and Lorentz factor.

The hadronic interaction of cosmic rays with solar atmosphere can produce high energy gamma rays. The gamma-ray luminosity is correlated both with the flux of primary cosmic rays and the intensity of the solar magnetic field. The gamma rays below 200 GeV have been observed by $Fermi$ without any evident energy cutoff. The bright gamma-ray flux above 100 GeV has been detected only during solar minimum. The only available data in TeV range come from the HAWC observations, however outside the solar minimum. The ARGO-YBJ dataset has been used to search for sub-TeV/TeV gamma rays from the Sun during the solar minimum from 2008 to 2010, the same time period covered by the Fermi data. A suitable model containing the Sun shadow, solar disk emission and inverse-Compton emission has been developed, and the chi-square minimization method was used to quantitatively estimate the disk gamma-ray signal. The result shows that no significant gamma-ray signal is detected and upper limits to the gamma-ray flux at 0.3$-$7 TeV are set at 95\% confidence level. In the low energy range these limits are consistent with the extrapolation of the Fermi-LAT measurements taken during solar minimum and are compatible with a softening of the gamma-ray spectrum below 1 TeV. They provide also an experimental upper bound to any solar disk emission at TeV energies. Models of dark matter annihilation via long-lived mediators predicting gamma-ray fluxes > $10^{-7}$ GeV $cm^{-2}$ $s^{-1}$ below 1 TeV are ruled out by the ARGO-YBJ limits.

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.

This paper reports on the observation of the sidereal large-scale anisotropy of cosmic rays using data collected by the ARGO-YBJ experiment over 5 years (2008$-$2012). This analysis extends previous work limited to the period from 2008 January to 2009 December,near the minimum of solar activity between cycles 23 and 24.With the new data sample the period of solar cycle 24 from near minimum to maximum is investigated. A new method is used to improve the energy reconstruction, allowing us to cover a much wider energy range, from 4 to 520 TeV. Below 100 TeV, the anisotropy is dominated by two wide regions, the so-called "tail-in" and "loss-cone" features. At higher energies, a dramatic change of the morphology is confirmed. The yearly time dependence of the anisotropy is investigated. Finally, no noticeable variation of cosmic-ray anisotropy with solar activity is observed for a median energy of 7 TeV.

A correlation between the secondary cosmic ray flux and the near-earth electric field intensity, measured during thunderstorms, has been found by analyzing the data of the ARGO-YBJ experiment, a full coverage air shower array located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China). The counting rates of showers with different particle multiplicities, have been found to be strongly dependent upon the intensity and polarity of the electric field measured during the course of 15 thunderstorms. In negative electric fields (i.e. accelerating negative charges downwards), the counting rates increase with increasing electric field strength. In positive fields, the rates decrease with field intensity until a certain value of the field EFmin (whose value depends on the event multiplicity), above which the rates begin increasing. By using Monte Carlo simulations, we found that this peculiar behavior can be well described by the presence of an electric field in a layer of thickness of a few hundred meters in the atmosphere above the detector, which accelerates/decelerates the secondary shower particles of opposite charge, modifying the number of particles with energy exceeding the detector threshold. These results, for the first time, give a consistent explanation for the origin of the variation of the electron/positron flux observed for decades by high altitude cosmic ray detectors during thunderstorms.

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.

The combined observation of GW170817 and its electromagnetic counterpart GRB170817A reveals that gravitational waves propagate at the speed of light in high precision. We apply the effective field theory approach to investigate the experimental consequences for the theory of $f(T)$ gravity. We find that the speed of gravitational waves within $f(T)$ gravity is exactly equal to the light speed, and hence the constraints from GW170817 and GRB170817A are trivially satisfied. The results are verified through the standard analysis of cosmological perturbations. Nevertheless, by examining the dispersion relation and the frequency of cosmological gravitational waves, we observe a deviation from the results of General Relativity, quantified by a new parameter. Although its value is relatively small in viable $f(T)$ models, its possible future measurement in advancing gravitational-wave astronomy would be the smoking gun of testing this type of modified gravity.

The hybrid Tibet AS array was successfully constructed in 2014. It has 4500 m$^{2}$ underground water Cherenkov pools used as the muon detector (MD) and 789 scintillator detectors covering 36900 m$^{2}$ as the surface array. At 100 TeV, cosmic-ray background events can be rejected by approximately 99.99\%, according to the full Monte Carlo (MC) simulation for $\gamma$-ray observations. In order to use the muon detector efficiently, we propose to extend the surface array area to 72900 m$^{2}$ by adding 120 scintillator detectors around the current array to increase the effective detection area. A new prototype scintillator detector is developed via optimizing the detector geometry and its optical surface, by selecting the reflective material and adopting dynode readout. This detector can meet our physics requirements with a positional non-uniformity of the output charge within 10\% (with reference to the center of the scintillator), time resolution FWHM of $\sim$2.2 ns, and dynamic range from 1 to 500 minimum ionization particles.

The ARGO-YBJ experiment, a full coverage extensive air shower (EAS) detector located at high altitude (4300 m a.s.l.) in Tibet, China, has smoothly taken data, with very high stability, since November 2007 to the beginning of 2013. The array consisted of a carpet of about 7000 m$^2$ Resistive Plate Chambers (RPCs) operated in streamer mode and equipped with both digital and analog readout, providing the measurement of particle densities up to few particles per cm$^2$. The unique detector features (full coverage, readout granularity, wide dynamic range, etc) and location (very high altitude) allowed a detailed study of the lateral density profile of charged particles at ground very close to the shower axis and its description by a proper lateral distribution function (LDF). In particular, the information collected in the first 10 m from the shower axis have been shown to provide a very effective tool for the determination of the shower development stage ("age") in the energy range 50 TeV - 10 PeV. The sensitivity of the age parameter to the mass composition of primary Cosmic Rays is also discussed.

The ARGO-YBJ detector, located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China), was a full coverage air shower array dedicated to gamma ray astronomy and cosmic ray studies. The wide field of view (~ 2 sr) and high duty cycle (> 86%), made ARGO-YBJ suitable to search for short and unexpected gamma ray emissions like gamma ray bursts (GRBs). Between 2007 November 6 and 2013 February 7, 156 satellite-triggered GRBs (24 of them with known redshift) occurred within the ARGO-YBJ field of view. A search for possible emission associated to these GRBs has been made in the two energy ranges 10-100 GeV and 10-1000 GeV. No significant excess has been found in time coincidence with the satellite detections nor in a time window of one hour after the bursts. Taking into account the EBL absorption, upper limits to the energy fluence at 99% of confidence level have been evaluated,with values ranging from ~ 10-5 erg cm-2 to ~10-1 erg cm-2.

Neutrino-dominated accretion flows (NDAFs) around rotating stellar-mass black holes (BHs) have been theorized as the central engine of relativistic jets launched in massive star core collapse events or compact star mergers. In this work, we calculate the electron neutrino/anti-neutrino spectra of NDAFs by fully taking into account the general relativistic effects, and investigate the effects of viewing angle, BH spin, and mass accretion rate on the results. We show that even though a typical NDAF has a neutrino luminosity lower than that of a typical supernova (SN), it can reach $10^{50}-10^{51}~{\rm erg~s^{-1}}$ peaking at $\sim 10$ MeV, making them potentially detectable with the upcoming sensitive MeV neutrino detectors if they are close enough to Earth. Based on the observed GRB event rate in the local universe and requiring that at least 3 neutrinos are detected to claim a detection, we estimate a detection rate up to $\sim$ (0.10-0.25) per century for GRB-related NDAFs by the Hyper-Kamiokande (Hyper-K) detector if one neglects neutrino oscillation. If one assumes that all Type Ib/c SNe have an engine-driven NDAF, the Hyper-K detection rate would be $\sim$ (1-3) per century. By considering neutrino oscillations, the detection rate may decrease by a factor of 2-3. Detecting one such event would establish the observational evidence of NDAFs in the universe.

We report on a measurement of thermal neutrons, generated by the hadronic component of extensive air showers (EAS), by means of a small array of EN-detectors developed for the PRISMA project (PRImary Spectrum Measurement Array), novel devices based on a compound alloy of ZnS(Ag) and $^{6}$LiF. This array has been operated within the ARGO-YBJ experiment at the high altitude Cosmic Ray Observatory in Yangbajing (Tibet, 4300 m a.s.l.). Due to the tight correlation between the air shower hadrons and thermal neutrons, this technique can be envisaged as a simple way to estimate the number of high energy hadrons in EAS. Coincident events generated by primary cosmic rays of energies greater than 100 TeV have been selected and analyzed. The EN-detectors have been used to record simultaneously thermal neutrons and the air shower electromagnetic component. The density distributions of both components and the total number of thermal neutrons have been measured. The correlation of these data with the measurements carried out by ARGO-YBJ confirms the excellent performance of the EN-detector.

Black holes (BHs) hide themselves behind various astronomical phenomena, and their properties, i.e., mass and spin, are usually difficult to constrain. One leading candidate for the central engine model of gamma-ray bursts (GRBs) invokes a stellar mass BH and a neutrino-dominated accretion flow (NDAF), with the relativistic jet launched due to neutrino-anti-neutrino annihilations. Such a model gives rise to a matter-dominated fireball, and is suitable to interpret GRBs with a dominant thermal component with a photospheric origin. We propose a method to constrain BH mass and spin within the framework of this model, and apply the method to a thermally-dominant GRB 101219B whose initial jet launching radius $r_0$ is constrained from the data. Using our numerical model of NDAF jets, we estimate the following constraints on the central BH: mass $M_{\rm BH} \sim 5-9~M_\odot$, spin parameter $a_* \gtrsim 0.6$, and disk mass $3~M_\odot \lesssim M_{\rm disk} \lesssim 4~M_\odot$. Our results also suggest that the NDAF model is a competitive candidate for the central engine of GRBs with a strong thermal component.

We report on the extensive multi-wavelength observations of the blazar Markarian 421 (Mrk 421) covering radio to gamma-rays, during the 4.5 year period of ARGO-YBJ and Fermi common operation time, from August 2008 to February 2013. In particular, thanks to the ARGO-YBJ and Fermi data, the whole energy range from 100 MeV to 10 TeV is covered without any gap. In the observation period, Mrk 421 showed both low and high activity states at all wavebands. The correlations among flux variations in different wavebands were analyzed. Seven large flares, including five X-ray flares and two GeV gamma-ray flares with variable durations (3-58 days), and one X-ray outburst phase were identified and used to investigate the variation of the spectral energy distribution with respect to a relative quiescent phase. During the outburst phase and the seven flaring episodes, the peak energy in X-rays is observed to increase from sub-keV to few keV. The TeV gamma-ray flux increases up to 0.9-7.2 times the flux of the Crab Nebula. The behavior of GeV gamma-rays is found to vary depending on the flare, a feature that leads us to classify flares into three groups according to the GeV flux variation. Finally, the one-zone synchrotron self-Compton model was adopted to describe the emission spectra. Two out of three groups can be satisfactorily described using injected electrons with a power-law spectral index around 2.2, as expected from relativistic diffuse shock acceleration, whereas the remaining group requires a harder injected spectrum. The underlying physical mechanisms responsible for different groups may be related to the acceleration process or to the environment properties.

The events recorded by ARGO-YBJ in more than five years of data collection have been analyzed to determine the diffuse gamma-ray emission in the Galactic plane at Galactic longitudes 25{\deg} < l < 100{\deg} and Galactic latitudes . The energy range covered by this analysis, from ~350 GeV to ~2 TeV, allows the connection of the region explored by Fermi with the multi-TeV measurements carried out by Milagro. Our analysis has been focused on two selected regions of the Galactic plane, i.e., 40{\deg} < l < 100{\deg} and 65{\deg} < l < 85{\deg} (the Cygnus region), where Milagro observed an excess with respect to the predictions of current models. Great care has been taken in order to mask the most intense gamma-ray sources, including the TeV counterpart of the Cygnus cocoon recently identified by ARGO-YBJ, and to remove residual contributions. The ARGO-YBJ results do not show any excess at sub-TeV energies corresponding to the excess found by Milagro, and are consistent with the predictions of the Fermi model for the diffuse Galactic emission. From the measured energy distribution we derive spectral indices and the differential flux at 1 TeV of the diffuse gamma-ray emission in the sky regions investigated.

By proposing a pure leptonic radiation model, we study the potential gamma-ray emissions from jets of the low-mass X-ray binaries. In this model, the relativistic electrons that are accelerated in the jets are responsible for radiative outputs. Nevertheless, dynamics of jets are dominated by the magnetic and proton-matter kinetic energies. The model involves all kinds of related radiative processes and considers the evolution of relativistic electrons along the jet by numerically solving the kinetic equation. Numerical results show that the spectral energy distributions can extend up to TeV bands, in which synchrotron radiation and synchrotron self-Compton scattering are dominant components. As an example, we apply the model to the low-mass X-ray binary GX 339-4. The results can not only reproduce the currently available observations from GX 339-4, but also predict detectable radiation at GeV and TeV bands by Fermi and CTA telescopes. The future observations with Fermi and CTA can be used to test our model, which could be employed to distinguish the origin of X-ray emissions.

A neutrino-dominated accretion flow (NDAF) around a rotating stellar-mass black hole (BH) is one of the plausible candidates for the central engine of gamma-ray bursts (GRBs). Two mechanisms, i.e., Blandford-Znajek (BZ) mechanism and neutrino annihilation process, are generally considered to power GRBs. Using the analytic solutions from Xue et al. (2013) and ignoring the effects of the magnetic field configuration, we estimate the BZ and neutrino annihilation luminosities as the functions of the disk masses and BH spin parameters to contrast the observational jet luminosities of GRBs. The results show that, although the neutrino annihilation processes could account for most of GRBs, the BZ mechanism is more effective, especially for long-duration GRBs. Actually, if the energy of afterglows and flares of GRBs is included, the distinction between these two mechanisms is more significant. Furthermore, massive disk mass and high BH spin are beneficial to power high luminosities of GRBs. Finally, we discuss possible physical mechanisms to enhance the disk mass or the neutrino emission rate of NDAFs and relevant difference between these two mechanisms.

The ARGO-YBJ experiment has been in stable data taking from November 2007 till February 2013 at the YangBaJing Cosmic Ray Observatory (4300 m a.s.l.). The detector consists of a single layer of Resistive Plate Chambers (RPCs) ( about 6700 m^2}) operated in streamer mode. The signal pick-up is obtained by means of strips facing one side of the gas volume. The digital readout of the signals, while allows a high space-time resolution in the shower front reconstruction, limits the measurable energy to a few hundred TeV. In order to fully investigate the 1-10 PeV region, an analog readout has been implemented by instrumenting each RPC with two large size electrodes facing the other side of the gas volume. Since December 2009 the RPC charge readout has been in operation on the entire central carpet (about 5800 m^2). In this configuration the detector is able to measure the particle density at the core position where it ranges from tens to many thousands of particles per m^2. Thus ARGO-YBJ provides a highly detailed image of the charge component at the core of air showers. In this paper we describe the analog readout of RPCs in ARGO-YBJ and discuss both the performance of the system and the physical impact on the EAS measurements.

The ARGO-YBJ experiment is a full-coverage air shower detector located at the Yangbajing Cosmic Ray Observatory (Tibet, People's Republic of China, 4300 m a.s.l.). The high altitude, combined with the full-coverage technique, allows the detection of extensive air showers in a wide energy range and offer the possibility of measuring the cosmic ray proton plus helium spectrum down to the TeV region, where direct balloon/space-borne measurements are available. The detector has been in stable data taking in its full configuration from November 2007 to February 2013. In this paper the measurement of the cosmic ray proton plus helium energy spectrum is presented in the region 3-300 TeV by analyzing the full collected data sample. The resulting spectral index is $\gamma = -2.64 \pm 0.01$. These results demonstrate the possibility of performing an accurate measurement of the spectrum of light elements with a ground based air shower detector.

The search for Gamma Ray Burst (GRB) emission in the energy range 1-100 GeV in coincidence with the satellite detection has been carried out using the Astrophysical Radiation with Ground-based Observatory at YangBaJing (ARGO-YBJ) experiment. The high altitude location (4300 m a.s.l.), the large active surface ($\sim$ 6700 m$^2$ of Resistive Plate Chambers), the wide field of view ($\sim 2~$sr, limited only by the atmospheric absorption) and the high duty cycle ($>$ 86 %) make the ARGO-YBJ experiment particularly suitable to detect short and unexpected events like GRBs. With the scaler mode technique, i.e., counting all the particles hitting the detector with no measurement of the primary energy and arrival direction, the minimum threshold of $\sim$ 1 GeV can be reached, overlapping the direct measurements carried out by satellites. During the experiment lifetime, from December 17, 2004 to February 7, 2013, a total of 206 GRBs occurring within the ARGO-YBJ field of view (zenith angle $\theta$ $\le$ 45$^{\circ}$) have been analyzed. This is the largest sample of GRBs investigated with a ground-based detector. Two lightcurve models have been assumed and since in both cases no significant excess has been found, the corresponding fluence upper limits in the 1-100 GeV energy region have been derived, with values as low as 10$^{-5}~$erg cm$^{-2}$. The analysis of a subset of 24 GRBs with known redshift has been used to constrain the fluence extrapolation to the GeV region together with possible cutoffs under different assumptions on the spectrum.

The measurement of cosmic ray energy spectra, in particular for individual species, is an essential approach in finding their origin. Locating the "knees" of the spectra is an important part of the approach and has yet to be achieved. Here we report a measurement of the mixed Hydrogen and Helium spectrum using the combination of the ARGO-YBJ experiment and of a prototype Cherenkov telescope for the LHAASO experiment. A knee feature at 640+/-87 TeV, with a clear steepening of the spectrum, is observed. This gives fundamental inputs to galactic cosmic ray acceleration models.

We numerically investigate the thermally unstable accretion disks around black holes. We adopt an evolutionary viscous stress equation to replace the standard alpha-prescription based on the results of two MHD simulations. We find a kind of interesting oscillations on some running models in limit-cycle outburst state. The oscillations arise near the inner boundary and propagate radially outwards. We deem that they are the trapped $p$-mode oscillations excited by sonic-point instability. We directly integrate the local radiation cooling fluxes to construct the mimic bolometric light-curve. We find a series of overtones beside the fundamental harmonic on the power spectra of mimic light-curves. The frequency of the fundamental harmonic is very close to the maximum epicyclic frequency of the disk and the frequency ratio of the fundamental harmonic and overtones is a regular integer series. We suggest that the code for ray-tracing calculation must be time-dependent in virtual observation and point out the robustness of the black hole spin measurement with high frequency QPOs.

The extended TeV gamma-ray source ARGO J2031+4157 (or MGRO J2031+41) is positionally consistent with the Cygnus Cocoon discovered by $Fermi$-LAT at GeV energies in the Cygnus superbubble. Reanalyzing the ARGO-YBJ data collected from November 2007 to January 2013, the angular extension and energy spectrum of ARGO J2031+4157 are evaluated. After subtracting the contribution of the overlapping TeV sources, the ARGO-YBJ excess map is fitted with a two-dimensional Gaussian function in a square region of $10^{\circ}\times 10^{\circ}$, finding a source extension $\sigma_{ext}$= 1$^{\circ}$.8$\pm$0$^{\circ}$.5. The observed differential energy spectrum is $dN/dE =(2.5\pm0.4) \times 10^{-11}(E/1 TeV)^{-2.6\pm0.3}$ photons cm$^{-2}$ s$^{-1}$ TeV$^{-1}$, in the energy range 0.2-10 TeV. The angular extension is consistent with that of the Cygnus Cocoon as measured by $Fermi$-LAT, and the spectrum also shows a good connection with the one measured in the 1-100 GeV energy range. These features suggest to identify ARGO J2031+4157 as the counterpart of the Cygnus Cocoon at TeV energies. The Cygnus Cocoon, located in the star-forming region of Cygnus X, is interpreted as a cocoon of freshly accelerated cosmic rays related to the Cygnus superbubble. The spectral similarity with Supernova Remnants indicates that the particle acceleration inside a superbubble is similar to that in a SNR. The spectral measurements from 1 GeV to 10 TeV allows for the first time to determine the possible spectrum slope of the underlying particle distribution. A hadronic model is adopted to explain the spectral energy distribution.

The energy spectrum of cosmic Hydrogen and Helium nuclei has been measured, below the so-called "knee", by using a hybrid experiment with a wide field-of-view Cherenkov telescope and the Resistive Plate Chamber (RPC) array of the ARGO-YBJ experiment at 4300 m above sea level. The Hydrogen and Helium nuclei have been well separated from other cosmic ray components by using a multi-parameter technique. A highly uniform energy resolution of about 25% is achieved throughout the whole energy range (100 TeV - 700 TeV). The observed energy spectrum is compatible with a single power law with index gamma=-2.63+/-0.06.

The ARGO-YBJ detector is an extensive air shower array that has been used to monitor the northern $\gamma$-ray sky at energies above 0.3 TeV from 2007 November to 2013 January. In this paper, we present the results of a sky survey in the declination band from $-10^{\circ}$ to $70^{\circ}$, using data recorded over the past five years. With an integrated sensitivity ranging from 0.24 to $\sim$1 Crab units depending on the declination, six sources have been detected with a statistical significance greater than 5 standard deviations. Several excesses are also reported as potential $\gamma$-ray emitters. The features of each source are presented and discussed. Additionally, $95\%$ confidence level upper limits of the flux from the investigated sky region are shown. Specific upper limits for 663 GeV $\gamma$-ray AGNs inside the ARGO-YBJ field of view are reported. The effect of the absorption of $\gamma$-rays due to the interaction with extragalactic background light is estimated.

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.

Neutrino-dominated accretion flows (NDAFs) around rotating stellar-mass black holes are plausible candidates for the central engines of gamma-ray bursts (GRBs). We investigate one-dimensional global solutions of NDAFs, taking account of general relativity in Kerr metric, neutrino physics and nucleosynthesis more precisely than previous works. We calculate sixteen solutions with different characterized accretion rates and black hole spins to exhibit the radial distributions of various physical properties in NDAFs. \iffalse We find that the gas pressure and the neutrino cooling always become dominated in the inner region for large accretion rate, $\gtrsim 0.1 M_\odot$ $\rm s^{-1}$. \fi We confirm that the electron degeneracy has important effects in NDAFs and we find that the electron fraction is about 0.46 in the outer region for all the sixteen solutions. From the perspective of the mass fraction, free nucleons, $\rm ^4He$, and $\rm ^{56}Fe$ dominate in the inner, middle, and outer region, respectively. The influence of neutrino trapping on the annihilation is of importance for the superhigh accretion ($\dot M=10M_{\odot} \rm s^{-1}$) and most of the sixteen solutions have an adequate annihilation luminosity for GRBs.

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 the observation of a very high energy \gamma-ray source, whose position is coincident with HESS J1841-055. This source has been observed for 4.5 years by the ARGO-YBJ experiment from November 2007 to July 2012. Its emission is detected with a statistical significance of 5.3 standard deviations. Parameterizing the source shape with a two-dimensional Gaussian function we estimate an extension \sigma=(0.40(+0.32,-0.22}) degree, consistent with the HESS measurement. The observed energy spectrum is dN/dE =(9.0-+1.6) x 10^{-13}(E/5 TeV)^{-2.32-+0.23} photons cm^{-2} s^{-1} TeV^{-1}, in the energy range 0.9-50 TeV. The integral \gamma-ray flux above 1 TeV is 1.3-+0.4 Crab units, which is 3.2-+1.0 times the flux derived by HESS. The differences in the flux determination between HESS and ARGO-YBJ, and possible counterparts at other wavelengths are discussed.

We investigate the vertical structure and elements distribution of neutrino-dominated accretion flows around black holes in spherical coordinates with the reasonable nuclear statistical equilibrium. According our calculations, heavy nuclei tend to be produced in a thin region near the disk surface, whose mass fractions are primarily determined by the accretion rate and the vertical distribution of temperature and density. In this thin region, we find that $^{56}\rm Ni$ is dominant for the flow with low accretion rate (e.g., $0.05 $$M_{\odot}$ $\rm s^{-1}$) but $^{56}\rm Fe$ is dominant for the high counterpart (e.g., $1 M_{\odot}$ $\rm s^{-1}$). The dominant $^{56}\rm Ni$ in the special region may provide a clue to understand the bumps in the optical light curve of core-collapse supernovae.

As one of the brightest active blazars in both X-ray and very high energy $\gamma$-ray bands, Mrk 501 is very useful for physics associated with jets from AGNs. The ARGO-YBJ experiment is monitoring it for $\gamma$-rays above 0.3 TeV since November 2007. Starting from October 2011 the largest flare since 2005 is observed, which lasts to about April 2012. In this paper, a detailed analysis is reported. During the brightest $\gamma$-rays flaring episodes from October 17 to November 22, 2011, an excess of the event rate over 6 $\sigma$ is detected by ARGO-YBJ in the direction of Mrk 501, corresponding to an increase of the $\gamma$-ray flux above 1 TeV by a factor of 6.6$\pm$2.2 from its steady emission. In particular, the $\gamma$-ray flux above 8 TeV is detected with a significance better than 4 $\sigma$. Based on time-dependent synchrotron self-Compton (SSC) processes, the broad-band energy spectrum is interpreted as the emission from an electron energy distribution parameterized with a single power-law function with an exponential cutoff at its high energy end. The average spectral energy distribution for the steady emission is well described by this simple one-zone SSC model. However, the detection of $\gamma$-rays above 8 TeV during the flare challenges this model due to the hardness of the spectra. Correlations between X-rays and $\gamma$-rays are also investigated.

The anisotropy of cosmic rays (CRs) in the solar vicinity is generally at- tributed to the CR streaming due to the discrete distribution of CR sources or local magnetic field modulation. Recently, the two dimensional large scale CR anisotropy has been measured by many experiments in TeV-PeV energy range in both hemispheres. The tail-in excess along the tangential direction of the local spiral arm and the loss cone deficit pointing to the north Galactic pole direction agree with what have been obtained in tens to hundreds of GeV. The persistence of the two large scale anisotropy structures in such a wide range of energy suggests that the anisotropy might be due to a global streaming of the Galactic CRs (GCRs). This work tries to extend the observed CR anisotropy picture from solar system to the whole galaxy. In such a case, we can find a new interesting signature, a loop of GCR streaming, of the GCR propagation. We further calculate the overall GCR streaming induced magnetic field, and find a qualitative consistence with the observed structure of the halo magnetic field.

We report the observation of TeV gamma-rays from the Cygnus region using the ARGO-YBJ data collected from 2007 November to 2011 August. Several TeV sources are located in this region including the two bright extended MGRO J2019+37 and MGRO J2031+41. According to the Milagro data set, at 20 TeV MGRO J2019+37 is the most significant source apart from the Crab Nebula. No signal from MGRO J2019+37 is detected by the ARGO-YBJ experiment, and the derived flux upper limits at 90% confidence level for all the events above 600 GeV with medium energy of 3 TeV are lower than the Milagro flux, implying that the source might be variable and hard to be identified as a pulsar wind nebula. The only statistically significant (6.4 standard deviations) gamma-ray signal is found from MGRO J2031+41, with a flux consistent with the measurement by Milagro.

We revisit the vertical structure of neutrino dominated accretion flows (NDAFs) in spherical coordinates with a new boundary condition based on the mechanical equilibrium. The solutions show that NDAF is significantly thick. The Bernoulli parameter and neutrino trapping are determined by the mass accretion rate and the viscosity parameter. According to the distribution of the Bernoulli parameter, the possible outflow may appear in the outer region of the disk. The neutrino trapping can essentially affect the neutrino radiation luminosity. The vertical structure of NDAF is like a "sandwich", and the multilayer accretion may account for the flares in gamma-ray bursts.

ARGO-YBJ is an air shower detector array with a fully covered layer of resistive plate chambers. It is operated with a high duty cycle and a large field of view. It continuously monitors the northern sky at energies above 0.3 TeV. In this paper, we report a long-term monitoring of Mrk 421 over the period from 2007 November to 2010 February. This source was observed by the satellite-borne experiments Rossi X-ray Timing Explorer and Swift in the X-ray band. Mrk 421 was especially active in the first half of 2008. Many flares are observed in both X-ray and gamma-ray bands simultaneously. The gamma-ray flux observed by ARGO-YBJ has a clear correlation with the X-ray flux. No lag between the X-ray and gamma-ray photons longer than 1 day is found. The evolution of the spectral energy distribution is investigated by measuring spectral indices at four different flux levels. Hardening of the spectra is observed in both X-ray and gamma-ray bands. The gamma-ray flux increases quadratically with the simultaneously measured X-ray flux. All these observational results strongly favor the synchrotron self-Compton process as the underlying radiative mechanism.

For the first time ever, we derive equations governing the time-evolution of fully relativistic slim accretion disks in the Kerr metric, and numerically construct their detailed non-stationary models. We discuss applications of these general results to a possible limit-cycle behavior of thermally unstable disks. Our equations and numerical method are applicable in a wide class of possible viscosity prescriptions, but in this paper we use a diffusive form of the "standard alpha prescription" that assumes the viscous torque is proportional to the total pressure. In this particular case, we find that the parameters which dominate the limit-cycle properties are the mass-supply rate and the value of the alpha-viscosity parameter. Although the duration of the cycle (or the outburst) does not exhibit any clear dependence on the black hole spin, the maximal outburst luminosity (in the Eddington units) is positively correlated with the spin value. We suggest a simple method for a rough estimate of the black hole spin based on the maximal luminosity and the ratio of outburst to cycle durations. We also discuss a temperature-luminosity relation for the Kerr black hole accretion discs limit-cycle. Based on these results we discuss the limit-cycle behavior observed in microquasar GRS 1915+105. We also extend this study to several non-standard viscosity prescriptions, including a "delayed heating" prescription recently stimulated by the recent MHD simulations of accretion disks.

The sun blocks cosmic ray particles from outside the solar system, forming a detectable shadow in the sky map of cosmic rays detected by the ARGO-YBJ experiment in Tibet. Because the cosmic ray particles are positive charged, the magnetic field between the sun and the earth deflects them from straight trajectories and results in a shift of the shadow from the true location of the sun. Here we show that the shift measures the intensity of the field which is transported by the solar wind from the sun to the earth.

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 revisit the vertical structure of black hole accretion disks in spherical coordinates. By comparing the advective cooling with the viscous heating, we show that advection-dominated disks are geometrically thick, i.e., with the half-opening angle > 2\pi/5, rather than slim as supposed previously in the literature.

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 take the vertical distribution of the radial and azimuthal velocity into account in spherical coordinates, and find that the analytic relation c_{s0}/(v_K \Theta) = [(\gamma -1)/(2\gamma)]^{1/2} is valid for both geometrically thin and thick accretion flows, where c_{s0} is the sound speed on the equatorial plane, v_K is the Keplerian velocity, \Theta is the half-opening angle of the flow, and \gamma is the adiabatic index.

We show that there exists a maximal possible accretion rate, beyond which global slim disk solutions cannot be constructed because in the vertical direction the gravitational force would be unable to balance the pressure force to gather the accreted matter. The principle for this restriction is the same as that for the Eddington luminosity and the corresponding critical accretion rate, which were derived for spherical accretion by considering the same force balance in the radial direction. If the assumption of hydrostatic equilibrium is waived and vertical motion is included, this restriction may become even more serious as the value of the maximal possible accretion rate becomes smaller. Previous understanding in the literature that global slim disk solutions could stand for any large accretion rates is due to the overestimation of the vertical gravitational force by using an approximate potential. For accretion flows with large accretion rates at large radii, outflows seem unavoidable in order for the accretion flow to reduce the accretion rate and follow a global solution till the central black hole.

Using the data taken from Tibet II High Density (HD) Array (1997 February-1999 September) and Tibet-III array (1999 November-2005 November), our previous northern sky survey for TeV $\gamma-$ray point sources has now been updated by a factor of 2.8 improved statistics. From $0.0^{\circ}$ to $60.0^{\circ}$ in declination (Dec) range, no new TeV $\gamma-$ray point sources with sufficiently high significance were identified while the well-known Crab Nebula and Mrk421 remain to be the brightest TeV $\gamma-$ray sources within the field of view of the Tibet air shower array. Based on the currently available data and at the 90% confidence level (C.L.), the flux upper limits for different power law index assumption are re-derived, which are approximately improved by 1.7 times as compared with our previous reported limits.

We present a unified description of all the three known classes of optically thick accretion disks around black holes, namely Shakura-Sunyaev disks, slim disks, and neutrino-dominated accretion flows (NDAFs). It is found that NDAFs have both a maximal and a minimal possible mass accretion rate at their each radius. This may be suggestive of an interpretation for the origin of X-ray flares observed in gamma-ray bursts.

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 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.

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.

We present a numerical method for spatially 1.5-dimensional and time-dependent studies of accretion disks around black holes, that is originated from a combination of the standard pseudo-spectral method and the adaptive domain decomposition method existing in the literature, but with a number of improvements in both the numerical and physical senses. In particular, we introduce a new treatment for the connection at the interfaces of decomposed subdomains, construct an adaptive function for the mapping between the Chebyshev-Gauss-Lobatto collocation points and the physical collocation points in each subdomain, and modify the over-simplified 1-dimensional basic equations of accretion flows to account for the effects of viscous stresses in both the azimuthal and radial directions. Our method is verified by reproducing the best results obtained previously by Szuszkiewicz & Miller on the limit-cycle behavior of thermally unstable accretion disks with moderate viscosity. A new finding is that, according to our computations, the Bernoulli function of the matter in such disks is always and everywhere negative, so that outflows are unlikely to originate from these disks. We are encouraged to study the more difficult case of thermally unstable accretion disks with strong viscosity, and wish to report our results in a subsequent paper.

Neutrino-cooled hyperaccretion disks around stellar mass black holes are plausible candidates for the central engine of gamma-ray bursts. We calculate the one-dimensional structure and the annihilation luminosity of such disks. The neutrino optical depth is of crucial importance in determining the neutrino cooling rate and is in turn dependent on the electron fraction, the free nucleon fraction, and the electron degeneracy, with given density and temperature of the disk matter. We construct a bridging formula for the electron fraction that works for various neutrino optical depths, and give exact definitions for the free proton fraction and free neutron fraction. We show that the electron degeneracy has important effects in the sense that it enlarges the absorption optical depth for neutrinos, and it along with the neutronization processes favored by high temperature cause the electron fraction to drop to be below 0.1 in the inner region of the disk. The resulting neutrino annihilation luminosity is considerably reduced comparing with that obtained in previous works where the electron degeneracy was not considered and the electron fraction was simply taken to be 0.5, but it is still likely to be adequate for gamma-ray bursts, and it is ejected mainly from the inner region of the disk and has an anisotropic distribution.

We study the massless scalar wave propagation in the time-dependent Schwarzschild black hole background. We find that the Kruskal coordinate is an appropriate framework to investigate the time-dependent spacetime. A time-dependent scattering potential is derived by considering dynamical black hole with parameters changing with time. It is shown that in the quasinormal ringing both the decay time-scale and oscillation are modified in the time-dependent background.