We report studies of the relationships between the total bolometric luminosity ($L_{\rm bol}$ or $L_{\rm TIR}$) and the molecular line luminosities of $J=1-0$ transitions of H$^{13}$CN, H$^{13}$CO$^+$, HCN, and HCO$^+$ with data obtained from ACA observations in the "ATOMS" survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities $L'_{\rm mol}$ of the four transitions all appear to be approximately linear. Line emission of isotopologues shows as large scatters in $L_{\rm bol}$-$L'_{\rm mol}$ relations as their main line emission. The log($L_{\rm bol}$/$L'_{\rm mol}$) for different molecular line tracers have similar distributions. The $L_{\rm bol}$-to-$L'_{\rm mol}$ ratios do not change with galactocentric distances ($R_{\rm GC}$) and clump masses ($M_{\rm clump}$). The molecular line luminosity ratios (HCN-to-HCO$^+$, H$^{13}$CN-to-H$^{13}$CO$^+$, HCN-to-H$^{13}$CN and HCO$^+$-to-H$^{13}$CO$^+$) all appear constant against $L_{\rm bol}$, dust temperature ($T_{\rm d}$), $M_{\rm clump}$ and $R_{\rm GC}$. Our studies suggest that both the main lines and isotopologue lines are good tracers of the total masses of dense gas in Galactic molecular clumps. The large optical depths of main lines do not affect the interpretation of the slopes in star formation relations. We find that the mean star formation efficiency (SFE) of massive Galactic clumps in the "ATOMS" survey is reasonably consistent with other measures of the SFE for dense gas, even those using very different tracers or examining very different spatial scales.

The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the "ATOMS" survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS $J = 2-1$, HCO$^+$ $J = 1-0$ and HCN $J = 1-0$, are found to show extended gas emission in low density regions within the clump; less than 25\% of their emission is from dense cores. SO, CH$_3$OH, H$^{13}$CN and HC$_3$N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over $\sim$1 pc), which may be caused by slow shocks from large--scale colliding flows or H{\sc ii} regions. Stellar feedback from an expanding H{\sc ii} region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analyzed with other survey data, e.g., "MALT90", "Orion B", "EMPIRE", "ALMA\_IMF", and "ALMAGAL", to deepen our understandings of "dense gas" star formation scaling relations and massive proto-cluster formation.

We have introduced a new method of estimating the electron temperature and density of H II regions by using single dish observations. In this method, multiple hydrogen radio recombination lines of different bands are computed under the assumption of low optical depth. We use evolutionary hydrodynamical models of H II regions to model hydrogen recombination line emission from a variety of H II regions and assess the reliability of the method. According to the simulated results, the error of the estimated temperature is commonly < 13%, and that of the estimated density is < 25% for a < 1% uncertainty of the observed line fluxes. A reasonable estimated value of electron density can be achieved if the uncertainty of the line fluxes are lower than 3%. In addition, the estimated values are more representative of the properties in the relatively high-density region if the gas density gradient is present in the H II region. Our method can be independent of the radio continuum observations. But the accuracy will be improved if a line-to-continuum ratio at millimeter wavelengths is added to the estimation. Our method provides a way to measure the temperature and density in ionized regions without interferometers.

The Large High Altitude Air Shower Observatory (LHAASO) project is a new generation multi-component instrument, to be built at 4410 meters of altitude in the Sichuan province of China, with the aim to study with unprecedented sensitivity the spec trum, the composition and the anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$ eV, as well as to act simultaneously as a wide aperture (one stereoradiant), continuously-operated gamma ray telescope in the energy range between 10$^{11}$ and $10^{15}$ eV. The experiment will be able of continuously surveying the TeV sky for steady and transient sources from 100 GeV to 1 PeV, t hus opening for the first time the 100-1000 TeV range to the direct observations of the high energy cosmic ray sources. In addition, the different observables (electronic, muonic and Cherenkov/fluorescence components) that will be measured in LHAASO will allow to investigate origin, acceleration and propagation of the radiation through a measurement of energy spec trum, elemental composition and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO in cosmic rays physics and gamma astronomy would play a key-role in the comprehensive general program to explore the High Energy Universe. LHAASO will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. In this document we introduce the concept of LHAASO and the main science goals, providing an overview of the project.

