The size of the dust torus in Active Galactic Nuclei (AGN) and their high-luminosity counterparts, quasars, can be inferred from the time delay between UV/optical accretion disk continuum variability and the response in the mid-infrared (MIR) torus emission. This dust reverberation mapping (RM) technique has been successfully applied to $\sim 70$ $z\lesssim 0.3$ AGN and quasars. Here we present first results of our dust RM program for distant quasars covered in the SDSS Stripe 82 region combining $\sim 20$-yr ground-based optical light curves with 10-yr MIR light curves from the WISE satellite. We measure a high-fidelity lag between W1-band (3.4 $\mu$m) and $g$ band for 587 quasars over $0.3\lesssim z\lesssim 2$ ($\left<z\right>\sim 0.8$) and two orders of magnitude in quasar luminosity. They tightly follow (intrinsic scatter $\sim 0.17$ dex in lag) the IR lag-luminosity relation observed for $z<0.3$ AGN, revealing a remarkable size-luminosity relation for the dust torus over more than four decades in AGN luminosity, with little dependence on additional quasar properties such as Eddington ratio and variability amplitude. This study motivates further investigations in the utility of dust RM for cosmology, and strongly endorses a compelling science case for the combined 10-yr Vera C. Rubin Observatory Legacy Survey of Space and Time (optical) and 5-yr Nancy Grace Roman Space Telescope 2$\mu$m light curves in a deep survey for low-redshift AGN dust RM with much lower luminosities and shorter, measurable IR lags. The compiled optical and MIR light curves for 7,384 quasars in our parent sample are made public with this work.

We present optical light curves from the Transiting Exoplanet Survey Satellite (TESS) for the archetypical dwarf AGN in the nearby galaxy NGC 4395 hosting a $\sim 10^5\,M_\odot$ supermassive black hole. Significant variability is detected on timescales from weeks to hours before reaching the background noise level. The $\sim$month-long, 30-minute-cadence, high-precision TESS light curve can be well fit by a simple damped random walk model, with the damping timescale $\tau_{\rm DRW}$ constrained to be $2.3_{-0.7}^{+1.8}$~days ($1\sigma$). \obj\ lies almost exactly on the extrapolation of the $\tau_{\rm DRW}-M_{\rm BH}$ relation measured for AGN with BH masses that are more than three orders of magnitude larger. The optical variability periodogram can be well fit by a broken power-law with the high-frequency slope ($-1.88\pm0.15$) and the characteristic timescale ($\tau_{\rm br}\equiv 1/(2\pi f_{\rm br})=1.4_{-0.5}^{+1.9}\,$days) consistent with the DRW model within 1$\sigma$. This work demonstrates the power of TESS light curves in identifying low-mass accreting SMBHs with optical variability, and a potential global $\tau_{\rm DRW}-M_{\rm BH}$ relation that can be used to estimate SMBH masses with optical variability measurements.

We study the optical light curves - primarily probing the variable emission from the accretion disk - of ~ 900 extreme variability quasars (EVQs, with maximum flux variations more than 1 mag) over an observed-frame baseline of ~ 16 years using public data from the SDSS Stripe 82, PanSTARRS-1 and the Dark Energy Survey. We classify the multi-year long-term light curves of EVQs into three categories roughly in the order of decreasing smoothness: monotonic decreasing or increasing (3.7%), single broad peak and dip (56.8%), and more complex patterns (39.5%). The rareness of monotonic cases suggests that the major mechanisms driving the extreme optical variability do not operate over timescales much longer than a few years. Simulated light curves with a damped random walk model generally under-predict the first two categories with smoother long-term trends. Despite the different long-term behaviors of these EVQs, there is little dependence of the long-term trend on the physical properties of quasars, such as their luminosity, BH mass, and Eddington ratio. The large dynamic range of optical flux variability over multi-year timescales of these EVQs allows us to explore the ensemble correlation between the short-term (< 6 months) variability and the seasonal-average flux across the decade-long baseline (the rms-mean flux relation). We find that unlike the results for X-ray variability studies, the linear short-term flux variations do not scale with the seasonal-average flux, indicating different mechanisms that drive the short-term flickering and long-term extreme variability of accretion disk emission. Finally, we present a sample of 16 EVQs, where the approximately bell-shaped large amplitude variation in the light curve can be reasonably well fit by a simple microlensing model.

Transition Radiation (TR) plays an important role in particle identification in high-energy physics and its characteristics provide a feasible method of energy calibration in the energy range up to 10 TeV, which is of interest for dark matter searches in cosmic rays. In a Transition Radiation Detector (TRD), the TR signal is superimposed onto the ionization energy loss signal induced by incident charged particles. In order to make the TR signal stand out from the background of ionization energy loss in a significant way, we optimized both the radiators and the detector. We have designed a new prototype of regular radiator optimized for a maximal TR photon yield, combined with the Side-On TRD which is supposed to improve the detection efficiency of TR. We started a test beam experiment with the Side-On TRD at Conseil Europ\'{e}en pour la Recherche Nucl\'{e}aire (CERN), and found that the experimental data is consistent with the simulation results.

