Using multi-band photometric images of M 51 and its companion NGC 5195 from ultraviolet to optical and infrared, we investigate spatially resolved stellar population properties of this interacting system with stellar population synthesis models. The observed IRX is used to constrain dust extinction. Stellar mass is also inferred from the model fitting. By fitting observed spectral energy distributions (SEDs) with synthetical ones, we derive two-dimensional distributions of stellar age, metallicity, dust extinction, and stellar mass. In M 51, two grand-designed spiral arms extending from the bulge show young age, rich metallicity, and abundant dust. The inter-arm regions are filled with older, metal-poorer, and less dusty stellar populations. Except for the spiral arm extending from M 51 into NGC 5195, the stellar population properties of NGC 5195 are quite featureless. NGC 5195 is much older than M 51, and its core is very dusty with $A_V$ up to 1.67 mag and dense in stellar mass surface density. The close encounters might drive the dust in the spiral arm of M51 into the center of NGC 5195.

M51 and NGC 5195 is an interacting system that can be explored in great details with ground-based telescopes. The H II regions in M51 were observed using the 2.16 m telescope of the National Astronomical Observatories of the Chinese Academy of Sciences and the 6.5 m Multiple Mirror Telescope with spatial resolution of less than $\sim100$ pc. We obtain a total of 113 spectra across the galaxy and combine the literature data of Croxall et al. to derive a series of physical properties, including the gas-phase extinction, stellar population age, star formation rate (SFR) surface density, and oxygen abundance. The spatial distributions and radial profiles of these properties are investigated in order to study the characteristics of M51 and the clues to the formation and evolution of this galaxy. M51 presents a mild radial extinction gradient. The lower gas-phase extinction in the north spiral arms compared to the south arms are possibly caused by the past encounters with the companion galaxy of NGC 5195. A number of H II regions have the stellar age between 50 and 500 Myr, consistent with the recent interaction history by simulations in the literatures. The SFR surface density presents a mild radial gradient, which is ubiquitous in spiral galaxies. There is a negative metallicity gradient of $-0.08$ dex $R_{e}^{-1}$ in the disk region, which is also commonly found in many spiral galaxies. It is supported by the "inside-out" scenario of galaxy formation. We find a positive abundance gradient of 0.26 dex $R_{e}^{-1}$ in the inner region. There are possible reasons causing the positive gradient, including the freezing of the chemical enrichment due to the star-forming quenching in the bulge and the gas infall and dilution due to the pseudobulge growth and/or galactic interaction.

Site-testing is crucial for achieving the goal of scientific research and analysis of meteorological and optical observing conditions is one of the basic tasks of it. As one of three potential sites to host 12-meter Large Optical/infrared Telescope (LOT), Muztagh-ata site which is located on the Pamirs Plateau in west China's Xinjiang began its site-testing task in the spring of 2017. In this paper, we firstly start with an introduction to the site and then present a statistical analysis of the ground-level meteorological properties such as air temperature, barometric pressure, relative humidity, wind speed and direction, recorded by automatic weather station with standard meteorological sensors for two-year long. We also show the monitoring results of sky brightness during this period.

In this article, we present a detailed analysis of the statistical properties of seeing for the Muztagh-ata site which is the candidate site for hosting future Chinese Large Optical/infrared Telescope (LOT) project. The measurement was obtained with Differential Image Motion Monitor (DIMM) from April 2017 to November 2018 at different heights during different periods. The median seeing at 11 meters and 6 meters are very close but different significantly from that on the ground. We mainly analyzed the seeing at 11 meters monthly and hourly, having found that the best season for observing was from late autumn to early winter and seeing tended to improve during the night only in autumn. The analysis of the dependence on temperature inversion, wind speed, direction also was made and the best meteorological conditions for seeing is given.