Observations in the local universe show a tight correlation between the masses of supermassive black holes (SMBHs; $M_{\rm BH}$) and host-galaxy bulges ($M_{\rm bulge}$), suggesting a strong connection between SMBH and bulge growth. However, direct evidence for such a connection in the distant universe remains elusive. We have studied sample-averaged SMBH accretion rate ($\overline{\rm BHAR}$) for bulge-dominated galaxies at $z=0.5-3$. While previous observations found $\overline{\rm BHAR}$ is strongly related to host-galaxy stellar mass ($M_\star$) for the overall galaxy population, our analyses show that, for the bulge-dominated population, $\overline{\rm BHAR}$ is mainly related to SFR rather than $M_\star$. This ${\overline{\rm BHAR}}$-SFR relation is highly significant, e.g. $9.0\sigma$ (Pearson statistic) at $z=0.5-1.5$. Such a $\overline{\rm BHAR}$-SFR connection does not exist among our comparison sample of galaxies that are not bulge-dominated, for which $M_\star$ appears to be the main determinant of SMBH accretion. This difference between the bulge-dominated and comparison samples indicates that SMBHs only coevolve with bulges rather than the entire galaxies, explaining the tightness of the local $M_{\rm BH}$-$M_{\rm bulge}$ correlation. Our best-fit $\overline{\rm BHAR}$-SFR relation for the bulge-dominated sample is $\log\overline{\rm BHAR} = \log\mathrm{SFR} - (2.48\pm0.05)$ (solar units). The best-fit $\overline{\rm BHAR}/\mathrm{SFR}$ ratio ($10^{-2.48}$) for bulge-dominated galaxies is similar to the observed $M_{\rm BH}/M_{\rm bulge}$ values in the local universe. Our results reveal that SMBH and bulge growth are in lockstep, and thus non-causal scenarios of merger averaging are unlikely the origin of the $M_{\rm BH}$-$M_{\rm bulge}$ correlation. This lockstep growth also predicts that the $M_{\rm BH}$-$M_{\rm bulge}$ relation should not have strong redshift dependence.

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.

