We present 5.5 GHz observations with the VLA of a sample of nearby galaxies with energetic nuclear outbursts at mid-infrared (MIR) bands. These observations reach a uniform depth down to a median rms of ~10 uJy, representing one of most sensitive searches for radio emission associated with nuclear transients. We detect radio emission in 12 out of 16 galaxies at a level of >5sigma, corresponding to a detection rate of 75%. Such a high detection is remarkably different from previous similar searches in stellar tidal disruption events. The radio emission is compact and not resolved for the majority of sources on scales of ~<0.5" (<0.9 kpc at z<0.1). We find the possibility of the star-formation contributing to the radio emission is low, but an AGN origin remains a plausible scenario, especially for sources that show evidence of AGN activity in their optical spectra. If the detections could represent radio emission associated with nuclear transient phenomenon such as jet or outflow, we use the blast wave model by analogy with the GRB afterglows to describe the evolution of radio light curves. In this context, the observations are consistent with a decelerating jet with an energy of ~10^{51-52} erg viewed at 30\degree-60\degree off-axis at later times, suggesting that powerful jets may be ubiquitous among MIR-burst galaxies. Future continuous monitoring observations will be crucial to decipher the origin of radio emission through detections of potential flux and spectral evolution. Our results highlight the importance of radio observations to constrain the nature of nuclear MIR outbursts in galaxies.

The acceleration of electrons in 3C 279 is investigated through analyzing the injected electron energy distribution (EED) in a time-dependent synchrotron self-Compton + external Compton emission model. In this model, it is assumed that relativistic electrons are continuously injected into the emission region, and the injected EED [$Q_e^\prime(\gamma^\prime)$] follows a single power-law form with low- and high-energy cutoffs $\rm \gamma_{min}'$ and $\rm \gamma_{max}'$, respectively, and the spectral index $n$, i.e, $Q_e^\prime(\gamma^\prime)\propto\gamma^{\prime-n}$. This model is applied to 14 quasi-simultaneous spectral energy distributions (SEDs) of 3C 279. The Markov Chain Monte Carlo fitting technique is performed to obtain the best-fitting parameters and the uncertainties on the parameters. The results show that the injected EED is well constrained in each state. The value of $n$ is in the range of 2.5 to 3.8, which is larger than that expected by the classic non-relativistic shock acceleration. However, the large value of $n$ can be explained by the relativistic oblique shock acceleration. The flaring activity seems to be related to an increased acceleration efficiency, reflected in an increased $\gamma'_{\rm min}$ and electron injection power.

The statistical properties of X-ray flares from two separate locations (nucleus and HST-1) in the M87 jet are investigated to reveal the physical origin of the flares. We analyse the archival \textit{Chandra} data for M87, and identify 14 flares in the nucleus and 9 flares in HST-1. The peak intensity ($I_{\rm{P}}$) and the flaring duration time ($T_{\rm{fl}}$) for each flare are obtained. It is found that the distributions of both $I_{\rm{P}}$ and $T_{\rm{fl}}$ for the nucleus obey a power-law form with a similar index. A similar result is also obtained for HST-1, and no significant inconsistency between the nucleus and HST-1 is found for the indices. Similar to solar X-ray flares, the power-law distributions of the flare event parameters can be well explained by a self-organized criticality (SOC) system, which are triggered by magnetic reconnection. Our results suggest that the flares from nucleus and HST-1 are possibly triggered by magnetic reconnection process. The consistent indices for the distributions of $I_{\rm{P}}$ and $T_{\rm{fl}}$ in the CORE and HST-1 indicate that the dimensions of the energy dissipation of the magnetic reconnection is identical in the two regions. A strong correlation between the flares in the two regions also suggests a similar physical origin for the flares.

Generating mocks for future sky surveys requires large volumes and high resolutions, which is computationally expensive even for fast simulations. In this work we try to develop numerical schemes to calibrate various halo and matter statistics in fast low resolution simulations compared to high resolution N-body and hydrodynamic simulations. For the halos, we improve the initial condition accuracy and develop a halo finder "relaxed-FOF", where we allow different linking length for different halo mass and velocity dispersions. We show that our relaxed-FoF halo finder improves the common statistics, such as halo bias, halo mass function, halo auto power spectrum in real space and in redshift space, cross correlation coefficient with the reference halo catalog, and halo-matter cross power spectrum. We also incorporate the potential gradient descent (PGD) method into fast simulations to improve the matter distribution at nonlinear scale. By building a lightcone output, we show that the PGD method significantly improves the weak lensing convergence tomographic power spectrum. With these improvements FastPM is comparable to the high resolution full N-body simulation of the same mass resolution, with two orders of magnitude fewer time steps. These techniques can be used to improve the halo and matter statistics of FastPM simulations for mock catalogs of future surveys such as DESI and LSST.

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.