We have assembled the historical light curves of the BL Lac Object 3FGL J0449.4-4350 at optical and gamma-ray bands, the time spanning about 10 years, analyzed the periodic variability of the light curves by using four different methods (Lomb-Scargle periodogram, REDFIT38, Jurkevich and DACF). We detected a marginally possible quasi-periodic oscillation (QPO) of ~450 days. Assuming it originates from the helical motion jet in a supermassive binary black hole (SMBBH) system undergoing major merger, we estimate the primary black hole mass M~7.7*10^{9} M_sun. To explore the origin of the gamma-ray, we investigated the optical-gamma-ray correlations using discrete correlation function (DCF) method, and found that the correlation between the two bands is very significant. This strong correlation tends to imply lepton self-synchro-Compton (LSSC) model to produce the gamma-ray.

In order to measure the Universal Time (UT1) in real time, National Time Service Center (NTSC) has built a VGOS-like (VLBI Global Observing System) broadband VLBI network, which includes three 13-m radio telescopes located in Jilin, Sanya and Kashi, and a data analysis center in Xi'an. Each station is equipped with a highly stable hydrogen atomic clock and a self-developed VLBI backend, and is co-located with two GPS receivers. This VGOS-like VLBI network may play an important role in improving the Chinese broadband VLBI technology and making valuable contributions to domestic VLBI measurements of UT1. In this paper, we introduce the specifications of this VLBI network, and present the UT1 measurements at C-band conducted in 2018 using the Jilin-Kashi baseline of this network. The comparisons between our UT1 estimates and those provided by IERS suggest that the NTSC VLBI network is capable to determine UT1 accurate at the level of 58.8 microseconds.

The Large Optical/infrared Telescope (LOT) is a ground-based 12m diameter optical/infrared telescope which is proposed to be built in the western part of China in the next decade. Based on satellite remote sensing data, along with geographical, logistical and political considerations, three candidate sites were chosen for ground-based astronomical performance monitoring. These sites include: Ali in Tibet, Daocheng in Sichuan, and Muztagh Ata in Xinjiang. Up until now, all three sites have continuously collected data for two years. In this paper, we will introduce this site testing campaign, and present its monitoring results obtained during the period between March 2017 and March 2019.

The proposed mission "Space Atomic Gravity Explorer" (SAGE) has the scientific objective to investigate gravitational waves, dark matter, and other fundamental aspects of gravity as well as the connection between gravitational physics and quantum physics using new quantum sensors, namely, optical atomic clocks and atom interferometers based on ultracold strontium atoms.

The internal rotational dynamics of massive stars are poorly understood. If angular momentum (AM) transport between the core and the envelope is inefficient, the large core AM upon core-collapse will produce rapidly rotating neutron stars (NSs). However, observations of low-mass stars suggest an efficient AM transport mechanism is at work, which could drastically reduce NS spin rates. Here we study the effects of the baroclinic instability and the magnetic Tayler instability in differentially rotating radiative zones. Although the baroclinic instability may occur, the Tayler instability is likely to be more effective for AM transport. We implement Tayler torques as prescribed by Fuller et al. 2019 into models of massive stars, finding they remove the vast majority of the core's AM as it contracts between the main sequence and helium-burning phases of evolution. If core AM is conserved during core-collapse, we predict natal NS rotation periods of $P_{\rm NS} \approx 50-200 \, {\rm ms}$, suggesting these torques help explain the relatively slow rotation rates of most young NSs, and the rarity of rapidly rotating engine-driven supernovae. Stochastic spin-up via waves just before core-collapse, asymmetric explosions, and various binary evolution scenarios may increase the initial rotation rates of many NSs.

The age of gravitational wave (GW) astronomy has begun, and black hole (BH) mergers detected by LIGO are providing novel constraints on massive star evolution. A major uncertainty in stellar theory is the angular momentum (AM) transport within the star that determines its core rotation rate and the resulting BH's spin. Internal rotation rates of low-mass stars measured from asteroseismology prove that AM transport is efficient, suggesting that massive stellar cores may rotate slower than prior expectations. We investigate AM transport via the magnetic Tayler instability, which can largely explain the rotation rates of low-mass stars and white dwarfs. Implementing an updated AM transport prescription into models of high-mass stars, we compute the spins of their BH remnants. We predict that BHs born from single stars rotate very slowly, with $a \sim 10^{-2}$, regardless of initial rotation rate, possibly explaining the low $\chi_{\rm eff}$ of most BH binaries detected by LIGO thus far. A limited set of binary models suggests slow rotation for many binary scenarios as well, although homogeneous evolution and tidal spin-up of post-common envelope helium stars can create moderate or high BH spins. We make predictions for the values of $\chi_{\rm eff}$ in future LIGO events, and we discuss implications for engine-powered transients.