Jet launching in radio loud (RL) quasars is one of the fundamental problems in astrophysics. Exploring the differences in the inner accretion disk properties between RL and radio quiet (RQ) quasars might yield helpful clues to this puzzle. We previously discovered that the shorter term UV/optical variations of quasars are bluer than the longer term ones, i.e., the so-called timescale-dependent color variation. This is consistent with the scheme that the faster variations come from the inner and hotter disk regions, thus providing a useful tool to map the accretion disk which is otherwise unresolvable. In this work we compare the UV/optical variations of RL quasars in SDSS Stripe 82 to those of several RQ samples, including those matched in redshift-luminosity-black hole mass and/or color-magnitude. We find that while both RL and RQ populations appear bluer when they brighten, RL quasars potentially show a weaker/flatter dependence on timescale in their color variation. We further find that while both RL and RQ populations on average show similar variation amplitudes at long timescales, fast variations of RL sources appear weaker/smaller (at timescales of ~ 25 -- 300 days in the observer's frame), and the difference is more prominent in the g-band than in the r-band. Inhomogeneous disk simulations can qualitatively reproduce these observed differences if the inner accretion disk of RL quasars fluctuates less based on simple toy models. Though the implications are likely model dependent, the discovery points to an interesting diagram that magnetic fields in RL quasars may be prospectively stronger and play a key role in both jet launching and the stabilization of the inner accretion disk.

We develop a new method, which is based on the optimal redshift weighting scheme, to extract the maximal tomographic information of baryonic acoustic oscillations (BAO) and redshift space distortions (RSD) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 quasar (DR14Q) survey. We validate our method using the EZ mocks, and apply our pipeline to the eBOSS DR14Q sample in the redshift range of $0.8<z<2.2$. We report a joint measurement of $f\sigma_8$ and two-dimensional BAO parameters $D_{\rm A}$ and $H$ at four effective redshifts of $z_{\rm eff}=0.98, 1.23, 1.52$ and $1.94$, and provide the full data covariance matrix. Using our measurement combined with BOSS DR12, MGS and 6dFGS BAO measurements, we find that the existence of dark energy is supported by observations at a $7.4\sigma$ significance level. Combining our measurement with BOSS DR12 and Planck observations, we constrain the gravitational growth index to be $\gamma=0.580\pm0.082$, which is fully consistent with the prediction of general relativity. This paper is part of a set that analyses the eBOSS DR14 quasar sample.

Observations have shown that UV/optical variation amplitude of quasars depend on several physi- cal parameters including luminosity, Eddington ratio, and likely also black hole mass. Identifying new factors which correlate with the variation is essential to probe the underlying physical processes. Combining ~ten years long quasar light curves from SDSS stripe 82 and X-ray data from Stripe 82X, we build a sample of X-ray detected quasars to investigate the relation between UV/optical variation amplitude ($\sigma_{rms}$) and X-ray loudness. We find that quasars with more intense X-ray radiation (com- pared to bolometric luminosity) are more variable in UV/optical. Such correlation remains highly significant after excluding the effect of other parameters including luminosity, black hole mass, Ed- dington ratio, redshift, rest-frame wavelength (i.e., through partial correlation analyses). We further find the intrinsic link between X-ray loudness and UV/optical variation is gradually more prominent on longer timescales (up to 10 years in the observed frame), but tends to disappear at timescales < 100 days. This suggests a slow and long-term underlying physical process. The X-ray reprocessing paradigm, in which UV/optical variation is produced by a variable central X-ray emission illuminating the accretion disk, is thus disfavored. The discovery points to an interesting scheme that both the X-ray corona heating and UV/optical variation is quasars are closely associated with magnetic disc turbulence, and the innermost disc turbulence (where corona heating occurs) correlates with the slow turbulence at larger radii (where UV/optical emission is produced).

The variation in the high energy cut-off E_c in active galactic nuclei uniquely probes the corona physics. In this work we show that the ratio of two NuSTAR spectra (in analogy to difference imaging technique widely used in astronomy) is uniquely useful in studying E_c variations. The spectra ratio could directly illustrate potential E_c variation between two spectra. By comparing with the ratio of two spectral fitting models, it also examines the reliability of the spectral-fitting measured E_c variation. Assisted with this technique, we revisit the 5 AGNs in literature (MCG -5-23-16, 3C 382, NGC 4593, NGC 5548 and Mrk 335) for which E_c (kT_e) variations have been claimed with NuSTAR observations. We show the claimed E_c variations appear inconsistent with the spectra ratios in three of them, thus need to be revised, demonstrating the striking usefulness of spectra ratio. We present thereby improved spectral fitting results and E_c variations. We also report a new source with E_c variations based on NuSTAR observations (radio galaxy 4C +74.26). We find the corona tends to be hotter when it brightens (hotter-when-brighter) in 3C 382, NGC 5548, Mrk 335 and 4C +74.27, but MCG -5-23-16 and NGC 4593 show no evidence of significant E_c variations. Meanwhile all 6 sources in this small sample appear softer-when-brighter. Changes in corona geometry are required to explain the observed hotter-when-brighter trends.

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.

The reprocessing scenario is widely adopted in literature to explain the observed tight inter-band correlation and short lags in the UV/optical variations of active galactic nuclei (AGNs). In this work we look into the color variability of the famous Seyfert galaxy NGC 5548 with high quality Swift multi-band UV/optical light curves. We find the color variation of NGC 5548 is clearly timescale-dependent, in a way that it is more prominent at shorter timescales. This is similar to that previously detected in quasar samples, but for the first time in an individual AGN. We show that while a reprocessing model with strict assumptions on the driving source and the disk size can apparently match the observed light curves and inter-band lags, it fails to reproduce the observed timescale dependency in the color variation. Such discrepancy raises a severe challenge to, and can hardly be reconciled under the widely accepted reprocessing diagram. It also demonstrates that the timescale dependency of the color variation is uniquely powerful in probing the physics behind AGN UV/optical variations.

The tight inter-band correlation and the lag-wavelength relation among UV/optical continua of active galactic nuclei have been firmly established. They are usually understood within the widespread reprocessing scenario, however, the implied inter-band lags are generally too small. Furthermore, it is challenged by new evidences, such as the X-ray reprocessing yields too much high frequency UV/optical variations as well as it fails to reproduce the observed timescale-dependent color variations among {\it Swift} lightcurves of NGC 5548. In a different manner, we demonstrate that an upgraded inhomogeneous accretion disk model, whose local {\it independent} temperature fluctuations are subject to a speculated {\it common} large-scale temperature fluctuation, can intrinsically generate the tight inter-band correlation and lag across UV/optical, and be in nice agreement with several observational properties of NGC 5548, including the timescale-dependent color variation. The emergent lag is a result of the {\it differential regression capability} of local temperature fluctuations when responding to the large-scale fluctuation. An average speed of propagations as large as $\gtrsim 15\%$ of the speed of light may be required by this common fluctuation. Several potential physical mechanisms for such propagations are discussed. Our interesting phenomenological scenario may shed new light on comprehending the UV/optical continuum variations of active galactic nuclei.

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 present a measurement of the anisotropic and isotropic Baryon Acoustic Oscillations (BAO) from the extended Baryon Oscillation Spectroscopic Survey Data Release 14 quasar sample with optimal redshift weights. Applying the redshift weights improves the constraint on the BAO dilation parameter $\alpha(z_{\rm eff})$ by 17\%. We reconstruct the evolution history of the BAO distance indicators in the redshift range of $0.8<z<2.2$. This paper is part of a set that analyses the eBOSS DR14 quasar sample.

We measure the growth rate and its evolution using the anisotropic clustering of the extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 (DR14) quasar sample, which includes $148\,659$ quasars covering the wide redshift range of $0.8 < z < 2.2$ and a sky area of $2112.90$ $\rm deg^2$. To optimise measurements we deploy a redshift-dependent weighting scheme, which allows us to avoid binning, and perform the data analysis consistently including the redshift evolution across the sample. We perform the analysis in Fourier space, and use the redshift evolving power spectrum multipoles to measure the redshift space distortion parameter $f\sigma_8$ and parameters controlling the anisotropic projection of the cosmological perturbations. We measure $f \sigma_8(z=1.52)=0.43 \pm 0.05 $ and $df\sigma_8/dz (z=1.52)= - 0.16 \pm 0.08$, consistent with the expectation for a $\Lambda$CDM cosmology as constrained by the Planck experiment.

We present an anisotropic analysis of Baryon Acoustic Oscillation (BAO) signal from the SDSS-IV extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 (DR14) quasar sample. The sample consists of 147,000 quasars distributed over a redshift range of $0.8 < z < 2.2$. We apply the redshift weights technique to the clustering of quasars in this sample and achieve a 4.6 per cent measurement of the angular distance measurement $D_M$ at $z = 2.2$ and Hubble parameter $H$ at $z=0.8$. We parameterize the distance-redshift relation, relative to a fiducial model, as a quadratic expansion. The coefficients of this expansion are used to reconstruct the distance-redshift relation and obtain distance and Hubble parameter measurements at all redshifts within the redshift range of the sample. Reporting the result at two characteristic redshifts, we determine $D_M(z=1) = 3405\pm305 \ (r_{\rm d} / r_{\rm d, fid}) \ {\rm Mpc}$, $H(z=1) = 120.7\pm 7.3 \ (r_{\rm d,fid} / r_{\rm d}) \ {\rm km} \ {\rm s}^{-1}{\rm Mpc}^{-1}$ and $D_M(z=2) = 5325\pm249 \ (r_{\rm d} / r_{\rm d, fid}) \ {\rm Mpc}$, $H(z=2) = 189.9\pm 32.9 \ (r_{\rm d,fid} / r_{\rm d}) \ {\rm km} \ {\rm s}^{-1}{\rm Mpc}^{-1}$. These measurements are highly correlated. We assess the outlook of BAO analysis from the final quasar sample by testing the method on a set of mocks that mimic the noise level in the final sample. We demonstrate on these mocks that redshift weighting shrinks the measurement error by over 25 per cent on average. We conclude redshift weighting can bring us closer to the cosmological goal of the final quasar sample.

The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is one of the first of a new generation of galaxy redshift surveys that will cover a large range in redshift with sufficient resolution to measure the baryon acoustic oscillations (BAO) signal. For surveys covering a large redshift range we can no longer ignore cosmological evolution, meaning that either the redshift shells analysed have to be significantly narrower than the survey, or we have to allow for the averaging over evolving quantities. Both of these have the potential to remove signal: analysing small volumes increases the size of the Fourier window function, reducing the large-scale information, while averaging over evolving quantities can, if not performed carefully, remove differential information. It will be important to measure cosmological evolution from these surveys to explore and discriminate between models. We apply a method to optimally extract this differential information to mock catalogues designed to mimic the eBOSS quasar sample. By applying a set of weights to extract redshift space distortion measurements as a function of redshift, we demonstrate an analysis that does not invoke the problems discussed above. We show that our estimator gives unbiased constraints.

We present measurements of the Baryon Acoustic Oscillation (BAO) scale in redshift-space using the clustering of quasars. We consider a sample of 147,000 quasars from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) distributed over 2044 square degrees with redshifts $0.8 < z < 2.2$ and measure their spherically-averaged clustering in both configuration and Fourier space. Our observational dataset and the 1400 simulated realizations of the dataset allow us to detect a preference for BAO that is greater than 2.8$\sigma$. We determine the spherically averaged BAO distance to $z = 1.52$ to 3.8 per cent precision: $D_V(z=1.52)=3843\pm147 \left(r_{\rm d}/r_{\rm d, fid}\right)\ $Mpc. This is the first time the location of the BAO feature has been measured between redshifts 1 and 2. Our result is fully consistent with the prediction obtained by extrapolating the Planck flat $\Lambda$CDM best-fit cosmology. All of our results are consistent with basic large-scale structure (LSS) theory, confirming quasars to be a reliable tracer of LSS, and provide a starting point for numerous cosmological tests to be performed with eBOSS quasar samples. We combine our result with previous, independent, BAO distance measurements to construct an updated BAO distance-ladder. Using these BAO data alone and marginalizing over the length of the standard ruler, we find $\Omega_{\Lambda} > 0$ at 6.6$\sigma$ significance when testing a $\Lambda$CDM model with free curvature.

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 perform a tomographic baryon acoustic oscillations analysis using the two-point galaxy correlation function measured from the combined sample of BOSS DR12, which covers the redshift range of $0.2<z<0.75$. Splitting the sample into multiple overlapping redshift slices to extract the redshift information of galaxy clustering, we obtain a measurement of $D_A(z)/r_d$ and $H(z)r_d$ at nine effective redshifts with the full covariance matrix calibrated using MultiDark-Patchy mock catalogues. Using the reconstructed galaxy catalogues, we obtain the precision of $1.3\%-2.2\%$ for $D_A(z)/r_d$ and $2.1\%-6.0\%$ for $H(z)r_d$. To quantify the gain from the tomographic information, we compare the constraints on the cosmological parameters using our 9-bin BAO measurements, the consensus 3-bin BAO and RSD measurements at three effective redshifts in \citet{Alam2016}, and the non-tomographic (1-bin) BAO measurement at a single effective redshift. Comparing the 9-bin with 1-bin constraint result, it can improve the dark energy Figure of Merit by a factor of 1.24 for the Chevallier-Polarski-Linder parametrisation for equation of state parameter $w_{\rm DE}$. The errors of $w_0$ and $w_a$ from 9-bin constraints are slightly improved when compared to the 3-bin constraint result.

