We describe system verification tests and early science results from the pulsar processor (PTUSE) developed for the newly-commissioned 64-dish SARAO MeerKAT radio telescope in South Africa. MeerKAT is a high-gain (~2.8 K/Jy) low-system temperature (~18 K at 20cm) radio array that currently operates from 580-1670 MHz and can produce tied-array beams suitable for pulsar observations. This paper presents results from the MeerTime Large Survey Project and commissioning tests with PTUSE. Highlights include observations of the double pulsar J0737-3039A, pulse profiles from 34 millisecond pulsars from a single 2.5h observation of the Globular cluster Terzan 5, the rotation measure of Ter5O, a 420-sigma giant pulse from the Large Magellanic Cloud pulsar PSR J0540-6919, and nulling identified in the slow pulsar PSR J0633-2015. One of the key design specifications for MeerKAT was absolute timing errors of less than 5 ns using their novel precise time system. Our timing of two bright millisecond pulsars confirm that MeerKAT delivers exceptional timing. PSR J2241-5236 exhibits a jitter limit of <4 ns per hour whilst timing of PSR J1909-3744 over almost 11 months yields an rms residual of 66 ns with only 4 min integrations. Our results confirm that the MeerKAT is an exceptional pulsar telescope. The array can be split into four separate sub-arrays to time over 1000 pulsars per day and the future deployment of S-band (1750-3500 MHz) receivers will further enhance its capabilities.

Conventional cosmic-ray propagation models usually assume an isotropic diffusion coefficient to account for the random deflection of cosmic rays by the turbulent interstellar magnetic field. Such a picture is very successful in explaining many observational phenomena related to the propagation of Galactic cosmic rays, such as broken power-law energy spectra, secondary-to-primary ratios, etc. However, the isotropic diffusion presupposition is facing severe challenges from recent observations. In particular, such observations on the large-scale anisotropy of TeV cosmic rays show that the dipole direction differs from the prediction of the conventional model. One possible reason is that the large-scale regular magnetic field, which leads to an anisotropic diffusion of cosmic rays, has not been included in the model provided by the public numerical packages. In this work, we propose two numerical schemes to solve the $3$-dimensional anisotropic transport equation: the pseudo source method and Hundsdorfer-Verwer scheme. Both methods are verified by reproducing the measured B/C and proton spectrum and the radial variation of spectral index expected by former 2D simulation. As a demonstration of the prediction capability, dipole anisotropy is also calculated by a toy simulation with a rough magnetic field.

