The LIGO-Virgo and Fermi collaborations recently reported a possible joint detection of a sub-threshold gravitational wave (GW) event and a sub-threshold gamma-ray burst (GRB), GBM-190816, that occurred 1.57 s after the merger. We perform an independent analysis of the publicly available data and investigate the physical implications of this potential association. By carefully studying the following properties of GBM-190816 using Fermi/GBM data, including signal-to-noise ratio, duration, f-parameter, spectral properties, energetic properties, and its compliance with some GRB statistical correlations, we confirm that this event is likely a typical short GRB. Assuming its association with the sub-threshold GW event, the inferred luminosity is $1.47_{-1.04}^{+3.40} \times 10^{49}$ erg s$^{-1}$. Based on the available information of the sub-threshold GW event, we infer the mass ratio q of the compact binary as $q=2.26_{-1.43}^{+2.75}$ according to the reported range of luminosity distance. If the heavier compact object has a mass > 3 solar masses, q can be further constrained to $q=2.26_{-0.12}^{+2.75}$. The leading physical scenario invokes an NS-BH merger system with the NS tidally disrupted. Within this scenario, we constrain the physical properties of such a system to produce a GRB. The GW data may also allow an NS-BH system with no tidal disruption of the NS or a BH-BH merger. We apply the charged compact binary coalescence (cCBC) theory (for both a constant charge and an increasing charge for the merging members) to derive the model parameters to account for GBM-190816 and found that the required parameters are extreme. Finally, we argue that the fact that the observed GW-GRB delay time scale is comparable to that of GW170817/GRB 170817A suggests that the GW-GRB time delay of these two cases is mainly defined by the time scale for the jet to propagate to the GRB emission site.

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

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

In this paper, we estimate the eccentricity of the 10 BBHs in GWTC-1 by using the inspiral-only BBH waveform template EccentricFD. Firstly, we test our method with simulated eccentric BBHs. Afterwards we apply the method to the real BBH gravitational wave data. We find that the BBHs in GWTC-1 except GW151226, GW170608 and GW170729 admit almost zero eccentricity. Their upper limits on eccentricity rang from 0.033 to 0.084 with 90% credible interval at the reference frequency 10Hz. And their median eccentricities are 0.000. For GW151226, GW170608 and GW170729, the upper limits are higher than 0.1. But the median eccentricities are respectively 0.000, 0.064 and 0.000. We also point out the limitations of the inspiral-only waveform template in eccentricity measurement. The current eccentricity measurement results of GWTC-1 is not enough to constrain the formation mechanisms of BBHs, and more observational data are needed in the future.

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

The assumptions of large-scale homogeneity and isotropy underly the familiar Friedmann-Lema\^{\i}tre-Robertson-Walker (FLRW) metric that appears to be an accurate description of our Universe. In this paper, we propose a new strategy of testing the validity of the FLRW metric, based on the galactic-scale lensing systems where strongly lensed gravitational waves and their electromagnetic counterparts can be simultaneously detected. Each strong lensing system creates opportunity to infer the curvature parameter of the Universe. Consequently, combined analysis of many such systems will provide a model-independent tool to test the validity of the FLRW metric. Our study demonstrates that the third-generation ground based GW detectors, like the Einstein Telescope (ET) and space-based detectors, like the Big Bang Observer (BBO), are promising concerning determination of the curvature parameter or possible detection of deviation from the FLRW metric. Such accurate measurements of the FLRW metric can become a milestone in precision GW cosmology.

The equivalence principle is important in fundamental physics. The fifth force, as a describing formalism of the equivalence principle, may indicate the property of an unknown theory. Dark matter is one of the most mysterious objects in the current natural science. It is interesting to constrain the fifth force of dark matter. We propose a new method to use perihelion precession of planets to constrain the long-range fifth force of dark matter. Due to the high accuracy of perihelion precession observation, and the large difference of matter composition between the Sun and planets, we get one of the strongest constraints on the fifth force of dark matter. In the near future, the BepiColombo mission will be capable to improve the test by another factor of ten.

Effective-one-body (EOB) numerical-relativity (NR) waveform models for spin-aligned binary black holes (BBHs), known as the SEOBNR waveform models, are based on the EOB theoretical framework and NR simulations. SEOBNR models have played an important role in the LIGO scientific collaboration (LSC) gravitational wave (GW) data analysis for both signal search and parameter estimation. SEOBNR models for quasi-circular orbits have evolved through version 1 to version 4 by extending their validity domain and including more NR results. Along another direction, we recently extended SEOBNRv1 model to SEOBNRE model which is valid for spin-aligned BBH coalescence along eccentric orbits. In this paper we validate this theoretical waveform model by comparing them against the numerical relativity simulation bank, Simulating eXtreme Spacetimes (SXS) catalog. In total, 278 NR waveforms are investigated which include binaries with large eccentricity; large spin and large mass ratio. Our SEOBNRE can model the NR waveforms quite well. The fitting factor for most of the 278 waveforms is larger than 99\%. It indicates that the SEOBNRE model could be used as template waveforms for eccentric spin-aligned BBH coalescence.

Deep learning method develops very fast as a tool for data analysis these years. Such a technique is quite promising to treat gravitational wave detection data. There are many works already in the literature which used deep learning technique to process simulated gravitational wave data. In this paper we apply deep learning to LIGO O1 data. In order to improve the weak signal recognition we adjust the convolutional neural network (CNN) a little bit. Our adjusted convolutional neural network admits comparable accuracy and efficiency of signal recognition as other deep learning works published in the literature. Based on our adjusted CNN, we can clearly recognize the eleven confirmed gravitational wave events included in O1 and O2. And more we find about 2000 gravitational wave triggers in O1 data.

We investigate the scalar-tensor gravity of Damour and Esposito-Far\`ese (DEF), which predicts non-trivial phenomena in the nonperturbative strong-field regime for neutron stars (NSs). Instead of solving the modified Tolman-Oppenheimer-Volkoff equations, we construct reduced-order surrogate models, coded in the pySTGROM package, to predict the relations of a NS radius, mass, and effective scalar coupling to its central density. Our models are accurate at $\sim1\%$ level and speed up large-scale calculations by two orders of magnitude. As an application, we use pySTGROM and Markov-chain Monte Carlo techniques to constrain parameters in the DEF theory, with five well-timed binary pulsars, the binary NS (BNS) inspiral GW170817, and a hypothetical BNS inspiral in the Advanced LIGO and next-generation GW detectors. In the future, as more binary pulsars and BNS mergers are detected, our surrogate models will be helpful in constraining strong-field gravity with essential speed and accuracy.

When a stellar massive compact object, such as a black hole (BH), inspirals into an intermediate massive black hole (IMBH), an intermediate-mass-ratio-inspiral (IMRI) system forms. Such kind of systems are important sources for space-based gravitational wave detectors including LISA, Taiji and Tianqin. Dark matter (DM) minispikes may form around IMBHs. We study the effect of dynamical friction against DM minispike on the evolution of eccentric IMRI. For such investigation we construct the dynamical equations which describes the evolution of eccentric IMRI under the effect of dynamical friction. As dynamical friction is large for small velocity, the dissipation of energy near apastron is more than that near periastron. This will greatly enhance the eccentricity. For instance, with an initial semi-latus rectum of $1\rm AU$, even a moderate DM minispike can make the eccentricity grow from $0.3$ to $0.85$. In the extremal case the eccentricity could be enhanced to near $1$. We also study a specific case which corresponds to an IMRI in the center of a globular cluster (GC) and find the eccentricity can keep its value above $0.95$ until the IMRI enters LISA band. These gravitational wave with enhanced eccentricity by DM minispikes can be easly distinguished from that without DM at $10^{-3}\rm Hz$ due to the eccentricity difference. These anticipations can be tested by future space-based GW detectors such as LISA.

For non-spinning, charged (Reissner-Nordstr\"om) black holes, the particles with an opposite sign of charge with respect to that of the black hole will be pulled into the black hole by the extra electromagnetic force. Such a hole will be quickly neutralized so that there should not exist significantly charged, non-spinning black holes in the universe. The case of spinning, charged (Kerr-Newmann) black holes is more complicated. For a given initial position and initial velocity of the particle, an oppositely charged particle does not always more easily fall into the black hole than a neutral particle. The possible existence of a magnetosphere further complicate the picture. One therefore cannot straightforwardly conclude that a charged spinning black hole will be neutralized. In this paper, we make the first step to investigate the neutralization of Kerr-Newmann black holes without introducing a magnetosphere. We track the particle trajectories under the influence of the curved spacetime and the electromagnetic field carried by the spinning, charged black hole. A statistical method is used to investigate the neutralization problem. We find a universal dependence of the falling probability into the black hole on the charge of the test particle, with the oppositely charged particles having a higher probability of falling. We therefore conclude that charged, spinning black holes without a magnetosphere should be quickly neutralized, consistent with people's intuition. The neutralization problem of Kerr-Newmann black holes with a co-rotating force-free magnetosphere is subject to further studies.

With a selected sample of neutron star (NS) equation-of-states (EOSs) that are consistent with the current observations and have a range of maximum masses, we investigate the relations between NS gravitational mass $M_g$ and baryonic mass $M_b$, and the relations between the maximum NS mass supported through uniform rotation ($M_{\rm max}$) and that of nonrotating NSs ($M_{\rm TOV}$). We find that if one intends to apply an EOS-independent quadratic, universal transformation formula ($M_b=M_g+A\times M_{g}^2$) to all EOSs, the best fit $A$ value is 0.080 for non-rotating NSs only and 0.073 when different spin periods are considered. The residual error of the transformation is as large as $\sim0.1M_{\odot}$. For different EOSs, we find that the parameter $A$ for non-rotating NSs is proportional to $R_{1.4}^{-1}$ (where $R_{1.4}$ is NS radius for 1.4$M_\odot$ in unit of km). For a particular EOS, if one adopts the best-fit parameters for different spin periods, the residual error of the transformation is smaller, which is of the order of 0.01$M_\odot$ for the quadratic form and less than 0.01$M_\odot$ for the cubic form ($M_b=M_g+A_1\times M_{g}^2+A_2\times M_{g}^3$). We also find a very tight and general correlation between the normalized mass gain due to spin $\Delta m\equiv(M_{\rm max}-M_{\rm TOV})/M_{\rm TOV}$ and the spin period normalized to the Keplerian period ${\cal P}$, i.e. ${\rm log_{10}}\Delta m = (-2.74\pm0.05){\rm log_{10}}{\cal P}+{\rm log_{10}}(0.20\pm 0.01)$, which is independent of EOS models. Applications of our results to GW170817 is discussed.

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

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

The models currently used in the detection of gravitational waves (GWs) either do not consider a relative motion between the center-of-mass of the source and the observer, or usually only consider its effect on the frequencies of GWs. However, it is known for light waves that a relative motion not only changes the frequencies but also the brightness of the source, the latter of which is called the `beaming effect'. Here we investigate such an effect for GWs and find that the observed amplitude of a GW signal, unlike the behavior of light, is not a monotonic function of the relative velocity and responds differently to the two GW polarizations. We attribute the difference to a rotation of the wave-vector, as well as a reorientation of the GW polarizations. We find that even for velocities as small as $0.25\%$ of the speed of light, ignoring the aforementioned beaming effect could induce a systematic error that is larger than the designated calibration accuracy of LIGO. This error could lead to an incorrect estimation of the distance and orbital inclination of a GW source, or result in a spurious signal that appears to be incompatible with general relativity.

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.