In this paper we present three new extrasolar planets from the Qatar Exoplanet Survey (QES). Qatar-8b is a hot Saturn, with Mpl = 0.37 Mjup and Rpl = 1.3 Rjup, orbiting a solar-like star every Porb = 3.7 days. Qatar-9b is a hot Jupiter with a mass of Mpl = 1.2 Mjup and a radius of Rpl = 1 Rjup, in a Porb = 1.5 days orbit around a low mass, Mstar = 0.7 Msun, mid-K main-sequence star. Finally, Qatar-10b is a hot, Teq ~ 2000 K, sub-Jupiter mass planet, Mpl = 0.7 Mjup, with a radius of Rpl = 1.54 Rjup and an orbital period of Porb = 1.6 days, placing it on the edge of the sub-Jupiter desert.

Large High Altitude Air Shower Observatory(LHAASO) is a composite cosmic ray observatory consisting of three detector arrays: kilometer square array (KM2A) which includes the electromagnetic detector array and muon detector array, water Cherenkov detector array (WCDA) and wide field of view Cherenkov telescope array (WFCTA). One of the main scientific objectives of LHAASO is to precisely measure the cosmic rays energy spectrum of individual components from 1014 eV to 1018 eV. The hybrid observation will be employed by LHAASO experiment, in which the lateral and longitudinal distributions of the extensive air shower can be observed simultaneously. Thus many kinds of parameters can be used for primary nuclei identification. In this paper, high purity cosmic ray simulation samples of light nuclei component are obtained through Multi-Variable Analysis. The apertures of 1/4 LHAASO array for pure proton and mixed proton and helium (H&He) samples are 900 m2Sr and 1800 m2Sr respectively. A prospect of proton and H&He spectra from 100 TeV to 4 PeV is discussed.

An increasing number of white dwarf stars show atmospheric chemical composition polluted by heavy elements accreted from debris disk material. The existence of such debris disks strongly suggests the presence of one or more planet(s) whose gravitational interaction with rocky planetesimals is responsible for their disruption by tidal effect. The ZZ Ceti pulsator and polluted DAZ white dwarf GD 133 is a good candidate for searching for such a potential planet. We started in 2011 a photometric follow-up of its pulsations. As a result of this work in progress, we used the data gathered from 2011 to 2015 to make an asteroseismological analysis of GD 133, providing the star parameters from a best fit model with $M$/$M_{\odot}$ = 0.630 $\pm$ 0.002, $T_{\rm eff}$ = 12400 K $\pm$ 70 K, log($M_{\rm He}/M$) = -2.00 $\pm$ 0.02, log($M_{\rm H}/M$) = -4.50 $\pm$ 0.02 and determining a rotation period of $\approx$ 7 days.

A dense neutrino medium can support flavor oscillation waves which are coherent among different momentum modes of the neutrinos. The dispersion relation (DR) branches of such a wave with complex frequencies and/or wave numbers can lead to the exponential growth of the wave amplitude which in turn will engender a collective flavor transformation in the neutrino medium. In this work we propose that the complex DR branches of the neutrino oscillation wave should be bound by the critical points of the DR. We demonstrate how this theory can be applied to the neutrino medium with an (approximate) axial symmetry about the propagation direction of the neutrino oscillation wave. We also show how the flavor instabilities in this medium can be identified by tracing the critical points of the DR as the electron lepton number distribution of the neutrino medium is changed continuously.

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