The Water Cherenkov Detector Array (WCDA) is a major component of the Large High Altitude Air Shower Array Observatory (LHAASO), a new generation cosmic-ray experiment with unprecedented sensitivity, currently under construction. The WCDA is aimed at the study of TeV $\gamma$-rays. In order to evaluate the prospects of searching for TeV $\gamma$-ray sources with the WCDA, we present in this paper a projection for the one-year sensitivity of the WCDA to TeV $\gamma$-ray sources from TeVCat (footnote: http://tevcat.uchicago.edu) using an all-sky approach. Out of 128 TeVCat sources observable to the WCDA up to a zenith angle of $45^\circ$, we estimate that 42 would be detectable for one year of observations at a median energy of 1 TeV. Most of them are Galactic sources, and the extragalactic sources are Active Galactic Nuclei (AGN).

In the standard cosmological model, dark matter drives the structure formation and constructs potential wells within which galaxies may form. The baryon fraction in dark halos can reach the universal value (15.7%) in massive clusters and decreases rapidly as the mass of the system decreases. The formation of dwarf galaxies is sensitive both to baryonic processes and the properties of dark matter owing to the shallow potential wells in which they form. In dwarf galaxies in the Local Group, dark matter dominates the mass content even within their optical-light half-radii (r_e ~ 1 kpc). However, recently it has been argued that not all dwarf galaxies are dominated by dark matter. Here we report 19 dwarf galaxies that could consist mainly of baryons up to radii well beyond r_e, at which point they are expected to be dominated by dark matter. Of these, 14 are isolated dwarf galaxies, free from the influence of nearby bright galaxies and high dense environments. This result provides observational evidence that could challenge the formation theory of low-mass galaxies within the framework of standard cosmology. Further observations, in particular deep imaging and spatially-resolved kinematics, are needed to constrain the baryon fraction better in such galaxies.

In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which are regularly observed by these groups. A basic noise analysis, including the processes which are both correlated and uncorrelated in time, provides noise models and timing ephemerides for the pulsars. We find that the timing precisions of pulsars are generally improved compared to the previous data release, mainly due to the addition of new data in the combination. The main purpose of this work is to create the most up-to-date IPTA data release. These data are publicly available for searches for low-frequency gravitational waves and other pulsar science.

Recent timing observation provides an intriguing result for the braking index of the X-ray pulsar PSR J1640-4631, which has a measured braking index $n=3.15\pm0.03$. The decrease of the inclination angle between between the spin axis and the magnetic axis can be responsible for such a high braking index. However, the physical mechanisms causing the change of the magnetic inclination angle have not been fully understood. In this Letter, we apply a two-dipole model given by Hamil et al (2016) to explain the decrease of the magnetic inclination angle of PSR J1640-4631. The rotation effect of a charged sphere and the magnetization of ferromagnetically ordered material produce magnetic moments $M_{1}$ and $M_{2}$, respectively. There exist a minimum of the potential energy for the magnetic moment $M_{2}$ in the magnetic field of $M_{1}$, hence the $M_{2}$ will freely rotate around the minimum energy position (i. e equilibrium position), similar to a simple pendulum. Our calculation indicate the magnetic moment $M_{2}$ would evolve towards alignment with the spin axis for PSR J1640-4631, and cause the magnetic inclination angle to decrease. The single peak in the pulse profile favors a relatively low change rate of the magnetic inclination angle.

An unexpected strong geomagnetic storm occurred on 2018 August 26, which was caused by a slow coronal mass ejection (CME) from a gradual eruption of a large quiet-region filament. We investigate the eruption and propagation characteristics of this CME in relation to the strong geomagnetic storm with remote sensing and in situ observations. Coronal magnetic fields around the filament are extrapolated and compared with EUV observations. We determine the propagation direction and tilt angle of the CME flux rope near the Sun using a graduated cylindrical shell (GCS) model and the Sun-to-Earth kinematics of the CME with wide-angle imaging observations from STEREO A. We reconstruct the flux-rope structure using a Grad-Shafranov technique based on the in situ measurements at the Earth and compare it with those from solar observations and the GCS results. Our conclusions are as follows: (1) the eruption of the filament was unusually slow and occurred in the regions with relatively low critical heights of the coronal field decay index; (2) the axis of the CME flux rope rotated in the corona as well as in interplanetary space, which tended to be aligned with the local heliospheric current sheet; (3) the CME was bracketed between slow and fast solar winds, which enhanced the magnetic field inside the CME at 1 AU; (4) the geomagnetic storm was caused by the enhanced magnetic field and a southward orientation of the flux rope at 1 AU from the rotation of the flux rope.

In order to have a comprehensive view of the propagation and evolution of coronal mass ejections (CMEs) from the Sun to deep interplanetary space beyond 1 au, we carry out a kinematic analysis of 7 CMEs in solar cycle 23. The events are required to have coordinated coronagraph observations, interplanetary type II radio bursts, and multi-point in-situ measurements at the Earth and Ulysses. A graduated cylindrical shell model, an analytical model without free parameters and a magnetohydrodynamic model are used to derive CME kinematics near the Sun, to quantify the CME/shock propagation in the Sun-Earth space, and to connect in-situ signatures at the Earth and Ulysses, respectively. We find that each of the 7 CME-driven shocks experienced a major deceleration before reaching 1 au and thereafter propagated with a gradual deceleration from the Earth to larger distances. The resulting CME/shock propagation profile for each case is roughly consistent with all the data, which verifies the usefulness of the simple analytical model for CME/shock propagation in the heliosphere. The statistical analysis of CME kinematics indicates a tendency that the faster the CME, the larger the deceleration, and the shorter the deceleration time period within 1 au. For several of these events, the associated geomagnetic storms were mainly caused by the southward magnetic fields in the sheath region. In particular, the interaction between a CME-driven shock and a preceding ejecta significantly enhanced the preexisting southward magnetic fields and gave rise to a severe complex geomagnetic storm.

Context. Solar gravity modes (g modes) are buoyancy waves trapped in the solar radiative zone that have been very difficult to detect at the surface. Solar g modes would complement solar pressure modes (p modes) in probing the central regions of the Sun, for example the core rotation rate. Aims. A detection of g modes using changes in the large frequency separation of p modes has recently been reported. However, it is unclear how p and g modes interact. The aim of this study is to evaluate to what extent g modes can perturb the frequencies of p modes. Methods. We computed the first-order perturbation to global p-mode frequencies due to a flow field and perturbations to solar structure caused by a g mode. We focused on long-period g modes and assumed that the g-mode perturbations are constant in time. The surface amplitude of g modes is assumed to be $1$ mm s$^{-1}$, which is close to the observational limit. Results. Gravity modes do perturb p-mode frequencies to first order if the harmonic degree of the g mode is even and if its azimuthal order is zero. The effect is extremely small. For dipole and quadrupole p modes, all frequency shifts are smaller than $0.1$ nHz, or $2\times10^{-8}$ in relative numbers. This is because the relative perturbation to solar structure quantities caused by a g mode of realistic amplitude is of the order of $10^{-6}$ to $10^{-5}$. We find that structural changes dominate over advection. Surprisingly, the interaction of g and p modes takes place to a large part near the surface, where p modes spend most of their propagation times and g modes generate the largest relative changes to solar structure. Conclusions. It appears to be impossible to detect g modes solely through their signature in p-mode frequency shifts. Whether g modes leave a detectable signature in p-mode travel times under a given observational setup remains an open question.

The Giant Radio Array for Neutrino Detection (GRAND) is a planned large-scale observatory of ultra-high-energy (UHE) cosmic particles, with energies exceeding 10^8 GeV. Its goal is to solve the long-standing mystery of the origin of UHE cosmic rays. To do this, GRAND will detect an unprecedented number of UHE cosmic rays and search for the undiscovered UHE neutrinos and gamma rays associated to them with unmatched sensitivity. GRAND will use large arrays of antennas to detect the radio emission coming from extensive air showers initiated by UHE particles in the atmosphere. Its design is modular: 20 separate, independent sub-arrays, each of 10 000 radio antennas deployed over 10 000 km^2. A staged construction plan will validate key detection techniques while achieving important science goals early. Here we present the science goals, detection strategy, preliminary design, performance goals, and construction plans for GRAND.

We investigate a global extreme-ultraviolet (EUV) wave associated with a coronal mass ejection (CME)-driven shock on 2017 September 10. The EUV wave is transmitted by north- and south-polar coronal holes (CHs), which is observed by the Solar Dynamics Observatory (SDO) and Solar Terrestrial Relations Observatory A (STEREO-A) from opposite sides of the Sun. We obtain key findings on how the EUV wave interacts with multiple coronal structures, and on its connection with the CME-driven shock: (1) the transmitted EUV wave is still connected with the shock that is incurvated to the Sun, after the shock has reached the opposite side of the eruption; (2) the south CH transmitted EUV wave is accelerated inside an on-disk, low-density region with closed magnetic fields, which implies that an EUV wave can be accelerated in both open and closed magnetic field regions; (3) part of the primary EUV wavefront turns around a bright point (BP) with a bipolar magnetic structure when it approaches a dim, low-density filament channel near the BP; (4) the primary EUV wave is diffused and apparently halted near the boundaries of remote active regions (ARs) that are far from the eruption, and no obvious AR related secondary waves are detected; (5) the EUV wave extends to an unprecedented scale of ~360{\deg} in latitudes, which is attributed to the polar CH transmission. These results provide insights into the effects of coronal density and magnetic field distributions on the evolution of an EUV wave, and into the connection between the EUV wave and the associated CME-driven shock.

We report on the highest energy photons from the Crab Nebula observed by the Tibet air shower array with the underground water-Cherenkov-type muon detector array. Based on the criterion of muon number measured in an air shower, we successfully suppress 99.92% of the cosmic-ray background events with energies $E>100$ TeV. As a result, we observed 24 photon-like events with $E>100$ TeV against 5.5 background events, which corresponds to 5.6$\sigma$ statistical significance. This is the first detection of photons with $E>100$ TeV from an astrophysical source.

We report on the highest energy photons from the Crab Nebula observed by the Tibet air shower array with the underground water-Cherenkov-type muon detector array. Based on the criterion of muon number measured in an air shower, we successfully suppress 99.92% of the cosmic-ray background events with energies $E>100$ TeV. As a result, we observed 24 photon-like events with $E>100$ TeV against 5.5 background events, which corresponds to 5.6$\sigma$ statistical significance. This is the first detection of photons with $E>100$ TeV from an astrophysical source.

TREND50 is a radio detection setup of 50 self-triggered antennas working in the 50-100MHz frequency range and deployed in a radio-quiet valley of the Tianshan mountains (China). TREND50 achieved its goal: the autonomous radiodetection and identification of air showers. Thanks to a dedicated offine selection algorithm, 564 air shower candidates were indeed selected out of $7\cdot10^8$ transient radio signals recorded during the 314 live days of data taken during the first two years of operation of this setup (2011 and 2012). This event rate, as well as the distribution of the candidate directions of arrival, is consistent with what is expected from cosmic-ray-induced air showers according to simulations, assuming an additional $\sim$20% contamination of the final sample by background events. This result is obtained at the cost of a reduced air shower detection efficiency, estimated to be $\sim$3%. This low efficiency is mostly due to the large amount of dead time of the setup. This result paves the way for the GRANDProto35 experiment, the first stage of the GRAND project.

It is paramount from both scientific and societal perspectives to understand the generation of extreme space weather. We discuss the formation of solar superstorms based on a comparative study of the 2012 July 23 and 2017 July 23 eruptions. The first one is Carrington-class, and the second could rival the 1989 March event that caused the most intense geomagnetic storm of the space age. Observations of these events in the historically weak solar cycle 24 indicate that a solar superstorm can occur in any solar cycle and at any phase of the cycle. Recurrent patterns are identified in both cases, including the long-lived eruptive nature of the active region, a complex event composed of successive eruptions from the same active region, and in-transit interaction between the successive eruptions resulting in exceptionally strong ejecta magnetic fields at 1 AU. Each case also shows unique characteristics. Preconditioning of the upstream solar wind leading to unusually high solar wind speeds at 1 AU is observed in the first case whereas not in the latter. This may suggest that the concept of "preconditioning" appears to be necessary for making a Carrington-class storm. We find a considerable deflection by nearby coronal holes in the second case but not in the first. On the basis of these results, we propose a hypothesis for further investigation that superstorms are "perfect storms" in nature, i.e., a combination of circumstances that results in an event of unusual magnitude. Historical records of some extreme events seem to support our hypothesis.

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

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

(Abridged) This is the Maunakea Spectroscopic Explorer 2018 book. It is intended as a concise reference guide to all aspects of the scientific and technical design of MSE, for the international astronomy and engineering communities, and related agencies. The current version is a status report of MSE's science goals and their practical implementation, following the System Conceptual Design Review, held in January 2018. MSE is a planned 10-m class, wide-field, optical and near-infrared facility, designed to enable transformative science, while filling a critical missing gap in the emerging international network of large-scale astronomical facilities. MSE is completely dedicated to multi-object spectroscopy of samples of between thousands and millions of astrophysical objects. It will lead the world in this arena, due to its unique design capabilities: it will boast a large (11.25 m) aperture and wide (1.52 sq. degree) field of view; it will have the capabilities to observe at a wide range of spectral resolutions, from R2500 to R40,000, with massive multiplexing (4332 spectra per exposure, with all spectral resolutions available at all times), and an on-target observing efficiency of more than 80%. MSE will unveil the composition and dynamics of the faint Universe and is designed to excel at precision studies of faint astrophysical phenomena. It will also provide critical follow-up for multi-wavelength imaging surveys, such as those of the Large Synoptic Survey Telescope, Gaia, Euclid, the Wide Field Infrared Survey Telescope, the Square Kilometre Array, and the Next Generation Very Large Array.

We report the observations of a moderate but relatively intense geo-effective solar eruption on 2015 November 4 from the peripheral diffusive polarities of active region 12443. We use space-borne Solar Dynamics Observatory and ACE observations. EUV images identified helical pattern along a filament channel and we regard this channel as flux-rope structure. Flow velocity derived from tracked magnetograms infers converging motion along the polarity inversion line beneath the filament channel. An associated magnetic cancellation process was detected in the converging region. Further, the pre-eruptive EUV brightening was observed in the converging region, the most intense part of which appeared in the magnetic cancellation region. These observations imply that the converging and cancelling flux probably contributed to the formation of the helical magnetic fields associated with the flux rope. A filament-height estimation method suggests that the middle part of the filament probably lies at a low altitude and was consistent with the initial place of the eruption. A thick current channel associated with the flux rope is also determined. For a expanding thick current channel, the critical height of the decay index for torus instability lies in the range of 37 - 47 Mm. Southward magnetic fields in the sheath and the ejecta induced a geomagnetic storm with a Dst global minimum of ~-90 nT.

We examine the possible influence of Earth-directed coronal mass ejections (ECMEs) on the Sun's shadow in the 3~TeV cosmic-ray intensity observed by the Tibet-III air shower (AS) array. We confirm a clear solar-cycle variation of the intensity deficit in the Sun's shadow during ten years between 2000 and 2009. This solar-cycle variation is overall reproduced by our Monte Carlo (MC) simulations of the Sun's shadow based on the potential field model of the solar magnetic field averaged over each solar rotation period. We find, however, that the magnitude of the observed intensity deficit in the Sun's shadow is significantly less than that predicted by MC simulations, particularly during the period around solar maximum when a significant number of ECMEs is recorded. The $\chi^2$ tests of the agreement between the observations and the MC simulations show that the difference is larger during the periods when the ECMEs occur, and the difference is reduced if the periods of ECMEs are excluded from the analysis. This suggests the first experimental evidence of the ECMEs affecting the Sun's shadow observed in the 3~TeV cosmic-ray intensity.

We examine the possible influence of Earth-directed coronal mass ejections (ECMEs) on the Sun's shadow in the 3~TeV cosmic-ray intensity observed by the Tibet-III air shower (AS) array. We confirm a clear solar-cycle variation of the intensity deficit in the Sun's shadow during ten years between 2000 and 2009. This solar-cycle variation is overall reproduced by our Monte Carlo (MC) simulations of the Sun's shadow based on the potential field model of the solar magnetic field averaged over each solar rotation period. We find, however, that the magnitude of the observed intensity deficit in the Sun's shadow is significantly less than that predicted by MC simulations, particularly during the period around solar maximum when a significant number of ECMEs is recorded. The $\chi^2$ tests of the agreement between the observations and the MC simulations show that the difference is larger during the periods when the ECMEs occur, and the difference is reduced if the periods of ECMEs are excluded from the analysis. This suggests the first experimental evidence of the ECMEs affecting the Sun's shadow observed in the 3~TeV cosmic-ray intensity.

We investigate the kinetic properties of a typical fast-mode shock inside an interplanetary coronal mass ejection (ICME) observed on 1998 August 6 at 1 AU, including particle distributions and wave analysis with the in situ measurements from Wind. Key results are obtained concerning the shock and the shock-ICME interaction at kinetic scales: (1) gyrating ions, which may provide energy dissipation at the shock in addition to wave-particle interactions, are observed around the shock ramp; (2) despite the enhanced proton temperature anisotropy of the shocked plasma, the low plasma beta inside the ICME constrains the shocked plasma under the thresholds of the ion cyclotron and mirror-mode instabilities; (3) whistler heat flux instabilities, which can pitch--angle scatter halo electrons through a cyclotron resonance, are observed around the shock, and can explain the disappearance of bidirectional electrons inside the ICME together with normal betatron acceleration; (4) whistler waves near the shock are likely associated with the whistler heat flux instabilities excited at the shock ramp, which is consistent with the result that the waves may originate from the shock ramp; (5) the whistlers share a similar characteristic with the shocklet whistlers observed by Wilson et al, providing possible evidence that the shock is decaying because of the strong magnetic field inside the ICME.

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.

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.

We investigate how a weak coronal mass ejection (CME) launched on 2016 October 8 without obvious signatures in the low corona produced a relatively intense geomagnetic storm. Remote sensing observations from SDO, STEREO and SOHO and in situ measurements from WIND are employed to track the CME from the Sun to the Earth. Using a graduated cylindrical shell (GCS) model, we estimate the propagation direction and the morphology of the CME near the Sun. CME kinematics are determined from the wide-angle imaging observations of STEREO A and are used to predict the CME arrival time and speed at the Earth. We compare ENLIL MHD simulation results with in situ measurements to illustrate the background solar wind where the CME was propagating. We also apply a Grad--Shafranov technique to reconstruct the flux rope structure from in situ measurements in order to understand the geo-effectiveness associated with the CME magnetic field structure. Key results are obtained concerning how a weak CME can generate a relatively intense geomagnetic storm: (1) there were coronal holes at low latitudes, which could produce high speed streams (HSSs) to interact with the CME in interplanetary space; (2) the CME was bracketed between a slow wind ahead and a HSS behind, which enhanced the southward magnetic field inside the CME and gave rise to the unexpected geomagnetic storm.

Spectra of Cosmic Rays (CRs), particularly their features, may hold essential information about astroparticle and fundamental physics. Realizing that both $4$ PeV CR knee and electron's TeV break have nearly identical Lorentz factor, i.e. $\gamma \sim 10^6$, we propose that both CR nuclei and electron experience a threshold interaction with unknown light particle $\rm X$, widespread in the diffusive halo. The interaction occurs only when the Lorentz factor of CRs is larger than $10^6$. Via this process, the energy spectrum of each component has deformation above the threshold Lorentz factor, so both knee and spectral break of electrons can be naturally reproduced. Given the large uncertainties of the available measurements and simplified assumption, our model accommodates a wide mass range of the particle X, which is allowable from ultralight value to $\sim 1$ eV.

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.

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.

The active reflector is one of the three main innovations of the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The deformation of such a huge spherically shaped reflector into different transient parabolic shapes is achieved by using 2225 hydraulic actuators which change the position of the 2225 nodes through the connected down tied cables. For each different tracking process of the telescope, more than 1/3 of these 2225 actuators must be in operation to tune the parabolic aperture accurately to meet the surface error restriction. It means that some of these actuators are inevitably located within the main beam of the receiver, and the Electromagnetic Interference (EMI) from the actuators must be mitigated to ensure the scientific output of the telescope. Based on the threshold level of interference detrimental to radio astronomy presented in ITU-R Recommendation RA.769 and EMI measurements, the shielding efficiency (SE) requirement of each actuator is set to be 80dB in the frequency range from 70MHz to 3GHz. Therefore, Electromagnetic Compatibility (EMC) was taken into account in the actuator design by measures such as power line filters, optical fibers, shielding enclosures and other structural measures. In 2015, all the actuators had been installed at the FAST site. Till now, no apparent EMI from the actuators has been detected by the receiver, which proves the effectiveness of these EMC measures.

Diffusion of elements in the atmosphere and envelope of a star can drastically alter its surface composition, leading to extreme chemical peculiarities. We consider the case of hot subdwarfs, where surface helium abundances range from practically zero to almost 100 percent. Since hot subdwarfs can form via a number of different evolution channels, a key question concerns how the formation mechanism is connected to the present surface chemistry. A sequence of extreme horizontal branch star models was generated by producing post-common envelope stars from red giants. Evolution was computed with MESA from envelope ejection up to core-helium ignition. Surface abundances were calculated at the zero-age horizontal branch for models with and without diffusion. A number of simulations also included radiative levitation. The goal was to study surface chemistry during evolution from cool giant to hot subdwarf and determine when the characteristic subdwarf surface is established. Only stars leaving the giant branch close to core-helium ignition become hydrogen-rich subdwarfs at the zero-age horizontal branch. Diffusion, including radiative levitation, depletes the initial surface helium in all cases. All subdwarf models rapidly become more depleted than observations allow. Surface abundances of other elements follow observed trends in general, but not in detail. Additional physics is required.

The remnant of supernova explosion is widely believed to be the acceleration site of high-energy cosmic ray particles. The acceleration timescale is, however, typically very long. Here we report the detection of a variable $\gamma$-ray source with the Fermi Large Area Telescope, which is positionally and temporally consistent with a peculiar supernova, iPTF14hls. A quasi-stellar object SDSS J092054.04+504251.5, which is probably a blazar according to the infrared data, is found in the error circle of the $\gamma$-ray source. More data about the $\gamma$-ray source and SDSS J092054.04+504251.5 are needed to confirm their association. On the other hand, if the association between the $\gamma$-ray source and the supernova is confirmed, this would be the first time to detect high-energy $\gamma$-ray emission from a supernova, suggesting very fast particle acceleration by supernova explosions.

We present an investigation of the rotation and non-radial motion of a coronal mass ejection (CME) from AR 12468 on 2015 December 16 using observations from SDO, SOHO, STEREO A and Wind. The EUV and HMI observations of the source region show that the associated magnetic flux rope (MFR) axis pointed to the east before the eruption. We use a nonlinear fore-free field (NLFFF) extrapolation to determine the configuration of the coronal magnetic field and calculate the magnetic energy density distributions at different heights. The distribution of the magnetic energy density shows a strong gradient toward the northeast. The propagation direction of the CME from a Graduated Cylindrical Shell (GCS) modeling deviates from the radial direction of the source region by about 45 degr in longitude and about 30 degr in latitude, which is consistent with the gradient of the magnetic energy distribution around the AR.The MFR axis determined by the GCS modeling points southward, which has rotated counterclockwise by about 95 degr compared with the orientation of the MFR in the low corona.The MFR reconstructed by a Grad-Shafranov (GS) method at 1 AU has almost the same orientation as the MFR from the GCS modeling, which indicates that the MFR rotation occurred in the low corona. It is the rotation of the MFR that caused the intense geomagnetic storm with the minimum Dst of -155 nT. These results suggest that the coronal magnetic field surrounding the MFR plays a crucial role in the MFR rotation and propagation direction.

We analyze the Sun's shadow observed with the Tibet-III air shower array and find that the shadow's center deviates northward (southward) from the optical solar disc center in the "Away" ("Toward") IMF sector. By comparing with numerical simulations based on the solar magnetic field model, we find that the average IMF strength in the "Away" ("Toward") sector is $1.54 \pm 0.21_{\rm stat} \pm 0.20_{\rm syst}$ ($1.62 \pm 0.15_{\rm stat} \pm 0.22_{\rm syst}$) times larger than the model prediction. These demonstrate that the observed Sun's shadow is a useful tool for the quantitative evaluation of the average solar magnetic field.

We analyze the Sun's shadow observed with the Tibet-III air shower array and find that the shadow's center deviates northward (southward) from the optical solar disc center in the "Away" ("Toward") IMF sector. By comparing with numerical simulations based on the solar magnetic field model, we find that the average IMF strength in the "Away" ("Toward") sector is $1.54 \pm 0.21_{\rm stat} \pm 0.20_{\rm syst}$ ($1.62 \pm 0.15_{\rm stat} \pm 0.22_{\rm syst}$) times larger than the model prediction. These demonstrate that the observed Sun's shadow is a useful tool for the quantitative evaluation of the average solar magnetic field.

The discrepancy between the theoretical prediction of primordial lithium abundances and astronomical observations is called the Lithium Problem. We find that extra contributions from non-thermal hydrogen and helium during Big Bang nucleosynthesis can explain the discrepancy, for both Li-7 and Li-6, and will change the deuterium abundance only little. The allowed parameter space of such an amount of non-thermal particles and the energy range is shown. The hypothesis is stable regardless of the cross-section uncertainty of relevant reactions and the explicit shape of the energy spectrum.

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

The DArk Matter Particle Explorer (DAMPE), a high energy cosmic ray and $\gamma$-ray detector in space, has recently reported the new measurement of the total electron plus positron flux between 25 GeV and 4.6 TeV. A spectral softening at $\sim0.9$ TeV and a tentative peak at $\sim1.4$ TeV have been reported. We study the physical implications of the DAMPE data in this work. The presence of the spectral break significantly tightens the constraints on the model parameters to explain the electron/positron excesses. The spectral softening can either be explained by the maximum acceleration limits of electrons by astrophysical sources, or a breakdown of the common assumption of continuous distribution of electron sources at TeV energies in space and time. The tentive peak at $\sim1.4$ TeV implies local sources of electrons/positrons with quasi-monochromatic injection spectrum. We find that the cold, ultra-relativistic $e^+e^-$ winds from pulsars may give rise to such a structure. The pulsar is requird to be middle-aged, relatively slowly-rotated, mildly magnetized, and isolated in a density cavity. The annihilation of DM particles ($m_{\chi}\sim1.5$ TeV) into $e^+e^-$ pairs in a nearby clump or an over-density region may also explain the data. In the DM scenario, the inferred clump mass (or density enhancement) is about $10^7-10^8$ M$_\odot$ (or $17-35$ times of the canonical local density) assuming a thermal production cross section, which is relatively extreme compared with the expectation from numerical simulations. A moderate enhancement of the annihilation cross section via, e.g., the Sommerfeld mechanism or non-thermal production, is thus needed.

By the end of 2014, a cosmic ray muon telescope was installed at Zhongshan Station in Antarctic and has been continuously collecting data since then. It is the first surface muon telescope to be built in Antarctic. In June 2015, five CMEs were ejected towards the Earth initiating a big large Forbush decrease (FD) event. We conduct a comprehensive study of the galactic cosmic ray intensity fluctuations during the FD using the data from cosmic ray detectors of multiple stations (Zhongshan, McMurdo, South Polar and Nagoya) and he solar wind measurements from ACE and WIND. A pre-increase before the shock arrival was observed. Distinct differences exist in the timelines of the galactic cosmic ray recorded by the neutron monitors and the muon telescopes. FD onset for Zhongshan muon telescope is delayed (2.5h) with respect to SSC onset. This FD had a profile of four-step decrease. The traditional one- or two-step classification of FDs was inadequate to explain this FD.

The Giant Radio Array for Neutrino Detection (GRAND) aims at detecting ultra-high energy extraterrestrial neutrinos via the extensive air showers induced by the decay of tau leptons created in the interaction of neutrinos under the Earth's surface. Consisting of an array of $\sim10^5$ radio antennas deployed over $\sim 2\times10^5\,\rm {km}^2$, GRAND plans to reach, for the first time, an all-flavor sensitivity of $\sim1.5\times10^{-10} \,\rm GeV\, cm^{-2} \,s^{-1}\, sr^{-1}$ above $5\times10^{17}$ eV and a sub-degree angular resolution, beyond the reach of other planned detectors. We describe here preliminary designs and simulation results, plans for the ongoing, staged approach to the construction of GRAND, and the rich research program made possible by GRAND's design sensitivity and angular resolution.

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

We investigate the coronal and interplanetary evolution of a coronal mass ejection (CME) launched on 2010 September 4 from a source region linking two active regions (ARs) 11101 and 11103, using extreme ultraviolet imaging, magnetogram, white-light and in situ observations from SDO, STEREO, SOHO, VEX and Wind. A potential-field source-surface model is employed to examine the configuration of the coronal magnetic field surrounding the source region. The graduated cylindrical shell model and a triangulation method are applied to determine the kinematics of the CME in the corona and interplanetary space. From the remote sensing and in situ observations we obtain some key results: (1) the CME was deflected in both the eastward and southward directions in the low corona by the magnetic pressure from the two ARs and possibly interacted with another ejection, which caused that the CME arrived at VEX that was longitudinally distant from the source region; (2) although VEX was closer to the Sun, the observed and derived CME arrival times at VEX are not earlier than those at Wind, which suggests the importance of determining both the frontal shape and propagation direction of the CME in interplanetary space; (3) the ICME was compressed in the radial direction while the longitudinal transverse size was extended.

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.

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.

Propagation of coronal mass ejections (CMEs) from the Sun far into interplanetary space is not well understood due to limited observations. In this study we examine the propagation characteristics of two geo-effective CMEs, which occurred on 2005 May 6 and 13, respectively. Significant heliospheric consequences associated with the two CMEs are observed, including interplanetary CMEs (ICMEs) at the Earth and Ulysses, interplanetary shocks, a long-duration type II radio burst, and intense geomagnetic storms. We use coronagraph observations from SOHO/LASCO, frequency drift of the long-duration type II burst, in situ measurements at the Earth and Ulysses, and magnetohydrodynamic (MHD) propagation of the observed solar wind disturbances at 1 AU to track the CMEs from the Sun far into interplanetary space. We find that both of the two CMEs underwent a major deceleration within 1 AU and thereafter a gradual deceleration when they propagated from the Earth to deep interplanetary space due to interactions with the ambient solar wind. The results also reveal that the two CMEs interacted with each other in the distant interplanetary space even though their launch times on the Sun were well separated. The intense geomagnetic storm for each case was caused by the southward magnetic fields ahead of the CME, stressing the critical role of the sheath region in geomagnetic storm generation, although for the first case there is a corotating interaction region involved.

The Giant Radio Array for Neutrino Detection (GRAND) is a planned array of ~200 000 radio antennas deployed over ~200 000 km2 in a mountainous site. It aims primarly at detecting high-energy neutrinos via the observation of extensive air showers induced by the decay in the atmosphere of taus produced by the interaction of cosmic neutrinos under the Earth surface. GRAND aims at reaching a neutrino sensitivity of 5.10$^{11}$ E$^{-2}$ GeV$^{-1}$cm$^{-2}$s$^{-1}$sr$^{-1}$ above 3.10$^{16}$ eV. This ensures the detection of cosmogenic neutrinos in the most pessimistic source models, and ~50 events per year are expected for the standard models. The instrument will also detect UHECRs and possibly FRBs. Here we show how our preliminary design should enable us to reach our sensitivity goals, and discuss the steps to be taken to achieve GRAND.

We examine the structure, propagation and expansion of the shock associated with the 2012 July 23 extreme coronal mass ejection (CME). Characteristics of the shock determined from multi-point imaging observations are compared to in situ measurements at different locations and a complex radio type II burst, which according to our definition has multiple branches that may not all be fundamental-harmonic related. The white-light shock signature can be modeled reasonably well by a spherical structure and was expanding backward even on the opposite side of the Sun. The expansion of the shock, which was roughly self-similar after the first $\sim$1.5 hours from launch, largely dominated over the translation of the shock center for the time period of interest. Our study also suggests a bow-shock morphology around the nose at later times due to the outward motion in combination with the expansion of the ejecta. The shock decayed and failed to reach Mercury in the backward direction and STEREO B and Venus in the lateral directions, as indicated by the imaging and in situ observations. The shock in the nose direction, however, may persist to the far outer heliosphere, with predicted impact on Dawn around 06 UT on July 25 and on Jupiter around 23:30 UT on July 27 by an MHD model. The type II burst shows properties generally consistent with the spatial/temporal variations of the shock deduced from imaging and in situ observations. In particular, the low-frequency bands agree well with the in situ measurements of a very low density ahead of the shock at STEREO A.

We substantially update the capabilities of the open-source software instrument Modules for Experiments in Stellar Astrophysics (MESA). MESA can now simultaneously evolve an interacting pair of differentially rotating stars undergoing transfer and loss of mass and angular momentum, greatly enhancing the prior ability to model binary evolution. New MESA capabilities in fully coupled calculation of nuclear networks with hundreds of isotopes now allow MESA to accurately simulate advanced burning stages needed to construct supernova progenitor models. Implicit hydrodynamics with shocks can now be treated with MESA, enabling modeling of the entire massive star lifecycle, from pre-main sequence evolution to the onset of core collapse and nucleosynthesis from the resulting explosion. Coupling of the GYRE non-adiabatic pulsation instrument with MESA allows for new explorations of the instability strips for massive stars while also accelerating the astrophysical use of asteroseismology data. We improve treatment of mass accretion, giving more accurate and robust near-surface profiles. A new MESA capability to calculate weak reaction rates "on-the-fly" from input nuclear data allows better simulation of accretion induced collapse of massive white dwarfs and the fate of some massive stars. We discuss the ongoing challenge of chemical diffusion in the strongly coupled plasma regime, and exhibit improvements in MESA that now allow for the simulation of radiative levitation of heavy elements in hot stars. We close by noting that the MESA software infrastructure provides bit-for-bit consistency for all results across all the supported platforms, a profound enabling capability for accelerating MESA's development.

The 2015 March 15 coronal mass ejection as one of the two that together drove the largest geomagnetic storm of solar cycle 24 so far was associated with sympathetic filament eruptions. We investigate the relations between the different filaments involved in the eruption. A surge-like small-scale filament motion is confirmed as the trigger that initiated the erupting filament with multi-wavelength observations and using a forced magnetic field extrapolation method. When the erupting filament moved to an open magnetic field region, it experienced an obvious acceleration process and was accompanied by a C-class flare and the rise of another larger filament that eventually failed to erupt. We measure the decay index of the background magnetic field, which presents a critical height of 118 Mm. Combining with a potential field source surface extrapolation method, we analyze the distributions of the large-scale magnetic field, which indicates that the open magnetic field region may provide a favorable condition for F2 rapid acceleration and have some relation with the largest solar storm. The comparison between the successful and failed filament eruptions suggests that the confining magnetic field plays an important role in the preconditions for an eruption.

We analyze multi-spacecraft observations associated with the 2012 July 12 Coronal Mass Ejection (CME), covering the source region on the Sun from SDO, stereoscopic imaging observations from STEREO, magnetic field characteristics at MESSENGER, and type II radio burst and in situ measurements from Wind. A triangulation method based on STEREO stereoscopic observations is employed to determine the kinematics of the CME, and the outcome is compared with the result derived from the type II radio burst with a solar wind electron density model. A Grad-Shafranov technique is applied to Wind in situ data to reconstruct the flux-rope structure and compare it with the observation of the solar source region, which helps understand the geo-effectiveness associated with the CME structure. Conclusions are as follows: (1) the CME undergoes an impulsive acceleration, a rapid deceleration before reaching MESSENGER, and then a gradual deceleration out to 1 AU, which should be noticed in CME kinematics models; (2) the type II radio burst was probably produced from a high-density interaction region between the CME-driven shock and a nearby streamer or from the shock flank with lower heights, which implies uncertainties in the determination of CME kinematics using solely type II radio bursts; (3) the flux-rope orientation and chirality deduced from in situ reconstruction at Wind agree with those obtained from solar source observations; (4) the prolonged southward magnetic field near the Earth is mainly from the axial component of the largely southward inclined flux rope, which indicates the importance of predicting both the flux-rope orientation and magnetic field components in geomagnetic activity forecasting.

We carry out two-dimensional global particle-in-cell simulations of the interaction between the solar wind and a dipole field to study the formation of the bow shock and magnetosphere. A self-reforming bow shock ahead of a dipole field is presented by using relatively high temporal-spatial resolutions. We find that (1) the bow shock and the magnetosphere are formed and reach a quasi-stable state after several ion cyclotron periods, and (2) under the Bz southward solar wind condition the bow shock undergoes a self-reformation for low \b{eta}i and high MA. Simultaneously, a magnetic reconnection in the magnetotail is found. For high \b{eta}i and low MA, the shock becomes quasi-stationary, and the magnetotail reconnection disappears. In addition, (3) the magnetopause deflects the magnetosheath plasmas. The sheath particles injected at the quasi-perpendicular region of the bow shock can be convected to downstream of an oblique shock region. A fraction of these sheath particles can leak out from the magnetosheath at the wings of the bow shock. Hence, the downstream situation is more complicated than that for a planar shock produced in local simulations.

We study the solar energetic particle (SEP) event associated with the 2012 July 23 extreme solar storm, for which STEREO and the spacecraft at L1 provide multi-point remote sensing and in situ observations. The extreme solar storm, with a superfast shock and extremely enhanced ejecta magnetic fields observed near 1 AU at STEREO A, was caused by the combination of successive coronal mass ejections (CMEs). Meanwhile, energetic particles were observed by STEREO and near-Earth spacecraft such as ACE and SOHO, suggestive of a wide longitudinal spread of the particles at 1 AU. Combining the SEP observations with in situ plasma and magnetic field measurements we investigate the longitudinal distribution of the SEP event in connection with the associated shock and CMEs. Our results underscore the complex magnetic configuration of the inner heliosphere formed by solar eruptions. The examinations of particle intensities, proton anisotropy distributions, element abundance ratios, magnetic connectivity and spectra also give important clues for the particle acceleration, transport and distribution.

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

As a follow-up study on Sun-to-Earth propagation of fast coronal mass ejections (CMEs), we examine the Sun-to-Earth characteristics of slow CMEs combining heliospheric imaging and in situ observations. Three events of particular interest, the 2010 June 16, 2011 March 25 and 2012 September 25 CMEs, are selected for this study. We compare slow CMEs with fast and intermediate-speed events, and obtain key results complementing the attempt of \citet{liu13} to create a general picture of CME Sun-to-Earth propagation: (1) the Sun-to-Earth propagation of a typical slow CME can be approximately described by two phases, a gradual acceleration out to about 20-30 solar radii, followed by a nearly invariant speed around the average solar wind level, (2) comparison between different types of CMEs indicates that faster CMEs tend to accelerate and decelerate more rapidly and have shorter cessation distances for the acceleration and deceleration, (3) both intermediate-speed and slow CMEs would have a speed comparable to the average solar wind level before reaching 1 AU, (4) slow CMEs have a high potential to interact with other solar wind structures in the Sun-Earth space due to their slow motion, providing critical ingredients to enhance space weather, and (5) the slow CMEs studied here lack strong magnetic fields at the Earth but tend to preserve a flux-rope structure with axis generally perpendicular to the radial direction from the Sun. We also suggest a "best" strategy for the application of a triangulation concept in determining CME Sun-to-Earth kinematics, which helps to clarify confusions about CME geometry assumptions in the triangulation and to improve CME analysis and observations.

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.

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 precise spectra of Cosmic Ray (CR) electrons and positrons have been published by the measurement of AMS-02. It is reasonable to regard the difference between the electrons and positrons spectra ( $\triangle \Phi= \Phi_{e^-}-\Phi_{e^+}$ ) as being dominated by primary electrons. Noticing that the resulting electron spectrum shows no sign of spectral softening above 20 GeV, which is in contrast with the prediction of standard model. In this work, we generalize the analytic one dimensional two-halo model of diffusion to a three dimensional realistic calculation by implementing a spatial variant diffusion coefficients in DRAGON package. As a result, we can reproduce the spectral hardening of protons observed by several experiments, and predict an excess of high energy primary electrons which agrees with the measurement reasonably well. Unlike the break spectrum obtained for protons, the model calculation predicts a smooth electron excess and thus slightly over predicts the flux from tens of GeV to 100GeV. To understand this issue, further experimental and theoretical studies are necessary.

We examine two successive flare eruptions (X5.4 and X1.3) on 2012 March 7 in the NOAA active region 11429 and investigate the magnetic field reconfiguration associated with the two eruptions. Using an advanced non-linear force-free field (NLFFF) extrapolation method based on the SDO/HMI vector magnetograms, we obtain a stepwise decrease in the magnetic free energy during the eruptions, which is roughly 20%-30% of the energy of the pre-flare phase. We also calculate the magnetic helicity, and suggest that the changes of the sign of the helicity injection rate might be associated with the eruptions. Through the investigation of the magnetic field evolution, we find that the appearance of the "implosion" phenomenon has a strong relationship with the occurrence of the first X-class flare. Meanwhile, the magnetic field changes of the successive eruptions with implosion and without implosion were well observed.

We study the role of the coronal magnetic field configuration of an active region in determining the propagation direction of a coronal mass ejection (CME). The CME occurred in the active region 11944 (S09W01) near the disk center on 2014 January 7 and was associated with an X1.2 flare. A new CME reconstruction procedure based on a polarimetric technique is adopted, which shows that the CME changed its propagation direction by around 28$^\circ$ in latitude within 2.5 R$_\odot$ and 43$^\circ$ in longitude within 6.5 R$_\odot$ with respect to the CME source region. This significant non-radial motion is consistent with the finding of M$\ddot{o}$stl et al. (2015). We use nonlinear force-free field (NLFFF) and potential field source surface (PFSS) extrapolation methods to determine the configurations of the coronal magnetic field. We also calculate the magnetic energy density distributions at different heights based on the extrapolations. Our results show that the active region coronal magnetic field has a strong influence on the CME propagation direction. This is consistent with the "channelling" by the active region coronal magnetic field itself, rather than deflection by nearby structures. These results indicate that the active region coronal magnetic field configuration has to be taken into account in order to determine CME propagation direction correctly.

In gamma ray astronomy, the energy range from sub-100GeV to TeV is crucial due to where there is a gap between space experiments and ground-based ones. In addition, observations in this energy range are expected to provide more details about the high energy emission from GRBs,and thus to understand EBL. Based on the observation results and the related knowledge, scientists may be able to unveil the mysteries of galaxy formation and the evolution of early universe. One of the principal issues for next generation Imaging Atmospheric Cherenkov Telescopes (IACT) is to achieve larger field of view (FoV). In this work, we report a refractive water convex lens as light collector to test the feasibility of a new generation of IACT, and some preliminary test results on the optical properties (the focal length, spot size, transmittance, etc.) of a 0.9 m diameter water lens, the photodetectors and DAQ system of a prototype are presented and discussed.

High-energy neutrino astronomy will probe the working of the most violent phenomena in the Universe. The Giant Radio Array for Neutrino Detection (GRAND) project consists of an array of $\sim10^5$ radio antennas deployed over $\sim$200000km$^2$ in a mountainous site. It aims at detecting high-energy neutrinos via the measurement of air showers induced by the decay in the atmosphere of $\tau$ leptons produced by the interaction of the cosmic neutrinos under the Earth surface. Our objective with GRAND is to reach a neutrino sensitivity of $3\times10^{-11}E^{-2}$GeV$^{-1}$cm$^{-2}$s$^{-1}$sr$^{-1}$ above $3 \times10^{16}$eV. This sensitivity ensures the detection of cosmogenic neutrinos in the most pessimistic source models, and about 100 events per year are expected for the standard models. GRAND would also probe the neutrino signals produced at the potential sources of UHECRs. We show how our preliminary design should enable us to reach our sensitivity goals, and present the experimental characteristics. We assess the possibility to adapt GRAND to other astrophysical radio measurements. We discuss in this token the technological options for the detector and the steps to be taken to achieve the GRAND project.

The largest geomagnetic storms of solar cycle 24 so far occurred on 2015 March 17 and June 22 with $D_{\rm st}$ minima of $-223$ and $-195$ nT, respectively. Both of the geomagnetic storms show a multi-step development. We examine the plasma and magnetic field characteristics of the driving coronal mass ejections (CMEs) in connection with the development of the geomagnetic storms. A particular effort is to reconstruct the in situ structure using a Grad-Shafranov technique and compare the reconstruction results with solar observations, which gives a larger spatial perspective of the source conditions than one-dimensional in situ measurements. Key results are obtained concerning how the plasma and magnetic field characteristics of CMEs control the geomagnetic storm intensity and variability: (1) a sheath-ejecta-ejecta mechanism and a sheath-sheath-ejecta scenario are proposed for the multi-step development of the 2015 March 17 and June 22 geomagnetic storms, respectively; (2) two contrasting cases of how the CME flux-rope characteristics generate intense geomagnetic storms are found, which indicates that a southward flux-rope orientation is not a necessity for a strong geomagnetic storm; and (3) the unexpected 2015 March 17 intense geomagnetic storm resulted from the interaction between two successive CMEs plus the compression by a high-speed stream from behind, which is essentially the "perfect storm" scenario proposed by \citet[][i.e., a combination of circumstances results in an event of unusual magnitude]{liu14a}, so the "perfect storm" scenario may not be as rare as the phrase implies.

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 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 Large sky Area Multi-Object Spectroscopic Telescope (LAMOST) General Survey is a spectroscopic survey that will eventually cover approximately half of the celestial sphere and collect 10 million spectra of stars, galaxies and QSOs. Objects both in the pilot survey and the first year general survey are included in the LAMOST First Data Release (DR1). The pilot survey started in October 2011 and ended in June 2012, and the data have been released to the public as the LAMOST Pilot Data Release in August 2012. The general survey started in September 2012, and completed its first year of operation in June 2013. The LAMOST DR1 includes a total of 1202 plates containing 2,955,336 spectra, of which 1,790,879 spectra have observed signal-to-noise S/N >10. All data with S/N>2 are formally released as LAMOST DR1 under the LAMOST data policy. This data release contains a total of 2,204,696 spectra, of which 1,944,329 are stellar spectra, 12,082 are galaxy spectra and 5,017 are quasars. The DR1 includes not only spectra, but also three stellar catalogues with measured parameters: AFGK-type stars with high quality spectra (1,061,918 entries), A-type stars (100,073 entries), and M stars (121,522 entries). This paper introduces the survey design, the observational and instrumental limitations, data reduction and analysis, and some caveats. Description of the FITS structure of spectral files and parameter catalogues is also provided.

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 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 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 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 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 standard model of cosmic ray propagation has been very successful in explaining all kinds of the Galactic cosmic ray spectra. However, high precision measurement recently revealed the appreciable discrepancy between data and model expectation, from spectrum observations of $\gamma$-rays, $e^+/e^-$ and probably the $B/C$ ratio starting from $\sim$10 GeV energy. In this work, we propose that the fresh cosmic rays, which are supplied by the young accelerators and detained by local magnetic field, can contribute additional secondary particles interacting with local materials. As this early cosmic ray has a hard spectrum, the model calculation results in a two-component $\gamma$-ray spectrum, which agree very well with the observation. Simultaneously, the expected neutrino number from the galactic plane could contribute $\sim60\%$ of IceCube observation neutrino number below a few hundreds of TeV. The same pp-collision process can account for a significant amount of the positron excesses. Under this model, it is expected that the excesses in $\overline p/p$ and $B/C$ ratio will show up when energy is above $\sim$10 GeV. We look forward that the model will be tested in the near future by new observations from AMS02, IceCube, AS$\gamma$, HAWC and future experiments such as LHASSO, HiSCORE and CTA.

The origin of the knee in cosmic ray spectrum remains to be an unsolved fundamental problem. There are various kinds of models which predict different break positions and the compositions of the knee. In this work, we suggest to use diffuse $\gamma$-rays and neutrinos as probes to test these models. Based on several typical types of the composition models, the diffuse $\gamma$-ray and neutrino spectra are calculated, which show distinctive cutoff behaviours at energies from tens of TeV to multi-PeV. The expected flux will be observable by the newly upgraded Tibet-AS$\gamma$+MD (muon detector) experiment as well as more sensitive future projects, such as LHAASO and HiSCORE. By comparing the neutrino spectrum with the recent observations by IceCube experiment, we find that the diffuse neutrinos from interactions between the cosmic rays and the interstellar medium may not be responsible to the majority of the IceCube events. Future measurements of the neutrinos may be able to identify the Galactic diffuse component and further shed light on the problem of the knee of cosmic rays.

We have computed a new grid of evolutionary subdwarf B star (sdB) models from the start of central He burning, taking into account atomic diffusion due to radiative levitation, gravitational settling, concentration diffusion, and thermal diffusion. We have computed the non-adiabatic pulsation properties of the models and present the predicted p-mode and g-mode instability strips. In previous studies of the sdB instability strips, artificial abundance enhancements of Fe and Ni were introduced in the pulsation driving layers. In our models, the abundance enhancements of Fe and Ni occur naturally, eradicating the need to use artificial enhancements. We find that the abundance increases of Fe and Ni were previously underestimated and show that the instability strip predicted by our simulations solves the so-called blue edge problem of the subdwarf B star g-mode instability strip. The hottest known g-mode pulsator, KIC 10139564, now resides well within the instability strip {even when only modes with low spherical degrees (l<=2) are considered.

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

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

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

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

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

We report on a clear solar-cycle variation of the Sun's shadow in the 10 TeV cosmic-ray flux observed by the Tibet air shower array during a full solar cycle from 1996 to 2009. In order to clarify the physical implications of the observed solar cycle variation, we develop numerical simulations of the Sun's shadow, using the Potential Field Source Surface (PFSS) model and the Current Sheet Source Surface (CSSS) model for the coronal magnetic field. We find that the intensity deficit in the simulated Sun's shadow is very sensitive to the coronal magnetic field structure, and the observed variation of the Sun's shadow is better reproduced by the CSSS model. This is the first successful attempt to evaluate the coronal magnetic field models by using the Sun's shadow observed in the TeV cosmic-ray flux.

We report on a clear solar-cycle variation of the Sun's shadow in the 10 TeV cosmic-ray flux observed by the Tibet air shower array during a full solar cycle from 1996 to 2009. In order to clarify the physical implications of the observed solar cycle variation, we develop numerical simulations of the Sun's shadow, using the Potential Field Source Surface (PFSS) model and the Current Sheet Source Surface (CSSS) model for the coronal magnetic field. We find that the intensity deficit in the simulated Sun's shadow is very sensitive to the coronal magnetic field structure, and the observed variation of the Sun's shadow is better reproduced by the CSSS model. This is the first successful attempt to evaluate the coronal magnetic field models by using the Sun's shadow observed in the TeV cosmic-ray flux.

We present new quasars discovered in the vicinity of the Andromeda and Triangulum galaxies with the LAMOST during the 2010 and 2011 observational seasons. Quasar candidates are selected based on the available SDSS, KPNO 4 m telescope, XSTPS optical, and WISE near infrared photometric data. We present 509 new quasars discovered in a stripe of ~135 sq. deg from M31 to M33 along the Giant Stellar Stream in the 2011 pilot survey datasets, and also 17 new quasars discovered in an area of ~100 sq. deg that covers the central region and the southeastern halo of M31 in the 2010 commissioning datasets. These 526 new quasars have i magnitudes ranging from 15.5 to 20.0, redshifts from 0.1 to 3.2. They represent a significant increase of the number of identified quasars in the vicinity of M31 and M33. There are now 26, 62 and 139 known quasars in this region of the sky with i magnitudes brighter than 17.0, 17.5 and 18.0 respectively, of which 5, 20 and 75 are newly-discovered. These bright quasars provide an invaluable collection with which to probe the kinematics and chemistry of the ISM/IGM in the Local Group of galaxies. A total of 93 quasars are now known with locations within 2.5 deg of M31, of which 73 are newly discovered. Tens of quasars are now known to be located behind the Giant Stellar Stream, and hundreds behind the extended halo and its associated substructures of M31. The much enlarged sample of known quasars in the vicinity of M31 and M33 can potentially be utilized to construct a perfect astrometric reference frame to measure the minute PMs of M31 and M33, along with the PMs of substructures associated with the Local Group of galaxies. Those PMs are some of the most fundamental properties of the Local Group.

A new hybrid experiment has been started by AS{\gamma} experiment at Tibet, China, since August 2011, which consists of a low threshold burst-detector-grid (YAC-II, Yangbajing Air shower Core array), the Tibet air-shower array (Tibet-III) and a large underground water Cherenkov muon detector (MD). In this paper, the capability of the measurement of the chemical components (proton, helium and iron) with use of the (Tibet-III+YAC-II) is investigated by means of an extensive Monte Carlo simulation in which the secondary particles are propagated through the (Tibet-III+YAC-II) array and an artificial neural network (ANN) method is applied for the primary mass separation. Our simulation shows that the new installation is powerful to study the chemical compositions, in particular, to obtain the primary energy spectrum of the major component at the knee.

A new hybrid experiment has been started by AS{\gamma} experiment at Tibet, China, since August 2011, which consists of a low threshold burst-detector-grid (YAC-II, Yangbajing Air shower Core array), the Tibet air-shower array (Tibet-III) and a large underground water Cherenkov muon detector (MD). In this paper, the capability of the measurement of the chemical components (proton, helium and iron) with use of the (Tibet-III+YAC-II) is investigated by means of an extensive Monte Carlo simulation in which the secondary particles are propagated through the (Tibet-III+YAC-II) array and an artificial neural network (ANN) method is applied for the primary mass separation. Our simulation shows that the new installation is powerful to study the chemical compositions, in particular, to obtain the primary energy spectrum of the major component at the knee.

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

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

We revisit the idea that the Galactic center (GC) is the dominant source of Galactic cosmic rays (GCRs), based on a series of new observational evidence. A unified model is proposed to explain the new phenomena of GCRs and $\gamma$-rays simultaneously. The GCRs are thought to be accelerated during past activities of the GC. The pair production process of GCRs in the strong radiation field due to the GC activity is responsible for the knee structure of the cosmic ray spectra. A f raction of $e^+e^-$ produced by pair production interactions, can be reaccelerated in the induced bipolar jets and be transported into the halo, leaving the Fermi $\gamma$-ray bubbles and WMAP microwave haze as the remnant signal. Finally, the CRs diffuse in the bulge could further interact with the interstellar medium to produce low energy $e^+e^-$. After cooling down, these positrons may annihilate to produce the 511 keV line emission as discovered by INTEGRAL.

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.

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 have used 605 days of photometric data from the Kepler spacecraft to study KIC 6614501, a close binary system with an orbital period of 0.15749747(25) days (3.779939 hours), that consists of a low-mass subdwarf B (sdB) star and a white dwarf. As seen in many other similar systems, the gravitational field of the white dwarf produces an ellipsoidal deformation of the sdB which appears in the light curve as a modulation at two times the orbital frequency. The ellipsoidal deformation of the sdB implies that the system has a maximum inclination of ~40 degrees, with i \approx 20 degrees being the most likely. The orbital radial velocity of the sdB star is high enough to produce a Doppler beaming effect with an amplitude of 432 \pm 5 ppm, clearly visible in the folded light curve. The photometric amplitude that we obtain, K1 = 85.8 km/s, is ~12 per cent less than the spectroscopic RV amplitude of 97.2 \pm 2.0 km/s. The discrepancy is due to the photometric contamination from a close object at about 5 arcsec North West of KIC 6614501, which is difficult to remove. The atmospheric parameters of the sdB star, Teff = 23 700 \pm 500 K and log g = 5.70 \pm 0.10, imply that it is a rare object below the Extreme Horizontal Branch (EHB), similar to HD 188112 (Heber et al. 2003). The comparison with different evolutionary tracks suggests a mass between ~0.18 and ~0.25 Msun, too low to sustain core helium burning. If the mass was close to 0.18-0.19 Msun, the star could be already on the final He-core WD cooling track. A higher mass, up to ~0.25 Msun, would be compatible with a He-core WD progenitor undergoing a cooling phase in a H-shell flash loop. A third possibility, with a mass between ~0.32 and ~0.40 Msun, can not be excluded and would imply that the sdB is a "normal" (but with an unusually low mass) EHB star burning He...

We report on the properties of strong pulses from PSR B0656+14 by analyzing the data obtained using Urumqi 25-m radio telescope at 1540 MHz from August 2007 to September 2010. In 44 hrs of observational data, a total of 67 pulses with signal-to-noise ratios above a 5-{\sigma} threshold were detected. The peak flux densities of these pulses are 58 to 194 times that of the average profile, and the pulse energies of them are 3 to 68 times that of the average pulse. These pulses are clustered around phases about 5 degrees ahead of the peak of the average profile. Comparing with the width of the average profile, they are relatively narrow, with the full widths at half-maximum range from 0.28 to 1.78 degrees. The distribution of pulse-energies of the pulses follows a lognormal distribution. These sporadic strong pulses detected from PSR B0656+14 are different in character from the typical giant pulses, and from its regular pulses.

We present our analyses of 15 months of Kepler data on KIC 10139564. We detected 57 periodicities with a variety of properties not previously observed all together in one pulsating subdwarf B star. Ten of the periodicities were found in the low-frequency region, and we associate them with nonradial g-modes. The other periodicities were found in the high-frequency region, which are likely p-modes. We discovered that most of the periodicities are components of multiplets with a common spacing. Assuming that multiplets are caused by rotation, we derive a rotation period of 25.6(1.8) days. The multiplets also allow us to identify the pulsations to an unprecedented extent for this class of pulsator. We also detect l<=2 multiplets, which are sensitive to the pulsation inclination and can constrain limb darkening via geometric cancellation factors. While most periodicities are stable, we detected several regions that show complex patterns. Detailed analyses showed these regions are complicated by several factors. Two are combination frequencies that originate in the superNyquist region and were found to be reflected below the Nyquist frequency. The Fourier peaks are clear in the superNyquist region, but the orbital motion of Kepler smears the Nyquist frequency in the barycentric reference frame and this effect is passed on to the subNyquist reflections. Others are likely multiplets but unstable in amplitudes and/or frequencies. The density of periodicities also make KIC 10139564 challenging to explain using published models. This menagerie of properties should provide tight constraints on structural models, making this subdwarf B star the most promising for applying asteroseismology.

A particle-in-cell code is used to examine contributions of the pickup ions (PIs) and the solar wind ions (SWs) to the cross shock electric field at the supercritical, perpendicular shocks. The code treats the pickup ions self-consistently as a third component. Herein, two different runs with relative pickup ion density of 25% and 55% are presented in this paper. Present preliminary results show that: (1) in the low percentage (25%) pickup ion case, the shock front is nonstationary. During the evolution of this perpendicular shock, a nonstationary foot resulting from the reflected solar wind ions is formed in front of the old ramp, and its amplitude becomes larger and larger. At last, the nonstationary foot grows up into a new ramp and exceeds the old one. Such a nonstationary process can be formed periodically. hen the new ramp begins to be formed in front of the old ramp, the Hall term mainly contributed by the solar wind ions becomes more and more important. The electric field Ex is dominated by the Hall term when the new ramp exceeds the old one. Furthermore, an extended and stationary foot in pickup ion gyro-scale is located upstream of the nonstationary/self-reforming region within the shock front, and is always dominated by the Lorentz term contributed by the pickup ions; (2) in the high percentage (55%) pickup ion case, the amplitude of the stationary foot is increased as expected. One striking point is that the nonstationary region of the shock front evidenced by the self-reformation disappears. Instead, a stationary extended foot dominated by Lorentz term contributed by the pickup ions, and a tationary ramp dominated by Hall term contributed by the solar wind ions are clearly evidenced. The significance of the cross electric field on ion dynamics is also discussed.

There may be non-thermal cosmic rays during big-bang nucleosynthesis (BBN) epoch (dubbed as BBNCRs). This paper investigated whether such BBNCRs can be the origin of Lithium problem or not. It can be expected that BBNCRs flux will be small in order to keep the success of standard BBN (SBBN). With favorable assumptions on the BBNCR spectrum between 0.09 -- 4 MeV, our numerical calculation showed that extra contributions from BBNCRs can account for the $^7$Li abundance successfully. However $^6$Li abundance is only lifted an order of magnitude, which is still much lower than the observed value. As the deuteron abundance is very sensitive to the spectrum choice of BBNCRs, the allowed parameter space for the spectrum is strictly constrained. We should emphasize that the acceleration mechanism for BBNCRs in the early universe is still an open question. For example, strong turbulent magnetic field is probably the solution to the problem. Whether such a mechanism can provide the required spectrum deserves further studies.

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

We report on extended photometry of two pulsating sdB stars in close binaries. For both cases, we use rotational splitting of the pulsation frequencies to show that the sdB component rotates much too slowly to be in synchronous rotation. We use a theory of tidal interaction in binary stars to place limits on the mass ratios that are independent of estimates based on the radial velocity curves. The companions have masses below 0.26 M\odot. The pulsation spectra show the signature of high-overtone g-mode pulsation. One star, KIC 11179657, has a clear sequence of g-modes with equal period spacings as well as several periodicities that depart from that trend. KIC 02991403 shows a similar sequence, but has many more modes that do not fit the simple pattern.

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

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

Context. Subdwarf B pulsators exhibit pressure (p) and/or gravity (g) modes. Their frequency spectra range from very simple, with few frequencies, to very rich, with more than fifty peaks in some cases. Balloon09 is a hybrid pulsating subdwarf B, showing a large number of p- and g-modes including a triplet and a quintuplet-like structures which are interpreted as p-mode frequency splittings due to stellar rotation. Photometric observations undertaken in two subsequent years revealed a change in these rotational splittings of 0.24 - 0.58 {\mu}Hz/yr. Aims. We analyse the possibility of angular momentum interchange between stellar rotation and internal gravity waves as a mechanism for the observed rotational splitting variations. Methods. The expected change in the rotational splitting of eigenmodes resulting from this mechanism are computed in the nonadiabatic linear approximation for a stellar structure model that is representative of the target star. The fact that g-modes are also observed in Balloon09 makes it a particularly suitable candidate for our study, because the change in the rotational splittings are proportional to the amplitude squared of the g-modes which, in this case, can be estimated from the observations. Results. We find that this mechanism is able to predict changes in the splittings that can be of the same order of magnitude as the observed variations.

Pulsations in subdwarf B stars are attributed to radiative levitation of iron-group elements in the stellar envelope. Until now, only iron diffusion is accounted for in stellar models used for sdB seismology. However, nickel has also been suggested as a contributor to the opacity bump that drives the pulsation modes. Stellar models including time-dependent atomic diffusion, as we compute here, are needed to evaluate the importance of different iron-group elements for mode driving. We perform detailed calculations of radiative accelerations of H, He, C, N, O, Ne, Mg, Fe and Ni and include these in Burgers' diffusion equations. We compute the evolution and non-adiabatic pulsations of a typical subdwarf B star. We show that, despite its lower initial abundance, nickel accumulates to comparable mass fractions as iron in the sdB envelope. For accurate determination of pulsation frequencies and mode instability, it is essential that diffusion of both metals are included in stellar models. The role of other iron-group elements remain to be evaluated.

Subdwarf B stars show chemical peculiarities that cannot be explained by diffusion theory alone. Both mass loss and turbulence have been invoked to slow down atomic diffusion in order to match observed abundances. The fact that some sdB stars show pulsations gives upper limits on the amount of mass loss and turbulent mixing allowed. Consequently, non-adiabatic asteroseismology has the potential to decide which process is responsible for the abundance anomalies. We compute for the first time seismic properties of sdB models with atomic diffusion included consistently during the stellar evolution. The diffusion equations with radiative forces are solved for H, He, C, N, O, Ne, Mg, Fe and Ni. We examine the effects of various mass-loss rates and mixed surface masses on the abundances and mode stability. It is shown that the mass-loss rates needed to simulate the observed He abundances (10^{-14}<=Mdot [Msun/yr]<=10^{-13}) are not consistent with observed pulsations. We find that for pulsations to be driven the rates should be Mdot<=10^{-15} Msun/yr. On the other hand, weak turbulent mixing of the outer 10^{-6} Msun can explain the He abundance anomalies while still allowing pulsations to be driven. The origin of the turbulence remains unknown but the presence of pulsations gives tight constraints on the underlying turbulence model.

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.

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.

We analyze individual pulses detected from RRAT J1819$-$1458. From April 2007 to April 2010, we carried out observations using the Nanshan 25-m radio telescope of Urumqi Observatory at a central frequency of 1541.25 MHz. We obtain a dispersion measure $DM=195.7\pm0.3$ pc cm^{-3} by analyzing all the 423 detected bursts. The tri-band pattern of arrival time residuals is confirmed by a single pulse timing analysis. Twenty-seven bimodal bursts located in the middle residual band are detected, and, profiles of two typical bimodal bursts and two individual single-peak pulses are presented. We determine the statistical properties of SNR and W$_{50}$ of bursts in different residual bands. The W$_{50}$ variation with SNR shows that the shapes of bursts are quite different from each other. The cumulative probability distribution of intensity for a possible power law with index $\alpha=1.6\pm0.2$ is inferred from the number of those bursts with $SNR\ge6$ and high intensities.

An antenna array devoted to the autonomous radio-detection of high energy cosmic rays is being deployed on the site of the 21 cm array radio telescope in XinJiang, China. Thanks in particular to the very good electromagnetic environment of this remote experimental site, self-triggering on extensive air showers induced by cosmic rays has been achieved with a small scale prototype of the foreseen antenna array. We give here a detailed description of the detector and present the first detection of extensive air showers with this prototype.

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 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 paper investigates the overall and detailed features of cosmic ray (CR) spectra in the knee region using the scenario of nuclei-photon interactions around the acceleration sources. Young supernova remnants can be the physical realities of such kind of CR acceleration sites. The results show that the model can well explain the following problems simultaneously with one set of source parameters: the knee of CR spectra and the sharpness of the knee, the detailed irregular structures of CR spectra, the so-called "component B" of Galactic CRs, and the electron/positron excesses reported by recent observations. The coherent explanation serves as evidence that at least a portion of CRs might be accelerated at the sources similar to young supernova remnants, and one set of source parameters indicates that this portion mainly comes from standard sources or from a single source.

We analyze the large-scale two-dimensional sidereal anisotropy of multi-TeV cosmic rays by Tibet Air Shower Array, with the data taken from 1999 November to 2008 December. To explore temporal variations of the anisotropy, the data set is divided into nine intervals, each in a time span of about one year. The sidereal anisotropy of magnitude about 0.1% appears fairly stable from year to year over the entire observation period of nine years. This indicates that the anisotropy of TeV Galactic cosmic rays remains insensitive to solar activities since the observation period covers more than a half of the 23rd solar cycle.

We analyze the large-scale two-dimensional sidereal anisotropy of multi-TeV cosmic rays by Tibet Air Shower Array, with the data taken from 1999 November to 2008 December. To explore temporal variations of the anisotropy, the data set is divided into nine intervals, each in a time span of about one year. The sidereal anisotropy of magnitude about 0.1% appears fairly stable from year to year over the entire observation period of nine years. This indicates that the anisotropy of TeV Galactic cosmic rays remains insensitive to solar activities since the observation period covers more than a half of the 23rd solar cycle.

Diffusion of atoms can be important during quiescent phases of stellar evolution. Particularly in the very thin inert envelopes of subdwarf B stars, diffusive movements will considerably change the envelope structure and the surface abundances on a short timescale. Also, the subdwarfs will inherit the effects of diffusion in their direct progenitors, namely giants near the tip of the red giant branch. This will influence the global evolution and the pulsational properties of subdwarf B stars. We investigate the impact of gravitational settling, thermal diffusion and concentration diffusion on the evolution and pulsations of subdwarf B stars. Our diffusive stellar models are compared with models evolved without diffusion. We constructed subdwarf B models with a mass of 0.465 Msun from a 1 and 3 Msun ZAMS progenitor. The low mass star ignited helium in an energetic flash, while the intermediate mass star started helium fusion gently. For each progenitor type we computed series with and without atomic diffusion. Atomic diffusion in red giants causes the helium core mass at the onset of helium ignition to be larger. We find an increase of 0.0015 Msun for the 1 Msun model and 0.0036 Msun for the 3 Msun model. The effects on the red giant surface abundances are small after the first dredge up. The evolutionary tracks of the diffusive subdwarf B models are shifted to lower surface gravities and effective temperatures due to outward diffusion of hydrogen. This affects both the frequencies of the excited modes and the overall frequency spectrum. Especially the structure and pulsations of the post-non-degenerate sdB star are drastically altered, proving that atomic diffusion cannot be ignored in these stars.

Using the Tibet-III air shower array, we search for TeV gamma-rays from 27 potential Galactic sources in the early list of bright sources obtained by the Fermi Large Area Telescope at energies above 100 MeV. Among them, we observe 7 sources instead of the expected 0.61 sources at a significance of 2 sigma or more excess. The chance probability from Poisson statistics would be estimated to be 3.8 x 10^-6. If the excess distribution observed by the Tibet-III array has a density gradient toward the Galactic plane, the expected number of sources may be enhanced in chance association. Then, the chance probability rises slightly, to 1.2 x 10^-5, based on a simple Monte Carlo simulation. These low chance probabilities clearly show that the Fermi bright Galactic sources have statistically significant correlations with TeV gamma-ray excesses. We also find that all 7 sources are associated with pulsars, and 6 of them are coincident with sources detected by the Milagro experiment at a significance of 3 sigma or more at the representative energy of 35 TeV. The significance maps observed by the Tibet-III air shower array around the Fermi sources, which are coincident with the Milagro >=3sigma sources, are consistent with the Milagro observations. This is the first result of the northern sky survey of the Fermi bright Galactic sources in the TeV region.

Using the Tibet-III air shower array, we search for TeV gamma-rays from 27 potential Galactic sources in the early list of bright sources obtained by the Fermi Large Area Telescope at energies above 100 MeV. Among them, we observe 7 sources instead of the expected 0.61 sources at a significance of 2 sigma or more excess. The chance probability from Poisson statistics would be estimated to be 3.8 x 10^-6. If the excess distribution observed by the Tibet-III array has a density gradient toward the Galactic plane, the expected number of sources may be enhanced in chance association. Then, the chance probability rises slightly, to 1.2 x 10^-5, based on a simple Monte Carlo simulation. These low chance probabilities clearly show that the Fermi bright Galactic sources have statistically significant correlations with TeV gamma-ray excesses. We also find that all 7 sources are associated with pulsars, and 6 of them are coincident with sources detected by the Milagro experiment at a significance of 3 sigma or more at the representative energy of 35 TeV. The significance maps observed by the Tibet-III air shower array around the Fermi sources, which are coincident with the Milagro >=3sigma sources, are consistent with the Milagro observations. This is the first result of the northern sky survey of the Fermi bright Galactic sources in the TeV region.

This paper briefly summarizes the status of the cosmic ray observations by EAS (Extended Air Shower) experiments with energy below 10**16eV and the related studies of the hadronic interaction models. Based on the observed sharp knee structure and the irregularities of the cosmic ray spectrum around knee energy, plus the newly discovered electron and positron excess, the origin of the galactic cosmic rays and the single source model interpretation are discussed, but convincing evidence is not yet available. High precision measurements of the mass composition of primary cosmic rays at knee energy will be very useful to disentangle the problem. To reach this goal, a better understanding of the hadronic interaction models is crucial. It is good to see that more dedicated accelerator and cosmic ray experiments will be conducted soon. As one EAS component, the muon distribution and muon charge ratio are important for testing the hadronic interaction models. In addition muons are an important background to neutrino experiments and all underground ultra-low background experiments. They are also a very useful tool for the meteorological studies.

Realistic stellar models are essential to the forward modelling approach in asteroseismology. For practicality however, certain model assumptions are also required. For example, in the case of subdwarf B stars, one usually starts with zero-age horizontal branch structures without following the progenitor evolution. We analyse the effects of common assumptions in subdwarf B models on the g-mode pulsational properties. We investigate if and how the pulsation periods are affected by the H-profile in the core-envelope transition zone. Furthermore, the effects of C-production and convective mixing during the core helium flash are evaluated. Finally, we reanalyse the effects of stellar opacities on the mode excitation in subdwarf B stars. We find that helium settling causes a shift in the theoretical blue edge of the g-mode instability domain to higher effective temperatures. This results in a closer match to the observed instability strip of long-period sdB pulsators, particularly for l<=3 modes. We show further that the g-mode spectrum is extremely sensitive to the H-profile in the core-envelope transition zone. If atomic diffusion is efficient, details of the initial shape of the profile become less important in the course of evolution. Diffusion broadens the chemical gradients, and results in less effective mode trapping and different pulsation periods. Furthermore, we report on the possible consequences of the He-flash for the g-modes. The outer edge of a flash-induced convective region introduces an additional chemical transition in the stellar models, and the corresponding spike in the Brunt-Vaisala frequency produces a complicated mode trapping signature in the period spacings.

We develop a model anisotropy best-fitting to the two-dimensional sky-map of multi-TeV galactic cosmic ray (GCR) intensity observed with the Tibet III air shower (AS) array. By incorporating a pair of intensity excesses in the hydrogen deflection plane (HDP) suggested by Gurnett et al., together with the uni-directional and bi-directional flows for reproducing the observed global feature, this model successfully reproduces the observed sky-map including the "skewed" feature of the excess intensity from the heliotail direction, whose physical origin has long remained unknown. These additional excesses are modeled by a pair of the northern and southern Gaussian distributions, each placed ~50 degree away from the heliotail direction. The amplitude of the southern excess is as large as ~0.2 %, more than twice the amplitude of the northern excess. This implies that the Tibet AS experiment discovered for the first time a clear evidence of the significant modulation of GCR intensity in the heliotail and the asymmetric heliosphere.

We develop a model anisotropy best-fitting to the two-dimensional sky-map of multi-TeV galactic cosmic ray (GCR) intensity observed with the Tibet III air shower (AS) array. By incorporating a pair of intensity excesses in the hydrogen deflection plane (HDP) suggested by Gurnett et al., together with the uni-directional and bi-directional flows for reproducing the observed global feature, this model successfully reproduces the observed sky-map including the "skewed" feature of the excess intensity from the heliotail direction, whose physical origin has long remained unknown. These additional excesses are modeled by a pair of the northern and southern Gaussian distributions, each placed ~50 degree away from the heliotail direction. The amplitude of the southern excess is as large as ~0.2 %, more than twice the amplitude of the northern excess. This implies that the Tibet AS experiment discovered for the first time a clear evidence of the significant modulation of GCR intensity in the heliotail and the asymmetric heliosphere.

The UV-Excess Survey of the Northern Galactic Plane images a 10x185 degree wide band, centered on the Galactic Equator using the 2.5m Isaac Newton Telescope in four bands (U,g,r,HeI5875) down to ~21st-22nd magnitude (~20th in HeI5875). The setup and data reduction procedures are described. Simulations of the colours of main-sequence stars, giant, supergiants, DA and DB white dwarfs and AM CVn stars are made, including the effects of reddening. A first look at the data of the survey (currently 30% complete) is given.

Supernova remnants have long been regarded as sources of the Galactic cosmic rays up to petaelectronvolts, but convincing evidence is still lacking. In this work we explore the common origin of the subtle features of the cosmic ray spectra, such as the knee of cosmic ray spectra and the excesses of electron/positron fluxes recently observed by ATIC, H.E.S.S., Fermi-LAT and PAMELA. Numerical calculation shows that those features of cosmic ray spectra can be well reproduced in a scenario with e$^+$e$^-$ pair production by interactions between high energy cosmic rays and background photons in an environment similar to the young supernova remnant. The success of such a coherent explanation serves in turn as an evidence that at least a portion of cosmic rays might be accelerated at young supernova remnants.

There are many unknowns in the formation of subdwarf B stars. Different formation channels are considered to be possible and to lead to a variety of helium-burning subdwarfs. All seismic models to date, however, assume that a subdwarf B star is a post-helium-flash-core surrounded by a thin inert layer of hydrogen. We examine an alternative formation channel, in which the subdwarf B star originates from a massive (>~2 Msun) red giant with a non-degenerate helium-core. Although these subdwarfs may evolve through the same region of the log g-Teff diagram as the canonical post-flash subdwarfs, their interior structure is rather different. We examine how this difference affects their pulsation modes and whether it can be observed. Using detailed stellar evolution calculations we construct subdwarf B models from both formation channels. The iron accumulation in the driving region due to diffusion, which causes the excitation of the modes, is approximated by a Gaussian function. The pulsation modes and frequencies are calculated with a non-adiabatic pulsation code. A detailed comparison of two subdwarf B models from different channels, but with the same log g and Teff, shows that their mode excitation is different. The excited frequencies are lower for the post-flash than for the post-non-degenerate subdwarf B star. This is mainly due to the differing chemical composition of the stellar envelope. A more general comparison between two grids of models shows that the excited frequencies of most post-non-degenerate subdwarfs cannot be well-matched with the frequencies of post-flash subdwarfs. In the rare event that an acceptable seismic match is found, additional information, such as mode identification and log g and Teff determinations, allows us to distinguish between the two formation channels.

The INT/WFC Photometric H-alpha Survey of the Northern Galactic Plane (IPHAS) is an imaging survey being carried out in H-alpha, r' and i' filters, with the Wide Field Camera (WFC) on the 2.5-metre Isaac Newton Telescope (INT) to a depth of r'=20 (10 sigma). The survey is aimed at revealing large scale structure in our local galaxy, and also the properties of key early and late populations making up the Milky Way. Mapping emission line objects enables a particular focus on objects in the young and old stages of stellar evolution ranging from early T-Tauri stars to late planetary nebulae. In this paper we present the IPHAS Initial Data Release, primarily a photometric catalogue of about 200 million unique objects, coupled with associated image data covering about 1600 square degrees in three passbands. We note how access to the primary data products has been implemented through use of standard virtual observatory publishing interfaces. Simple traditional web access is provided to the main IPHAS photometric catalogue, in addition to a number of common catalogues (such as 2MASS) which are of immediate relevance. Access through the AstroGrid VO Desktop opens up the full range of analysis options, and allows full integration with the wider range of data and services available through the Virtual Observatory. The IDR represents the largest dataset published primarily through VO interfaces to date, and so stands as an examplar of the future of survey data mining. Examples of data access are given, including a cross-matching of IPHAS photometry with sources in the UKIDSS Galactic Plane Survey that validates the existing calibration of the best data.

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

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

The amplitude of the Compton-Getting (CG) anisotropy contains the power-law index of the cosmic-ray energy spectrum. Based on this relation and using the Tibet air-shower array data, we measure the cosmic-ray spectral index to be $-3.03 \pm 0.55_{stat} \pm < 0.62_{syst}$ between 6 TeV and 40 TeV, consistent with $-$2.7 from direct energy spectrum measurements. Potentially, this CG anisotropy analysis can be utilized to confirm the astrophysical origin of the ``knee'' against models for non-standard hadronic interactions in the atmosphere.

The amplitude of the Compton-Getting (CG) anisotropy contains the power-law index of the cosmic-ray energy spectrum. Based on this relation and using the Tibet air-shower array data, we measure the cosmic-ray spectral index to be $-3.03 \pm 0.55_{stat} \pm < 0.62_{syst}$ between 6 TeV and 40 TeV, consistent with $-$2.7 from direct energy spectrum measurements. Potentially, this CG anisotropy analysis can be utilized to confirm the astrophysical origin of the ``knee'' against models for non-standard hadronic interactions in the atmosphere.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We are planning to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV (10 - 1000 TeV) energy region by discriminating them from cosmic-ray hadrons. We report on the possibility of detection of gamma rays in the 100 TeV energy region in our field of view, based on the improved sensitivity of our air shower array deduced from the full Monte Carlo simulation.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We are planning to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV (10 - 1000 TeV) energy region by discriminating them from cosmic-ray hadrons. We report on the possibility of detection of gamma rays in the 100 TeV energy region in our field of view, based on the improved sensitivity of our air shower array deduced from the full Monte Carlo simulation.