Retrieving the mass of a gravitational-wave (GW) source is a fundamental but difficult problem because the mass is degenerate with redshift. In astronomy, three types of redshift exist, namely cosmological, Doppler, and gravitational redshift, but the latter two are normally too weak to affect the observation. In this Letter, we show that the current astrophysical models allow binary black holes (BBHs) to merge within $10$ Schwarzschild radii of a supermassive black hole (SMBH). We find that in this case both the Doppler and gravitational redshift are significant, and in the most extreme condition they could increase the "apparent" black-hole mass and distance by a factor of $1.9-3.4$. We show that such a factor is consistent with the distribution in the distance-mass diagram of the ten BBHs detected so far by LIGO/Virgo. We also discuss the difficulties of this redshift scenario caused by the low event rate predicted by the current models, as well the potential solutions.

Gaia data release 2 (DR2) provides the best non-rotating optical frame aligned with the radio frame (ICRF) thanks to the inclusion of about half-million quasars in the 5-parameter astrometric solution. Given their crucial diagnostic role for characterizing the properties of the celestial reference frame, we aim to make an independent assessment of the astrometry of quasars in DR2. We cross-match with Gaia DR2 the quasars from LQAC3, SDSS and LAMOST, obtaining 208743 new sources (denominated as KQCG). With the quasars already identified in DR2, we estimate the global offsets of parallaxes and proper motions of different quasar subsets to check their astrometric consistency. The features of the proper motion field are analyzed by means of vectorial spherical harmonics (VSH); the scalar field of parallaxes is expanded into spherical harmonics to investigate their spatial correlation. We find a bias of $\sim$ $-0.030$ $mas$ in KQCG parallaxes, and a bias of $-9.1$ $\mu as/yr$ in $\mu_{\alpha\ast}$; a $\sim$ +10 $\mu as/yr$ bias in $\mu_{\delta}$ consistently found in the entire quasar sample. The results of the VSH analysis on different subsets indicate good agreement between them. The proper motion field exhibits a ((-6,-5,-5) $\pm$ 1) $\mu as/yr$ rotation in the northern hemisphere and a rotation of ((0,+1,+3) $\pm$ 1 ) $\mu as /yr$ in the southern one. The spherical harmonics expansion of the parallax field reveals an angular scale of systematics of $\approx$ 18 degrees with a RMS amplitude of 13 $\mu$as. The positional comparison shows that the axes of the current Gaia Celestial Reference Frame and the ICRF2 are aligned within 30 $\mu$as.

Context. Quasars are spatially stationary, and they are essential objects in astrometry when defining reference frames. However, the census of quasars is far from complete. Mid-infared colors can be used to find quasar candidates because AGNs show a peculiar appearance in mid-infrared color, but these methods are incapable of separating quasars from AGNs. Aims. The aim of our study is to use astrometric and mid-infrared methods to select quasars and get a reliable quasar candidates catalog. Methods. We used a near-zero proper motion criterion in conjuction with WISE (all-sky Wide-field Infrared Survey Explorer) [W1-W2] color to select quasar candidates. The [W1-W2] color criterion is defined by the linear boundary of two samples: LAMOST DR5 quasars, which serve as the quasar sample, and LAMOST DR5 stars and galaxies, which serve as the non-quasar sample. The contamination and completeness are evaluated. Results. We present a catalog of 662 753 quasar candidates, with a completeness of about 75% and a reliability of 77.2%.

Recently, two of us have found a family of singularity-free rotating black hole solutions in Einstein's conformal gravity. These spacetimes are characterized by three parameters: the black hole mass $M$, the black hole spin angular momentum $J$, and a parameter $L$ that is not specified by the theory but can be expected to be proportional to the black hole mass $M$. The Kerr black hole solution of Einstein's gravity is recovered for $L = 0$. In a previous paper, we showed that X-ray data of astrophysical black holes require $L/M < 1.2$. In the present paper, we report the results of a more sophisticated analysis. We apply the X-ray reflection model {\sc relxill\_nk} to \textsl{NuSTAR} and \textsl{Swift} data of the supermassive black hole in 1H0707-495. We find the constraint $L/M < 0.45$ (90% confidence level).

We study the metallicity distribution and kinematic properties of 4,680 A/F/G/K-type giant stars with $|z|>$ 5 kpc selected from the LAMOST spectroscopic survey. The metallicity distribution of giant stars with 5 $<|z|\leqslant$ 15 kpc can be described by a three-peak Gaussian model with peaks at [Fe/H] $\sim-0.6\pm0.1$, $-1.2\pm0.3$ and $-2.0\pm0.2$, corresponding to the ratio of 19$\%$, 74$\%$ and 7$\%$, respectively. The $\rm{[\alpha/Fe]}$ is used to associate the three peaks with the thick disk, inner-halo and outer-halo components of the Galaxy. The metallicity distribution of these giant stars, which is fit with Gaussians corresponding to the three components, show a growing fraction of inner-halo component and declining fraction of the thick-disk component with increasing distance from the Galactic plane. Adopting a galaxy potential model, we also derive the orbital parameters of the sample stars, such as orbit eccentricity and rotation velocity. The peak values of derived orbital eccentricity for stars covering different metallicity regions maintain $e\sim$ 0.75, independent of height above the plane, within the range 5$<|z|<$ 15 kpc. By comparing the MDFs of stars in different rotation velocity intervals, we find that the majority of the retrograde stars are more metal-poor than the prograde stars.

Einstein's theory of general relativity was proposed over 100 years ago and has successfully passed a large number of observational tests in the weak field regime. However, the strong field regime is largely unexplored, and there are many modified and alternative theories that have the same predictions as Einstein's gravity for weak fields and present deviations when gravity becomes strong. RELXILL_NK is the first relativistic reflection model for probing the spacetime metric in the vicinity of astrophysical black holes and testing Einstein's gravity in the strong field regime. Here we present our current constraints on possible deviations from Einstein's gravity obtained from the black holes in 1H0707-495, Ark 564, GX 339-4, and GS 1354-645.

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.

Quasars are essential for astrometric in the sense that they are spatial stationary because of their large distance from the Sun. The European Space Agency (ESA) space astrometric satellite Gaia is scanning the whole sky with unprecedented accuracy up to a few muas level. However, Gaia's two fields of view observations strategy may introduce a parallax bias in the Gaia catalog. Since it presents no significant parallax, quasar is perfect nature object to detect such bias. More importantly, quasars can be used to construct a Celestial Reference Frame in the optical wavelengths in Gaia mission. In this paper, we compile the most reliable quasars existing in literatures. The final compilation (designated as Known Quasars Catalog for Gaia mission, KQCG) contains 1843850 objects, among of them, 797632 objects are found in Gaia DR1 after cross-identifications. This catalog will be very useful in Gaia mission.

We introduce a new code for cosmological simulations, PHoToNs, which has features on performing massive cosmological simulations on heterogeneous high performance Computer (HPC) and threads oriented programming. PHoToNs adopts a hybrid scheme to compute gravity force, with the conventional PM to compute the long-range force, the Tree algorithm to compute the short range force, and the direct summation PP to compute the gravity from very close particles. A self-similar space filling Peano-Hilbert curve is used to decompose computing domain. Threads programming is highly used to more flexibly manage the domain communication, PM calculation and synchronization, as well as Dual Tree Traversal on the CPU+MIC platform. The scalability of the PHoToNs performs well and the efficiency of PP kernel achieves 68.6% of peak performance on MIC and 74.4% on CPU platforms. We also test the accuracy of the code against the much used Gadget-2 in the community and found excellent agreement.