The energy dissipation mechanism in blazar jet is unknown. Blazar's flares could provide insights on this problem. Here we report statistical results of XMM-Newton X-ray flares of Mrk 421. We analyze all public XMM-Newton X-ray observations for Mrk 421, and construct the light curves. Through fitting light curves, we obtain the parameters of flare-profiles, such as peak flux ($F_{\rm p}$) and flaring time duration ($T_{\rm fl}$). It is found that both the distributions of $F_{\rm p}$ and $T_{\rm fl}$ obey a power-law form, with the same index of $\alpha_{\rm F}=\alpha_{\rm T}\approx1$. The statistical properties are consistent with the predictions by a self-organized criticality (SOC) system with energy dissipation in one-dimensional space. This is similar to solar flare, but with different space dimensions of the energy dissipation domain. This suggests that X-ray flaers of Mrk 421 are possibly driven by a magnetic reconnection mechanism. Moreover, in the analysis, we find that variability on timescale of $\sim1000\ $s frequently appears. Such rapid variability indicates a magnetic field of $\geq 2.1\delta_{\rm D}^{-1/3}$ G ($\delta_{\rm D}$ is the Doppler factor) in emission region.

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.

Fast N-body PM simulations with a small number of time steps such as FastPM or COLA have been remarkably successful in modeling the galaxy statistics, but their lack of small scale force resolution and long time steps cannot give accurate halo matter profiles or matter power spectrum. High resolution N-body simulations can improve on this, but lack baryonic effects, which can only be properly included in hydro simulations. Here we present a scheme to calibrate the fast simulations to mimic the precision of the hydrodynamic simulations or high resolution N-body simulations. The scheme is based on a gradient descent of either effective gravitational potential, which mimics the short range force, or of effective enthalpy, which mimics gas hydrodynamics and feedback. The scheme is fast and differentiable, and can be incorporated as a post-processing step into any simulation. It gives very good results for the matter power spectrum for several of the baryonic feedback and dark matter simulations, and also gives improved dark matter halo profiles. The scheme is even able to find the large subhalos, and increase the correlation coefficient between the fast simulations and the high resolution N-body or hydro simulations. It can also be used to add baryonic effects to the high resolution N-body simulations. While the method has free parameters that can be calibrated on various simulations, they can also be viewed as astrophysical nuisance parameters describing baryonic effects that can be marginalized over during the data analysis. In this view these parameters can be viewed as an efficient parametrization of baryonic effects.

Recent analyses of the Pal 5 and GD-1 tidal streams suggest that the inner dark matter halo of the Milky Way is close to spherical, in tension with predictions from collisionless N-body simulations of cosmological structure formation. We use the Eris simulation to test whether the combination of dissipative physics and hierarchical structure formation can produce Milky Way-like galaxies whose dark matter halos match the tidal stream constraints from the GD-1 and Pal 5 clusters. We use a dynamical model of the simulated Eris galaxy to generate many realizations of the GD-1 and Pal 5 tidal streams, marginalize over observational uncertainties in the cluster galactocentric positions and velocities, and compare with the observational constraints. We find that the total density and potential of Eris contributed by baryons and dark matter satisfies constraints from the existing Milky Way stellar stream data, as the baryons both round and redistribute the dark matter during the dissipative formation of the galaxy, and provide a centrally-concentrated mass distribution that rounds the inner potential. The Eris dark matter halo or a spherical Navarro-Frenk-White dark matter work comparably well in modeling the stream data. In contrast, the equivalent dark matter-only ErisDark simulation produces a prolate halo that cannot reproduce the observed stream data. The on-going Gaia mission will provide decisive tests of the consistency between $\Lambda$CDM and Milky Way streams, and should distinguish between models like Eris and more spherical halos.

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.