We perform a tomographic baryon acoustic oscillations (BAO) analysis using the monopole, quadrupole and hexadecapole of the redshift-space galaxy power spectrum measured from the pre-reconstructed combined galaxy sample of the completed Sloan Digital Sky Survey (SDSS-III) Baryon Oscillation Spectroscopic Survey (BOSS) Data Release (DR)12 covering the redshift range of $0.20<z<0.75$. By allowing for overlap between neighbouring redshift slices, we successfully obtained the isotropic and anisotropic BAO distance measurements within nine redshift slices to a precision of $1.5\%-3.4\%$ for $D_V/r_d$, $1.8\% -4.2\%$ for $D_A/r_d$ and $3.7\% - 7.5\%$ for $H \ r_d$, depending on effective redshifts. We provide our BAO measurement of $D_A/r_d$ and $H \ r_d$ with the full covariance matrix, which can be used for cosmological implications. Our measurements are consistent with those presented in \citet{Acacia}, in which the BAO distances are measured at three effective redshifts. We constrain dark energy parameters using our measurements, and find an improvement of the Figure-of-Merit of dark energy in general due to the temporal BAO information resolved. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS.

The low statistical errors on cosmological parameters promised by future galaxy surveys will only be realised with the development of new, fast, analysis methods that reduce potential systematic problems to low levels. We present an efficient method for measuring the evolution of the growth of structure using Redshift Space Distortions (RSD), that removes the need to make measurements in redshift shells. We provide sets of galaxy-weights that cover a wide range in redshift, but are optimised to provide differential information about cosmological evolution. These are derived to optimally measure the coefficients of a parameterisation of the redshift-dependent matter density, which provides a framework to measure deviations from the concordance $\Lambda$CDM cosmology, allowing for deviations in both geometric and/or growth. We test the robustness of the weights by comparing with alternative schemes and investigate the impact of galaxy bias. We extend the results to measure the combined anisotropic Baryon Acoustic Oscillation (BAO) and RSD signals.