The Beijing-Arizona Sky Survey (BASS) is a wide and deep imaging survey to cover a 5400 deg$^2$ area in the Northern Galactic Cap with the 2.3m Bok telescope using two filters ($g$ and $r$ bands). The Mosaic $z$-band Legacy Survey (MzLS) covers the same area in $z$ band with the 4m Mayall telescope. These two surveys will be used for spectroscopic targeting of the Dark Energy Spectroscopic Instrument (DESI). The BASS survey observations were completed in 2019 March. This paper describes the third data release (DR3) of BASS, which contains the photometric data from all BASS and MzLS observations between 2015 January and 2019 March. The median astrometric precision relative to {\it Gaia} positions is about 17 mas and the median photometric offset relative to the PanSTARRS1 photometry is within 5 mmag. The median $5\sigma$ AB magnitude depths for point sources are 24.2, 23.6, and 23.0 mag for $g$, $r$, and $z$ bands, respectively. The photometric depth within the survey area is highly homogeneous, with the difference between the 20\% and 80\% depth less than 0.3 mag. The DR3 data, including raw data, calibrated single-epoch images, single-epoch photometric catalogs, stacked images, and co-added photometric catalogs, are publicly accessible at \url{http://batc.bao.ac.cn/BASS/doku.php?id=datarelease:home}.

The characteristic axion couplings could be generated via effective couplings between the Standard Model (SM) fermions to a pseudo-Goldstone from a high-scale $U(1)$ Peccei-Quinn (PQ) symmetry breaking. Assuming that the UV-induced effective operators generate necessary couplings before the PQ symmetry breaking, and any low-scale couplings to the SM are restricted to the Yukawa sector, three minimal natural scenarios can be formulated, which provides a connection between the QCD-axions and mediators at the GUT/string scales. We find that the PQ symmetry breaking scale could be about $10^{15}$ GeV, higher than the classical QCD dark matter axion window but possible if the anthropic window is considered. We also propose an experiment to probe such scenarios. If the dark matter axion is discovered, they might suggest that we live in an atypical Hubble volume.

The broad Mg II line in quasars has distinct variability properties compared with broad Balmer lines: it is less variable, and usually does not display a "breathing" mode, the increase in the average cloud distance when luminosity increases. We demonstrate that these variability properties of Mg II can be reasonably well explained by simple Locally Optimally Emitting Cloud (LOC) photoionization models, confirming earlier photoionization results. In the fiducial LOC model, the Mg II-emitting gas is on average more distant from the ionizing source than the H$\alpha$/H$\beta$ gas, and responds with a lower amplitude to continuum variations. If the broad-line region (BLR) is truncated at a physical radius of $\sim 0.3$ pc (for a $10^{8.5}M_{\odot}$ BH accreting at Eddington ratio of 0.1), most of the Mg II flux will always be emitted near this outer boundary and hence will not display breathing. These results indicate that reverberation mapping results on broad Mg II, while generally more difficult to obtain due to the lower line responsivity, can still be used to infer the Mg II BLR size and hence black hole mass. But it is possible that Mg II does not have a well defined intrinsic BLR size-luminosity relation for individual quasars, even though a global one for the general population may still exist. The dramatic changes in broad H$\alpha$/H$\beta$ emission in the observationally-rare changing-look quasars are fully consistent with photoionization responses to extreme continuum variability, and the LOC model provides natural explanations for the persistence of broad Mg II in changing-look quasars defined on H$\alpha$/H$\beta$, and the rare population of broad Mg II emitters in the spectra of massive inactive galaxies.

Observations by the Wilkinson Microwave Anisotropy Probe and the Planck mission suggest a hemispherical power amplitude asymmetry in the cosmic microwave background, with a correlation length on the order of the size of the observable Universe. We find that this anomaly can be naturally explained by an axion-like particle (ALP) cosmic string formed near our visible Universe. The field variation associated to this cosmic string creates particle density fluctuations after inflation, which consequently decay into radiation before the Big Bang Nucleosynthesis (BBN) era and resulted in the observed power asymmetry. We find in this scenario that the hemispherical power amplitude asymmetry is strongly scale dependent: $A(k)\propto {\rm exp}(-kl)/k$. Admittedly, typical inflation models predict a relic number density of topological defects of order one per observable Universe and so in our model the cosmic string must be tuned to have an impact factor of order $1/H_0$. Interestingly, the constraints based on purely cosmological considerations also give rise to a Peccei-Quinn scale $F_a$ of order $10^3$ larger then the Hubble scale of inflation $H_I$. Assuming $H_I\sim 10^{13}$GeV, we then have an ALP with $F_a\sim 10^{16}$GeV, which coincides with the presumed scale of grand unification. As we require ALP decays occur before the BBN era, which implies a relatively heavy mass or strong self-coupling, and considering that the associated potential should break the shift symmetry softly in order to protect the system from radiative corrections, we also conclude that the required ALP potential should be monodromic in nature.

The EDGES collaboration's observation of an anomalously strong 21 cm absorption feature around the cosmic dawn era has energised the cosmological community by suggesting a novel signature of dark matter in the cooling of cosmic hydrogen. In a recent letter we have argued that by virtue of the ability to mediate cooling processes whilst in the condensed phase, a small amount of axion dark matter can explain these observations within the context of standard models of axions and axion-like particles. These axions and axion-like particles (ALPs) can thermalize through gravitational self-interactions and so eventually form a Bose-Einstein condensate (BEC), whereupon large-scale long-range correlation can produce experimentally observable signals such as these. In this context the EDGES best-fit result favours an axion-like-particles mass in the (10, 450) meV range, which can be compressed for the QCD axion to (100, 450) meV in the absence of fine tuning. Future experiments and galaxy surveys, particularly the International Axion Observatory (IAXO) and EUCLID, should have the capability to directly test this scenario. In this paper, we will explore this mechanism in detail and give more thorough computational details of certain key points.

We report the discovery of an exceptional MIR flare in a Type 2 AGN, SDSS J165726.81+234528.1, at $z=0.059$. This object brightened by 3 mag in the Wide-field Infrared Survey Explorer (WISE) $W1$ and $W2$ bands between 2015 and 2017 (and is fading since 2018), without significant changes ($\lesssim$ 0.2 mag) in the optical over the same period of time. Based on the WISE light curves and near-IR imaging, the flare is more significant at longer wavelengths, suggesting an origin of hot dust emission. The estimated black hole mass ($\sim 10^{6.5}\,M_\odot$) from different methods places its peak bolometric luminosity around the Eddington limit. The high luminosity of the MIR flare and its multi-year timescale suggest that it most likely originated from reprocessed dust radiation in an extended torus surrounding the AGN, instead of from stellar explosions. The MIR color variability is consistent with known changing-look AGN and tidal disruption events (TDEs), but inconsistent with normal supernovae. We suggest that it is a turning-on Type 2 AGN or TDE, where the optical variability is obscured by the dust torus during the transition. This MIR flare event reveals a population of dramatic nuclear transients that are missed in the optical.

We present new Gemini/GMOS optical spectroscopy of 16 extreme variability quasars (EVQs) that dimmed by more than 1.5 mag in the $g$ band between the Sloan Digital Sky Survey (SDSS) and the Dark Energy Survey (DES) epochs (separated by a few years in the quasar rest frame). The quasar sample covers a redshift range of $0.5 < z < 2.1$. Nearly half of these EVQs brightened significantly (by more than 0.5 mag in the $g$ band) in a few years after reaching their previous faintest state, and some EVQs showed rapid (non-blazar) variations of greater than 1-2 mag on timescales of only months. Leveraging on the large dynamic range in continuum variability between the earlier SDSS and the new GMOS spectra, we explore the associated variations in the broad Mg II,$\lambda2798$ line, whose variability properties have not been well studied before. The broad Mg II flux varies in the same direction as the continuum flux, albeit with a smaller amplitude, which indicates at least some portion of Mg II is reverberating to continuum changes. However, the width (FWHM) of Mg II does not vary accordingly as continuum changes for most objects in the sample, in contrast to the case of the broad Balmer lines. Using the width of broad Mg II to estimate the black hole mass therefore introduces a luminosity-dependent bias.

The ubiquitous variability of quasars across a wide range of wavelengths and timescales encodes critical information about the structure and dynamics of the circumnuclear emitting regions that are too small to be directly resolved, as well as detailed underlying physics of accretion and feedback processes in these active supermassive black holes. We emphasize the importance of studying quasar variability with time-domain spectroscopy, focusing on two science cases: (1) reverberation mapping (RM) to measure the broad-line region sizes and black hole masses in distant quasars; (2) spectroscopic follow-up of extreme variability quasars that dramatically change their continuum and broad-line flux within several years. We highlight the need for dedicated optical-infrared spectroscopic survey facilities in the coming decades to accompany wide-area time-domain imaging surveys, including: (1) the next phase of the Sloan Digital Sky Survey (SDSS-V; ~2020-2025), an all-sky, time-domain multi-object spectroscopic survey with 2.5m-class telescopes; (2) the planned Maunakea Spectroscopic Explorer, a dedicated 10m-class spectroscopic survey telescope with a 1.5 sq. deg field-of-view and multiplex of thousands of fibers in both optical and near-IR (J+H) to begin operations in 2029; (3) the Time-domain Spectroscopic Observatory (TSO), a proposed Probe-class ~1.3m telescope at L2, with imaging and spectroscopy (R=200, 1800) in 4 bands (0.3 - 5 micron) and rapid slew capability to 90% of sky, which will extend the coverage of Hbeta to z=8.

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

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

We study the co-evolution of supermassive black holes (SMBHs) with galaxies by means of semi-analytic model (SAM) of galaxy formation based on sub-halo merger trees built from Millennium and Millennium-II simulation. We utilize the simulation results from Guo 2013 and Henriques 2015 to study two aspects of the co-evolution, \emph{i.e.} the stochastic gravitational wave (GW) background generated by SMBH merger and the SMBH/galaxy clustering. The characteristic strain amplitude of GW background predicted by Guo 2013 and Henriques 2015 models are $A_{yr^{-1}}=5.00\times10^{-16}$ and $A_{yr^{-1}}=9.42\times10^{-17}$, respectively. We find the GW amplitude is very sensitive to the galaxy merger rate. The difference in the galaxy merger rate between Guo 2013 and Henriques 2015, results in a factor $5$ deviation in the GW strain amplitude. For clusterings, we calculate the spatially isotropic two point auto- and cross-correlation functions (2PCFs) for both SMBHs and galaxies by using the mock catalogs generated from Guo 2013 model. We find that all 2PCFs have positive dependence on both SMBH and galaxy mass. And there exist a significant time evolution in 2PCFs, namely, the clustering effect is enhanced at lower redshifts. Interestingly, this result is not reported in the active galactic nuclei samples in SDSS. Our analysis also shows that, roughly, SMBHs and galaxies, with galaxy mass $10^2\sim10^3$ larger than SMBH mass, have similar pattern of clustering, which is a reflection of the co-evolution of SMBH and galaxy. Finally, we calculate the first ten multiples of the angular power spectrum of the energy density of GW background. We find the amplitude of angular power spectrum of the first ten multiples is about $10\%$ to $60\%$ of the monopole component in the whole frequency range.

We present the Data Release 4&5 quasar catalog from the quasar survey by Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST), which includes quasars observed between September 2015 and June 2017. There are a total of 19,253 quasars identified by visual inspections of the spectra. Among them, 11,458 are independently discovered by LAMOST, in which 3296 were reported by SDSS DR12 and DR14 quasar catalog after our survey began, while the rest 8162 are new discoveries of LAMOST. We provide the emission line measurements for the Halpha, Hbeta, MgII and/or CIV for 18100 quasars. Since LAMOST does not have absolute flux calibration information, we obtain the monochromatic continuum luminosities by fitting the SDSS photometric data using the quasar spectra, and then estimate the black hole masses. The catalog and spectra for these quasars are available online. This is the third installment in the series of LAMOST quasar survey which has released spectra for totally ~43,000 quasars hitherto. There are 24,772 independently discovered quasars, 17,128 of which are newly discovered. In addition to the great supplement to the new quasar discoveries, LAMOST has also provided a large database (overlapped with SDSS) for investigating the quasar spectral variability and discovering unusual quasars, including changing-look quasars, with ongoing and upcoming large surveys.