In the work, we explore the role of jet properties in the spectral evolution for a sample of Fermi-LAT bright blazars composed primarily by flat spectrum radio quasars (FSRQs). We introduce a near-equipartition log-parabola (NELP) model to fit the quasi-simultaneous multi-waveband spectral energy distributions (SEDs). The Markov Chain Monte Carlo (MCMC) method is employed to determine the best values of spectral parameters and its errors. The correlations of synchrotron peak luminosity $L_{pk}^*$ and its peak frequency $\nu_{pk}^*$, SSC-dominant factor $\xi_s$ and $\nu_{pk}^*$, and their implications as to the spectral evolution were studied. The statistical analysis and the comparison with theoretical approximation indicate that: (1)the spectral evolution in $L_{pk}^*-\nu_{pk}^*$ plane may not only ascribed to variation of the characteristic energy $\gamma_{br}^\prime$ of the employed electron energy distribution (EED), but that other quantities must be varying as well.(2) for the spectral evolution in $\xi_s-\nu_{pk}^*$ plane the magnetic field $B^\prime$ may play the dominant role, and $\gamma_{br}^\prime$ may be considered to be the subdominant role. (3) Variation of the beaming factor $\delta_b$ may be mainly responsible for the absence of a possible correlation, since $B^\prime$ is strongly anti-correlated with the size of radiation zone $R_b^\prime$. The relation is agreement with that derived from the analytical theory, and should be self-consistent with the underling acceleration mechanism of the employed EED. By assuming a conical jet with opening angle $\theta_{op}\sim1/\Gamma_b$, we find that the $\gamma$-ray emission site being located at a distance $\sim10^{17}-10^{19}$ cm.

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.

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.

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.

With the standard deviation for the logarithm of the re-scaled range $\langle |F(t+\tau)-F(t)|\rangle$ of simulated fractal Brownian motions $F(t)$ given in a previous paper \cite{q14}, the method of least squares is adopted to determine the slope, $S$, and intercept, $I$, of the log$(\langle |F(t+\tau)-F(t)|\rangle)$ vs $\rm{log}(\tau)$ plot to investigate the limitation of this procedure. It is found that the reduced $\chi^2$ of the fitting decreases with the increase of the Hurst index, $H$ (the expectation value of $S$), which may be attributed to the correlation among the re-scaled ranges. Similarly, it is found that the errors of the fitting parameters $S$ and $I$ are usually smaller than their corresponding standard deviations. These results show the limitation of using the simple least square method to determine the dimension of a fractal time series. Nevertheless, they may be used to reinterpret the fitting results of the least square method to determine the dimension of fractal Brownian motions more self-consistently. The currency exchange rate between Euro and Dollar is used as an example to demonstrate this procedure and a fractal dimension of 1.511 is obtained for spans greater than 30 transactions.

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 long-term optical and infrared color variability of blazars has been investigated with the SMARTS monitoring data. The sample in this study consists of 49 flat spectrum radio quasars and 22 BL Lac objects. The fractional variability amplitudes of each source have been calculated in both optical R band and infrared J band. Overall, the variability amplitudes of FSRQs are larger than those of BL Lac objects. The results also suggest that the variability amplitude of most FSRQs is larger at lower energy band (J band) than at higher one (R band), while the variability amplitude of BL Lacs are larger at higher energy band. Two types of blazars both display color variation along with the variability in brightness. However, they show different variation behaviors in general. With the whole data set, 35 FSRQs exhibit redder-when-brighter trends, and 11 FSRQs exhibit opposite trends; 11 BL Lacs follow bluer-when-brighter trends, and 7 BL Lacs follow opposite trends. The examination in detail shows that there are 10 blazars showing redder-when-brighter trend in their low state, and bluer-when-brighter or stable-when-brighter trends in their high state. Some more complicated color behaviors have also been detected in several blazars. The non-thermal jet emission and the thermal emission from the accretion disc are employed to explain the observed color behaviors.

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.

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.

PKS 2155-304 is a well studied BL Lac object in the southern sky. The historical optical data during different period have been collected and compiled. Light curves with a time span of 35 years have been constructed. The R-band light curve has been analyzed by means of three methods: epoch folding method, Jurkevich method and discrete correlation function (DCF) method. It is derived that there is an evident periodic component of 317 days (i.e. 0.87 yr) superposed on a long-term trend with large-amplitude variation in the light curve. The variability of this source is accompanied by a slight color variation, and the brightness and color index are correlated with each other. On the long time-scale, PKS 2155-304 exhibits a tendency of bluer-when-brighter, which means the spectrum becomes flatter when the source brightens.

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 investigate the electron energy distributions (EEDs) and the corresponding acceleration processes in the jet of PKS 0447$-$439 and estimate its redshift through modeling its observed spectral energy distribution (SED) in the frame of a one-zone synchrotron-self Compton (SSC) model. Three EEDs formed in different acceleration scenarios are assumed: the power-law with exponential cut-off (PLC) EED (shock-acceleration scenario or the case of the EED approaching equilibrium in the stochastic-acceleration scenario), the log-parabolic (LP) EED (stochastic-acceleration scenario and the acceleration dominating) and the broken power law (BPL) EED (no acceleration scenario), and then the corresponding fluxes of both synchrotron and SSC are calculated. The model is applied to PKS 0447-439 and modeling SEDs are compared to the observed SED of this object by using the Markov Chain Monte Carlo (MCMC) method. Calculating results show that PLC model fails to fit the observed SED well, while the LP and BPL models give comparably good fits for the observed SED. The results indicate that it is possible that stochastic acceleration process acts in the emitting region of PKS 0447-439 and the EED is far from equilibrium (acceleration dominating) or no acceleration process works (in the emitting region). The redshift of PKS 0447-439 is also estimated in our fitting, and $z=0.16\pm0.05$ for LP case and $z=0.17\pm0.04$ for BPL case.

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.