In recent work done by Sun et. al., the color variation of quasars, namely the bluer-when-brighter trend, was found to be timescale-dependent using SDSS $g/r$ band light curves in the Stripe 82. Such timescale dependence, i.e., bluer variation at shorter timescales, supports the thermal fluctuation origin of the UV/optical variation in quasars, and can be well modeled with the inhomogeneous accretion disk model. In this paper, we extend the study to much shorter wavelengths in the rest frame (down to extreme UV), using GALaxy Evolution eXplorer (GALEX) photometric data of quasars collected in two ultraviolet bands (near-UV and far-UV). We develop Monte-Carlo simulations to correct possible biases due to the considerably larger photometric uncertainties in GALEX light curves (particularly in far-UV, comparing with SDSS $g/r$ bands), which otherwise could produce artificial results. We securely confirm the previously discovered timescale dependence of the color variability with independent datasets and at shorter wavelengths. We further find the slope of the correlation between the amplitude of color variation and timescale however appears even steeper than that predicted by the inhomogeneous disk model, which assumes that disk fluctuations follow damped random walk process. In line with the much flatter structure function observed in far-UV comparing with that at longer wavelengths, this implies deviation from DRW process in the inner disk where rest frame extreme UV radiation is produced.

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.

The gamma-ray observation of dwarf spheroidal satellites (dSphs) is an ideal approach for probing the dark matter (DM) annihilation signature. The latest Fermi-LAT dSph searches have set stringent constraints on the velocity independent annihilation cross section in the small DM mass range, which gives very strong constraints on the scenario to explain the the AMS-02 positron excess by DM annihilation. However, the dSph constraints would change in the velocity dependent annihilation scenarios, because the velocity dispersion in the dSphs varies from that in the Milky Way. In this work, we use a likelihood map method to set constraints on the velocity dependent annihilation cross section from the Fermi-LAT observation of six dSphs. We consider three typical forms of the annihilation cross section, i.e. p-wave annihilation, Sommerfeld enhancement, and Breit-Wigner resonance. For the p-wave annihilation and Sommerfeld-enhancement, the dSph limits would become much weaker and stronger compared with those for the velocity independent annihilation, respectively. For the Breit-Wigner annihilation, the dSph limits would vary depending on the model parameters. We show that the scenario to explain the AMS-02 positron excess by DM annihilation is still viable in the velocity dependent cases.