The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) is the largest existing spectroscopic survey telescope, having 32 scientific charge-coupled-device (CCD) cameras for receiving spectra. Stability and automation of the camera-control software are essential, but cannot be provided by the existing system. The Remote Telescope System 2nd Version (RTS2) is an open-source and automatic observatory-control system. However, all previous RTS2 applications have concerned small telescopes. This paper focuses on implementation of an RTS2-based camera-control system for the 32 CCDs of LAMOST. A virtual camera module inherited from the RTS2 camera module is built as a device component working on the RTS2 framework. To improve the controllability and robustness, a virtualized layer is designed using the master-slave software paradigm, and the virtual camera module is mapped to the 32 real cameras of LAMOST. The new system is deployed in the actual environment and experimentally tested. Finally, multiple observations are conducted using this new RTS2-framework-based control system. The new camera-control system is found to satisfy the requirements for LAMOST automatic camera control. This is the first time that RTS2 has been applied to a large telescope, and provides a referential solution for full RTS2 introduction to the LAMOST observatory control system.

It is discussed that the formation of the spiral galaxies is driven by the cosmic background rotation, not a result of an isolated evolution proposed by the density wave theory. To analyze the motions of the galaxies, a simple double particle galaxy model is considered and the Coriolis force formed by the rotational background is introduced. The numerical analysis shows that not only the trajectory of the particle is the spiral shape, but also the relationship between the velocity and the radius reveals both the existence of spiral arm and the change of the arm number. In addition, the results of the three-dimensional simulation also give the warped structure of the spiral galaxies, and shows that the disc surface of the warped galaxy, like a spinning coins on the table, exists a whole overturning movement. Through the analysis, it can be concluded that the background environment of the spiral galaxies have a large-scale rotation, and both the formation and evolution of hurricane-like spiral galaxies are driven by this background rotation.

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.

Recently, the LIGO-Virgo collaboration reported their first detection of gravitational wave (GW) signals from a low mass compact binary merger GW170817, which is most likely due to a double neutron star (NS) merger. With the GW signals only, the chirp mass of the binary is precisely constrained to $1.188^{+0.004}_{-0.002}~\rm{M_{\odot}}$, but the mass ratio is loosely constrained in the range $0.4-1$, so that a very rough estimation of the individual NS masses ($0.86~{\rm M_{\odot}}<M_1<1.36~\rm{M_{\odot}}$ and $1.36~{\rm M_{\odot}}<M_2<2.26~\rm{M_{\odot}}$) was obtained. Here we propose that if one can constrain the dynamical ejecta mass through performing kilonova modeling of the optical/IR data, by utilizing an empirical relation between the dynamical ejecta mass and the mass ratio of NS binaries, one may place a more stringent constraint on the mass ratio of the system. For instance, considering that the red "kilonova" component is powered by the dynamical ejecta, we reach a tight constraint on the mass ratio in the range of $0.46-0.59$. Alternatively, if the blue "kilonova" component is powered by the dynamical ejecta, the mass ratio would be constrained in the range of $0.53-0.67$. Overall, such a multi-messenger approach could narrow down the mass ratio of GW170817 system to the range of $0.46-0.67$, which gives a more precise estimation of the individual NS mass than pure GW signal analysis, i.e. $0.90~{\rm M_{\odot}}<M_1<1.16~{\rm M_{\odot}}$ and $1.61~{\rm M_{\odot}}<M_2<2.11~{\rm M_{\odot}}$.

Recently, we have extended the X-ray reflection model relxill to test the spacetime metric in the strong gravitational field of astrophysical black holes. In the present Letter, we employ this extended model to analyze XMM-Newton, NuSTAR, and Swift data of the supermassive black hole in 1H0707-495 and test deviations from a Kerr metric parametrized by the Johannsen deformation parameter $\alpha_{13}$. Our results are consistent with the hypothesis that the spacetime metric around the black hole in 1H0707-495 is described by the Kerr solution.

We present the ultraviolet magnitudes for over three million stars in the LAMOST survey, in which 2,202,116 stars are detected by $GALEX$. For 889,235 undetected stars, we develop a method to estimate their upper limit magnitudes. The distribution of (FUV $-$ NUV) shows that the color declines with increasing effective temperature for stars hotter than 7000 K in our sample, while the trend disappears for the cooler stars due to upper atmosphere emission from the regions higher than their photospheres. For stars with valid stellar parameters, we calculate the UV excesses with synthetic model spectra, and find that the (FUV $-$ NUV) vs. $R'_{\mathrm{FUV}}$ can be fitted with a linear relation and late-type dwarfs tend to have high UV excesses. There are 87,178 and 1,498,103 stars detected more than once in the visit exposures of $GALEX$ in the FUV and NUV, respectively. We make use of the quantified photometric errors to determine statistical properties of the UV variation, including intrinsic variability and the structure function on the timescale of days. The overall occurrence of possible false positives is below 1.3\% in our sample. UV absolute magnitudes are calculated for stars with valid parallaxes, which could serve as a possible reference frame in the NUV. We conclude that the colors related to UV provide good criteria to distinguish between M giants and M dwarfs, and the variability of RR Lyrae stars in our sample is stronger than that of other A and F stars.

In this work, we present a catalog of 2651 carbon stars from the fourth Data Release (DR4) of the Large Sky Area Multi-Object Fiber Spectroscopy Telescope (LAMOST). Using an efficient machine-learning algorithm, we find out these stars from more than seven million spectra. As a by-product, 17 carbon-enhanced metal-poor (CEMP) turnoff star candidates are also reported in this paper, and they are preliminarily identified by their atmospheric parameters. Except for 176 stars that could not be given spectral types, we classify the other 2475 carbon stars into five subtypes including 864 C-H, 226 C-R, 400 C-J, 266 C-N, and 719 barium stars based on a series of spectral features. Furthermore, we divide the C-J stars into three subtypes of CJ( H), C-J(R), C-J(N), and about 90% of them are cool N-type stars as expected from previous literature. Beside spectroscopic classification, we also match these carbon stars to multiple broadband photometries. Using ultraviolet photometry data, we find that 25 carbon stars have FUV detections and they are likely to be in binary systems with compact white dwarf companions.

With many large science equipment constructing and putting into use, astronomy has stepped into the big data era. The new method and infrastructure of big data processing has become a new requirement of many astronomers. Cloud computing, Map/Reduce, Hadoop, Spark, etc. many new technology has sprung up in recent years. Comparing to the high performance computing(HPC), Data is the center of these new technology. So, a new computing architecture infrastructure is necessary, which can be shared by both HPC and big data processing. Based on Astronomy Cloud project of Chinese Virtual Observatory (China-VO), we have made much efforts to optimize the designation of the hybrid computing platform. which include the hardware architecture, cluster management, Job and Resource scheduling.

Rotating radio transients (RRATs) are sporadically emitting pulsars detectable only through searches for single pulses. While over 100 RRATs have been detected, only a small fraction (roughly 20\%) have phase-connected timing solutions, which are critical for determining how they relate to other neutron star populations. Detecting more pulses in order to achieve solutions is a key to understanding their physical nature. Astronomical signals collected by radio telescopes contain noise from many sources, making the detection of weak pulses difficult. Applying a denoising method to raw time series prior to performing a single-pulse search typically leads to a more accurate estimation of their times of arrival (TOAs). Taking into account some features of RRAT pulses and noise, we present a denoising method based on wavelet data analysis, an image-processing technique. Assuming that the spin period of an RRAT is known, we estimate the frequency spectrum components contributing to the composition of RRAT pulses. This allows us to suppress the noise, which contributes to other frequencies. We apply the wavelet denoising method including selective wavelet reconstruction and wavelet shrinkage to the de-dispersed time series of eight RRATs with existing timing solutions. The signal-to-noise ratio (S/N) of most pulses are improved after wavelet denoising. Compared to the conventional approach, we measure 12% to 69% more TOAs for the eight RRATs. The new timing solutions for the eight RRATs show 16% to 90% smaller estimation error of most parameters. Thus, we conclude that wavelet analysis is an effective tool for denoising RRATs signal.

We study the classification and composite spectra of galaxy in the fourth data release (DR4) of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST). We select 40,182 spectra of galaxies from LAMOST DR4, which have photometric in- formation but no spectroscopic observations in the Sloan Digital Sky Survey(SDSS). These newly observed spectra are re-calibrated and classified into six classes, i.e. pas- sive, H{\alpha}-weak, star-forming, composite, LINER and Seyfert using the line intensity (H\b{eta}, [OIII]{\lambda}5007, H{\alpha} and [NII]{\lambda}6585). We also study the correlation between spectral classes and morphological types through three parameters: concentration index, (u - r) color, and D4000n index. We calculate composite spectra of high signal-to-noise ra- tio(S/N) for six spectral classes, and using these composites we pick out some features that can differentiate the classes effectively, including H\b{eta}, Fe5015, H{\gamma}A, HK, and Mg2 band etc. In addition, we compare our composite spectra with the SDSS ones and analyse their difference. A galaxy catalogue of 40,182 newly observed spectra (36,601 targets) and the composite spectra of the six classes are available online.