The quasi-simultaneous multi-wavelength emission of TeV blazar 3C 66A is studied by using a one-zone multi-component leptonic jet model. It is found that the quasi-simultaneous spectral energy distribution (SED) of 3C 66A can be well reproduced, especially its {\it Fermi}-LAT first 3 months average spectrum can be well reproduced by the synchrotron-self Compton (SSC) component plus external Compton (EC) component of the broad line region (BLR). Clues on its redshift and gamma-ray emission location are obtained. The results indicate the following. (i) On the redshift; The theoretical intrinsic TeV spectra can be predicted by extrapolating the reproduced GeV spectra. Through comparing this extrapolated TeV spectra with the extragalactic background light (EBL) corrected observed TeV spectra, it is suggested that the redshift of 3C 66A could be between 0.1 and 0.3, the most likely value is $\sim$ 0.2. (ii) On the gamma-ray emission location; To well reproduce the GeV emission of 3C 66A under different assumptions on BLR, the gamma-ray emission region is always required to be beyond the inner zone of BLR. The BLR absorption effect on gamma-ray emission confirms this point.

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.

We present a detailed comparison between the 2-10 keV hard X-ray and infrared (IR) luminosity function (LF) of active galactic nuclei (AGN). The composite X-ray to IR spectral energy distributions (SEDs) of AGN used for connecting the hard X-ray LF (HXLF) and IR LF (IRLF) are modeled with a simple but well tested torus model based on the radiative transfer and photoionization code CLOUDY. Four observational determinations of the evolution of 2-10 keV HXLF and six evolution models of the obscured type-2 AGN fraction ($f_2$) have been considered. The 8.0 and 15 \micron LFs for the total, unobscured type-1 and obscured type-2 AGN are predicted from the HXLFs, and then compared with the measurements currently available. We find that the IRLFs predicted from HXLFs tend to underestimate the number of the most IR-luminous AGN. This is independent of the choices of HXLF and $f_2$, and even more obvious for the HXLFs recently measured. We show that the discrepancy between the HXLFs and IRLFs can be largely resolved when the anticorrelation between the UV to X-ray slope $\alpha_{\mathrm{ox}}$ and UV luminosity $L_{\rm UV}$ is appropriately considered. We also discuss other possible explanations for the discrepancy, such as the missing population of Compton-thick AGN and possible contribution of star-formation in the host to the mid-IR. Meanwhile, we find that the HXLFs and IRLFs of AGN can be more consistent with each other if the obscuration mechanisms of quasars and Seyferts are assumed to be different, corresponding to their different triggering and fueling mechanisms. More accurate measurements of the IRLFs of AGN, especially that determined at smaller redshift bins and more accurately separated to that for type-1 and type-2, are very helpful for clarifying these interesting issues.

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.

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.

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.

S5 0716+714 is a well-studied BL Lac object in the sky. Verifying the existence of correlations among the flux variations in different bands serves as an important tool to investigate the emission processes. To examine the possible existence of a lag between variations in different optical bands on this source, we employ a discrete correlation function (DCF) analysis on the light curves. In order to obtain statistically meaningful values for the cross-correlation time lags and their related uncertainties, we perform Monte Carlo simulations called "flux redistribution/random subset selection" (FR/RSS). Our analysis confirms that the variations in different optical light curves are strongly correlated. The time lags show a hint of the variations in high frequency band leading those in low frequency band of the order of a few minutes.