Large redshift surveys capable of measuring the Baryon Acoustic Oscillation (BAO) signal have proven to be an effective way of measuring the distance-redshift relation in cosmology. Building off the work in Zhu et al. (2015), we develop a technique to directly constrain the distance-redshift relation from BAO measurements without splitting the sample into redshift bins. We parametrize the distance-redshift relation, relative to a fiducial model, as a quadratic expansion. We measure its coefficients and reconstruct the distance-redshift relation from the expansion. We apply the redshift weighting technique in Zhu et al. (2015) to the clustering of galaxies from 1000 QuickPM (QPM) mock simulations after reconstruction and achieve a 0.75% measurement of the angular diameter distance $D_A$ at $z=0.64$ and the same precision for Hubble parameter H at $z=0.29$. These QPM mock catalogs are designed to mimic the clustering and noise level of the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12 (DR12). We compress the correlation functions in the redshift direction onto a set of weighted correlation functions. These estimators give unbiased $D_A$ and $H$ measurements at all redshifts within the range of the combined sample. We demonstrate the effectiveness of redshift weighting in improving the distance and Hubble parameter estimates. Instead of measuring at a single 'effective' redshift as in traditional analyses, we report our $D_A$ and $H$ measurements at all redshifts. The measured fractional error of $D_A$ ranges from 1.53% at $z=0.2$ to 0.75% at $z=0.64$. The fractional error of $H$ ranges from 0.75% at $z=0.29$ to 2.45% at $z = 0.7$. Our measurements are consistent with a Fisher forecast to within 10% to 20% depending on the pivot redshift. We further show the results are robust against the choice of fiducial cosmologies, galaxy bias models, and RSD streaming parameters.

We present a science forecast for the eBOSS survey, part of the SDSS-IV project, which is a spectroscopic survey using multiple tracers of large-scale structure, including luminous red galaxies (LRGs), emission line galaxies (ELGs) and quasars (both as a direct probe of structure and through the Ly-$\alpha$ forest). Focusing on discrete tracers, we forecast the expected accuracy of the baryonic acoustic oscillation (BAO), the redshift-space distortion (RSD) measurements, the $f_{\rm NL}$ parameter quantifying the primordial non-Gaussianity, the dark energy and modified gravity parameters. We also use the line-of-sight clustering in the Ly-$\alpha$ forest to constrain the total neutrino mass. We find that eBOSS LRGs ($0.6<z<1.0$) (combined with the BOSS LRGs at $z>0.6$), ELGs ($0.6<z<1.2$) and Clustering Quasars (CQs) ($0.6<z<2.2$) can achieve a precision of 1%, 2.2% and 1.6% precisions, respectively, for spherically averaged BAO distance measurements. Using the same samples, the constraint on $f\sigma_8$ is expected to be 2.5%, 3.3% and 2.8% respectively. For primordial non-Gaussianity, eBOSS alone can reach an accuracy of $\sigma(f_{\rm NL})\sim10-15$, depending on the external measurement of the galaxy bias and our ability to model large-scale systematic errors. eBOSS can at most improve the dark energy Figure of Merit (FoM) by a factor of $3$ for the Chevallier-Polarski-Linder (CPL) parametrisation, and can well constrain three eigenmodes for the general equation-of-state parameter (Abridged).

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 knee phenomenon of the cosmic ray spectrum, which plays an important role in studying the acceleration mechanism of cosmic rays, is still an unsolved mystery. We try to reconcile the knee spectra measured by ARGO-YBJ and Tibet-III. A simple broken power-law model fails to explain the experimental data. Therefore a modified broken power-law model with non-linear acceleration effects is adopted, which can describe the sharp knee structure. This model predicts that heavy elements dominate at the knee.