We present the final results from our survey of luminous $z \sim $ 5.5 quasars. This is the first systematic quasar survey focusing on quasars at $z \sim$ 5.5, during the post-reionization epoch. It has been challenging to select quasars at $5.3 < z < 5.7$ using conventional color selections, due to their similar optical colors to those of late-type stars, especially M dwarfs. We developed a new selection technique for $z \sim$ 5.5 quasars based on optical, near-IR, and mid-IR photometry, using data from the Sloan Digital Sky Survey (SDSS), PanSTARR1 (PS1), the UKIRT Infrared Deep Sky Surveys - Large Area Survey, the UKIRT Hemisphere Survey, the VISTA Hemisphere Survey, and the Wide Field Infrared Survey Explorer ($WISE$), covering $\sim$ 11000 deg$^2$ of high galactic latitude sky. In this paper, we present the discovery of 15 new quasars at $z\sim 5.5$. Together with results from Yang et al. (2017), our survey provides a complete, flux-limited sample of 31 quasars at $5.3 \le z \le 5.7$. We measure the quasar spatial density at $z \sim 5.5$ and $M_{1450} < -26.2$. Our result is consistent with the rapid decline of quasar spatial density from $z = 5$ to 6, with $k=-0.66 \pm 0.05$ ($\rho(z) \propto 10^{kz}$). In addition, we present a new survey using optical colors only from the full PS1 area for luminous quasars at $z = 5.0 - 5.5$, beyond the SDSS footprint, and report the preliminary results from this survey, including 51 new quasars discovered at $4.61\le z \le5.71$.

The EDGES Collaboration has reported an anomalously strong 21cm absorption feature corresponding to the era of first star formation, which may indirectly betray the influence of dark matter during this epoch. We demonstrate that, by virtue of the ability to mediate cooling processes whilst in the condensed phase, a small amount of axion dark matter can explain these observations within the context of standard models of axions and axion-like-particles. The EDGES best-fit result favours an axion-like-particles mass in the (10, 450) meV range, which can be compressed for the QCD axion to (100, 450) meV in the absence of fine tuning. Future experiments and large scale surveys, particularly the International Axion Observatory (IAXO) and EUCLID, should have the capability to directly test this scenario.

The Beijing-Arizona Sky Survey (BASS) is a new wide-field legacy imaging survey in the northern Galactic cap using the 2.3m Bok telescope. The survey will cover about 5400 deg$^2$ in the $g$ and $r$ bands, and the expected 5$\sigma$ depths (corrected for the Galactic extinction) in the two bands are 24.0 and 23.4 mag, respectively. BASS started observations in January 2015, and has completed about 41% of the whole area as of July 2016. The first data release contains both calibrated images and photometric catalogs obtained in 2015 and 2016. The depths of single-epoch images in the two bands are 23.4 and 22.9 mag, and the full depths of three epochs are about 24.1 and 23.5 mag, respectively.

The rare case of changing-look (CL) AGNs, with the appearance or disappearance of broad Balmer emission lines within a few years, challenges our understanding of the AGN unified model. We present a sample of 21 new CL AGNs at $0.08<z<0.58$, which doubles the number of such objects known to date. These new CL AGNs were discovered by several ways, from (1) repeat spectra in the SDSS, (2) repeat spectra in the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) and SDSS, and (3) photometric variability and new spectroscopic observations. We use the photometric data from surveys, including the SDSS imaging survey, the Pan-STARRS1, the DESI Legacy imaging survey, the Wide-field Infrared Survey Explorer (WISE), the Catalina Real-time Transient Survey, and the Palomar Transient Factory. The estimated upper limits of transition timescale of the CL AGNs in this sample spans from 0.9 to 13 years in the rest frame. The continuum flux in the optical and mid-infrared becomes brighter when the CL AGNs turn on, or vice versa. Variations of more than 0.2 mag in $W1$ band were detected in 15 CL AGNs during the transition. The optical and mid-infrared variability is not consistent with the scenario of variable obscuration in 10 CL AGNs at more than $3\sigma$ confidence level. We confirm a bluer-when-brighter trend in the optical. However, the mid-infrared WISE colors $W1-W2$ become redder when the objects become brighter in the $W1$ band, possibly due to a stronger hot dust contribution in the $W2$ band when the AGN activity becomes stronger. The physical mechanism of type transition is important for understanding the evolution of AGNs.

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

We present the ensemble variability analysis results of quasars using the Dark Energy Camera Legacy Survey (DECaLS) and the Sloan Digital Sky Survey (SDSS) quasar catalogs. Our dataset includes 119,305 quasars with redshifts up to 4.89. Combining the two datasets provides a 15-year baseline and permits analysis of the long timescale variability. Adopting a power-law form for the variability structure function, $V=A(t/1yr)^{\gamma}$, we use the multi-dimensional parametric fitting to explore the relationships between the quasar variability amplitude and a wide variety of quasar properties, including redshift (positive), bolometric luminosity (negative), rest-frame wavelength (negative), and black hole mass (uncertain). We also find that $\gamma$ can be also expressed as a function of redshift (negative), bolometric luminosity (positive), rest-frame wavelength (positive), and black hole mass (positive). Tests of the fitting significance with the bootstrap method show that, even with such a large quasar sample, some correlations are marginally significant. The typical value of $\gamma$ for the entire dataset is $\gtrsim 0.25$, consistent with the results in previous studies on both the quasar ensemble variability and the structure function. A significantly negative correlation between the variability amplitude and the Eddington ratio is found, which may be explained as an effect of accretion disk instability.

The hidden sector U(1) vector bosons created from inflationary fluctuations can be a substantial fraction of dark matter if their mass is around $10^{-5}$eV. The creation mechanism makes the vector bosons' energy spectral density $\rho_{cdm}/\Delta E$ very high. Therefore, the dark electric dipole transition rate in atoms is boosted if the energy gap between atomic states equals the mass of the vector bosons. By using the Zeeman effect, the energy gap between the 2S state and the 2P state in hydrogen atoms or hydrogen like ions can be tuned. The $2S$ state can be populated with electrons due to its relatively long life, which is about $1/7$s. When the energy gap between the semi-ground $2S$ state and the 2P state matches the mass of the cosmic vector bosons, induced transitions occur and the 2P state subsequently decays into the 1S state. The $2P\to1S$ decay emitted Lyman-$\alpha$ photons can then be registered. The choices of target atoms depend on the experimental facilities and the mass ranges of the vector bosons. Because the mass of the vector boson is connected to the inflation scale, the proposed experiment may provide a probe to inflation.

This is the second installment for the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST) Quasar Survey, which includes quasars observed from September 2013 to June 2015. There are 9024 confirmed quasars in DR2 and 10911 in DR3. After cross-match with the SDSS quasar catalogs and NED, 12126 quasars are discovered independently. Among them 2225 quasars were released by SDSS DR12 QSO catalogue in 2014 after we finalised the survey candidates. 1801 sources were identified by SDSS DR14 as QSOs. The remaining 8100 quasars are considered as newly founded, and among them 6887 quasars can be given reliable emission line measurements and the estimated black hole masses. Quasars found in LAMOST are mostly located at low-to-moderate redshifts, with a mean value of 1.5. The highest redshift observed in DR2 and DR3 is 5. We applied emission line measurements to H$\alpha$, H$\beta$, Mg{\sc ii} and C{\sc iv}. We deduced the monochromatic continuum luminosities using photometry data, and estimated the virial black hole masses for the newly discovered quasars. Results are compiled into a quasar catalog, which will be available online.

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