We present the discovery of Qatar-7b --- a very hot and inflated giant gas planet orbiting close its parent star. The host star is a relatively massive main sequence F-star with mass and radius Mstar = 1.41 +/- 0.03 Msun and Rstar = 1.56 +/- 0.02 Rsun, respectively, at a distance d = 726 +/- 26 pc, and an estimated age ~1 Gyr. With its orbital period of P = 2.032 days the planet is located less than 5 stellar radii from its host star and is heated to a high temperature Teq ~ 2100 K. From a global solution to the available photometric and radial velocity observations, we calculate the mass and radius of the planet to be Mpl = 1.88 +/- 0.25 Mjup and Rpl = 1.70 +/- 0.03 Rjup, respectively. The planet radius and equilibrium temperature put Qatar-7b in the top 6% of the hottest and largest known exoplanets. With its large radius and high temperature Qatar-7b is a valuable addition to the short list of targets that offer the best opportunity for studying their atmospheres through transmission spectroscopy.

We consider the flavor transformation of neutrinos through oscillatory matter profiles. We show that the neutrino oscillation Hamiltonian in this case describes a Rabi system with an infinite number of Rabi modes. We further show that, in a given physics problem, the majority of the Rabi modes have too small amplitudes to be relevant. We also go beyond the rotating wave approximation and derive the relative detuning of the Rabi resonance when multiple Rabi modes with small amplitudes are present. We provide an explicit criterion of whether an off-resonance Rabi mode can affect the parametric flavor transformation of the neutrino.

This paper presents the Charge-Coupled Device time-series photometric observations of the open cluster NGC 1582 and its surrounding field with Johnson B, V and R filters by using Nanshan 1m telescope of Xinjiang Astronomical Observatory. 19 variable stars and 3 variable candidates were detected in a 45 {\times} 48.75 field around the cluster. 12 of the variable stars are newly-discovered variable objects. The physical properties, classifications, and membership of the 22 objects are studied through their light curves, their positions on the color-magnitude diagram, and with the archival data from the Naval Observatory Merged Astrometric Dataset. Among these objects, 5 are eclipsing binary systems, 6 are pulsating variable stars including one known Delta Scuti star and one newly-discovered RR Lyrae star. The distance of the RR Lyrae star is estimated to be 7.9 {\pm} 0.3 kpc, indicates that the star locates far behind the cluster. 4 variable stars are probable members of the cluster, and 13 of the 22 objects are confirmed to be field stars.

In this paper, the high-precision light curves of $Kepler$ target KIC 5950759 are analyzed. The Fourier analysis of the long cadence light curve reveals 3 independent frequencies. Two of them are main pulsation modes: F0 = 14.221373(21) $\rm{d^{-1}}$ and F1 = 18.337249(44) $\rm{d^{-1}}$. The third independent frequency $f_m$ = 0.3193 d$^{-1}$ is found in long cadence data with a signal-to-noise ratio of 6.2. A weak modulation of $f_m$ to F0 and F1 modes (triplet structures centred on F0 and F1) are detected both in long and short cadence data. This is the first detection of the modulation effect in a double-mode HADS star. The most possible cause of the modulation effect in the light curves is amplitude modulation with the star's rotation frequency of 0.3193 d$^{-1}$. The preliminary analysis suggests that KIC 5950759 is in the bottom of the HADS instability strip and likely situated in the main sequence. Spectroscopic observations are necessary to verify the true nature of the modulation terms.

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

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

We present a study on the pulsations of the high-amplitude $\delta$ Scuti star YZ Boo based on photometric observations in Johnson $V$ and $R$ bands with both the Nanshan 1-m telescope of Xinjiang Astronomical Observatory (XAO) and the Xinglong 85-cm telescope of National Astronomical Observatories, Chinese Academy of Sciences (NAOC). The Fourier analysis of the light curves reveals the fundamental radial mode and its five harmonics, with the fourth and the fifth newly detected. 39 new times of light maximum are determined from the light curves, and combined with those in the literature, we construct the $O - C$ diagram and derive a new ephemeris and the determination of a new value of the updated period 0.104091579(2). Besides, the $O - C$ diagram reveals an increasing period change rate of YZ Boo. Theoretical models are calculated and constrained with the observationally determined parameters of YZ Boo. The mass and age of YZ Boo are hence derived as $M$ = 1.61 $\pm$ 0.05 $M_\odot$ and $age$ = (1.44 $\pm$ 0.14) $\times$ 10$^9$ yr, respectively. With both the frequency of the fundamental radial mode and the period change rate, YZ Boo is located at the post-main-sequence stage.

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.

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.

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

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