We have developed a full numerical method to study the gas dynamics of cometary ultra-compact (UC) H II regions, and associated photodissociation regions (PDRs). The bow-shock and champagne-flow models with a $40.9/21.9 M_\odot$ star are simulated. In the bow-shock models, the massive star is assumed to move through dense ($n=8000~cm^{-3}$) molecular material with a stellar velocity of $15~km~s^{-1}$. In the champagne-flow models, an exponential distribution of density with a scale height of 0.2 pc is assumed. The profiles of the [Ne II] 12.81\mum and $H_2~S(2)$ lines from the ionized regions and PDRs are compared for two sets of models. In champagne-flow models, emission lines from the ionized gas clearly show the effect of acceleration along the direction toward the tail due to the density gradient. The kinematics of the molecular gas inside the dense shell is mainly due to the expansion of the H II region. However, in bow-shock models the ionized gas mainly moves in the same direction as the stellar motion. The kinematics of the molecular gas inside the dense shell simply reflects the motion of the dense shell with respect to the star. These differences can be used to distinguish two sets of models.

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.

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 simulate evolution of cometary H II regions based on several champagne flow models and bow shock models, and calculate the profiles of the [Ne II] fine-structure line at $12.81\mu m$, the $H30\alpha$ recombination line and the [Ne III] fine-structure line at $15.55\mu m$ for these models at different inclinations of $0^o, 30^o \textrm{and} 60^o$. We find that the profiles in the bow shock models are generally different from those in the champagne flow models, but the profiles in the bow shock with lower stellar velocity ($\leq5km s^{-1}$) are similar to those in the champagne flow models. In champagne flow models, both the velocity of peak flux and the flux weighted central velocities of all three lines are pointing outward from molecular clouds. In bow shock models, the directions of these velocities rely on the speed of stars. They have the similar motion in high stellar speed case but opposite directions in low stellar speed case. We notice that the line profiles from the slit along the symmetrical axis of the projected 2D image of these models are useful for distinguishing bow shock models and champagne flow models. It is also confirmed by the calculation that the flux weighted central velocity and the line luminosity of the [Ne III] line can be estimated from the [Ne II] line and the $H30\alpha$ line.

Future baryon acoustic oscillation (BAO) surveys will survey very large volumes, covering wide ranges in redshift. We derive a set of redshift weights to compress the information in the redshift direction to a small number of modes. We suggest that such a compression preserves almost all of the signal for most cosmologies, while giving high signal-to-noise measurements for each combination. We present some toy models and simple worked examples. As an intermediate step, we give a precise meaning to the "effective redshift" of a BAO measurement.

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.

Using deep narrow-band $H_2S1$ and $K_{s}$-band imaging data obtained with CFHT/WIRCam, we identify a sample of 56 H$\alpha$ emission-line galaxies (ELGs) at $z=2.24$ with the 5$\sigma$ depths of $H_2S1=22.8$ and $K_{s}=24.8$ (AB) over 383 arcmin$^{2}$ area in the ECDFS. A detailed analysis is carried out with existing multi-wavelength data in this field. Three of the 56 H$\alpha$ ELGs are detected in Chandra 4 Ms X-ray observation and two of them are classified as AGNs. The rest-frame UV and optical morphologies revealed by HST/ACS and WFC3 deep images show that nearly half of the H$\alpha$ ELGs are either merging systems or with a close companion, indicating that the merging/interacting processes play a key role in regulating star formation at cosmic epoch z=2-3; About 14% are too faint to be resolved in the rest-frame UV morphology due to high dust extinction. We estimate dust extinction from SEDs. We find that dust extinction is generally correlated with H$\alpha$ luminosity and stellar mass (SM). Our results suggest that H$\alpha$ ELGs are representative of star-forming galaxies (SFGs). Applying extinction correction for individual objects, we examine the intrinsic H$\alpha$ luminosity function (LF) at $z=2.24$, obtaining a best-fit Schechter function characterized by a faint-end slope of $\alpha=-1.3$. This is shallower than the typical slope of $\alpha \sim -1.6$ in previous works based on constant extinction correction. We demonstrate that this difference is mainly due to the different extinction corrections. The proper extinction correction is thus key to recovering the intrinsic LF as the extinction globally increases with H$\alpha$ luminosity. Moreover, we find that our H$\alpha$ LF mirrors the SM function of SFGs at the same cosmic epoch. This finding indeed reflects the tight correlation between SFR and SM for the SFGs, i.e., the so-called main sequence.

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.