We report the result of a faint quasar survey in a one square degree field. The aim is to test the Y-K/g-z and J-K/i-Y color selection criteria for quasars at faint magnitude, to obtain a complete sample of quasars based on deep optical and near-infrared color-color selection, and to measure the faint end of quasar luminosity function (QLF) over a wide redshift range. We carried out a quasar survey based on the Y-K/g-z and J-K/i-Y quasar selection criteria, using the deep Y-band data obtained from our CFHT/WIRCam Y-band images in a two-degree field within the F22 field of the VIMOS VLT deep survey, optical co-added data from Sloan Digital Sky Survey Stripe 82 and deep near-infrared data from the UKIDSS Deep Extragalactic Survey in the same field. We discovered 25 new quasars at 0.5 < z < 4.5 and i < 22.5 mag within one square degree field. The survey significantly increases the number of faint quasars in this field, especially at z ~ 2-3. It confirms that our color selections are highly complete in a wide redshift range (z < 4.5), especially over the quasar number density peak at z ~ 2-3, even for faint quasars. Combining all previous known quasars and new discoveries, we construct a sample with 109 quasars, and measure the binned QLF and parametric QLF. Although the sample is small, our results agree with a pure luminosity evolution at lower redshift and luminosity evolution and density evolution model at redshift z > 2.5.

We modified the broadband photometric reverberation mapping (PRM) code, JAVELIN, and tested the availability to get broad line region (BLR) time delays that are consistent with the spectroscopic reverberation mapping (SRM) project SDSS-RM. The broadband light curves of SDSS-RM quasars produced by convolution with the system trans- mission curves were used in the test. We found that under similar sampling conditions (evenly and frequently sampled), the key factor determining whether the broadband PRM code can yield lags consistent with the SRM project is the flux ratio of the broad emission line to the reference continuum, which is in line with the previous findings. We further found a critical line-to-continuum flux ratio, about 6%, above which the mean of the ratios between the lags from PRM and SRM becomes closer to unity, and the scatter is pronouncedly reduced. We also tested our code on a subset of SDSS Stripe 82 quasars, and found that our program tends to give biased lag estimations due to the observation gaps when the R-L relation prior in Markov Chain Monte Carlo (MCMC) is discarded. The performance of the damped random walk (DRW) model and the power-law (PL) structure function model on broadband PRM were compared. We found that given both SDSS-RM-like or Stripe 82-like light curves, the DRW model performs better in carrying out broadband PRM than the PL model.

Studies of the most luminous quasars at high redshift directly probe the evolution of the most massive black holes in the early Universe and their connection to massive galaxy formation. However, extremely luminous quasars at high redshift are very rare objects. Only wide area surveys have a chance to constrain their population. The Sloan Digital Sky Survey (SDSS) has so far provided the most widely adopted measurements of the quasar luminosity function (QLF) at $z>3$. However, a careful re-examination of the SDSS quasar sample revealed that the SDSS quasar selection is in fact missing a significant fraction of $z\gtrsim3$ quasars at the brightest end. We have identified the purely optical color selection of SDSS, where quasars at these redshifts are strongly contaminated by late-type dwarfs, and the spectroscopic incompleteness of the SDSS footprint as the main reasons. Therefore we have designed the Extremely Luminous Quasar Survey (ELQS), based on a novel near-infrared JKW2 color cut using WISE AllWISE and 2MASS all-sky photometry, to yield high completeness for very bright ($m_{\rm{i}} < 18.0$) quasars in the redshift range of $3.0\leq z\leq5.0$. It effectively uses random forest machine-learning algorithms on SDSS and WISE photometry for quasar-star classification and photometric redshift estimation. The ELQS will spectroscopically follow-up $\sim 230$ new quasar candidates in an area of $\sim12000\,\rm{deg}^2$ in the SDSS footprint, to obtain a well-defined and complete quasars sample for an accurate measurement of the bright-end quasar luminosity function at $3.0\leq z\leq5.0$. In this paper we present the quasar selection algorithm and the quasar candidate catalog.

We present a new algorithm to estimate quasar photometric redshifts (photo-$z$s), by considering the asymmetries in the relative flux distributions of quasars. The relative flux models are built with multivariate Skew-t distributions in the multi-dimensional space of relative fluxes as a function of redshift and magnitude. For 151,392 quasars in the SDSS, we achieve a photo-$z$ accuracy, defined as the fraction of quasars with the difference between the photo-$z$ $z_p$ and the spectroscopic redshift $z_s$, $|\Delta z| = |z_s-z_p|/(1+z_s)$ within 0.1, of 74%. Combining the WISE W1 and W2 infrared data with the SDSS data, the photo-$z$ accuracy is enhanced to 87%. Using the Pan-STARRS1 or DECaLS photometry with WISE W1 and W2 data, the photo-$z$ accuracies are 79% and 72%, respectively. The prior probabilities as a function of magnitude for quasars, stars and galaxies are calculated respectively based on (1) the quasar luminosity function; (2) the Milky Way synthetic simulation with the Besan\c{c}on model; (3) the Bayesian Galaxy Photometric Redshift estimation. The relative fluxes of stars are obtained with the Padova isochrones, and the relative fluxes of galaxies are modeled through galaxy templates. We test our classification method to select quasars using the DECaLS $g$, $r$, $z$, and WISE W1 and W2 photometry. The quasar selection completeness is higher than 70% for a wide redshift range $0.5<z<4.5$, and a wide magnitude range $18<r<21.5$ mag. Our photo-$z$ regression and classification method has the potential to extend to future surveys. The photo-$z$ code will be publicly available.

Observations indicate that large-scale anomalies exist in the fluctuations of the cosmic microwave background. In these anomalies, the hemispherical power amplitude asymmetry has a correlation length comparable to that of the observable universe. We propose that a topological defect created by spontaneous breaking of the U(1) symmetry prior to inflation generated an initial phase variation, $\delta \theta$, across the observable region of the universe. The amplitude of this phase fluctuation is protected by topology if the defect is inside the horizon, and is frozen by causality if the defect exits the horizon. After inflation, the phase-corresponding boson field started to oscillate, when the Hubble rate decreased to a level comparable to the mass of the boson field. The energy density of the newly created boson particles varied across the observable universe. The bosons subsequently decayed into radiation prior to the BBN epoch, and the resulting fluctuations in the energy density produced the observed power asymmetry. This scenario predicts a scale-dependent modulation amplitude power asymmetry and in addition, as topological defects created by phase transitions are a very general phenomenon, the observed hemispherical asymmetry may be seen as an evidence for the cosmological inflation.

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 explore a scenario that the dark matter is a boson condensate created by the misalignment mechanism, in which a spin 0 boson (an axion-like particle) and a spin 1 boson (the dark photon) are considered, respectively. We find that although the sub-MeV dark matter boson is extremely stable, the huge number of dark matter particles in a galaxy halo makes the decaying signal detectable. A galaxy halo is a large structure bounded by gravity with a typical $\sim10^{12}$ solar mass, and the majority of its components are made of dark matter. For the axion-like particle case, it decays via $\phi\to \gamma\gamma$, therefore the photon spectrum is monochromatic. For the dark photon case, it is a three body decay $A'\to\gamma\gamma\gamma$. However, we find that the photon spectrum is heavily peaked at $M/2$ and thus can facilitate observation. We also suggest a physical explanation for the three body decay spectrum by comparing the physics in the decay of orthopositronium. In addition, for both cases, the decaying photon flux can be measured for some regions of parameter space using current technologies.

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 ARGO-YBJ detector, located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China), was a full coverage air shower array dedicated to gamma ray astronomy and cosmic ray studies. The wide field of view (~ 2 sr) and high duty cycle (> 86%), made ARGO-YBJ suitable to search for short and unexpected gamma ray emissions like gamma ray bursts (GRBs). Between 2007 November 6 and 2013 February 7, 156 satellite-triggered GRBs (24 of them with known redshift) occurred within the ARGO-YBJ field of view. A search for possible emission associated to these GRBs has been made in the two energy ranges 10-100 GeV and 10-1000 GeV. No significant excess has been found in time coincidence with the satellite detections nor in a time window of one hour after the bursts. Taking into account the EBL absorption, upper limits to the energy fluence at 99% of confidence level have been evaluated,with values ranging from ~ 10-5 erg cm-2 to ~10-1 erg cm-2.