It has been estimated that a significant proportion of binary neutron star merger events produce long-lived massive remnants supported by differential rotation and subject to rotational instabilities. To examine formation and oscillation of rapidly rotating neutron stars (NS) after merger, we present an exploratory study of fully general-relativistic hydrodynamic simulations using the public code Einstein Toolkit. The attention is focused on qualitative aspects of long-term postmerger evolution. As simplified test models, we use a moderately stiff Gamma=2 ideal-fluid equation of state and unmagnetized irrotational equal-mass binaries with three masses well below the threshold for prompt collapse. Our high resolution simulations generate postmerger "ringdown" gravitational wave (GW) signals of 170 ms, sustained by rotating massive NS remnants without collapsing to black holes. We observe that the high-density double-core structure inside the remnants gradually turns into a quasi-axisymmetric toroidal shape. It oscillates in a quasi-periodic manner and shrinks in size due to gravitational radiation. In the GW spectrograms, dominant double peaks persist throughout the postmerger simulations and slowly drift to higher frequencies. A new low-frequency peak emerges at about 100 ms after merger, owing to the growth of GW-driven unstable oscillation modes. The long-term effect of grid resolution is also investigated using the same initial model. Moreover, we comment on physical conditions that are favorable for the transient toroidal configuration to form, and discuss implication of our findings on future GW observation targeting rapidly rotating NSs.

The radial number density and flattening of the Milky Way's stellar halo is measured with $\mathrm{5351}$ metal-poor ([Fe/H]$<-1$) K giants from LAMOST DR3, using a nonparametric method which is model independent and largely avoids the influence of halo substucture. The number density profile is well described by a single power law with index $5.03^{+0.64}_{-0.64}$, and flattening that varies with radius. The stellar halo traced by LAMOST K giants is more flattened at smaller radii, and becomes nearly spherical at larger radii. The flattening, $q$, is about 0.64, 0.8, 0.96 at $r=15$, 20 and 30 kpc (where $r=\sqrt{R^2+\left[Z/q\left(r\right)\right]^2}$), respectively. Moreover, the leading arm of the Sagittarius dwarf galaxy tidal stream in the north, and the trailing arm in the south, are significant in the residual map of density distribution. In addition, an unknown overdensity is identified in the residual map at (R,Z)=(30,15) kpc.

We present the peculiar in-plane velocities derived from the LAMOST red clump stars, which are purified and separated by a novel approach into two groups with different ages. The samples are mostly contributed around the Galactic anti-centre direction such that we are able to map the radial profiles of the radial and azimuthal velocities in the outer disc. From the variations of the in-plane velocities with the Galactocentric radius for the younger and older populations, we find that both radial and azimuthal velocities are not axisymmetric at $8<R<14\,kpc$. The two red clump populations show that the mean radial velocity is negative within $R\sim9\,kpc$ and positive beyond. This is likely because of the perturbation induced by the rotating bar. The cross-zero radius, $R\sim9$\, kpc, essentially indicates the rough location of the outer Lindblad resonance (OLR) radius. Given the circular speed of 238\,km$\rm s^{-1}$, then the pattern speed of the bar can be approximated as $45$\,km$\rm s^{-1}\rm kpc^{-1}$. The young red clump stars show larger mean radial velocity than the old population by about 3$\,km\rm s^{-1}$ between $R\sim9$ and 12\,kpc. This is possibly because the younger population is more sensitive to the perturbation than the older one. The radial profiles of the mean azimuthal velocity for the two populations show an interesting U-shape, i.e. at $R<10.5\,kpc$, the azimuthal velocity declines with $R$ by about 10$\,km\rm s^{-1}$, while at $R>10.5$ it increases with $R$ to 240-245$\,km\rm s^{-1}$. It is not clear why the mean azimuthal velocity shows the U-shape along the Galactic anti-centre direction. Meanwhile, the azimuthal velocity for the younger population is slightly larger than the older one and the difference moderately declines with $R$. Beyond $R\sim12\,kpc$, the azimuthal velocities for the two populations are indistinguishable.

Neutron stars may sustain a non-axisymmetric deformation due to magnetic distortion and are potential sources of continuous gravitational waves (GWs) for ground-based interferometric detectors. With decades of searches using available GW detectors, no evidence of a GW signal from any pulsar has been observed. Progressively stringent upper limits of ellipticity have been placed on Galactic pulsars. In this work, we use the ellipticity inferred from the putative millisecond magnetars in short gamma-ray bursts (SGRBs) to estimate their detectability by current and future GW detectors. For $\sim 1$ ms magnetars inferred from the SGRB data, the detection horizon is $\sim 30$ Mpc and $\sim 600$ Mpc for advanced LIGO (aLIGO) and Einstein Telescope (ET), respectively. Using the ellipticity of SGRB millisecond magnetars as calibration, we estimate the ellipticity and gravitational wave strain of Galactic pulsars and magnetars assuming that the ellipticity is magnetic-distortion-induced. We find that the results are consistent with the null detection results of Galactic pulsars and magnetars with the aLIGO O1. We further predict that the GW signals from these pulsars/magnetars may not be detectable by the currently designed aLIGO detector. The ET detector may be able to detect some relatively low frequency signals ($<50$ Hz) from some of these pulsars. Limited by its design sensitivity, the eLISA detector seems not suitable for detecting the signals from Galactic pulsars and magnetars.

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 present a statistical method to derive the stellar density profiles of the Milky Way from spectroscopic survey data, taking into account selection effects. We assume that the selection function of the spectroscopic survey is based on photometric colors and magnitudes and possibly altered during observations and data reductions. Then the underlying selection function for a line-of-sight can be well recovered by comparing the distribution of the spectroscopic stars in a color-magnitude plane with that of the photometric dataset. Subsequently, the stellar density profile along a line-of-sight can be derived from the spectroscopically measured stellar density profile multiplied by the selection function. The method is validated using Galaxia mock data with two different selection functions. We demonstrate that the derived stellar density profiles well reconstruct the true ones not only for the full targets, but also for the sub-populations selected from the full dataset. Finally, the method is applied to map the density profiles for the Galactic disk and halo, respectively, using the LAMOST RGB stars. The Galactic disk extends to about R=19 kpc, where the disk still contributes about 10% to the total stellar surface density. Beyond this radius, the disk smoothly transitions to the halo without any truncation, bending, or broken. Moreover, no over-density corresponding to the Monoceros ring is found in the Galactic anti-center direction. The disk shows moderate north-south asymmetry at radii larger than 12 kpc. On the other hand, the R-Z tomographic map directly shows that the stellar halo is substantially oblate within a Galactocentric radius of 20 kpc and gradually becomes nearly spherical beyond 30 kpc.

Environment, such as the accretion disk, could modify the signal of the gravitational wave from the astrophysical black hole binaries. In this article, we model the matter field around the intermediate-mass binary black holes by means of an axion-like scalar field and investigate their joint evolution. In details, we consider the equal mass binary black holes surrounded by a shell of axion-like scalar field both in spherical symmetric and non-spherical symmetric cases, and with different strength of the scalar field. Our result shows that the environmental scalar field could essentially modify the dynamics. Firstly, in the spherical symmetric case, with increasing of the scalar field strength, the number of circular orbit of the binary black hole is reduced. It means that the scalar field could significantly accelerate the merger process. Secondly, once the scalar field strength exceeds certain critical value, the scalar field could collapse into a third black hole with its mass being larger than the binary. Consequently, the new black hole collapsed from the environmental scalar field could accrete the binary promptly and the binary collides head-on between each other. In this process, there is almost no any quadrupole signal produced, namely the gravitational wave is greatly suppressed. Thirdly, when the scalar field strength is relatively smaller than the critical value, the black hole orbit could develop eccentricity through the accretion of the scalar field. Fourthly, during the initial stage of the inspire, the gravitational attractive force from the axion-like scalar field could induce a sudden turn in the binary orbits, hence result in a transient wiggle in the gravitational waveform. Finally, in the non-spherical case, the scalar field could gravitationally attract the binary moving toward the mass center of the scalar field and slow down the merger process.

The direct detection of gravitational wave by Laser Interferometer Gravitational-Wave Observatory indicates the coming of the era of gravitational-wave astronomy and gravitational-wave cosmology. It is expected that more and more gravitational-wave events will be detected by currently existing and planned gravitational-wave detectors. The gravitational waves open a new window to explore the Universe and various mysteries will be disclosed through the gravitational-wave detection, combined with other cosmological probes. The gravitational-wave physics is not only related to gravitation theory, but also is closely tied to fundamental physics, cosmology and astrophysics. In this review article, three kinds of sources of gravitational waves and relevant physics will be discussed, namely gravitational waves produced during the inflation and preheating phases of the Universe, the gravitational waves produced during the first-order phase transition as the Universe cools down and the gravitational waves from the three phases: inspiral, merger and ringdown of a compact binary system, respectively. We will also discuss the gravitational waves as a standard siren to explore the evolution of the Universe.

Single-exposure spectra in large spectral surveys are valuable for time domain studies such as stellar variability, but there is no available method to eliminate cosmic rays for single-exposure, multi-fiber spectral images. In this paper, we describe a new method to detect and remove cosmic rays in multi-fiber spectroscopic single exposures. Through the use of two-dimensional profile fitting and a noise model that considers the position-dependent errors, we successfully detect as many as 80% of the cosmic rays and correct the cosmic ray polluted pixels to an average accuracy of 97.8%. Multiple tests and comparisons with both simulated data and real LAMOST data show that the method works properly in detection rate, false detection rate, and validity of cosmic ray correction.

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.