The synchrotron self-Compton (SSC) models and External Compton (EC) models of AGN jets with continually longitudinal and transverse bulk velocity structures are constructed. The observed spectra show complex and interesting patterns in different velocity structures and viewing angles. These models are used to calculate the synchrotron and inverse Compton spectra of two typical BL Lac objects (BLO) (Mrk 421 and 0716+714) and one Flat Spectrum Radio Quasars (FSRQs) (3c 279), and to discuss the implications of jet bulk velocity structures in unification of the BLO and FR I radio galaxies (FRI). By calculating the synchrotron spectra and SSC spectra of BL Lac object jets with continually bulk velocity structures, we find that the spectra are much different from ones in jets with uniform velocity structure under the increase of viewing angles. The unification of BLO and FRI is less constrained by viewing angles and would be imprinted by velocity structures intrinsic to the jet themselves. By considering the jets with bulk velocity structures constrained by apparent speed, we discuss the velocity structures imprinted on the observed spectra for different viewing angles. We find that the spectra are greatly impacted by longitudinal velocity structures, becasue the volume elements are compressed or expanded. Finally, we present the EC spectra of FSRQs and FR II radio galaxies (FRII) and find that they are weakly affected by velocity structures compared to synchrotron and SSC spectra.

Some studies suggested that a correlation between locations of BL Lacertae objects (BL Lacs) and the arrival directions of the ultrahigh energy cosmic rays (UHECRs) exists. Especially by assuming the primary particles charged +1 and using a galactic magnetic field (GMF) model to calculate the deflections of the UHECRs, the significance of correlation is improved. We construct a new GMF model by incorporating all progresses in the GMF measurements in recent years. Based on a thorough study of the deflections of the UHECRs measured by the AGASA experiment, we study the GFM model dependence of the correlation between the UHECRs and the selected BL Lacs using the new model together with others. It turns out that only specific one of those GMF models makes the correlation significant, even if neither GMF models themselves nor deflections of the UHECRs are not significantly different. It suggests that the significance of the correlation, calculated using a method suggested in those studies, is intensively depending on the GMF model. Great improvement in statistics may help to suppress the sensitivity to the GMF models.

We present the optical light curve for the gamma-ray blazar PKS 1510-089 with a significant variability of 2.0 mag within about 41 minutes in the R band which is the most violently rapid variability detected in our optical monitoring campaign operating since 1982. In the framework of the relativistic beaming model, the central supermassive black hole mass of PKS 1510-089 is of $(1.16-1.98)\times10^8M_\odot$. This is larger than some previous estimates for the central mass in blazars that have been based on models for nonthermal emission, but is consistent with the estimated result of the $M_\bullet-\sigma_\star$ correlation. Some other physical parameters, including the emission region size, Doppler factor $\delta$, the efficiency $\eta$ for conversion of accreted matter into energy and luminosity are calculated. The value of $\eta$ = 62.2 strongly implies that the relativistic beaming is responsible for the rapid variability of the gamma-ray loud blazar PKS 1510-089.

Using the sample of radio selected BL Lacertae objects (RBLs) and X-ray selected BL Lacertae objects (XBLs) presented by Sambruna et al. (1996), we calculated the luminosities of radio, optical and X-ray of each source and made the statistical analysis among the luminosities at different wave-bands, broad-band spectral indices from radio to X-ray ($\alpha_{\rm rx}$) and peak frequencies ($\nu_p$). Our results are as follows: (i) there is a positive correlation between radio luminosity $L_{\rm r}$ and $\alpha_{\rm rx}$ and a negative correlation between $L_{\rm r}$ and $\nu_p$. High-energy peak BL Lacs (HBLs) and low-energy peak BL Lacs (LBLs) can be distinguished very well, the dividing lines are probably those of $\log {L_{\rm r}}=43.25$ (erg/sec) and $\alpha_{\rm rx}>$(or $\leq $)0.75 for $L_{\rm r}$ - $\alpha_{\rm rx}$ plot and those of $\log {L_{\rm r}}\leq 43.25$ (erg/sec) and $\log {\nu_p}>14.7$ for the $L_{\rm r}$ - $\nu_p$ plot; (ii) there is a weak positive correlation between optical luminosity $L_o$ and $\alpha_{\rm rx}$ and a negatively weak correlation between $L_{\rm o}$ and $\nu_p$; (iii) there is no correlation between X-ray luminosity $L_X$ and $\alpha_{\rm rx}$ or between $L_X$ and $\nu_p$. From our analysis, we find that synchrotron radiation is the main X-ray radiation mechanism for HBLs while inverse Compton scattering for LBLs.