We present the first discoveries from a survey of $z\gtrsim6$ quasars using imaging data from the DECam Legacy Survey (DECaLS) in the optical, the UKIRT Deep Infrared Sky Survey (UKIDSS) and a preliminary version of the UKIRT Hemisphere Survey (UHS) in the near-IR, and ALLWISE in the mid-IR. DECaLS will image 9000 deg$^2$ of sky down to $z_{\rm AB}\sim23.0$, and UKIDSS and UHS, which will map the northern sky at $0<DEC<+60^{\circ}$, reaching $J_{\rm VEGA}\sim19.6$ (5-$\sigma$). The combination of these datasets allows us to discover quasars at redshift $z\gtrsim7$ and to conduct a complete census of the faint quasar population at $z\gtrsim6$. In this paper, we report on the selection method of our search, and on the initial discoveries of two new, faint $z\gtrsim6$ quasars and one new $z=6.63$ quasar in our pilot spectroscopic observations. The two new $z\sim6$ quasars are at $z=6.07$ and $z=6.17$ with absolute magnitudes at rest-frame wavelength 1450 \AA\ being $M_{1450}=-25.83$ and $M_{1450}=-25.76$, respectively. These discoveries suggest that we can find quasars close to or fainter than the break magnitude of the Quasar Luminosity Function (QLF) at $z\gtrsim6$. The new $z=6.63$ quasar has an absolute magnitude of $M_{1450}=-25.95$. This demonstrates the potential of using the combined DECaLS and UKIDSS/UHS datasets to find $z\gtrsim7$ quasars. Extrapolating from previous QLF measurements, we predict that these combined datasets will yield $\sim200$ $z\sim6$ quasars to $z_{\rm AB} < 21.5$, $\sim1{,}000$ $z\sim6$ quasars to $z_{\rm AB}<23$, and $\sim 30$ quasars at $z>6.5$ to $J_{\rm VEGA}<19.5$.

We present initial results from the first systematic survey of luminous $z\sim 5.5$ quasars. Quasars at $z \sim$ 5.5, the post-reionization epoch, are crucial tools to explore the evolution of intergalactic medium, quasar evolution and the early super-massive black hole growth. However, it has been very challenging to select quasars at redshifts 5.3 $\le z \le$ 5.7 using conventional color selections, due to their similar optical colors to late-type stars, especially M dwarfs, resulting in a glaring redshift gap in quasar redshift distributions. We develop a new selection technique for $z \sim$ 5.5 quasars based on optical, near-IR and mid-IR photometric data from Sloan Digital Sky Survey (SDSS), UKIRT InfraRed Deep Sky Surveys - Large Area Survey (ULAS), VISTA Hemisphere Survey (VHS) and Wide field Infrared Survey Explorer (WISE). From our pilot observations in SDSS-ULAS/VHS area, we have discovered 15 new quasars at 5.3 $\le z \le$ 5.7 and 6 new lower redshift quasars, with SDSS z band magnitude brighter than 20.5. Including other two $z \sim$ 5.5 quasars already published in our previous work, we now construct an uniform quasar sample at 5.3 $\le z \le$ 5.7 with 17 quasars in a $\sim$ 4800 square degree survey area. For further application in a larger survey area, we apply our selection pipeline to do a test selection by using the new wide field J band photometric data from a preliminary version of the UKIRT Hemisphere Survey (UHS). We successfully discover the first UHS selected $z \sim$ 5.5 quasar.

Very few low-ionization broad absorption line (LoBAL) QSOs have been found at high redshifts to date. One high-redshift LoBAL QSO, J0122+1216, was recently discovered at the Lijiang 2.4-m Telescope with an initial redshift determination of 4.76. Aiming to investigate its physical properties, we carried out follow-up observations in the optical and near-IR spectroscopy. Near-IR spectra from UKIRT and P200 confirms that it is a LoBAL, with a new redshift determination of $4.82\pm0.01$ based on the \mgii~ emission-line. The new \mgii~ redshift determination reveals strong blueshifts and asymmetry of the high-ionization emission lines. We estimated a black hole mass of $\sim 2.3\times 10^9 M_\odot$ and Eddington ratio of $\sim 1.0$ according to the empirical \mgii-based single-epoch relation and bolometric correction factor. It is possible that strong outflows are the result of an extreme quasar environment driven by the high Eddington ratio. A lower limit on the outflowing kinetic power ($>0.9\% L_{Edd}$) was derived from both emission and absorption lines, indicating these outflows play a significant role in the feedback process to regulate the growth of its black hole as well as host galaxy evolution.

Dark matter constitutes about $23\%$ of the total energy density of the universe but its properties are still little known besides that it should be composed by cold and weakly interacting particles. Many beyond standard model theories can provide proper candidates to serve as dark matter and the axion introduced to solve the strong CP problem turns out to be an attractive one. In this paper, we briefly review several important features of the axion and the axion dark matter.

If axion-like particles (ALPs) constitute a major part of dark matter (DM), due to the bosonic nature and a relative small mass, they could behave differently from the point like dark matter particles on the formation of the cosmic structures. When studying the structure formation, it is often useful to consider DM as a special fluid with a given density and a given velocity. ALP fluid obeys a same continuity equation comparing to the point-like collisionless DM but has a different first-order velocity equation. In the linear regime of structure formation, the resulted observational differences are negligible for the QCD axions but can be interesting for very light ALPs.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

We present the discovery of nine quasars at $z\sim6$ identified in the Sloan Digital Sky Survey (SDSS) imaging data. This completes our survey of $z\sim6$ quasars in the SDSS footprint. Our final sample consists of 52 quasars at $5.7<z\le6.4$, including 29 quasars with $z_{\rm AB}\le20$ mag selected from 11,240 deg$^2$ of the SDSS single-epoch imaging survey (the main survey), 10 quasars with $20\le z_{\rm AB}\le20.5$ selected from 4223 deg$^2$ of the SDSS overlap regions (regions with two or more imaging scans), and 13 quasars down to $z_{\rm AB}\approx22$ mag from the 277 deg$^2$ in Stripe 82. They span a wide luminosity range of $-29.0\le M_{1450}\le-24.5$. This well-defined sample is used to derive the quasar luminosity function (QLF) at $z\sim6$. After combining our SDSS sample with two faint ($M_{1450}\ge-23$ mag) quasars from the literature, we obtain the parameters for a double power-law fit to the QLF. The bright-end slope $\beta$ of the QLF is well constrained to be $\beta=-2.8\pm0.2$. Due to the small number of low-luminosity quasars, the faint-end slope $\alpha$ and the characteristic magnitude $M_{1450}^{\ast}$ are less well constrained, with $\alpha=-1.90_{-0.44}^{+0.58}$ and $M^{\ast}=-25.2_{-3.8}^{+1.2}$ mag. The spatial density of luminous quasars, parametrized as $\rho(M_{1450}<-26,z)=\rho(z=6)\,10^{k(z-6)}$, drops rapidly from $z\sim5$ to 6, with $k=-0.72\pm0.11$. Based on our fitted QLF and assuming an IGM clumping factor of $C=3$, we find that the observed quasar population cannot provide enough photons to ionize the $z\sim6$ IGM at $\sim90$\% confidence. Quasars may still provide a significant fraction of the required photons, although much larger samples of faint quasars are needed for more stringent constraints on the quasar contribution to reionization.