Weyl conformal symmetry can solve the problem the spacetime singularities present in Einstein's gravity. In a recent paper, two of us have found a singularity-free rotating black hole solution in conformal gravity. In addition to the mass $M$ and the spin angular momentum $J$ of the black hole, the new solution has a new parameter, $L$, which here we consider to be proportional to the black hole mass. Since the solution is conformally equivalent to the Kerr metric, photon trajectories are unchanged, while the structure of an accretion disk around a black hole is affected by the value of the parameter $L$. In this paper, we show that X-ray data of astrophysical black holes require $L/M < 1.2$.

We present a sample of quasars discovered in the area of Galactic Anti-Center (GAC) of $150^{\circ} \leq l \leq 210^{\circ}$ and $|b| \leq 30^{\circ}$, based on the LAMOST Data Release 3 (DR3). The sample contains 151 spectroscopically confirmed quasars. Among them 80 are newly discovered with the LAMOST. All those quasars are very bright, with $i$ magnitudes peaking around 17.5\, mag. All the newly quasars are discovered serendipitously, targeted originally with the LAMOST as stars of bluer colours, except for a few targeted as variable, young stellar object candidates. This bright quasar sample at low Galactic latitudes will help fill the gap in the spatial distribution of known quasars near the Galactic disk that are used to construct astrometric reference frame for the purpose of accurate proper motion measurements, for example, Gaia. They are also excellent tracers to probe the kinematics and chemistry of the interstellar medium of the Milky Way disk and halo via absorption line spectroscopy.

We present the second release of value-added catalogues of the LAMOST Spectroscopic Survey of the Galactic Anticentre (LSS-GAC DR2). The catalogues present values of radial velocity $V_{\rm r}$, atmospheric parameters --- effective temperature $T_{\rm eff}$, surface gravity log$g$, metallicity [Fe/H], $\alpha$-element to iron (metal) abundance ratio [$\alpha$/Fe] ([$\alpha$/M]), elemental abundances [C/H] and [N/H], and absolute magnitudes ${\rm M}_V$ and ${\rm M}_{K_{\rm s}}$ deduced from 1.8 million spectra of 1.4 million unique stars targeted by the LSS-GAC since September 2011 until June 2014. The catalogues also give values of interstellar reddening, distance and orbital parameters determined with a variety of techniques, as well as proper motions and multi-band photometry from the far-UV to the mid-IR collected from the literature and various surveys. Accuracies of radial velocities reach 5kms$^{-1}$ for late-type stars, and those of distance estimates range between 10 -- 30 per cent, depending on the spectral signal-to-noise ratios. Precisions of [Fe/H], [C/H] and [N/H] estimates reach 0.1dex, and those of [$\alpha$/Fe] and [$\alpha$/M] reach 0.05dex. The large number of stars, the contiguous sky coverage, the simple yet non-trivial target selection function and the robust estimates of stellar radial velocities and atmospheric parameters, distances and elemental abundances, make the catalogues a valuable data set to study the structure and evolution of the Galaxy, especially the solar-neighbourhood and the outer disk.

Virtual Observatory (VO) is a data-intensively online astronomical research and education environment, which takes advantages of advanced information technologies to achieve seamless and global access to astronomical information. AstroCloud is a cyber-infrastructure for astronomy research initiated by Chinese Virtual Observatory (China-VO) project, and also a kind of physical distributed platform which integrates lots of tasks such as telescope access proposal management, data archiving, data quality control, data release and open access, cloud based data processing and analysis. It consists of five application channels, i.e. observation, data, tools, cloud and public and is acting as a full lifecycle management system and gateway for astronomical data and telescopes. Physically, the platform is hosted in six cities currently, i.e. Beijing, Nanjing, Shanghai, Kunming, Lijiang and Urumqi, and serving more than 17 thousand users. Achievements from international Virtual Observatories and Cloud Computing are adopted heavily. In the paper, backgrounds of the project, architecture, Cloud Computing environment, key features of the system, current status and future plans are introduced.

Using the LAMOST-Gaia common stars, we demonstrate that the in-plane velocity field for the nearby young stars are significantly different from that for the old ones. For the young stars, the probably perturbed velocities similar to the old population are mostly removed from the velocity maps in the $X$--$Y$ plane. The residual velocity field shows that the young stars consistently move along $Y$ with faster $v_\phi$ at the trailing side of the local arm, while at the leading side, they move slower in azimuth direction. At both sides, the young stars averagely move inward with $v_R$ of $-5\sim-3$ km s$^{-1}$. The divergence of the velocity in $Y$ direction implies that the young stars are associated with a density wave nearby the local arm. We therefore suggest that the young stars may reflect the formation of the local spiral arm by correlating themselves with a density wave. The range of the age for the young stars is around 2 Gyr, which is sensible since the transient spiral arm can sustain that long. We also point out that alternative explanations of the peculiar velocity field for the young population cannot be ruled out only from this observed data.

We present a sample of 48 metal-poor galaxies at $ z < 0.14$ selected from 92,510 galaxies in the LAMOST survey. These galaxies are identified for their detection of the auroral emission line \oiii$\lambda$4363 above $3\sigma$ level, which allows a direct measurement of the electron temperature and the oxygen abundance. The emission line fluxes are corrected for internal dust extinction using Balmer decrement method. With electron temperature derived from \oiii$\lambda\lambda4959,5007/\oiii\lambda4363$ and electron density from $\sii\lambda6731/\sii\lambda6717$, we obtain the oxygen abundances in our sample which range from $\zoh= 7.63$ (0.09 $\Zsun$) to $8.46$ (0.6 $\Zsun$). We find an extremely metal-poor galaxy with $\zoh=7.63 \pm 0.01$. With multiband photometric data from FUV to NIR and $\ha$ measurements, we also determine the stellar masses and star formation rates, based on the spectral energy distribution fitting and $\ha$ luminosity, respectively. We find that our galaxies have low and intermediate stellar masses with $\rm 6.39 \le log(M/M_{\sun})\le 9.27$, and high star formation rates (SFRs) with $\rm -2.18 \le log(SFR/M_{\sun} yr^{-1}) \le 1.95$. We also find that the metallicities of our galaxies are consistent with the local $T_e$-based mass-metallicity relation, while the scatter is about 0.28 dex. Additionally, assuming the coefficient of $\rm \alpha =0.66$, we find most of our galaxies follow the local mass-metallicity-SFR relation, while a scatter about 0.24 dex exists, suggesting the mass-metallicity relation is weakly dependent on SFR for those metal-poor galaxies.

A sample of the LAMOST spectra of the early type M0-M3 dwarfs is compared with the Kepler observations. It is found that M dwarfs with strong chromospheric emission in $H_{\alpha}$ have large flare activity in general. The rotational periods derived from the Kepler measurements have close correlations with the sizes of the flares, the power-law distribution index and the equivalent widths of the $H_{\alpha}$ emission. A clear trend exists for higher magnetic activities being detected in faster rotating M dwarfs (rotation periods $<$ 20 day).

The present paper is a sequel to our previous work [Y. Ni et al., JCAP 1607, 049 (2016)] in which we studied the iron K$\alpha$ line expected in the reflection spectrum of Kerr black holes with scalar hair. These metrics are solutions of Einstein's gravity minimally coupled to a massive, complex scalar field. They form a continuous bridge between a subset of Kerr black holes and a family of rotating boson stars depending on one extra parameter, the dimensionless scalar hair parameter $q$, ranging from 0 (Kerr black holes) to 1 (boson stars). Here we study the limiting case $q=1$, corresponding to rotating boson stars. For comparison, spherical boson stars are also considered. We simulate observations with XIS/Suzaku. Using the fact that current observations are well fit by the Kerr solution and thus requiring that acceptable alternative compact objects must be compatible with a Kerr fit, we find that some boson star solutions are relatively easy to rule out as potential candidates to explain astrophysical black holes, while other solutions, which are neither too dilute nor too compact are more elusive and we argue that they cannot be distinguished from Kerr black holes by the analysis of the iron line with current X-ray facilities.

The nearly circular (mean eccentricity <e>~0.06) and coplanar (mean mutual inclination <i>~3 deg) orbits of the Solar System planets motivated Kant and Laplace to put forth the hypothesis that planets are formed in disks, which has developed into the widely accepted theory of planet formation. Surprisingly, the first several hundred extrasolar planets (mostly Jovian) discovered using the Radial Velocity (RV) technique are commonly on eccentric orbits (<e> ~ 0.3). This raises a fundamental question: Are the Solar System and its formation special? The Kepler mission has found thousands of transiting planets dominated by sub-Neptunes, but most of their orbital eccentricities remain unknown. By using the precise spectroscopic host star parameters from the LAMOST observations, we measure the eccentricity distributions for a large (698) and homogeneous Kepler planet sample with transit duration statistics. Nearly half of the planets are in systems with single transiting planets (singles), while the other half are multiple-transiting planets (multiples). We find an eccentricity dichotomy: on average, Kepler singles are on eccentric orbits with <e>~0.3, while the multiples are on nearly circular (<e> = 0.04^{+0.03}_{-0.04}) and coplanar (<i> = 1.4^{+0.8}_{-1.1} deg) orbits similar to the Solar System planets. Our results are consistent with previous studies of smaller samples and individual systems. We also show that Kepler multiples and solar system objects follow a common relation <e>~(1-2)x<i> between mean eccentricities and mutual inclinations. The prevalence of circular orbits and the common relation may imply that the solar system is not so atypical in the galaxy after all.