Kink instability is a possible mechanism for solar filament eruption. However, the twist of a solar filament is very difficult to directly measure from observation. In this paper, we carried out the measurement of the twist of a solar filament by analyzing its leg rotation. An inverse S-shaped filament in active region NOAA 11485 was observed by the Atmospheric Imaging Assembly (AIA) of Solar Dynamics Observatory (SDO) on 2012 May 22. During its eruption, the leg of the filament exhibited a significant rotation motion. The 304 \AA\ images were used to uncurl along the circles, whose centers are the axis of the filament's leg. The result shows that the leg of the filament rotated up to about 510 degrees (about 2.83$\pi$) around the axis of the filament within twenty-three minutes. The maximal rotation speed reached 100 degrees per minute (about 379.9 km/s at radius 18$^\prime$$^\prime$), which is the fastest rotation speed that has been reported. We also calculated the decay index along the polarity inversion line in this active region and found that the decline of the overlying field with height is not so fast enough to trigger the torus instability. According to the condition of kink instability, it is indicating that the kink instability is the trigger mechanism for the solar filament eruption.

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.

With the data from the Atmospheric Imaging Assembly (AIA) and the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO), we present a magnetic interaction between an isolated coronal hole (CH) and an emerging active region (AR). The AR emerged nearby the CH and interacted with it. Bright loops constantly formed between them, which led to a continuous retreat of the CH boundaries (CHBs). Meanwhile, two coronal dimmings respectively appeared at the negative polarity of the AR and the east boundary of the bright loops, and the AR was partly disturbed. Loop eruptions followed by a flare occurred in the AR. The interaction was also accompanied by many jets and an arc-shaped brightening that appeared to be observational signatures of magnetic reconnection at the CHBs. By comparing the observations with the derived coronal magnetic configuration, it is suggested that the interaction between the CH and the AR excellently fitted in with the model of interchange reconnection. It appears that our observations provide obvious evidences for interchange reconnection.

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 Chinese Area Positioning System (CAPS), a navigation system based on GEO communication satellites, was developed in 2002 by astronomers at Chinese Academy of Sciences. Extensive positioning experiments of CAPS have been performed since 2005. On the basis of CAPS, this paper studies the principle of navigation constellation composed of Slightly Inclined Geostationary Orbit (SIGSO) communication satellites. SIGSO satellites are derived from end-of-life Geostationary Orbit (GEO) satellites under inclined orbit operation. Considering the abundant frequency resources of SIGSO satellites, multi-frequency observations could be conducted to enhance the precision of pseudorange measurements and ameliorate the positioning performence. The constellation composed of two GEO satellites and four SIGSO satellites with inclination of 5 degrees can provide the most territory of China with 24-hour maximum PDOP less than 42. With synthetic utilization of the truncated precise (TP) code and physical augmentation factor in four frequencies, navigation system with this constellation is expected to obtain comparable positioning performance with that of coarse acquisition code of GPS. When the new approach of code-carrier phase combinations is adopted, the system has potential to possess commensurate accuracy of precise code in GPS. Additionally, the copious frequency resources can also be used to develop new anti-interference techniques and integrate navigation and communication.

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.