The active galactic nucleus (AGN) NGC 7213 shows a complex correlation between the monochromatic radio luminosity $L_R$ and the 2--10 keV X-ray luminosity $L_X$, i.e. the correlation is unusually weak with $p\sim 0$ (in the form $L_R\propto L_X^p$) when $L_X$ is below a critical luminosity, and steep with $p>1$ when $L_X$ is above that luminosity. Such a hybrid correlation in individual AGNs is unexpected as it deviates from the Fundamental Plane of AGN activity. Interestingly, a similar correlation pattern is observed in the black hole X-ray binary H1743--322, where it has been modelled by switching between different modes of accretion. We propose that the flat $L_R$--$L_X$ correlation of NGC 7213 is due to the presence of a luminous hot accretion flow, an accretion model whose radiative efficiency is sensitive to the accretion rate. Given the low luminosity of the source, $L_X\sim 10^{-4}$ of the Eddington luminosity, the viscosity parameter is determined to be small, $\alpha\approx 0.01$. We also modelled the broad-band spectrum from radio to $\gamma$-rays, the time lag between the radio and X-ray light curves, and the implied size and the Lorentz factor of the radio jet. We predict that NGC 7213 will enter into a two-phase accretion regime when $L_X > 1.5 \times 10^{42}\, {\rm erg\,s^{-1}}$. When this happens, we predict a softening of the X-ray spectrum with the increasing flux and a steep radio/X-ray correlation.

Luminous quasars at z>5.6 can be studied in detail with the current generation of telescopes and provide us with unique information on the first gigayear of the universe. Thus far these studies have been statistically limited by the number of quasars known at these redshifts. Such quasars are rare and therefore wide-field surveys are required to identify them and multiwavelength data are needed to separate them efficiently from their main contaminants, the far more numerous cool dwarfs. In this paper, we update and extend the selection for z~6 quasars presented in Banados et al. (2014) using the Pan-STARRS1 (PS1) survey. We present the PS1 distant quasar sample, which currently consists of 124 quasars in the redshift range 5.6<z<6.7 that satisfy our selection criteria. Seventy-seven of these quasars have been discovered with PS1, and 63 of them are newly identified in this paper. We present composite spectra of the PS1 distant quasar sample. This sample spans a factor of ~20 in luminosity and shows a variety of emission line properties. The number of quasars at z>5.6 presented in this work almost double the quasars previously known at these redshifts, marking a transition phase from studies of individual sources to statistical studies of the high-redshift quasar population, which was impossible with earlier, smaller samples.

This is the second paper in a series on a new luminous z ~ 5 quasar survey using optical and near-infrared colors. Here we present a new determination of the bright end of the quasar luminosity function (QLF) at z ~ 5. Combined our 45 new quasars with previously known quasars that satisfy our selections, we construct the largest uniform luminous z ~ 5 quasar sample to date, with 99 quasars in the range 4.7 <= z < 5.4 and -29 < M1450 <= -26.8, within the Sloan Digital Sky Survey (SDSS) footprint. We use a modified 1/Va method including flux limit correction to derive a binned QLF, and we model the parametric QLF using maximum likelihood estimation. With the faint-end slope of the QLF fixed as alpha = -2.03 from previous deeper samples, the best fit of our QLF gives a flatter bright end slope beta = -3.58+/-0.24 and a fainter break magnitude M*1450 = -26.98+/-0.23 than previous studies at similar redshift. Combined with previous work at lower and higher redshifts, our result is consistent with a luminosity evolution and density evolution (LEDE) model. Using the best fit QLF, the contribution of quasars to the ionizing background at z ~ 5 is found to be 18% - 45% with a clumping factor C of 2 - 5. Our sample suggests an evolution of radio loud fraction with optical luminosity but no obvious evolution with redshift.

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.

High-redshift quasars are important tracers of structure and evolution in the early universe. However, they are very rare and difficult to find when using color selection because of contamination from late-type dwarfs. High-redshift quasar surveys based on only optical colors suffer from incompleteness and low identification efficiency, especially at $z\gtrsim4.5$. We have developed a new method to select $4.7\lesssim z \lesssim 5.4$ quasars with both high efficiency and completeness by combining optical and mid-IR Wide-field Infrared Survey Explorer (WISE) photometric data, and are conducting a luminous $z\sim5$ quasar survey in the whole Sloan Digital Sky Survey (SDSS) footprint. We have spectroscopically observed 99 out of 110 candidates with $z$-band magnitudes brighter than 19.5 and 64 (64.6\%) of them are quasars with redshifts of $4.4\lesssim z \lesssim 5.5$ and absolute magnitudes of $-29\lesssim M_{1450} \lesssim -26.4$. In addition, we also observed 14 fainter candidates selected with the same criteria and identified 8 (57.1\%) of them as quasars with $4.7<z<5.4$ . Among 72 newly identified quasars, 12 of them are at $5.2 < z < 5.7$, which leads to an increase of $\sim$36\% of the number of known quasars at this redshift range. More importantly, our identifications doubled the number of quasars with $M_{1450}<-27.5$ at $z>4.5$, which will set strong constraints on the bright end of the quasar luminosity function. We also expand our method to select quasars at $z\gtrsim5.7$. In this paper we report the discovery of four new luminous $z\gtrsim5.7$ quasars based on SDSS-WISE selection.

The ultraviolet broad absorption lines have been seen in the spectra of quasars at high redshift, and are generally considered to be caused by outflows with velocities from thousands kilometers per second to one tenth of the speed of light. They provide crucial implications for the cosmological structures and physical evolutions related to the feedback of active galactic nuclei (AGNs). Recently, through a dedicated program of optically spectroscopic identifications of selected quasar candidates at redshift 5 by using the Lijiang 2.4 m telescope, we discovered two luminous broad absorption line quasars (BALQSOs) at redshift about 4.75. One of them may even have the potentially highest absorption Balnicity Index (BI) ever found to date, which is remarkably characterized by its deep, broad absorption lines and sub-relativistic outflows. Further physical properties, including the metal abundances, variabilities, evolutions of the supermassive black holes (SMBH) and accretion disks associated with the feedback process, can be investigated with multi-wavelength follow-up observations in the future.

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

We present preliminary results of the quasar survey in Large Sky Area Multi- Object Fiber Spectroscopic Telescope (LAMOST) first data release (DR1), which includes pilot survey and the first year regular survey. There are 3921 quasars identified with reliability, among which 1180 are new quasars discovered in the survey. These quasars are at low to median redshifts, with highest z of 4.83. We compile emission line measurements around the H{\alpha}, H{\beta}, Mg II, and C IV regions for the new quasars. The continuum luminosities are inferred from SDSS photo- metric data with model fitting as the spectra in DR1 are non-flux-calibrated. We also compile the virial black hole mass estimates, and flags indicating the selec- tion methods, broad absorption line quasars. The catalog and spectra for these quasars are available online. 28% of the 3921 quasars are selected with optical- infrared colours independently, indicating that the method is quite promising in completeness of quasar survey. LAMOST DR1 and the on-going quasar survey will provide valuable data in the studies of quasars.

The radiative mechanism of the black hole X-ray transients (BHXTs) in their quiescent states (defined as the 2-10 keV X-ray luminosity $\le 10^{34}{\rm \,erg\,s^{-1}}$) remains unclear. In this work, we investigate the quasi-simultaneous quiescent state spectrum (including radio, infrared, optical, ultraviolet and X-ray) of two BHXTs, A0620--00 and XTE J1118+480. We find that, these two sources can be well described by a coupled accretion -- jet model. More specifically, most of the emission (radio up-to infrared, and the X-ray waveband) comes from the collimated relativistic jet. Emission from hot accretion flow is totally insignificant, and it can only be observed in mid-infrared (the synchrotron peak). Emission from the outer cold disc is only evident in UV band. These results are consistent with our previous investigation on the quiescent state of V404 Cyg, and confirm that the quiescent state is jet-dominated.

Axions and axion like particles are very attractive dark matter candidates. In this review, we briefly investigate how the cosmological observations reveal the existence of dark matter and some unique properties of axions/axion like particles which make them more interesting.

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.

We report the discovery of an ultra-luminous quasar J030642.51+185315.8 (hereafter J0306+1853) at redshift 5.363, which hosts a super-massive black hole (SMBH) with $M_{BH} = (1.07 \pm 0.27) \times10^{10}~M_\odot$. With an absolute magnitude $M_{1450}=-28.92$ and bolometric luminosity $L_{bol}\sim3.4\times10^{14} L_{\odot}$, J0306+1853 is one of the most luminous objects in the early Universe. It is not likely to be a beamed source based on its small flux variability, low radio loudness and normal broad emission lines. In addition, a $z=4.986$ Damped Ly$\alpha$ system (DLA) with $\rm [M/H]=-1.3\pm0.1$, among the most metal rich DLAs at $z \gtrsim 5$, is detected in the absorption spectrum of this quasar. This ultra-luminous quasar puts strong constraint on the bright-end of quasar luminosity function and massive-end of black hole mass function. It will provide a unique laboratory to the study of BH growth and the co-evolution between BH and host galaxy with multi-wavelength follow-up observations. The future high resolution spectra will give more insights to the DLA and other absorption systems along the line-of-sight of J0306+1853.

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.