Based on the LAMOST survey and Sloan Digital Sky Survey (SDSS), we use low-resolution spectra of 130,043 F/G-type dwarf stars to study the kinematics and metallicity properties of the Galactic disk. Our study shows that the stars with poorer metallicity and larger vertical distance from Galactic plane tend to have larger eccentricity and velocity dispersion. After separating the sample stars into likely thin-disk and thick-disk sub-sample, we find that there exits a negative gradient of rotation velocity $V_{\phi}$ with metallicity [Fe/H] for the likely thin-disk sub-sample, and the thick-disk sub-sample exhibit a larger positive gradient of rotation velocity with metallicity. By comparing with model prediction, we consider the radial migration of stars appears to have influenced on the thin-disk formation. In addition, our results shows that the observed thick-disk stellar orbital eccentricity distribution peaks at low eccentricity ($e \sim 0.2$) and extends to a high eccentricity ($e \sim 0.8$). We compare this result with four thick-disk formation simulated models, and it imply that our result is consistent with gas-rich merger model.

The rotation curve (RC) of the Milky Way out to $\sim$ 100 kpc has been constructed using $\sim$ 16,000 primary red clump giants (PRCGs) in the outer disk selected from the LSS-GAC and the SDSS-III/APOGEE survey, combined with $\sim$ 5700 halo K giants (HKGs) selected from the SDSS/SEGUE survey. To derive the RC, the PRCG sample of the warm disc population and the HKG sample of halo stellar population are respectively analyzed using a kinematical model allowing for the asymmetric drift corrections and re-analyzed using the spherical Jeans equation along with measurements of the anisotropic parameter $\beta$ currently available. The typical uncertainties of RC derived from the PRCG and HKG samples are respectively 5-7 km/s and several tens km/s. We determine a circular velocity at the solar position, $V_c (R_0)$ = 240 $\pm$ 6 km/s and an azimuthal peculiar speed of the Sun, $V_{\odot}$ = 12.1 $\pm$ 7.6 km/s, both in good agreement with the previous determinations. The newly constructed RC has a generally flat value of 240 km/s within a Galactocentric distance $r$ of 25 kpc and then decreases steadily to 150 km/s at $r$ $\sim$ 100 kpc. On top of this overall trend, the RC exhibits two prominent localized dips, one at $r$ $\sim$ 11 kpc and another at $r$ $\sim$ 19 kpc. From the newly constructed RC, combined with other constraints, we have built a parametrized mass model for the Galaxy, yielding a virial mass of the Milky Way's dark matter halo of $0.90^{+0.07}_{-0.08} \times 10^{12}$ ${\rm M}_{\odot}$ and a local dark matter density, $\rho_{\rm \odot, dm} = 0.32^{+0.02}_{-0.02}$ GeV cm$^{-3}$.

All of the 14 subfields of the Kepler field have been observed at least once with the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST, Xinglong Observatory, China) during the 2012-2014 observation seasons. There are 88,628 reduced spectra with SNR$_g$ (signal-to-noise ratio in g band) $\geq$ 6 after the first round (2012-2014) of observations for the LAMOST-Kepler project (LK-project). By adopting the upgraded version of the LAMOST Stellar Parameter pipeline (LASP), we have determined the atmospheric parameters ($T_{\rm eff}$ , $\log g$, and $\rm [Fe/H]$) and heliocentric radial velocity $v_{\rm rad}$ for 51,406 stars with 61,226 spectra. Compared with atmospheric parameters derived from both high-resolution spectroscopy and asteroseismology method for common stars in Huber et al. (2014), an external calibration of LASP atmospheric parameters was made, leading to the determination of external errors for the giants and dwarfs, respectively. Multiple spectroscopic observations for the same objects of the LK-project were used to estimate the internal uncertainties of the atmospheric parameters as a function of SNR$_g$ with the unbiased estimation method. The LASP atmospheric parameters were calibrated based on both the external and internal uncertainties for the giants and dwarfs, respectively. A general statistical analysis of the stellar parameters leads to discovery of 106 candidate metal-poor stars, 9 candidate very metal-poor stars, and 18 candidate high-velocity stars. Fitting formulae were obtained segmentally for both the calibrated atmospheric parameters of the LK-project and the KIC parameters with the common stars. The calibrated atmospheric parameters and radial velocities of the LK-project will be useful for studying stars in the Kepler field.

Supernova remnants (SNRs) are believed to be the most important acceleration sites for cosmic rays (CRs) below $\sim10^{15}$ eV in the Galaxy. High energy photons, either directly from the shocks of the SNRs or indirectly from the interaction between SNRs and the nearby clouds, are crucial probes for the CR acceleration. Big progresses on observations of SNRs have been achieved by space- and ground-based $\gamma$-ray facilities. However, whether $\gamma$-rays come from accelerated hadrons or not, as well as their connection with the CRs observed at Earth, remains in debate. Large High Altitude Air Shower Observatory (LHAASO), the next generation experiment, is designed to survey the northern part of the very high energy $\gamma$-ray sky from $\sim 0.3$ TeV to PeV with the sensitivity of $\lesssim1\%$ of the Crab nebula flux. In this paper, we indicate that LHAASO will be dedicated to enlarging the $\gamma$-ray SNR samples and improving the spectral and morphological measurements. These measurements, especially at energies above 30 TeV, will be important for us to finally understand the CR acceleration in SNRs.

Air radon concentration measurement is useful for research on geophysical effects, but it is strongly sensitive to site geology and many geophysical and microclimatic processes such as wind, ventilation, air humidity and so on that induce very big fluctuations on the concentration of radon in air. On the contrary, monitoring the radon concentration in soil by measuring the thermal neutron flux reduces environmental effects. In this paper we report some experimental results on the natural thermal neutron flux as well as the concentration of air radon and its variations at 4300 m a.s.l. These results were obtained with unshielded thermal neutron scintillation detectors (en-detectors) and radon monitors located inside the ARGO-YBJ experimental hall. The correlation of these variations with the lunar month and 1-year period is undoubtedly confirmed. A method for earthquakes prediction provided by a global net of the en-detectors is currently under study.

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 present a catalogue including 11,204 spectra for 10,436 early-type emission-line stars from LAMOST DR2, among which 9,752 early-type emission-line spectra are newly discovered. For these early-type emission-line stars, we discuss the morphological and physical properties from their low-resolution spectra. In this spectral sample, the H$\alpha$ emission profiles display a wide variety of shapes. Based on the H$\alpha$ line profiles, these spectra are categorized into five distinct classes: single-peak emission, single-peak emission in absorption, double-peak emission, double-peak emission in absorption, and P-Cygni profiles. To better understand what causes the H$\alpha$ line profiles, we divide these objects into four types from the view of physical classification, which include classical Be stars, Herbig Ae/Be stars, close binaries and spectra contaminated by HII regions. The majority of Herbig Ae/Be stars and classical Be stars are identified and separated using the (H-K, K-W1) color-color diagram. We also discuss thirty one binary systems as listed in SIMBAD on-line catalogue and identify 3,600 spectra contaminated by HII regions after cross matching with positions in the Dubout-Crillon catalogue. A statistical analysis of line profiles versus classifications is then conducted in order to understand the distribution of H$\alpha$ profiles for each type in our sample. Finally, we also provide a table of 172 spectra with FeII emission lines and roughly calculate stellar wind velocities for seven spectra with P-Cygni profiles.

We present the peculiar in-plane velocities derived from the LAMOST red clump stars. From the variations of the in-plane velocity with the Galactocentric radius for the young and old red clump stars, we are able to identify two types of peculiar velocities: 1) both the two red clump populations show that the radial velocity is negative within $R=9.0$\,kpc and becomes positive beyond (denoted as the \emph{long-wave} mode); and 2) the young red clump stars show larger mean radial velocity than the old population by about 3\,km$\rm s^{-1}$ between $R=9$ and 12\,kpc (denoted as the \emph{short-wave} mode). We find that the elliptical disk induced by the rotating bar can well explain the \emph{long-wave} mode peculiar velocity. The axis ratio of the elliptical disk is around 0.8-0.95 and the disk keeps circular at $R=9.24\pm0.2$\,kpc, which should be the location of the outer Lindblad resonance radius (OLR). Adopting the circular speed of 238\,km$\rm s^{-1}$, the pattern speed of the bar is then derived as $48\pm3$\,km$\rm s^{-1}$kpc$^{-1}$ from the location of OLR. On the other hand, the \emph{short-wave} mode is likely the perturbation of the spiral arms as density waves.

We derive the fraction of substructure in the Galactic halo using a sample of over 10,000 spectroscopically-confirmed halo giant stars from the LAMOST spectroscopic survey. By observing 100 synthetic models along each line of sight with the LAMOST selection function in that sky area, we statistically characterize the expected halo populations. We define as SHARDS (Stellar Halo Accretion Related Debris Structures) any stars in >3-sigma excesses above the model predictions. We find that at least 10% of the Milky Way halo stars from LAMOST are part of SHARDS. By running our algorithm on smooth halos observed with the LAMOST selection function, we show that the LAMOST data contain excess substructure over all Galactocentric radii R_GC < 40 kpc, beyond what is expected due to statistical fluctuations and incomplete sampling of a smooth halo. The level of substructure is consistent with the fraction of stars in SHARDS in model halos created entirely from accreted satellites. This work illustrates the potential of vast spectroscopic surveys with high filling factors over large sky areas to recreate the merging history of the Milky Way.