Previous studies have found that the width of gamma-ray burst (GRB) pulse is energy dependent and that it decreases as a power-law function with increasing photon energy. In this work we have investigated the relation between the energy dependence of pulse and the so-called Band spectrum by using a sample including 51 well-separated fast rise and exponential decay long-duration GRB pulses observed by BATSE (Burst and Transient Source Experiment on the Compton Gamma Ray Observatory). We first decompose these pulses into rise, and decay phases and find the rise widths, and the decay widths also behavior as a power-law function with photon energy. Then we investigate statistically the relations between the three power-law indices of the rise, decay and total width of pulse (denoted as $\delta_r$, $\delta_d$ and $\delta_w$, respectively) and the three Band spectral parameters, high-energy index ($\alpha$), low-energy index ($\beta$) and peak energy ($E_p$). It is found that (1)$\alpha$ is strongly correlated with $\delta_w$ and $\delta_d$ but seems uncorrelated with $\delta_r$; (2)$\beta$ is weakly correlated with the three power-law indices and (3)$E_p$ does not show evident correlations with the three power-law indices. We further investigate the origin of $\delta_d-\alpha$ and $\delta_w-\alpha$. We show that the curvature effect and the intrinsic Band spectrum could naturally lead to the energy dependence of GRB pulse width and also the $\delta_d-\alpha$ and $\delta_w-\alpha$ correlations. Our results would hold so long as the shell emitting gamma rays has a curve surface and the intrinsic spectrum is a Band spectrum or broken power law. The strong $\delta_d-\alpha$ correlation and inapparent correlations between $\delta_r$ and three Band spectral parameters also suggest that the rise and decay phases of GRB pulses have different origins.

Time variable cosmological constant (TVCC) of holographic origin with interaction in Brans-Dicke theory is discussed in this paper. We investigate some characters for this model, and show the evolutions of deceleration parameter and equation of state (EOS) for dark energy. It is shown that in this scenario an accelerating universe can be obtained and the evolution of EOS for dark energy can cross over the boundary of phantom divide. In addition, a geometrical diagnostic method, jerk parameter is applied to this model to distinguish it with cosmological constant.

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.

A wide field of view Cherenkov/fluorescence telescope array is one of the main components of the Large High Altitude Air Shower Observatory project. To serve as Cherenkov and fluorescence detectors, a flexible and mobile design is adopted for easy reconfiguring of the telescope array. Two prototype telescopes have been constructed and successfully run at the site of the ARGO-YBJ experiment in Tibet. The features and performance of the telescopes are presented.