So far, roughly 40 quasars with redshifts greater than z=6 have been discovered. Each quasar contains a black hole with a mass of about one billion solar masses ($10^9 M_\odot$). The existence of such black holes when the Universe was less than 1 billion years old presents substantial challenges to theories of the formation and growth of black holes and the coevolution of black holes and galaxies. Here we report the discovery of an ultra-luminous quasar, SDSS J010013.02+280225.8, at redshift z=6.30. It has an optical and near-infrared luminosity a few times greater than those of previously known z>6 quasars. On the basis of the deep absorption trough on the blue side of the Ly $\alpha$ emission line in the spectrum, we estimate the proper size of the ionized proximity zone associated with the quasar to be 26 million light years, larger than found with other z>6.1 quasars with lower luminosities. We estimate (on the basis of a near-infrared spectrum) that the black hole has a mass of $\sim 1.2 \times 10^{10} M_\odot$, which is consistent with the $1.3 \times 10^{10} M_\odot$ derived by assuming an Eddington-limited accretion rate.

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

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

We describe directed searches for continuous gravitational waves in data from the sixth LIGO science data run. The targets were nine young supernova remnants not associated with pulsars; eight of the remnants are associated with non-pulsing suspected neutron stars. One target's parameters are uncertain enough to warrant two searches, for a total of ten. Each search covered a broad band of frequencies and first and second frequency derivatives for a fixed sky direction. The searches coherently integrated data from the two LIGO interferometers over time spans from 5.3-25.3 days using the matched-filtering F-statistic. We found no credible gravitational-wave signals. We set 95% confidence upper limits as strong (low) as $4\times10^{-25}$ on intrinsic strain, $2\times10^{-7}$ on fiducial ellipticity, and $4\times10^{-5}$ on r-mode amplitude. These beat the indirect limits from energy conservation and are within the range of theoretical predictions for neutron-star ellipticities and r-mode amplitudes.

We investigate the observed correlation between the 2--10 keV X-ray luminosity (in unit of the Eddington luminosity; $l_X \equiv L_X/L_{Edd}$) and the photon index ($\Gamma$) of the X-ray spectrum for both black hole X-ray binaries (BHBs) and active galactic nuclei (AGNs). We construct a large sample, with $10^{-9} < l_X < 10^{-1}$. We find that $\Gamma$ is positively and negatively correlated with $l_X$ when $l_X > 10^{-3}$ and $10^{-6.5} < l_X < 10^{-3}$ respectively, while $\Gamma$ is nearly a constant when $l_X < 10^{-6.5}$. We explain the above correlation in the framework of a coupled hot accretion flow -- jet model. The radio emission always come from the jet while the X-ray emission comes from the accretion flow and jet when $l_X$ is above and below $10^{-6.5}$, respectively. More specifically, we assume that with the increase of mass accretion rate, the hot accretion flow develops into a clumpy and further a disc -- corona two-phase structure because of thermal instability. We argue that such kind of two-phase accretion flow can explain the observed positive correlation.

We present results of a search for continuously-emitted gravitational radiation, directed at the brightest low-mass X-ray binary, Scorpius X-1. Our semi-coherent analysis covers 10 days of LIGO S5 data ranging from 50-550 Hz, and performs an incoherent sum of coherent $\mathcal{F}$-statistic power distributed amongst frequency-modulated orbital sidebands. All candidates not removed at the veto stage were found to be consistent with noise at a 1% false alarm rate. We present Bayesian 95% confidence upper limits on gravitational-wave strain amplitude using two different prior distributions: a standard one, with no a priori assumptions about the orientation of Scorpius X-1; and an angle-restricted one, using a prior derived from electromagnetic observations. Median strain upper limits of 1.3e-24 and 8e-25 are reported at 150 Hz for the standard and angle-restricted searches respectively. This proof of principle analysis was limited to a short observation time by unknown effects of accretion on the intrinsic spin frequency of the neutron star, but improves upon previous upper limits by factors of ~1.4 for the standard, and 2.3 for the angle-restricted search at the sensitive region of the detector.

In this paper we present the results of a coherent narrow-band search for continuous gravitational-wave signals from the Crab and Vela pulsars conducted on Virgo VSR4 data. In order to take into account a possible small mismatch between the gravitational wave frequency and two times the star rotation frequency, inferred from measurement of the electromagnetic pulse rate, a range of 0.02 Hz around two times the star rotational frequency has been searched for both the pulsars. No evidence for a signal has been found and 95$\%$ confidence level upper limits have been computed both assuming polarization parameters are completely unknown and that they are known with some uncertainty, as derived from X-ray observations of the pulsar wind torii. For Vela the upper limits are comparable to the spin-down limit, computed assuming that all the observed spin-down is due to the emission of gravitational waves. For Crab the upper limits are about a factor of two below the spin-down limit, and represent a significant improvement with respect to past analysis. This is the first time the spin-down limit is significantly overcome in a narrow-band search.

Only very few z>5 quasars discovered to date are radio-loud, with a radio-to-optical flux ratio (radio-loudness parameter) higher than 10. Here we report the discovery of an optically luminous radio-loud quasar, SDSS J013127.34-032100.1 (J0131-0321 in short), at z=5.18+-0.01 using the Lijiang 2.4m and Magellan telescopes. J0131-0321 has a spectral energy distribution consistent with that of radio-loud quasars. With an i-band magnitude of 18.47 and radio flux density of 33 mJy, its radio-loudness parameter is ~100. The optical and near-infrared spectra taken by Magellan enable us to estimate its bolometric luminosity to be L_bol ~ 1.1E48 erg/s, approximately 4.5 times greater than that of the most distant quasar known to date. The black hole mass of J0131-0321 is estimated to be 2.7E9 solar masses, with an uncertainty up to 0.4 dex. Detailed physical properties of this high-redshift, radio-loud, potentially super-Eddington quasar can be probed in the future with more dedicated and intensive follow-up observations using multi-wavelength facilities.

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

The dynamical and radiative properties of the quiescent state (X-ray luminosity $<10^{34}\ {\rm erg s^{-1} cm^{-2}}$) of black hole X-ray transients (BHXTs) remains unclear, mainly because of low-luminosity and poor data quantity. We demonstrate that, the simultaneous multi-wavelength (including radio, optical, ultraviolet and X-ray bands) spectrum of V404 Cyg in its bright quiescent state can be well described by the radiation from the companion star and more importantly, the compact jet. Neither the outer thin disc nor the inner hot accretion flow is important in the total spectrum. Together with several additional recent findings, i.e. the power-law X-ray spectrum and the constant X-ray spectral shape (or constant photon index) in contrast to the dramatic change in the X-ray luminosity, we argue the quiescent state spectrum of BHXTs is actually jet-dominated. Observational features consistent with this jet model are also discussed as supporting evidences.

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.