In the centenary year of Einstein's General Theory of Relativity, this paper reviews the current status of gravitational wave astronomy across a spectrum which stretches from attohertz to kilohertz frequencies. Sect. 1 of this paper reviews the historical development of gravitational wave astronomy from Einstein's first prediction to our current understanding the spectrum. It is shown that detection of signals in the audio frequency spectrum can be expected very soon, and that a north-south pair of next generation detectors would provide large scientific benefits. Sect. 2 reviews the theory of gravitational waves and the principles of detection using laser interferometry. The state of the art Advanced LIGO detectors are then described. These detectors have a high chance of detecting the first events in the near future. Sect. 3 reviews the KAGRA detector currently under development in Japan, which will be the first laser interferometer detector to use cryogenic test masses. Sect. 4 of this paper reviews gravitational wave detection in the nanohertz frequency band using the technique of pulsar timing. Sect. 5 reviews the status of gravitational wave detection in the attohertz frequency band, detectable in the polarisation of the cosmic microwave background, and discusses the prospects for detection of primordial waves from the big bang. The techniques described in sects. 1-5 have already placed significant limits on the strength of gravitational wave sources. Sects. 6 and 7 review ambitious plans for future space based gravitational wave detectors in the millihertz frequency band. Sect. 6 presents a roadmap for development of space based gravitational wave detectors by China while sect. 7 discusses a key enabling technology for space interferometry known as time delay interferometry.

Xinglong Observatory of the National Astronomical Observatories, Chinese Academy of Sciences (NAOC), is one of the major optical observatories in China, which hosts nine optical telescopes including the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) and the 2.16 m reflector. Scientific research from these telescopes is focused on stars, galaxies, and exoplanets using multicolor photometry and spectroscopic observations. Therefore, it is important to provide the observing conditions of the site, in detail, to the astronomers for an efficient use of these facilities. In this article, we present the characterization of observing conditions at Xinglong Observatory based on the monitoring of meteorology, seeing and sky brightness during the period from 2007 to 2014. Results suggest that Xinglong Observatory is still a good site for astronomical observations. Our analysis of the observing conditions at Xinglong Observatory can be used as a reference to the observers on targets selection, observing strategy, and telescope operation.

Part of a review paper entitled "Gravitational wave astronomy: the current status.", appeared in " Science China Physics, Mechanics & Astronomy 58.12 (2015): 1-41.

As the cyber-infrastructure for Astronomical research from Chinese Virtual Observatory (China-VO) project, AstroCloud has been archived solid progresses during the last one year. Proposal management system and data access system are re-designed. Several new sub-systems are developed, including China-VO PaperData, AstroCloud Statics and Public channel. More data sets and application environments are integrated into the platform. LAMOST DR1, the largest astronomical spectrum archive was released to the public using the platform. The latest progresses will be introduced.

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

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

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

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

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

The present work reports on a feasibility study commissioned by the Chinese Academy of Sciences of China to explore various possible mission options to detect gravitational waves in space alternative to that of the eLISA/LISA mission concept. Based on the relative merits assigned to science and technological viability, a few representative mission options descoped from the ALIA mission are considered. A semi-analytic Monte Carlo simulation is carried out to understand the cosmic black hole merger histories starting from intermediate mass black holes at high redshift as well as the possible scientific merits of the mission options considered in probing the light seed black holes and their coevolution with galaxies in early Universe. The study indicates that, by choosing the armlength of the interferometer to be three million kilometers and shifting the sensitivity floor to around one-hundredth Hz, together with a very moderate improvement on the position noise budget, there are certain mission options capable of exploring light seed, intermediate mass black hole binaries at high redshift that are not readily accessible to eLISA/LISA, and yet the technological requirements seem to within reach in the next few decades for China.

Data access and interoperability module connects the observation proposals, data, virtual machines and software. According to the unique identifier of PI (principal investigator), an email address or an internal ID, data can be collected by PI's proposals, or by the search interfaces, e.g. conesearch. Files associated with the searched results could be easily transported to cloud storages, including the storage with virtual machines, or several commercial platforms like Dropbox. Benefitted from the standards of IVOA (International Observatories Alliance), VOTable formatted searching result could be sent to kinds of VO software. Latter endeavor will try to integrate more data and connect archives and some other astronomical resources.