Hard X-ray emission in radio-quiet active galactic nuclei (AGNs) is believed to be produced via inverse Compton scattering by hot and compact coronae near the super massive black hole. However the origin and physical properties of the coronae, including geometry, kinematics and dynamics, yet remain poorly known. In this work, taking [OIV] 25.89um emission line as an isotropic indicator of AGN's intrinsic luminosity, we compare the intrinsic corona X-ray emission between Seyfert 1 and Compton-thin Seyfert 2 galaxies, which are viewed at different inclinations according to the unification scheme. We compile a sample of 130 Compton-thin Seyfert galaxies with both [OIV] 25.89um line luminosities measured with Spitzer-IRS and X-ray spectra observed by XMM-Newton, Chandra, Suzaku or Swift. Known radio-loud sources are excluded. We fit the X-ray spectra to obtain the absorption-corrected 2-10 keV continuum luminosities. We find that Seyfert 1 galaxies are intrinsically brighter in intrinsic 2-10 keV emission by a factor of 2.8+0.5-0.4 (2.2+0.9-0.3 in Swift-BAT 14-195 keV emission), comparing with Compton-thin Seyfert 2 galaxies. The Seyfert 1 and Compton-thin Seyfert 2 galaxies follow a statistically identical correlation between the absorption-corrected 2-10 keV luminosity and the 14-195 keV luminosity, indicating that our absorption correction to the 2-10 keV flux is sufficient. The difference in X-ray emission between the two populations is thus unlikely due to X-ray absorption, and instead implies an intrinsic anisotropy in the corona X-ray emission. This striking anisotropy of X-ray emission can be explained by a bipolar outflowing corona with a bulk velocity of ~0.3-0.5c. This would provide a natural link between the so-called coronae and weak jets in these systems. Other consequences of outflowing coronae are also discussed.

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

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

In our previous work, we have proposed two methods for computing the luminosity distance d_{L}^{\Lambda} in LCDM model. In this paper, two effective quadrature algorithms, known as Romberg Integration and composite Gaussian Quadrature, are presented to calculate the luminosity distance d_{L}^{CPL} in the Chevallier-Polarski-Linder parametrization(CPL) model. By comparing the efficiency and accuracy of the two algorithms, we find that the second is more promising. Moreover, we develop another strategy adapted for approximating d_{L}^{\Lambda} in flat LCDM universe. To some extent, our methods can make contributions to the recent numerical stimulation for the investigation of dark energy cosmology.

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.

The purpose of this work is to search the Yangbajing neutron monitor data obtained between 1998 October and 2000 June for solar neutrons associated with solar flares. Using the onset times of 166 BATSE-detected flares with the GOES peak flux (1 -- 8 \AA) higher than $1.0 \times 10^{-5}$ $\mathrm{Wm^{-2}}$, we prepare for each flare a light curve of the Yangbajing neutron monitor, spanning $\pm$ 1.5 hours from the BATSE onset time. Based on the light curves, a systematic search for solar neutrons in energies above 100 MeV from the 166 flares was performed. No statistically significant signals due to solar neutrons were found in the present work. Therefore, we put upper limits on the $>$ 100 MeV solar-neutron flux for 18 events consisting of 2 X and 16 M class flares. The calculation assumed a power-law shaped neutron energy spectrum and three types of neutron emission profiles at the Sun. Compared with the other positive neutron detections associated with X-class flares, typical 95% confidence level upper limits for the two X-class flares are found to be comparable to the lowest and second lowest neutron fluxes at the top of the atmosphere.In addition, the upper limits for M-class flares scatter in the range of $10^{-2}$ to 1 neutrons $\mathrm{cm^{-2}s^{-1}}$. This provides the first upper limits on the solar-neutron flux from M-class solar flares, using space observatories as well as ground-based neutron monitors.