Assuming an intrinsic `Band' shape spectrum and an intrinsic energy-independent emission profile we have investigated the connection between the evolution of the rest-frame spectral parameters and the spectral lags measured in gamma-ray burst (GRB) pulses by using a pulse model. We first focus our attention on the evolution of the peak energy, $E_{0,p}$, and neglect the effect of the curvature effect. It is found that the evolution of $E_{0,p}$ alone can produce the observed lags. When $E_{0,p}$ varies from hard to soft only the positive lags can be observed. The negative lags would occur in the case of $E_{0,p}$ varying from soft to hard. When the evolution of $E_{0,p}$ and the low-energy spectral index $\alpha_{0}$ varying from soft to hard then to soft we can find the aforesaid two sorts of lags. We then examine the combined case of the spectral evolution and the curvature effect of fireball and find the observed spectral lags would increase. A sample including 15 single pulses whose spectral evolution follows hard to soft has been investigated. All the lags of these pulses are positive, which is in good agreement with our theoretical predictions. Our analysis shows that only the intrinsic spectral evolution can produce the spectral lags and the observed lags should be contributed by the intrinsic spectral evolution and the curvature effect. But it is still unclear what cause the spectral evolution.

The Large High Altitude Air Shower Observatory project is proposed to study high energy gamma ray astronomy ( 40 GeV-1 PeV ) and cosmic ray physics ( 20 TeV-1 EeV ). The wide field of view Cherenkov telescope array, as a component of the LHAASO project, will be used to study energy spectrum and compositions of cosmic ray by measuring the total Cherenkov light generated by air showers and shower maximum depth. Two prototype telescopes have been in operation since 2008. The pointing accuracy of each telescope is crucial to the direction reconstruction of the primary particles. On the other hand the primary energy reconstruction relies on the shape of the Cherenkov image on the camera and the unrecorded photons due to the imperfect connections between photomultiplier tubes. UV bright stars are used as point-like objects to calibrate the pointing and to study the optical properties of the camera, the spot size and the fractions of unrecorded photons in the insensitive areas of the camera.

In this paper we have analyzed the temporal and spectral behavior of 52 Fast Rise and Exponential Decay (FRED) pulses in 48 long-duration gamma-ray bursts (GRBs) observed by the CGRO/BATSE, using a pulse model with two shape parameters and the Band model with three shape parameters, respectively. It is found that these FRED pulses are distinguished both temporally and spectrally from those in long-lag pulses. Different from these long-lag pulses only one parameter pair indicates an evident correlation among the five parameters, which suggests that at least $\sim$4 parameters are needed to model burst temporal and spectral behavior. In addition, our studies reveal that these FRED pulses have correlated properties: (i) long-duration pulses have harder spectra and are less luminous than short-duration pulses; (ii) the more asymmetric the pulses are the steeper the evolutionary curves of the peak energy ($E_{p}$) in the $\nu f_{\nu}$ spectrum within pulse decay phase are. Our statistical results give some constrains on the current GRB models.

In this paper, we revisit gamma-ray--emitting region for 10 GeV--1 TeV gamma rays from 3C 279 through studying the photon-photon absorption optical depth due to the diffuse radiation of the broad-line region (BLR) and the extragalactic background light (EBL). Based on the power-law spectrum detected by MAGIC, the preabsorbed spectra are inferred by correcting the photon-photon absorption on the diffuse photons of the BLR (internal absorption) and the EBL (external absorption). Position of gamma-ray emitting region $R_{\rm{\gamma}}$ determines the relative contributions of this two diffuse radiation to the total absorption. Our results indicate that $R_{\rm{\gamma}}$ may be within the BLR shell for 3C 279, likely closer to the inner radius, which is consistent with our previous results. This is neither consistent with the suggestions of B\"ottcher et al. (2008b), that VHE gamma-ray emission is produced far outside the BLR, nor with the assumptions of Tavecchio & Mazin (2008), that VHE gamma-ray--emitting region is inside the BLR cavity. $R_{\rm{\gamma}}$ is a key physical quantity that could set some constraints on emission mechanisms that produce the VHE gamma rays from 3C 279. Observations of $\it Fermi$-LAT, MAGIC, HESS, and VERITAS in the near future could give more constraints on the position of gamma-ray emitting region relative to the BLR.

Employing two samples containing of 56 and 59 well-separated FRED (fast rise and exponential decay) gamma-ray burst (GRB) pulses whose spectra are fitted by the Band spectrum and Compton model, respectively, we have investigated the evolutionary slope of $E_{p}$ (where $E_{p}$ is the peak energy in the $\nu F\nu$ spectrum) with time during the pulse decay phase. The bursts in the samples were observed by the Burst and Transient Source Experiment (BATSE) on the Compton Gamma-Ray Observatory. We first test the $E_{p}$ evolutionary slope during the pulse decay phase predicted by Lu et al. (2007) based on the model of highly symmetric expanding fireballs in which the curvature effect of the expanding fireball surface is the key factor concerned. It is found that the evolutionary slopes are normally distributed for both samples and concentrated around the values of 0.73 and 0.76 for Band and Compton model, respectively, which is in good agreement with the theoretical expectation of Lu et al. (2007). However, the inconsistence with their results is that the intrinsic spectra of most of bursts may bear the Comptonized or thermal synchrotron spectrum, rather than the Band spectrum. The relationships between the evolutionary slope and the spectral parameters are also checked. We show the slope is correlated with $E_{p}$ of time-integrated spectra as well as the photon flux but anticorrelated with the lower energy index $\alpha$. In addition, a correlation between the slope and the intrinsic $E_{p}$ derived by using the pseudo-redshift is also identified. The mechanisms of these correlations are unclear currently and the theoretical interpretations are required.

Employing a sample presented by Kaneko et al. (2006) and Kocevski et al. (2003), we select 42 individual tracking pulses (here we defined tracking as the cases in which the hardness follows the same pattern as the flux or count rate time profile) within 36 Gamma-ray Bursts (GRBs) containing 527 time-resolved spectra and investigate the spectral hardness, $E_{peak}$ (where $E_{peak}$ is the maximum of the $\nu F_{\nu}$ spectrum), evolutionary characteristics. The evolution of these pulses follow soft-to-hard-to-soft (the phase of soft-to-hard and hard-to-soft are denoted by rise phase and decay phase, respectively) with time. It is found that the overall characteristics of $E_{peak}$ of our selected sample are: 1) the $E_{peak}$ evolution in the rise phase always start on the high state (the values of $E_{peak}$ are always higher than 50 keV); 2) the spectra of rise phase clearly start at higher energy (the median of $E_{peak}$ are about 300 keV), whereas the spectra of decay phase end at much lower energy (the median of $E_{peak}$ are about 200 keV); 3) the spectra of rise phase are harder than that of the decay phase and the duration of rise phase are much shorter than that of decay phase as well. In other words, for a complete pulse the initial $E_{peak}$ is higher than the final $E_{peak}$ and the duration of initial phase (rise phase) are much shorter than the final phase (decay phase). This results are in good agreement with the predictions of Lu et al. (2007) and current popular view on the production of GRBs. We argue that the spectral evolution of tracking pulses may be relate to both of kinematic and dynamic process even if we currently can not provide further evidences to distinguish which one is dominant. Moreover, our statistical results give some witnesses to constrain the current GRB model.

In this paper, we study the photon-photon pair production optical depth for 10 GeV--1 TeV gamma rays from 3C 279 due to the diffuse radiation of broad-line region (BLR). Assuming a power-law spectrum of $E_{\gamma}^{-a_2}$ for the photon intensity of very high energy (VHE) gamma rays, $a_1 \gtrsim 405$ and $a_2\gtrsim 6.4$ are inferred by the integrated photon fluxes measured by MAGIC and HESS. Based on this power-law spectrum, the pre-absorbed spectra are inferred by correcting the photon-photon absorption on the diffuse photons of the BLR (internal absorption) and the extragalactic background light (external absorption). Position of gamma-ray emitting region $R_{\rm{\gamma}}$ determines the relative contributions of this two diffuse radiation to the total absorption for 10 GeV--1 TeV gamma rays. The internal absorption could make spectral shape of gamma rays more complex than only corrected for the external absorption, and could lead to the formation of arbitrary softening and hardening gamma-ray spectra. It should be necessary for the internal absorption to be considered in studying 10 GeV--1 TeV gamma rays from powerful blazars. The energies of annihilated gamma-ray photons due to the internal absorption are likely to be mainly reradiated around GeV. Our results indicate that $R_{\rm{\gamma}}$ may be between the inner and outer radii of the BLR for 3C 279. This implies for powerful blazars that $R_{\rm{\gamma}}$ might be neither inside the BLR cavity nor outside the BLR, but be within the BLR shell. Observations of $\it GLAST$, MAGIC, HESS, and VERITAS in the near future could give more constraints on the position of the gamma-ray emitting region relative to the BLR.

In this paper, standard accretion disk models of AGNs are tested using light curves of 26 objects well observed for reverberation mapping. Time scales of variations are estimated by the most common definition of the variability time scale and the zero-crossing time of the autocorrelation function of the optical light curves for each source. The measured time scales of variations by the two methods are consistent with each other. If the typical value of the viscosity parameter $\alpha \sim 0.1$ is adopted, the measured optical variability time scales are most close to the thermal time scales of the standard disks. If $\alpha$ is allowed to range from $\sim 0.03$ to $\sim 0.2$, the measured time scales are consistent with the thermal time scales of the standard disks. There is a linear relation between the measured variability time scales and black hole masses, and this linear relation is qualitatively consistent with expectation of the standard accretion disk models. The time lags measured by the ZDCF between different bands are on the order of days. The measured time lags of NGC 4151 and NGC 7469 are marginally consistent with the time lags estimated in the case of continuum thermal reprocessing for the standard accretion disk models. However, the measured time lags of NGC 5548 and Fairall 9 are unlikely to be the case of continuum thermal reprocessing. Our results are unlikely to be inconsistent with or are likely to be conditionally in favor of the standard accretion disk models of AGNs.