AstroCloud is a cyber-Infrastructure for Astronomy Research initiated by Chinese Virtual Observatory (China-VO) under funding support from NDRC (National Development and Reform commission) and CAS (Chinese Academy of Sciences){\url{http://astrocloud.china-vo.org}}\citep{O8-5_Cui_adassxxiv}. To archive the astronomical data in China, we present the implementation of the astronomical data archiving system (ADAS). Data archiving and quality control are the infrastructure for the AstroCloud. Throughout the data of the entire life cycle, data archiving system standardized data, transferring data, logging observational data, archiving ambient data, And storing these data and metadata in database. Quality control covers the whole process and all aspects of data archiving.

AstroCloud is a cyber-Infrastructure for Astronomy Research initiated by Chinese Virtual Observatory (China-VO) under funding support from NDRC (National Development and Reform commission) and CAS (Chinese Academy of Sciences). The ultimate goal of this project is to provide a comprehensive end-to-end astronomy research environment where several independent systems seamlessly collaborate to support the full lifecycle of the modern observational astronomy based on big data, from proposal submission, to data archiving, data release, and to in-situ data analysis and processing. In this paper, the architecture and key designs of the AstroCloud platform are introduced, including data access middleware, access control and security framework, extendible proposal workflow, and system integration mechanism.

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.

MGRO J2019+37 within the Cygnus region is a bright and extended source revealed by Milagro at 12-35 TeV. This source is almost as bright as Crab Nebula in northern sky, while it is not confirmed by ARGO-YBJ around TeV. Up to now, no obvious counterpart at low energy wavelengths has been found. Hence, MGRO J2019+37 becomes mysterious and its VHE gamma-ray emission mechanism is attractive. In this paper, a brief summary of the multi-wavelength observations from Radio to gamma-ray is presented. All the available data from the XMM-Newton and INTEGRAL at X-ray, and the Fermi-LAT at gamma-ray bands were used to get constraint on its emission flux at low energy wavelengths. Then, its possible counterparts and the VHE emission mechanism are discussed.

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

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

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

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

Searches for radio signatures of ultra-high energy neutrinos and cosmic rays could benefit from improved efficiency by using real-time beamforming or correlation triggering. For missions with power limitations, such as the ANITA-3 Antarctic balloon experiment, full speed high resolution digitization of incoming signals is not practical. To this end, the University of Hawaii has developed the Realtime Independent Three-bit Converter (RITC), a 3-channel, 3-bit, streaming analog-to-digital converter implemented in the IBM-8RF 0.13 um process. RITC is primarily designed to digitize broadband radio signals produced by the Askaryan effect, and thus targets an analog bandwidth of >1 GHz, with a sample-and-hold architecture capable of storing up to 2.6 gigasamples-per-second. An array of flash analog-to-digital converters perform 3-bit conversion of sets of stored samples while acquisition continues elsewhere in the sampling array. A serial interface is provided to access an array of on chip digital-to-analog converters that control the digitization thresholds for each channel as well as the overall sampling rate. Demultiplexed conversion outputs are read out simultaneously for each channel via a set of 36 LVDS links, each running at 650 Mb/s. We briefly describe the design architecture of RITC. Evaluation of the RITC is currently under way, and we will report testing updates as they become available, including prospects for the use of this architecture as the analog half of a novel triggering system for the ANITA-3 ultra-high energy neutrino experiment.

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

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

The ARGO-YBJ experiment is an extensive air shower array with full coverage RPC detectors located at Yangbajing (4300 m asl, Tibet, China). It is operated with high duty cycle (>86%) and a large field of view ($\sim$ 2sr). It continuously monitors the entire overhead sky at $\gamma$-ray energies above 0.1 TeV. In the talk, we will present the result of the northern sky survey (between declinations of -10$^{\circ}$ and 70$^{\circ}$) from an analysis of ~4 years of the ARGO-YBJ data (between July 2006 and February 2011). There are four known TeV sources observed with significance greater than 5 S.D.. The significance from Crab Nebula is more than 16 S.D.. 90% confidence level upper limits to the flux from all directions in the sky are also presented, which vary from 0.09 to 0.53 Crab unit for Crab-like point sources.

Although the roles of data centers and computing centers are becoming more and more important, and on-line research is becoming the mainstream for astronomy, individual research based on locally hosted data is still very common. With the increase of personal storage capacity, it is easy to find hundreds to thousands of FITS files in the personal computer of an astrophysicist. Because Flexible Image Transport System (FITS) is a professional data format initiated by astronomers and used mainly in the small community, data management toolkits for FITS files are very few. Astronomers need a powerful tool to help them manage their local astronomical data. Although Virtual Observatory (VO) is a network oriented astronomical research environment, its applications and related technologies provide useful solutions to enhance the management and utilization of astronomical data hosted in an astronomer's personal computer. FITSManager is such a tool to provide astronomers an efficient management and utilization of their local data, bringing VO to astronomers in a seamless and transparent way. FITSManager provides fruitful functions for FITS file management, like thumbnail, preview, type dependent icons, header keyword indexing and search, collaborated working with other tools and online services, and so on. The development of the FITSManager is an effort to fill the gap between management and analysis of astronomical data.

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

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

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

Gamma ray source detection above 30TeV is an encouraging approach for finding galactic cosmic ray origins. All sky survey for gamma ray sources using wide field of view detector is essential for population accumulation for various types of sources above 100GeV. To target the goals, the ARGO-YBJ experiment has been established. Significant progresses have been made in the experiment. A large air shower detector array in an area of 1km2 is proposed to boost the sensitivity. Hybrid detection with multi-techniques will allow a good discrimination between different types of primary particles, including photons and protons, thus enable an energy spectrum measurement for individual specie. Fluorescence light detector array will extend the spectrum measurement above 100PeV where the second knee is located. An energy scale determined by balloon experiments at 10TeV will be propagated to ultra high energy cosmic ray experiments.

We report studies of ultra-high energy cosmic ray composition via analysis of depth of airshower maximum (Xmax), for airshower events collected by the High Resolution Fly's Eye (HiRes) observatory. The HiRes data are consistent with a constant elongation rate d<Xmax>/d(log(E)) of 47.9 +- 6.0 (stat.) +- 3.2 (syst.)g/cm^2/decade for energies between 1.6 EeV and 63 EeV, and are consistent with a predominantly protonic composition of cosmic rays when interpreted via the QGSJET01 and QGSJET-II high-energy hadronic interaction models. These measurements constrain models in which the galactic-to-extragalactic transition is the cause of the energy spectrum "ankle' at 4 EeV.

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

Stereo data collected by the HiRes experiment over a six year period are examined for large-scale anisotropy related to the inhomogeneous distribution of matter in the nearby Universe. We consider the generic case of small cosmic-ray deflections and a large number of sources tracing the matter distribution. In this matter tracer model the expected cosmic ray flux depends essentially on a single free parameter, the typical deflection angle theta. We find that the HiRes data with threshold energies of 40 EeV and 57 EeV are incompatible with the matter tracer model at a 95% confidence level unless theta is larger than 10 degrees and are compatible with an isotropic flux. The data set above 10 EeV is compatible with both the matter tracer model and an isotropic flux.

The GZK cutoff predicted at the Ultra High Energy Cosmic Ray (UHECR) spectrum as been observed by the HiRes and Auger experiments. The results put severe constraints on the effect of Lorentz Invariance Violation(LIV) which has been introduced to explain the absence of GZK cutoff indicated in the AGASA data. Assuming homogeneous source distribution with a single power law spectrum, we calculate the spectrum of UHECRs observed on Earth by taking the processes of photopion production, $e^+e^-$ pair production and adiabatic energy loss into account. The effect of LIV is also taken into account in the calculation. By fitting the HiRes monocular spectra and the Auger combined spectra, we show that the LIV parameter is constrained to $\xi=-0.8^{+3.2}_{-0.5}\times10^{-23}$ and $0.0^{+1.0}_{-0.4}\times10^{-23}$ respectively, which is well consistent with strict Lorentz Invariance up to the highest energy.

The High Resolution Fly's Eye experiment has measured the flux of ultrahigh energy cosmic rays using the stereoscopic air fluorescence technique. The HiRes experiment consists of two detectors that observe cosmic ray showers via the fluorescence light they emit. HiRes data can be analyzed in monocular mode, where each detector is treated separately, or in stereoscopic mode where they are considered together. Using the monocular mode the HiRes collaboration measured the cosmic ray spectrum and made the first observation of the Greisen-Zatsepin-Kuzmin cutoff. In this paper we present the cosmic ray spectrum measured by the stereoscopic technique. Good agreement is found with the monocular spectrum in all details.

We have searched for correlations between the pointing directions of ultrahigh energy cosmic rays observed by the High Resolution Fly's Eye experiment and Active Galactic Nuclei (AGN) visible from its northern hemisphere location. No correlations, other than random correlations, have been found. We report our results using search parameters prescribed by the Pierre Auger collaboration. Using these parameters, the Auger collaboration concludes that a positive correlation exists for sources visible to their southern hemisphere location. We also describe results using two methods for determining the chance probability of correlations: one in which a hypothesis is formed from scanning one half of the data and tested on the second half, and another which involves a scan over the entire data set. The most significant correlation found occurred with a chance probability of 24%.

Air-fluorescence detectors such as the High Resolution Fly's Eye (HiRes) detector are very sensitive to upward-going, Earth-skimming ultrahigh energy electron-neutrino-induced showers. This is due to the relatively large interaction cross sections of these high-energy neutrinos and to the Landau-Pomeranchuk-Migdal (LPM) effect. The LPM effect causes a significant decrease in the cross sections for bremsstrahlung and pair production, allowing charged-current electron-neutrino-induced showers occurring deep in the Earth's crust to be detectable as they exit the Earth into the atmosphere. A search for upward-going neutrino-induced showers in the HiRes-II monocular dataset has yielded a null result. From an LPM calculation of the energy spectrum of charged particles as a function of primary energy and depth for electron-induced showers in rock, we calculate the shape of the resulting profile of these showers in air. We describe a full detector Monte Carlo simulation to determine the detector response to upward-going electron-neutrino-induced cascades and present an upper limit on the flux of electron-neutrinos.

Measurements are reported of the yield and spectrum of fluorescence, excited by a 28.5 GeV electron beam, in air at a range of pressures of interest to ultra-high energy cosmic ray detectors. The wavelength range was 300 - 420 nm. System calibration has been performed using Rayleigh scattering of a nitrogen laser beam. In atmospheric pressure dry air at 304 K the yield is 20.8 +/- 1.6 photons per MeV.

Based on the standard transfer matrix, a formally exact quantization condition for arbitrary potentials, which outflanks and unifies the historical approaches, is derived. It can be used to find the exact bound-state energy eigenvalues of the quantum system without solving an equation of motion for the system wave functions.

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

Cosmic Ray Tau Neutrino Telescopes (CRTNT) using for sub-EeV cosmic ray measurement is discussed. Performances of a stereoscope configuration with a tower of those telescopes plus two side-triggers are studied. This is done by using a detailed detector simulation driven by Corsika. Detector aperture as a function of shower energy above 10^17 eV is calculated. Event rate of about 20k per year for the second knee measurement is estimated. Event rate for cross calibration with detectors working on higher energy range is also estimated. Different configurations of the detectors are tried for optimization.

Measurements are reported on the fluorescence of air as a function of depth in electromagnetic showers initiated by bunches of 28.5 GeV electrons. The light yield is compared with the expected and observed depth profiles of ionization in the showers. It validates the use of atmospheric fluorescence profiles in measuring ultra high energy cosmic rays.

We have investigated the possibility of $\nu_\tau$ detection using Cosmic Ray Tau Neutrino Telescope (CRTNT) based on air shower fluorescence/Cerenkov light detector techniques. This approach requires an interaction of a $\nu_\tau$ with material such as a mountain or the earth's crust. $\tau$ lepton produced in the charged current interaction must escape from the earth and then decay and initiate a shower in the air. The probability for the conversion from $\nu_\tau$ to air shower has been calculated for an energy range from 1 PeV to 10 EeV. An air shower simulation program has been developed using the simulation package Corsika. The trigger efficiency has been estimated for a CRTNT detector similar to the HiRes/Dice detector in the shadow of Mt. Wheeler in Nevada, USA. A rate of about 8 triggered events per year is expected for the AGN neutrino source model with an optimized configuration and duty cycle of the detector.

The observation of a small change in spectral slope, or 'knee' in the fluxes of cosmic rays near energies 10^15 eV has caused much speculation since its discovery over 40 years ago. The origin of this feature remains unknown. A small workshop to review some modern experimental measurements of this region was held at the Adler Planetarium in Chicago, USA in June 2000. This paper summarizes the results presented at this workshop and the discussion of their interpretation in the context of hadronic models of atmospheric airshowers.

We study the spectrum and average mass composition of cosmic rays with primary energies between 10^{17} eV and 10^{18} eV using a hybrid detector consisting of the High Resolution Fly's Eye (HiRes) prototype and the MIA muon array. Measurements have been made of the change in the depth of shower maximum as a function of energy. A complete Monte Carlo simulation of the detector response and comparisons with shower simulations leads to the conclusion that the cosmic ray intensity is changing f rom a heavier to a lighter composition in this energy range. The spectrum is consistent with earlier Fly's Eye measurements and supports the previously found steepening near 4 \times 10^{17} eV .

We use the CORSIKA air shower simulation program to review the method for assigning energies to ultra-high energy (UHE) cosmic rays viewed with the air fluorescence technique. This technique uses the atmosphere as a calorimeter, and we determine the corrections that must be made to the calorimetric energy to yield the primary cosmic ray energy.

The average mass composition of cosmic rays with primary energies between $10^{17}$eV and $10^{18}$eV has been studied using a hybrid detector consisting of the High Resolution Fly's Eye (HiRes) prototype and the MIA muon array. Measurements have been made of the change in the depth of shower maximum, $X_{max}$, and in the change in the muon density at a fixed core location, $\rho_\mu(600m)$, as a function of energy. The composition has also been evaluated in terms of the combination of $X_{max}$ and $\rho_\mu(600m)$. The results show that the composition is changing from a heavy to lighter mix as the energy increases.