Identifications of quasars at intermediate redshifts (2.2<z<3.5) are inefficient in most previous quasar surveys as their optical colors are similar to those of stars. The near-IR K-band excess technique has been suggested to overcome this difficulty. Our study also proposed to use optical/near-IR colors for selecting z<4 quasars. To this method, we selected 105 unidentified bright targets with i<18.5 from the quasar candidates of SDSS DR6 with both SDSS ugriz optical and UKIDSS YJHK near-IR photometric data, which satisfy our proposed Y-K/g-z criterion and have photometric redshifts between 2.2 and 3.5 estimated from the 9-band SDSS-UKIDSS data. 43 of them were observed with the 2.16m telescope of NAOC in 2012. 36 of them were identified as quasars at 2.1<z<3.4. High success rate of discovering these quasars in the SDSS spectroscopic surveyed area demonstrates the robustness of both the Y-K/g-z selection criterion and the photometric redshift estimation technique. We also used the above criterion to investigate the possible star contamination rate to the quasar candidates of SDSS DR6, and found that it is much higher in selecting 3<z<3.5 quasar candidates than selecting lower redshift ones (z<2.2). The significant improvement in the photometric redshift estimation by using the 9-band SDSS-UKIDSS data than using the 5-band SDSS data is demonstrated and a catalog of 7,727 unidentified quasar candidates with photometric redshifts between 2.2 and 3.5 is provided. We also tested the Y-K/g-z selection criterion with the SDSS-III/DR9 quasar catalog, and found 96.2% of 17,999 DR9 quasars with UKIDSS Y and K-band data satisfy our criterion. With some samples of red and type II quasars, we found that 88% and 96.5% of them can be selected by the Y-K/g-z criterion respectively, which supports that using the Y-K/g-z criterion we can efficiently select both unobscured and obscured quasars. (abridged)

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

We propose an interferometry experiment for the detection of axion-like particles (ALPs). As in ordinary photon-regeneration (light shining through a wall) experiments, a laser beam traverses a region permeated by a magnetic field, where photons are converted to ALPs via the Primakoff process, resulting in a slight power loss and phase shift. The beam is then combined with a reference beam that originates from the same source. The detection of a change in the output intensity would signal the presence of ALPs (or possibly other particles that couple to the photon in a similar way). Because only one stage of conversion is needed, the signal is of ${\cal O}(g^2_{a\gamma \gamma})$, as opposed to ${\cal O}(g^4_{a\gamma \gamma})$ for photon-regeneration experiments, where $g_{a\gamma \gamma}$ is the coupling between ALPs and photons. This improvement over photon-regeneration is nullified by the presence of shot noise, which however can be reduced by the use of squeezed light, resulting in an improvement in the sensitivity to $g_{a\gamma\gamma}$ over ordinary photon-regeneration experiments by an order of $10^{1/2}$ assuming $10{\rm dB}$ noise suppression. Additionally, our setup can incorporate straightforwardly optical delay lines or Fabry-Perot cavities, boosting the signal by a factor of $n\sim10^3$, where $n$ is the number of times the laser beam is folded. This way, we can constrain $g_{a\gamma\gamma}$ better by yet another factor of $n^{1/2}\sim 10^{1.5}$, as compared to the $n^{1/4}$ boost that would be achieved in photon-regeneration experiments.

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.

Quasars with redshifts greater than 4 are rare, and can be used to probe the structure and evolution of the early universe. Here we report the discovery of six new quasars with $i$-band magnitudes brighter than 19.5 and redshifts between 2.4 and 4.6 from the YFOSC spectroscopy of the Lijiang 2.4m telescope in February, 2012. These quasars are in the list of $z>3.6$ quasar candidates selected by using our proposed $J-K/i-Y$ criterion and the photometric redshift estimations from the SDSS optical and UKIDSS near-IR photometric data. Nine candidates were observed by YFOSC, and five among six new quasars were identified as $z>3.6$ quasars. One of the other three objects was identified as a star and the other two were unidentified due to the lower signal-to-noise ratio of their spectra. This is the first time that $z>4$ quasars have been discovered using a telescope in China. Thanks to the Chinese Telescope Access Program (TAP), the redshift of 4.6 for one of these quasars was confirmed by the Multiple Mirror Telescope (MMT) Red Channel spectroscopy. The continuum and emission line properties of these six quasars, as well as their central black hole masses and Eddington ratios, were obtained.

We observe that photon cooling after big bang nucleosynthesis (BBN) but before recombination can remove the conflict between the observed and theoretically predicted value of the primordial abundance of $^7$Li. Such cooling is ordinarily difficult to achieve. However, the recent realization that dark matter axions form a Bose-Einstein condensate (BEC) provides a possible mechanism, because the much colder axions may reach thermal contact with the photons. This proposal predicts a high effective number of neutrinos as measured by the cosmic microwave anisotropy spectrum.

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.

Cold dark matter axions thermalize through gravitational self-interactions and form a Bose-Einstein condensate when the photon temperature reaches approximately 500 eV. Axion Bose-Einstein condensation provides an opportunity to distinguish axions from the other dark matter candidates on the basis of observation. The rethermalization of axions that are about to fall in a galactic potential well causes them to acquire net overall rotation, whereas ordinary cold dark matter falls in with an irrotational velocity field. The inner caustics of galactic halos are different in the two cases.

Axions differ from the other cold dark matter candidates in that they form a degenerate Bose gas. It is shown that their huge quantum degeneracy and large correlation length cause cold dark matter axions to thermalize through gravitational self-interactions when the photon temperature reaches approximately 500 eV. When they thermalize, the axions form a Bose-Einstein condensate. Their thermalization occurs in a regime, herein called the `condensed regime', where the Boltzmann equation is not valid because the energy dispersion of the particles is smaller than their interaction rate. We derive analytical expressions for the thermalization rate of particles in the condensed regime, and check the validity of these expressions by numerical simulation of a toy model. We revisit axion cosmology in light of axion Bose-Einstein condensation. It is shown that axions are indistinguishable from ordinary cold dark matter on all scales of observational interest, except when they thermalize or rethermalize. The rethermalization of axions that are about to fall in a galactic potential well causes them to acquire net overall rotation as they go to the lowest energy state consistent with the total angular momentum they acquired by tidal torquing. This phenomenon explains the occurrence of caustic rings of dark matter in galactic halos. We find that photons may reach thermal contact with axions and investigate the implications of this possibility for the measurements of cosmological parameters.

Cosmic ray protons propagating in a spatially-homogeneous but time-dependent field of axions or axion-like particles (ALPs) emit photons in a way that is reminiscent of Cherenkov radiation by charged particles in a preferred background. We compute the emission rate and energy spectrum of the photons, and discuss the possibility of their detection using the Square Kilometre Array which is currently under construction. In the case of a non-detection, constraints can be placed on the parameter space of ALPs whose mass lie between $10^{-7}$eV and $10^{-5}$ eV under the assumption that they are the primary constituent of dark matter.

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.

An optical photometric observation in 15 bands was carried out for nearby galaxy cluster Abell 1775 by the Beijing-Arizona-Taiwan-Connecticut (BATC) multi-color system. Over 5000 sources' spectral energy distributions (SEDs) were obtained. Since this cluster has also been observed by the Sloan Digital Sky Survey (SDSS), the BATC SEDs were combined with the SDSS five-band photometric data. Using the combined SEDs, 146 faint galaxies were selected as new member galaxies by the photometric redshift technique. Based on the positions, redshifts and 20-band SEDs of member galaxies, dynamical substructures and luminosity functions (LFs) of A1775 were investigated. The previous reported bimodal structure of A1775 has been confirmed: a poor subcluster with lower redshift, A1775A, is located ~14' southeast to the main concentration A1775B. After taking into account the new supplemented member galaxies, a new subcluster A1775C was found along the aligned direction of A1775A and A1775B. The different LF faint ends of the two subclusters indicate that A1775B is a more dynamically evolved system, while A1775A is still dynamically young. By the STARLIGHT spectral synthesis code, the star- formation histories of the member galaxies were studied. The dependence of the mean stellar ages upon the Hubble type was confirmed, and the environmental effect on star-formation activities for galaxies in A1775B has been explored.

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.

Cold dark matter axions form a Bose-Einstein condensate if the axions thermalize. Recently, it was found that they do thermalize when the photon temperature reaches T ~ 100 eV(f/10^12GeV)^1/2 and that they continue to do so thereafter. We discuss the differences between axion BEC and CDM in the linear regime and the non-linear regime of evolution of density perturbations. We find that axion BEC provides a mechanism for the production of net overall rotation in dark matter halos, and for the alignment of cosmic microwave anisotropy multi-poles.

We show that cold dark matter axions thermalize and form a Bose-Einstein condensate. We obtain the axion state in a homogeneous and isotropic universe, and derive the equations governing small axion perturbations. Because they form a BEC, axions differ from ordinary cold dark matter in the non-linear regime of structure formation and upon entering the horizon. Axion BEC provides a mechanism for the production of net overall rotation in dark matter halos, and for the alignment of cosmic microwave anisotropy multipoles.