The phantom dark energy remarkably boosts our prehension of the accelerating Universe. Various models are widely discussed in the phantom Universe without bulk viscosity. From the hydrodynamics' point of view, it is natural to introduce the nonperfect fluid in the study of the Universe, as an ideal fluid just an approximation to the real world after all and using the generalized equation of state (EoS) with bulk viscosity, the early inflationary universe and the accelerated expansion of the late-time universe are described by many authors. In this paper, in order to investigate how the viscosity will influence the evolution of the Universe, we study a class of phantom dark energy models with bulk viscosity by the method of dynamical analysis technique. We show that there are different cosmic late-time behaviors and the stability also brings some constraints on the models. We also plot the evolutionary trajectories of this model in the statefinder parameter-planes to see the different behaviors of the models from the statefinder viewpoint.

We report a measurement of the E-mode polarization power spectrum of the cosmic microwave background (CMB) using 150 GHz data taken from July 2014 to December 2016 with the POLARBEAR experiment. We reach an effective polarization map noise level of $32\,\mu\mathrm{K}$-$\mathrm{arcmin}$ across an observation area of 670 square degrees. We measure the EE power spectrum over the angular multipole range $500 \leq \ell <3000$, tracing the third to seventh acoustic peaks with high sensitivity. The statistical uncertainty on E-mode bandpowers is $\sim 2.3 \mu {\rm K}^2$ at $\ell \sim 1000$ with a systematic uncertainty of 0.5$\mu {\rm K}^2$. The data are consistent with the standard $\Lambda$CDM cosmological model with a probability-to-exceed of 0.38. We combine recent CMB E-mode measurements and make inferences about cosmological parameters in $\Lambda$CDM as well as in extensions to $\Lambda$CDM. Adding the ground-based CMB polarization measurements to the Planck dataset reduces the uncertainty on the Hubble constant by a factor of 1.2 to $H_0 = 67.20 \pm 0.57 {\rm km\,s^{-1} \,Mpc^{-1}}$. When allowing the number of relativistic species ($N_{eff}$) to vary, we find $N_{eff} = 2.94 \pm 0.16$, which is in good agreement with the standard value of 3.046. Instead allowing the primordial helium abundance ($Y_{He}$) to vary, the data favor $Y_{He} = 0.248 \pm 0.012$. This is very close to the expectation of 0.2467 from Big Bang Nucleosynthesis. When varying both $Y_{He}$ and $N_{eff}$, we find $N_{eff} = 2.70 \pm 0.26$ and $Y_{He} = 0.262 \pm 0.015$.

Recently, a new kind of $f(z)$ theory is proposed to provide a different perspective for the development of reliable alternative models of gravity in which the $f(R)$ Lagrangian terms are reformulated as polynomial parameterizations $f(z)$. In the previous study, the parameters in the $f(z)$ models have been constrained by using cosmological data. In this paper, these models will be tested by the observations in the solar system. After solving the Ricci scalar as a function of the redshift, one could obtain $f(R)$ that could be used to calculate the standard Parameterized-Post-Newtonian (PPN) parameters. First, we fit the parametric models with the latest cosmological observational data. Then the tests are performed by solar system observations. And last we combine the constraints of solar system and cosmology together and reconstruct the $f(R)$ actions of the $f(z)$ parametric models.

The spectral energy distribution (SED) sequence for type Ia supernovae (SN Ia) is modeled by an artificial neural network. The SN Ia luminosity is characterized as a function of phase, wavelength, a color parameter and a decline rate parameter. After training and testing the neural network, the SED sequence could give both the spectrum with wavelength range from 3000\AA~to 8000\AA~ and the light curve with phase from 20 days before to 50 days after the maximum luminosity for the supernovae with different colors and decline rates. Therefore, we call this the Artificial Neural Network Spectral Light Curve Template (ANNSLCT) model. We retrain the Joint Light-curve Analysis (JLA) supernova sample by using the ANNSLCT model and obtain the parameters for each supernova to make a constraint on the cosmological $\Lambda$CDM model. We find that the best fitting values of these parameters are almost the same as those from the JLA sample trained with the Spectral Adaptive Lightcurve Template 2 (SALT2) model. So we believe that the ANNSLCT model could be used to analyze a large number of SN Ia multi-color light curves measured in the current and future observational projects.

Using only cosmic microwave background polarization data from the POLARBEAR experiment, we measure $B$-mode polarization delensing on subdegree scales at more than $5\sigma$ significance. We achieve a 14% $B$-mode power variance reduction, the highest to date for internal delensing, and improve this result to 2% by applying for the first time an iterative maximum a posteriori delensing method. Our analysis demonstrates the capability of internal delensing as a means of improving constraints on inflationary models, paving the way for the optimal analysis of next-generation primordial $B$-mode experiments.

We present a measurement of the gravitational lensing deflection power spectrum reconstructed with two seasons cosmic microwave background polarization data from the POLARBEAR experiment. Observations were taken at 150 GHz from 2012 to 2014 which survey three patches of sky totaling 30 square degrees. We test the consistency of the lensing spectrum with a Cold Dark Matter (CDM) cosmology and reject the no-lensing hypothesis at a confidence of 10.9 sigma including statistical and systematic uncertainties. We observe a value of A_L = 1.33 +/- 0.32 (statistical) +/- 0.02 (systematic) +/- 0.07 (foreground) using all polarization lensing estimators, which corresponds to a 24% accurate measurement of the lensing amplitude. Compared to the analysis of the first year data, we have improved the breadth of both the suite of null tests and the error terms included in the estimation of systematic contamination.

We study the polarisation properties of extragalactic sources at 95 and 150 GHz in the SPTpol 500 deg$^2$ survey. We estimate the polarised power by stacking maps at known source positions, and correct for noise bias by subtracting the mean polarised power at random positions in the maps. We show that the method is unbiased using a set of simulated maps with similar noise properties to the real SPTpol maps. We find a flux-weighted mean-squared polarisation fraction $\langle p^2 \rangle= [8.9\pm1.1] \times 10^{-4}$ at 95 GHz and $[6.9\pm1.1] \times 10^{-4}$ at 150~GHz for the full sample. This is consistent with the values obtained for a sub-sample of active galactic nuclei. For dusty sources, we find 95 per cent upper limits of $\langle p^2 \rangle_{\rm 95}<16.9 \times 10^{-3}$ and $\langle p^2 \rangle_{\rm 150}<2.6 \times 10^{-3}$. We find no evidence that the polarisation fraction depends on the source flux or observing frequency. The 1-$\sigma$ upper limit on measured mean squared polarisation fraction at 150 GHz implies that extragalactic foregrounds will be subdominant to the CMB E and B mode polarisation power spectra out to at least $\ell\lesssim5700$ ($\ell\lesssim4700$) and $\ell\lesssim5300$ ($\ell\lesssim3600$), respectively at 95 (150) GHz.

The cosmic infrared background (CIB) anisotropies and cosmic microwave background (CMB) lensing are powerful measurements for exploring cosmological and astrophysical problems. In this work, we measure the auto-correlation power spectrum of the CIB anisotropies in the Herschel-SPIRE HerMES Large Mode Survey (HeLMS) field, which covers more than 280 square degrees of the sky, and the cross power spectrum between HeLMS CIB and the CMB lensing from Planck. We use the Herschel Level 1 timeline data to merge the CIB maps at 250, 350 and 500 um bands, and mask the areas where the flux is greater than 3-sigma (~50 mJy/beam) or no measured data. We obtain the final CIB power spectra at 100<ell<20,000 by considering several effects, such as beam function, mode coupling, transfer function, and so on. We also calculate the theoretical CIB auto- and cross-power spectra of CIB and CMB lensing by assuming that the CIB emissivity follows Gaussian distribution in redshift. The Markov chain Monte Carlo (MCMC) method is used to perform a joint fit of the CIB auto-power spectra and cross-power spectra of CIB and CMB lensing data. We find that our model can fit the power spectrum data very well, and obtain basically consistent results with higher accuracy comparing to previous studies.

If a fraction $f_{\rm dcdm}$ of the Dark Matter decays into invisible and massless particles (so-called "dark radiation") with the decay rate (or inverse lifetime) $\Gamma_{\rm dcdm}$, such decay will leave distinctive imprints on cosmological observables. With a full consideration of the Boltzmann hierarchy, we calculate the decay-induced impacts not only on the CMB but also on the redshift distortion and the kinetic Sunyaev-Zel'dovich effect, while providing detailed physical interpretations based on evaluating the evolution of gravitational potential. By using the current cosmological data with a combination of Planck 2015, Baryon Acoustic Oscillation and redshift distortion measurements which can improve the constraints, we update the $1\sigma$ bound on the fraction of decaying DM from $f_{\rm dcdm}\lesssim5.26\%$ to $f_{\rm dcdm}\lesssim2.73\%$ for the short-lived DM (assuming $\Gamma_{\rm dcdm}/H_0\gtrsim10^4$). However, no constraints are improved from RSD data ($f_{\rm dcdm}\lesssim0.94\%$) for the long-lived DM (i.e., $\Gamma_{\rm dcdm}/H_0\lesssim10^4$). We also find the fractional DM decay can only slightly reduce the $H_0$ and $\sigma_8$ tensions, which is consistent with other previous works. Furthermore, our calculations show that the kSZ effect in future would provide a further constraining power on the decaying DM.

We report a 4.8$\sigma$ measurement of the cross-correlation signal between the cosmic microwave background (CMB) lensing convergence reconstructed from measurements of the CMB polarization made by the POLARBEAR experiment and the infrared-selected galaxies of the Herschel-ATLAS survey. This is the first measurement of its kind. We infer a best-fit galaxy bias of $b = 5.76 \pm 1.25$, corresponding to a host halo mass of $\log_{10}(M_h/M_\odot) =13.5^{+0.2}_{-0.3}$ at an effective redshift of $z \sim 2$ from the cross-correlation power spectrum. Residual uncertainties in the redshift distribution of the sub-mm galaxies are subdominant with respect to the statistical precision. We perform a suite of systematic tests, finding that instrumental and astrophysical contaminations are small compared to the statistical error. This cross-correlation measurement only relies on CMB polarization information that, differently from CMB temperature maps, is less contaminated by galactic and extra-galactic foregrounds, providing a clearer view of the projected matter distribution. This result demonstrates the feasibility and robustness of this approach for future high-sensitivity CMB polarization experiments.

We present the first measurement of cross-correlation between the lensing potential, reconstructed from cosmic microwave background (CMB) {\it polarization} data, and the cosmic shear field from galaxy shapes. This measurement is made using data from the POLARBEAR CMB experiment and the Subaru Hyper Suprime-Cam (HSC) survey. By analyzing an 11~deg$^2$ overlapping region, we reject the null hypothesis at 3.5$\sigma$\ and constrain the amplitude of the {\bf cross-spectrum} to $\widehat{A}_{\rm lens}=1.70\pm 0.48$, where $\widehat{A}_{\rm lens}$ is the amplitude normalized with respect to the Planck~2018{} prediction, based on the flat $\Lambda$ cold dark matter cosmology. The first measurement of this {\bf cross-spectrum} without relying on CMB temperature measurements is possible due to the deep POLARBEAR map with a noise level of ${\sim}$6\,$\mu$K-arcmin, as well as the deep HSC data with a high galaxy number density of $n_g=23\,{\rm arcmin^{-2}}$. We present a detailed study of the systematics budget to show that residual systematics in our results are negligibly small, which demonstrates the future potential of this cross-correlation technique.

We propose in this White Paper a concept for a space experiment using cold atoms to search for ultra-light dark matter, and to detect gravitational waves in the frequency range between the most sensitive ranges of LISA and the terrestrial LIGO/Virgo/KAGRA/INDIGO experiments. This interdisciplinary experiment, called Atomic Experiment for Dark Matter and Gravity Exploration (AEDGE), will also complement other planned searches for dark matter, and exploit synergies with other gravitational wave detectors. We give examples of the extended range of sensitivity to ultra-light dark matter offered by AEDGE, and how its gravitational-wave measurements could explore the assembly of super-massive black holes, first-order phase transitions in the early universe and cosmic strings. AEDGE will be based upon technologies now being developed for terrestrial experiments using cold atoms, and will benefit from the space experience obtained with, e.g., LISA and cold atom experiments in microgravity. This paper is based on a submission (v1) in response to the Call for White Papers for the Voyage 2050 long-term plan in the ESA Science Programme. ESA limited the number of White Paper authors to 30. However, in this version (v2) we have welcomed as supporting authors participants in the Workshop on Atomic Experiments for Dark Matter and Gravity Exploration held at CERN: ({\tt https://indico.cern.ch/event/830432/}), as well as other interested scientists, and have incorporated additional material.

We present a measurement of the $B$-mode polarization power spectrum of the cosmic microwave background (CMB) using taken from July 2014 to December 2016 with the POLARBEAR experiment. The CMB power spectra are measured using observations at 150 GHz with an instantaneous array sensitivity of $\mathrm{NET}_\mathrm{array}=23\, \mu \mathrm{K} \sqrt{\mathrm{s}}$ on a 670 square degree patch of sky centered at (RA, Dec)=($+0^\mathrm{h}12^\mathrm{m}0^\mathrm{s},-59^\circ18^\prime$). A continuously rotating half-wave plate is used to modulate polarization and to suppress low-frequency noise. We achieve $32\,\mu\mathrm{K}$-$\mathrm{arcmin}$ effective polarization map noise with a knee in sensitivity of $\ell = 90$, where the inflationary gravitational wave signal is expected to peak. The measured $B$-mode power spectrum is consistent with a $\Lambda$CDM lensing and single dust component foreground model over a range of multipoles $50 \leq \ell \leq 600$. The data disfavor zero $C_\ell^{BB}$ at $2.2\sigma$ using this $\ell$ range of POLARBEAR data alone. We cross-correlate our data with Planck high frequency maps and find the low-$\ell$ $B$-mode power in the combined dataset to be consistent with thermal dust emission. We place an upper limit on the tensor-to-scalar ratio $r < 0.90$ at 95% confidence level after marginalizing over foregrounds.

The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

We report the first detection of gravitational lensing due to galaxy clusters using only the polarization of the cosmic microwave background (CMB). The lensing signal is obtained using a new estimator that extracts the lensing dipole signature from stacked images formed by rotating the cluster-centered Stokes $Q/U$ map cutouts along the direction of the locally measured background CMB polarization gradient. Using data from the SPTpol 500 deg$^{2}$ survey at the locations of roughly 18,000 clusters with richness $\lambda \ge 10$ from the Dark Energy Survey (DES) Year-3 full galaxy cluster catalog, we detect lensing at $4.8\sigma$. The mean stacked mass of the selected sample is found to be $(1.43 \pm 0.4)\ \times 10^{14}\ {\rm M_{\odot}}$ which is in good agreement with optical weak lensing based estimates using DES data and CMB-lensing based estimates using SPTpol temperature data. This measurement is a key first step for cluster cosmology with future low-noise CMB surveys, like CMB-S4, for which CMB polarization will be the primary channel for cluster lensing measurements.

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

The cosmic infrared background (CIB) is slightly polarized. Polarization directions of individual galaxies could be aligned with tidal fields around galaxies, resulting in nonzero CIB polarization. We use a linear intrinsic alignment model to theoretically predict angular correlations of the CIB polarization fluctuations and find that electriclike and curl-like ($B$-mode) polarization modes are equally generated with power four orders of magnitude less than its intensity. The CIB $B$-mode signal is negligible and not a concerning foreground for the inflationary $B$-mode searches at nominal frequencies for cosmic microwave background measurements, but could be detected at submillimetre wavelengths by future space missions.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

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

The cosmic infrared background (CIB) is gravitationally lensed. A quadratic-estimator technique that is inherited from lensing analyses of the cosmic microwave background (CMB) can be applied to detect the CIB lensing effects. However, the CIB fluctuations are intrinsically strongly non-Gaussian, making CIB lensing reconstruction highly biased. We perform numerical simulations to estimate the intrinsic non-Gaussianity and establish a cross-correlation approach to precisely extract the CIB lensing signal from raw data. We apply this technique to CIB data from the Planck satellite and cross-correlate the resulting lensing estimate with the CIB data, galaxy number counts and the CMB lensing potential. We detect an excess that is consistent with a lensing contribution at $>4\sigma$.

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

The DArk Matter Particle Explorer (DAMPE) developed in China was designed to search for evidence of dark matter particles by observing primary cosmic rays and gamma rays in the energy range from 5 GeV to 10 TeV. Since its launch in December 2015, a large quantity of data has been recorded. With the data set acquired during more than a year of operation in space, a precise time-dependent calibration for the energy measured by the BGO ECAL has been developed. In this report, the instrumentation and development of the BGO Electromagnetic Calorimeter (BGO ECAL) are briefly described. The calibration on orbit, including that of the pedestal, attenuation length, minimum ionizing particle peak, and dynode ratio, is discussed, and additional details about the calibration methods and performance in space are presented.

The polarization of the atmosphere has been a long-standing concern for ground-based experiments targeting cosmic microwave background (CMB) polarization. Ice crystals in upper tropospheric clouds scatter thermal radiation from the ground and produce a horizontally-polarized signal. We report the detailed analysis of the cloud signal using a ground-based CMB experiment, POLARBEAR, located at the Atacama desert in Chile and observing at 150 GHz. We observe horizontally-polarized temporal increases of low-frequency fluctuations ("polarized bursts," hereafter) of $\lesssim$0.1 K when clouds appear in a webcam monitoring the telescope and the sky. The hypothesis of no correlation between polarized bursts and clouds is rejected with $>$24$\sigma$ statistical significance using three years of data. We consider many other possibilities including instrumental and environmental effects, and find no other reasons other than clouds that can explain the data better. We also discuss the impact of the cloud polarization on future ground-based CMB polarization experiments.

The near-infrared background between 0.5 $\mu$m to 2 $\mu$m contains a wealth of information related to radiative processes in the universe. Infrared background anisotropies encode the redshift-weighted total emission over cosmic history, including any spatially diffuse and extended contributions. The anisotropy power spectrum is dominated by undetected galaxies at small angular scales and diffuse background of Galactic emission at large angular scales. In addition to these known sources, the infrared background also arises from intra-halo light (IHL) at $z < 3$ associated with tidally-stripped stars during galaxy mergers. Moreover, it contains information on the very first galaxies from the epoch of reionization (EoR). The EoR signal has a spectral energy distribution (SED) that goes to zero near optical wavelengths due to Lyman absorption, while other signals have spectra that vary smoothly with frequency. Due to differences in SEDs and spatial clustering, these components may be separated in a multi-wavelength-fluctuation experiment. To study the extent to which EoR fluctuations can be separated in the presence of IHL, extra-galactic and Galactic foregrounds, we develop a maximum likelihood technique that incorporates a full covariance matrix among all the frequencies at different angular scales. We apply this technique to simulated deep imaging data over a 2$\times$100 deg$^2$ sky area from 0.75 $\mu$m to 5 $\mu$m in 9 bands and find that such a "frequency tomography" can successfully reconstruct both the amplitude and spectral shape for representative EoR, IHL and the foreground signals.

We propose a hybrid quadratic estimator to measure cross correlations between gravitational lensing of the cosmic microwave background (CMB) and differential screening effects arising from fluctuations in the electron column density, such as could arise from patchy reionization. The hybrid quadratic estimators are validated by simulated data sets with both Planck and CMB-Stage 4 (CMB-S4) instrumental properties and found to be able to recover the cross-power spectra with almost no biases. We apply this technique to Planck 2015 temperature data and obtain cross-power spectra between gravitational lensing and differential screening effects. Planck data alone cannot detect the patchy-reionization-induced cross-power spectrum but future experiment like CMB-S4 will be able to robustly measure the expected signal and deliver new insights on reionization.

We report a measurement of the B-mode polarization power spectrum in the cosmic microwave background (CMB) using the POLARBEAR experiment in Chile. The faint B-mode polarization signature carries information about the Universe's entire history of gravitational structure formation, and the cosmic inflation that may have occurred in the very early Universe. Our measurement covers the angular multipole range 500 < l < 2100 and is based on observations of an effective sky area of 25 square degrees with 3.5 arcmin resolution at 150 GHz. On these angular scales, gravitational lensing of the CMB by intervening structure in the Universe is expected to be the dominant source of B-mode polarization. Including both systematic and statistical uncertainties, the hypothesis of no B-mode polarization power from gravitational lensing is rejected at 97.1% confidence. The band powers are consistent with the standard cosmological model. Fitting a single lensing amplitude parameter A_BB to the measured band powers, A_BB = 1.12 +/- 0.61 (stat) +0.04/-0.12 (sys) +/- 0.07 (multi), where A_BB = 1 is the fiducial WMAP-9 LCDM value. In this expression, "stat" refers to the statistical uncertainty, "sys" to the systematic uncertainty associated with possible biases from the instrument and astrophysical foregrounds, and "multi" to the calibration uncertainties that have a multiplicative effect on the measured amplitude A_BB.

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

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 revisit the global 21cm signal calculation incorporating a possible radio background at early times, and find that the global 21cm signal shows a much stronger absorption feature, which could enhance detection prospects for future 21 cm experiments. In light of recent reports of a possible low-frequency excess radio background, we propose that detailed 21 cm calculations should include a possible early radio background.

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 constrain interacting dark matter and dark energy (IDMDE) models using a 450-degree-square cosmic shear data from the Kilo Degree Survey (KiDS) and the angular power spectra from Planck's latest cosmic microwave background measurements. We revisit the discordance problem in the standard Lambda cold dark matter ($\Lambda$CDM) model between weak lensing and Planck datasets and extend the discussion by introducing interacting dark sectors. The IDMDE models are found to be able to alleviate the discordance between KiDS and Planck as previously inferred from the $\Lambda$CDM model, and moderately favored by a combination of the two datasets.

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

Recent progress in high sensitivity Cosmic Microwave Background (CMB) polarization experiments opens up a window on large scale structure (LSS), as CMB polarization fluctuations on small angular scales can arise from a combination of LSS and ionization fluctuations in the late universe. Gravitational lensing effects can be extracted from CMB datasets with quadratic estimators but reconstructions of electron density fluctuations (EDFs) with quadratic estimators are found to be significantly biased by the much larger lensing effects in the secondary CMB fluctuations. In this paper we establish a bispectrum formalism using tracers of LSS to extract the subdominant EDFs from CMB polarization data. We find that this bispectrum can effectively reconstruct angular band-powers of cross correlation between EDFs and LSS tracers. Next generation CMB polarization experiments in conjunction with galaxy surveys and cosmic infrared background experiments can detect signatures of EDFs with high significance.

We construct and train an artificial neural network called the back-propagation neural network to describe the evolution of the type Ia supernova spectrum by using the data from the CfA Supernova Program. This network method has many attractive features, and one of them is that the constructed model is differentiable. Benefitting from this, we calculate the absorption velocity and its variation. The model we constructed can well describe not only the spectrum of SNe Ia with wavelength range from $3500\AA$ to $8000\AA$, but also the light-curve evolution with phase time from $-15$ to $50$ with different colors. Moreover, the number of parameters needed during the training process is much less than the usual methods.

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

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

We report an improved measurement of the cosmic microwave background (CMB) $B$-mode polarization power spectrum with the POLARBEAR experiment at 150 GHz. By adding new data collected during the second season of observations (2013-2014) to re-analyzed data from the first season (2012-2013), we have reduced twofold the band-power uncertainties. The band powers are reported over angular multipoles $500 \leq \ell \leq 2100$, where the dominant $B$-mode signal is expected to be due to the gravitational lensing of $E$-modes. We reject the null hypothesis of no $B$-mode polarization at a confidence of 3.1$\sigma$ including both statistical and systematic uncertainties. We test the consistency of the measured $B$-modes with the $\Lambda$ Cold Dark Matter ($\Lambda$CDM) framework by fitting for a single lensing amplitude parameter $A_L$ relative to the Planck best-fit model prediction. We obtain $A_L = 0.60 ^{+0.26} _{-0.24} ({\rm stat}) ^{+0.00} _{-0.04}({\rm inst}) \pm 0.14 ({\rm foreground}) \pm 0.04 ({\rm multi})$, where $A_{L}=1$ is the fiducial $\Lambda$CDM value, and the details of the reported uncertainties are explained later in the manuscript.

We study 21cm and Ly$\mathrm{\alpha}$ fluctuations, as well as H$\mathrm{\alpha}$, while distinguishing between Ly$\mathrm{\alpha}$ emission of galactic, diffuse, and scattered intergalactic medium (IGM) origin. Cross-correlation information about the state of the IGM is obtained, testing neutral versus ionized medium cases with different tracers in a seminumerical simulation setup. In order to pave the way toward constraints on reionization history and modeling beyond power spectrum information, we explore parameter dependencies of the cross-power signal between 21$\,$cm and Ly$\mathrm{\alpha}$, which displays a characteristic morphology and a turnover from negative to positive correlation at scales of a couple Mpc$^{-1}$. In a proof of concept for the extraction of further information on the state of the IGM using different tracers, we demonstrate the use of the 21$\,$cm and H$\mathrm{\alpha}$ cross-correlation signal to determine the relative strength of galactic and IGM emission in Ly$\mathrm{\alpha}$. We conclude by showing the detectability of the 21$\,$cm and Ly$\mathrm{\alpha}$ cross-correlation signal over more than one decade in scale at high signal-to-noise ratio for upcoming probes like SKA and the proposed all-sky intensity mapping satellites SPHEREx and CDIM, while also including the Ly$\mathrm{\alpha}$ damping tail and 21cm foreground avoidance in the modeling.

We estimate uncertainties of cosmological parameters for phenomenological interacting dark energy models using weak lensing convergence power spectrum and bispectrum. We focus on the bispectrum tomography and examine how well the weak lensing bispectrum with tomography can constrain the interactions between dark sectors, as well as other cosmological parameters. Employing the Fisher matrix analysis, we forecast parameter uncertainties derived from weak lensing bispectra with a two-bin tomography and place upper bounds on strength of the interactions between the dark sectors. The cosmic shear will be measured from upcoming weak lensing surveys with high sensitivity, thus it enables us to use the higher order correlation functions of weak lensing to constrain the interaction between dark sectors and will potentially provide more stringent results with other observations combined.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

We present a halo-model-based approach to calculate the cross-correlation between $\rm{21\,cm}$ HI intensity fluctuations and $\rm{Ly}\alpha$ emitters (LAE) during the epoch of reionization (EoR). Ionizing radiation around dark matter halos are modeled as bubbles with the size and growth determined based on the reionization photon production, among other physical parameters. The cross-correlation shows a clear negative-to-positive transition, associated with transition from ionized to neutral hydrogen in the intergalactic medium during EoR. The cross-correlation is subject to several foreground contaminants, including foreground radio point sources important for $\rm{21\,cm}$ experiments and low-$z$ interloper emission lines, such as $\rm{H}\alpha$, OIII, and OII, for $\rm{Ly}\alpha$ experiments. Our calculations show that by masking out high fluxes in the $\rm{Ly}\alpha$ measurement, the correlated foreground contamination on the $\rm{21\,cm}$-$\rm{Ly}\alpha$ cross-correlation can be dramatically reduced. We forecast the detectability of $\rm{21\,cm}$-$\rm{Ly}\alpha$ cross-correlation at different redshifts and adopt a Fisher matrix approach to estimate uncertainties on the key EoR parameters that have not been well constrained by other observations of reionization. This halo-model-based approach enables us to explore the EoR parameter space rapidly for different $\rm{21\,cm}$ and $\rm{Ly}\alpha$ experiments.

We present new limits on the amplitude of potential primordial magnetic fields (PMFs) using temperature and polarization measurements of the cosmic microwave background (CMB) from Planck, BICEP2/Keck Array, POLARBEAR, and SPTpol. We reduce twofold the 95% CL upper limit on the CMB anisotropy power due to a nearly-scale-invariant PMF, with an allowed B-mode power at $\ell=1500$ of $D_{\ell=1500}^{BB} < 0.071 \mu K^2$ for Planck versus $D_{\ell=1500}^{BB} < 0.034 \mu K^2$ for the combined dataset. We also forecast the expected limits from soon-to-deploy CMB experiments (like SPT-3G, Adv. ACTpol, or the Simons Array) and the proposed CMB-S4 experiment. Future CMB experiments should dramatically reduce the current uncertainties, by one order of magnitude for the near-term experiments and two orders of magnitude for the CMB-S4 experiment. The constraints from CMB-S4 have the potential to rule out much of the parameter space for PMFs.

We measure the Cosmic Microwave Background (CMB) skewness power spectrum in $\textit{Planck}$, using frequency maps of the HFI instrument and the Sunyaev-Zel'dovich (SZ) component map. The two-to-one skewness power spectrum measures the cross-correlation between CMB lensing and the thermal SZ effect. We also directly measure the same cross-correlation using $\textit{Planck}$ CMB lensing map and the SZ map and compare it to the cross-correlation derived from the skewness power spectrum. We model fit the SZ power spectrum and CMB lensing-SZ cross power spectrum via the skewness power spectrum to constrain the gas pressure profile of dark matter halos. The gas pressure profile is compared to existing measurements in the literature including a direct estimate based on the stacking of SZ clusters in $\textit{Planck}$.

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.

It is known that the Cardassian universe is successful in describing the accelerated expansion of the universe, but its dynamical equations are hard to get from the action principle. In this paper, we establish the connection between the Cardassian universe and $f(T, \mathcal{T})$ gravity, where $T$ is the torsion scalar and $\mathcal{T}$ is the trace of the matter energy-momentum tensor. For dust matter, we find that the modified Friedmann equations from $f(T, \mathcal{T})$ gravity can correspond to those of Cardassian models, and thus, a possible origin of Cardassian universe is given. We obtain the original Cardassian model, the modified polytropic Cardassian model, and the exponential Cardassian model from the Lagrangians of $f(T,\mathcal{T})$ theory. Furthermore, by adding an additional term to the corresponding Lagrangians, we give three generalized Cardassian models from $f(T,\mathcal{T})$ theory. Using the observation data of type Ia supernovae, cosmic microwave background radiation, and baryon acoustic oscillations, we get the fitting results of the cosmological parameters and give constraints of model parameters for all of these models.

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

The extragalactic $\gamma$-ray background and its spatial anisotropy could potentially contain a signature of dark matter (DM) annihilation or particle decay. Astrophysical foregrounds, such as blazars and star-forming galaxies (SFGs), however, dominate the $\gamma$-ray background, precluding an easy detection of the signal associated with the DM annihilation or decay in the background intensity spectrum. The DM imprint on the $\gamma$-ray background is expected to be correlated with large-scale structure tracers. In some cases, such a cross-correlation is even expected to have a higher signal-to-noise ratio than the auto-correlation. One reliable tracer of the DM distribution in the large-scale structure is lensing of the cosmic microwave background (CMB), and the cosmic infrared background (CIB) is a reliable tracer of SFGs. We analyze Fermi-LAT data taken over 92 months and study the cross-correlation with Planck CMB lensing, Planck CIB, and Fermi-$\gamma$ maps. We put upper limits on the DM annihilation cross-section from the cross-power spectra with the $\gamma$-ray background anisotropies. The unbiased power spectrum estimation is validated with simulations that include cross-correlated signals. We also provide a set of systematic tests and show that no significant contaminations are found for the measurements presented here. Using $\gamma$-ray background map from data gathered over 92 months, we find the best constraint on the DM annihilation with a $1\sigma$ confidence level upper limit of $10^{-25}$-$10^{-24}$ cm$^{3}$ s$^{-1}$, when the mass of DM particles is between 20 and 100 GeV.

Intensity mapping is now becoming a useful tool to study the large-scale structure of the universe through spatial variations in the integrated emission from galaxies and the intergalactic medium. We study intensity mapping of the H-alpha 6563AA, [OIII]5007AA, [OII]3727AA and H-beta 4861AA lines at 0.8<z<5.2. The mean intensities of these four emission lines are estimated using the observed luminosity functions (LFs), cosmological simulations, and the star formation rate density (SFRD) derived from observations at z<5. We calculate the intensity power spectra and consider the foreground contamination of other lines at lower redshifts. We use the proposed NASA small explorer SPHEREx (the Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer) as a case study for the detectability of the intensity power spectra of the four emission lines. We also investigate the cross correlation with the 21-cm line probed by CHIME (the Canadian Hydrogen Intensity Mapping Experiment), Tianlai experiment and SKA (the Square Kilometer Array) at 0.8<z<2.4. We find both the auto and cross power spectra can be well measured for the H-alpha, [OIII] and [OII] lines at z<3, while it is more challenging for the H-beta line. Finally, we estimate the constraint on the SFRD from intensity mapping, and find we can reach accuracy higher than 7% at z<4, which is better than usual measurements using the LFs of galaxies.

POLARBEAR-2 (PB-2) is a cosmic microwave background (CMB) polarization experiment that will be located in the Atacama highland in Chile at an altitude of 5200 m. Its science goals are to measure the CMB polarization signals originating from both primordial gravitational waves and weak lensing. PB-2 is designed to measure the tensor to scalar ratio, r, with precision {\sigma}(r) < 0.01, and the sum of neutrino masses, {\Sigma}m{\nu}, with {\sigma}({\Sigma}m{\nu}) < 90 meV. To achieve these goals, PB-2 will employ 7588 transition-edge sensor bolometers at 95 GHz and 150 GHz, which will be operated at the base temperature of 250 mK. Science observations will begin in 2017.

SPHEREx is a proposed SMEX mission selected for Phase A. SPHEREx will carry out the first all-sky spectral survey and provide for every 6.2" pixel a spectra between 0.75 and 4.18 $\mu$m [with R$\sim$41.4] and 4.18 and 5.00 $\mu$m [with R$\sim$135]. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. It is readily apparent, however, that many other questions in astrophysics and planetary sciences could be addressed with the SPHEREx data. The SPHEREx team convened a community workshop in February 2016, with the intent of enlisting the aid of a larger group of scientists in defining these questions. This paper summarizes the rich and varied menu of investigations that was laid out. It includes studies of the composition of main belt and Trojan/Greek asteroids; mapping the zodiacal light with unprecedented spatial and spectral resolution; identifying and studying very low-metallicity stars; improving stellar parameters in order to better characterize transiting exoplanets; studying aliphatic and aromatic carbon-bearing molecules in the interstellar medium; mapping star formation rates in nearby galaxies; determining the redshift of clusters of galaxies; identifying high redshift quasars over the full sky; and providing a NIR spectrum for most eROSITA X-ray sources. All of these investigations, and others not listed here, can be carried out with the nominal all-sky spectra to be produced by SPHEREx. In addition, the workshop defined enhanced data products and user tools which would facilitate some of these scientific studies. Finally, the workshop noted the high degrees of synergy between SPHEREx and a number of other current or forthcoming programs, including JWST, WFIRST, Euclid, GAIA, K2/Kepler, TESS, eROSITA and LSST.

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.

To probe the late evolution history of the Universe, we adopt two kinds of optimal basis systems. One of them is constructed by performing the principle component analysis (PCA) and the other is build by taking the multidimensional scaling (MDS) approach. Cosmological observables such as the luminosity distance can be decomposed into these basis systems. These basis are optimized for different kinds of cosmological models that based on different physical assumptions, even for a mixture model of them. Therefore, the so-called feature space that projected from the basis systems is cosmological model independent, and it provide a parameterization for studying and reconstructing the Hubble expansion rate from the supernova luminosity distance and even gamma-ray bursts (GRBs) data with self-calibration. The circular problem when using GRBs as cosmological candles is naturally eliminated in this procedure. By using the Levenberg-Marquardt (LM) technique and the Markov Chain Monte Carlo (MCMC) method, we perform an observational constraint on this kind of parameterization. The data we used include the "joint light-curve analysis" (JLA) data set that consists of $740$ Type Ia supernovae (SNIa) as well as $109$ long gamma-ray bursts with the well-known Amati relation.

We constrain anisotropic cosmic birefringence using four-point correlations of even-parity $E$-mode and odd-parity $B$-mode polarization in the cosmic microwave background measurements made by the POLARization of the Background Radiation (POLARBEAR) experiment in its first season of observations. We find that the anisotropic cosmic birefringence signal from any parity-violating processes is consistent with zero. The Faraday rotation from anisotropic cosmic birefringence can be compared with the equivalent quantity generated by primordial magnetic fields if they existed. The POLARBEAR nondetection translates into a 95% confidence level (C.L.) upper limit of 93 nanogauss (nG) on the amplitude of an equivalent primordial magnetic field inclusive of systematic uncertainties. This four-point correlation constraint on Faraday rotation is about 15 times tighter than the upper limit of 1380 nG inferred from constraining the contribution of Faraday rotation to two-point correlations of $B$-modes measured by Planck in 2015. Metric perturbations sourced by primordial magnetic fields would also contribute to the $B$-mode power spectrum. Using the POLARBEAR measurements of the $B$-mode power spectrum (two-point correlation), we set a 95% C.L. upper limit of 3.9 nG on primordial magnetic fields assuming a flat prior on the field amplitude. This limit is comparable to what was found in the Planck 2015 two-point correlation analysis with both temperature and polarization. We perform a set of systematic error tests and find no evidence for contamination. This work marks the first time that anisotropic cosmic birefringence or primordial magnetic fields have been constrained from the ground at subdegree scales.

We present an overview of the design and status of the \Pb-2 and the Simons Array experiments. \Pb-2 is a Cosmic Microwave Background polarimetry experiment which aims to characterize the arc-minute angular scale B-mode signal from weak gravitational lensing and search for the degree angular scale B-mode signal from inflationary gravitational waves. The receiver has a 365~mm diameter focal plane cooled to 270~milli-Kelvin. The focal plane is filled with 7,588 dichroic lenslet-antenna coupled polarization sensitive Transition Edge Sensor (TES) bolometric pixels that are sensitive to 95~GHz and 150~GHz bands simultaneously. The TES bolometers are read-out by SQUIDs with 40 channel frequency domain multiplexing. Refractive optical elements are made with high purity alumina to achieve high optical throughput. The receiver is designed to achieve noise equivalent temperature of 5.8~$\mu$K$_{CMB}\sqrt{s}$ in each frequency band. \Pb-2 will deploy in 2016 in the Atacama desert in Chile. The Simons Array is a project to further increase sensitivity by deploying three \Pb-2 type receivers. The Simons Array will cover 95~GHz, 150~GHz and 220~GHz frequency bands for foreground control. The Simons Array will be able to constrain tensor-to-scalar ratio and sum of neutrino masses to $\sigma(r) = 6\times 10^{-3}$ at $r = 0.1$ and $\sum m_\nu (\sigma =1)$ to 40 meV.

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.

Atmosphere is one of the most important noise sources for ground-based cosmic microwave background (CMB) experiments. By increasing optical loading on the detectors, it amplifies their effective noise, while its fluctuations introduce spatial and temporal correlations between detected signals. We present a physically motivated 3d-model of the atmosphere total intensity emission in the millimeter and sub-millimeter wavelengths. We derive a new analytical estimate for the correlation between detectors time-ordered data as a function of the instrument and survey design, as well as several atmospheric parameters such as wind, relative humidity, temperature and turbulence characteristics. Using an original numerical computation, we examine the effect of each physical parameter on the correlations in the time series of a given experiment. We then use a parametric-likelihood approach to validate the modeling and estimate atmosphere parameters from the POLARBEAR-I project first season data set. We derive a new 1.0% upper limit on the linear polarization fraction of atmospheric emission. We also compare our results to previous studies and weather station measurements. The proposed model can be used for realistic simulations of future ground-based CMB observations.

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.

Aiming at the observation of cosmic-ray chemical composition at the "knee" energy region, we have been developinga new type air-shower core detector (YAC, Yangbajing Air shower Core detector array) to be set up at Yangbajing (90.522$^\circ$ E, 30.102$^\circ$ N, 4300 m above sea level, atmospheric depth: 606 g/m$^2$) in Tibet, China. YAC works together with the Tibet air-shower array (Tibet-III) and an underground water cherenkov muon detector array (MD) as a hybrid experiment. Each YAC detector unit consists of lead plates of 3.5 cm thick and a scintillation counter which detects the burst size induced by high energy particles in the air-shower cores. The burst size can be measured from 1 MIP (Minimum Ionization Particle) to $10^{6}$ MIPs. The first phase of this experiment, named "YAC-I", consists of 16 YAC detectors each having the size 40 cm $\times$ 50 cm and distributing in a grid with an effective area of 10 m$^{2}$. YAC-I is used to check hadronic interaction models. The second phase of the experiment, called "YAC-II", consists of 124 YAC detectors with coverage about 500 m$^2$. The inner 100 detectors of 80 cm $\times $ 50 cm each are deployed in a 10 $\times$ 10 matrix from with a 1.9 m separation and the outer 24 detectors of 100 cm $\times$ 50 cm each are distributed around them to reject non-core events whose shower cores are far from the YAC-II array. YAC-II is used to study the primary cosmic-ray composition, in particular, to obtain the energy spectra of proton, helium and iron nuclei between 5$\times$$10^{13}$ eV and $10^{16}$ eV covering the "knee" and also being connected with direct observations at energies around 100 TeV. We present the design and performance of YAC-II in this paper.

The ARGO-YBJ experiment has been in stable data taking from November 2007 till February 2013 at the YangBaJing Cosmic Ray Observatory (4300 m a.s.l.). The detector consists of a single layer of Resistive Plate Chambers (RPCs) ( about 6700 m^2}) operated in streamer mode. The signal pick-up is obtained by means of strips facing one side of the gas volume. The digital readout of the signals, while allows a high space-time resolution in the shower front reconstruction, limits the measurable energy to a few hundred TeV. In order to fully investigate the 1-10 PeV region, an analog readout has been implemented by instrumenting each RPC with two large size electrodes facing the other side of the gas volume. Since December 2009 the RPC charge readout has been in operation on the entire central carpet (about 5800 m^2). In this configuration the detector is able to measure the particle density at the core position where it ranges from tens to many thousands of particles per m^2. Thus ARGO-YBJ provides a highly detailed image of the charge component at the core of air showers. In this paper we describe the analog readout of RPCs in ARGO-YBJ and discuss both the performance of the system and the physical impact on the EAS measurements.

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

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

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

Non-Gaussianity of the primordial density perturbations provides an important measure to constrain models of inflation. At cubic order the non-Gaussianity is captured by two parameters $\tau_{\rm NL}$ and $g_{\rm NL}$ that determine the amplitude of the density perturbation trispectrum. Here we report measurements of the kurtosis power spectra of the cosmic microwave background (CMB) temperature as mapped by Planck by making use of correlations between square temperature-square temperature and cubic temperature-temperature anisotropies. In combination with noise simulations, we find the best joint estimates to be $\tau_{\rm{NL}}=0.3 \pm 0.9 \times 10^4$ and $g_{\rm{NL}}=-1.2 \pm 2.8 \times 10^5$. If $\tau_{\rm NL}=0$, we find $g_{\rm NL}= -1.3\pm 1.8 \times 10^5$.

Gas response to the underlying stellar spirals is explored for M81 using unmagnetized hydrodynamic simulations. Constrained within the uncertainty of observations, 18 simulations are carried out to study the effects of selfgravity and to cover the parameter space comprising three different sound speeds and three different arm strengths. The results are confronted with those data observed at wavelengths of 8 $\mu$m and 21 cm. In the outer disk, the ring-like structure observed in 8 $\mu$m image is consistent with the response of cold neutral medium with an effective sound speed 7 km s$^{-1}$, while for the inner disk, the presence of spiral shocks can be understood as a result of 4:1 resonances associated with the warm neutral medium with an effective sound speed 19 km s$^{-1}$. Simulations with single effective sound speed alone cannot simultaneously explain the structures in the outer and inner disks. This justifies the coexistence of cold and warm neutral media in M81. The anomalously high streaming motions observed in the northeast arm and the outward shifted turning points in the iso-velocity contours seen along the southwest arm are interpreted as signatures of interactions with companion galaxies. The level of simulated streaming motions narrows down the uncertainty of observed arm strength toward larger amplitudes.

POLARBEAR-2 is a next-generation receiver for precision measurements of the polarization of the cosmic microwave background (Cosmic Microwave Background (CMB)). Scheduled to deploy in early 2015, it will observe alongside the existing POLARBEAR-1 receiver, on a new telescope in the Simons Array on Cerro Toco in the Atacama desert of Chile. For increased sensitivity, it will feature a larger area focal plane, with a total of 7,588 polarization sensitive antenna-coupled Transition Edge Sensor (TES) bolometers, with a design sensitivity of 4.1 uKrt(s). The focal plane will be cooled to 250 milliKelvin, and the bolometers will be read-out with 40x frequency domain multiplexing, with 36 optical bolometers on a single SQUID amplifier, along with 2 dark bolometers and 2 calibration resistors. To increase the multiplexing factor from 8x for POLARBEAR-1 to 40x for POLARBEAR-2 requires additional bandwidth for SQUID readout and well-defined frequency channel spacing. Extending to these higher frequencies requires new components and design for the LC filters which define channel spacing. The LC filters are cold resonant circuits with an inductor and capacitor in series with each bolometer, and stray inductance in the wiring and equivalent series resistance from the capacitors can affect bolometer operation. We present results from characterizing these new readout components. Integration of the readout system is being done first on a small scale, to ensure that the readout system does not affect bolometer sensitivity or stability, and to validate the overall system before expansion into the full receiver. We present the status of readout integration, and the initial results and status of components for the full array.

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 DAMPE (DArk Matter Particle Explorer) is a scientific satellite being developed in China, aimed at cosmic ray study, gamma ray astronomy, and searching for the clue of dark matter particles, with a planned mission period of more than 3 years and an orbit altitude of about 500 km. The BGO Calorimeter, which consists of 308 BGO (Bismuth Germanate Oxid) crystal bars, 616 PMTs (photomultiplier tubes) and 1848 dynode signals, has approximately 32 radiation lengths. It is a crucial sub-detector of the DAMPE payload, with the functions of precisely measuring the energy of cosmic particles from 5 GeV to 10TeV, distinguishing positrons/electrons and gamma rays from hadron background, and providing trigger information for the whole DAMPE payload. The dynamic range for a single BGO crystal is about 2?105 and there are 1848 detector signals in total. To build such an instrument in space, the major design challenges for the readout electronics come from the large dynamic range, the high integrity inside the very compact structure, the strict power supply budget and the long term reliability to survive the hush environment during launch and in orbit. Currently the DAMPE mission is in the end of QM (Qualification Model) stage. This paper presents a detailed description of the readout electronics for the BGO calorimeter.

Why the cosmological constant $\Lambda$ observed today is so much smaller than the Planck scale or why the universe is accelerating at present? This is so-called the cosmological constant fine-tuning problem. In this paper, we find that this problem is solved with the help of Higgs inflation by simply assuming a variable cosmological "constant" during the inflation epoch. In the meanwhile, it could predict a large tensor-to-scalar ratio $r\approx 0.20$ and a large running of spectral index $n'_s \approx -0.028$ with a red-tilt spectrum $n_s \approx 0.96$, as well as a big enough number of e-folds $N\approx 40$ that required to solve the problems in the Big Bang cosmology with the help of $\Lambda$.

The measurement of the tensor-to-scalar ratio $r$ shows a very powerful constraint to theoretical inflation models through the detection of B-mode. In this paper, we propose a single inflation model with infinity pows series potential called the Infinity Power Series (IPS) inflation model, which is well consistent with latest observations from \textit{Planck} and BICEP2. Furthermore, it is found that in the IPS model, the absolute value the running of the spectral index increases while the tensor-to-scalar ratio becomes large, namely both large $r \approx 0.20$ and $|n_s'|\gtrsim \mathcal{O}(10^{-2})$ are realized in the IPS model. In the meanwhile, the number of e-folds is also big enough $N\approx 35\sim45$ to solve the problems in Big Bang cosmology.

Gravitational lensing due to the large-scale distribution of matter in the cosmos distorts the primordial Cosmic Microwave Background (CMB) and thereby induces new, small-scale $B$-mode polarization. This signal carries detailed information about the distribution of all the gravitating matter between the observer and CMB last scattering surface. We report the first direct evidence for polarization lensing based on purely CMB information, from using the four-point correlations of even- and odd-parity $E$- and $B$-mode polarization mapped over $\sim30$ square degrees of the sky measured by the POLARBEAR experiment. These data were analyzed using a blind analysis framework and checked for spurious systematic contamination using null tests and simulations. Evidence for the signal of polarization lensing and lensing $B$-modes is found at 4.2$\sigma$ (stat.+sys.) significance. The amplitude of matter fluctuations is measured with a precision of $27\%$, and is found to be consistent with the Lambda Cold Dark Matter ($\Lambda$CDM) cosmological model. This measurement demonstrates a new technique, capable of mapping all gravitating matter in the Universe, sensitive to the sum of neutrino masses, and essential for cleaning the lensing $B$-mode signal in searches for primordial gravitational waves.

The detection of B-mode shows a very powerful constraint to theoretical inflation models through the measurement of the tensor-to-scalar ratio $r$. Higgs boson is the most likely candidate of the inflaton field. But usually, Higgs inflation models predict a small value of $r$, which is not quite consistent with the recent results from BICEP2. In this paper, we explored whether a cosmological constant energy component is needed to improve the situation. And we found the answer is yes. For the so-called Higgs chaotic inflation model with a quadratic potential, it predicts $r\approx 0.2$, $n_s\approx0.96$ with e-folds number $N\approx 56$, which is large enough to overcome the problems such as the horizon problem in the Big Bang cosmology. The required energy scale of the cosmological constant is roughly $\Lambda \sim (10^{14} \text{GeV})^2 $, which means a mechanism is still needed to solve the fine-tuning problem in the later time evolution of the universe, e.g. by introducing some dark energy component.

The power spectra of the scalar and tensor perturbations in the noncommutative k-inflation model are calculated in this paper. In this model, all the modes created when the stringy space-time uncertainty relation is satisfied are generated inside the sound/Hubble horizon during inflation for the scalar/tensor perturbations. It turns out that a linear term describing the noncommutative space-time effect contributes to the power spectra of the scalar and tensor perturbations. Confronting the general noncommutative k-inflation model with latest results from \textit{Planck} and BICEP2, and taking $c_S$ and $\lambda$ as free parameters, we find that it is well consistent with observations. However, for the two specific models, i.e. the tachyon and DBI inflation models, it is found that the DBI model is not favored, while the tachyon model lies inside the $1\sigma$ contour, if the e-folds number is assumed to be around $50\sim60$.

In this paper, we propose a new method to calculate the mode functions in the noncommutative power-law inflation model. In this model, all the modes created when the stringy space-time uncertainty relation is satisfied are generated inside the Hubble horizon during inflation. It turns out that a linear term describing the noncommutative space-time effect contributes to the power spectra of the scalar and tensor perturbations. Confronting this model with latest results from \textit{Planck} and BICEP2, we constrain the parameters in this model and we find it is well consistent with observations.

Usually, in order to investigate the evolution of a theory, one may find the critical points of the system and then perform perturbations around these critical points to see whether they are stable or not. This local method is very useful when the initial values of the dynamical variables are not far away from the critical points. Essentially, the nonlinear effects are totally neglected in such kind of approach. Therefore, one can not tell whether the dynamical system will evolute to the stable critical points or not when the initial values of the variables do not close enough to these critical points. Furthermore, when there are two or more stable critical points in the system, local analysis can not provide the informations that which one the system will finally evolute to. In this paper, we have further developed the nullcline method to study the bifurcation phenomenon and global dynamical behaviour of the $f(T)$ theory. We overcome the shortcoming of local analysis. And it is very clear to see the evolution of the system under any initial conditions.

We reconstruct the gravitational lensing convergence signal from Cosmic Microwave Background (CMB) polarization data taken by the POLARBEAR experiment and cross-correlate it with Cosmic Infrared Background (CIB) maps from the Herschel satellite. From the cross-spectra, we obtain evidence for gravitational lensing of the CMB polarization at a statistical significance of 4.0$\sigma$ and evidence for the presence of a lensing $B$-mode signal at a significance of 2.3$\sigma$. We demonstrate that our results are not biased by instrumental and astrophysical systematic errors by performing null-tests, checks with simulated and real data, and analytical calculations. This measurement of polarization lensing, made via the robust cross-correlation channel, not only reinforces POLARBEAR auto-correlation measurements, but also represents one of the early steps towards establishing CMB polarization lensing as a powerful new probe of cosmology and astrophysics.

Aside from the grand-design stellar spirals appearing in the disk of M81, a pair of stellar spiral arms situated well inside the bright bulge of M81 has been recently discovered by Kendall et al. (2008). The seemingly unrelated pairs of spirals pose a challenge to the theory of spiral density waves. To address this problem, we have constructed a three component model for M81, including the contributions from a stellar disk, a bulge, and a dark matter halo subject to observational constraints. Given this basic state for M81, a modal approach is applied to search for the discrete unstable spiral modes that may provide an understanding for the existence of both spiral arms. It is found that the apparently separated inner and outer spirals can be interpreted as a single trailing spiral mode. In particular, these spirals share the same pattern speed 25.5 km s$^{-1}$ kpc$^{-1}$ with a corotation radius of 9.03 kpc. In addition to the good agreement between the calculated and the observed spiral pattern, the variation of the spiral amplitude can also be naturally reproduced.

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.

Some experiments have found multi-core events in cosmic rays extensive air showers which should be interpreted by hadronic interaction theory. In this paper, the multi-core events are reproduced by Monte Carlo simulation with CORSIKA. The origin of each sub-cores is tracked back from the observation level. The interaction mechanism and original particles of sub-core are studied in this paper.

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

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

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

Recently, the vacuum energy of the QCD ghost in a time-dependent background is proposed as a kind of dark energy candidate to explain the acceleration of the Universe. In this model, the energy density of the dark energy is proportional to the Hubble parameter $H$, which is the Hawking temperature on the Hubble horizon of the Friedmann-Robertson-Walker (FRW) Universe. In this paper, we generalized this model and choice the Hawking temperature on the so-called trapping horizon, which will coincides with the Hubble temperature in the context of flat FRW Universe dominated by the dark energy component. We study the thermodynamics of Universe with this kind of dark energy and find that the entropy-area relation is modified, namely, there is an another new term besides the area term.

We cross correlate the gravitational lensing map extracted from cosmic microwave background measurements by the Wilkinson Microwave Anisotropy Probe (WMAP) with the radio galaxy distribution from the NRAO VLA Sky Survey (NVSS) by using a quadratic estimator technique. We use the full covariance matrix to filter the data, and calculate the cross-power spectra for the lensing-galaxy correlation. We explore the impact of changing the values of cosmological parameters on the lensing reconstruction, and obtain statistical detection significances at $>3\sigma$. The results of all cross correlations pass the curl null test as well as a complementary diagnostic test using the NVSS data in equatorial coordinates. We forecast the potential for Planck and NVSS to constrain the lensing-galaxy cross correlation as well as the galaxy bias. The lensing-galaxy cross-power spectra are found to be Gaussian distributed.

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.

In this paper, we propose a new class of parametrization of the equation of state of dark energy. In contrast with the famous CPL parametrization, these new parametrization of the equation of state does not divergent during the evolution of the Universe even in the future. Also, we perform a observational constraint on two simplest dark energy models belonging to this new class of parametrization, by using the Markov Chain Monte Carlo (MCMC) method and the combined latest observational data from the type Ia supernova compilations including Union2(557), cosmic microwave background, and baryon acoustic oscillation.

In this paper, we shall study the dynamical behavior of the universe accelerated by the so called Veneziano ghost dark energy component locally and globally by using the linearization and nullcline method developed in this paper. The energy density is generalized to be proportional to the Hawking temperature defined on the trapping horizon instead of Hubble horizon of the Friedmann-Robertson-Walker (FRW) universe. We also give a prediction of the fate of the universe and present the bifurcation phenomenon of the dynamical system of the universe. It seems that the universe could be dominated by dark energy at present in some region of the parameter space.

Recently, the vacuum energy of the QCD ghost in a time-dependent background is proposed as a kind of dark energy candidate to explain the acceleration of the universe. In this model, the energy density of the dark energy is proportional to the Hubble parameter $H$, which is the Hawking temperature on the Hubble horizon of the Friedmann-Robertson-Walker (FRW) universe. In this paper, we perform a constraint on the ghost dark energy model with and without bulk viscosity, by using the Markov Chain Monte Carlo (MCMC) method and the combined latest observational data from the type Ia supernova compilations including Union2.1(580) and Union2(557), cosmic microwave background, baryon acoustic oscillation, and the observational Hubble parameter data.

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.

Gravitational lensing by large scale structure introduces non-Gaussianity into the Cosmic Microwave Background and imprints a new observable, which can be used as a cosmological probe. We apply a four-point estimator to the Wilkinson Microwave Anisotropy Probe (WMAP) 7-year coadded temperature maps alone to reconstruct the gravitational lensing signal. The Gaussian bias is simulated and subtracted, and the higher order bias is investigated. We measure a gravitational lensing signal with a statistical amplitude of $\mathcal {C}$ = $1.27\pm 0.98$ using all the correlations of the W- and V-band Differencing Assemblies (DAs). We therefore conclude that WMAP 7-year data alone, can not detect lensing.

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

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

In this paper, we will use $\delta \mathcal{N}$-formalism to calculate the primordial curvature perturbation for the curvaton model with a Lagrange multiplier field. We calculate the non-linearity parameters $f_{NL}$ and $g_{NL}$ in the sudden-decay approximation in this kind of model, and we find that one could get a large non-Gaussinity even if the curvaton dominates the total energy density before it decays, and this property will make the curvaton model much richer. We also calculate the probability density function of the primordial curvature perturbation in the sudden-decay approximation, as well as some moments of it.

We have investigated the necessary conditions that prevent phantom inflation from being eternal. Allowing additionally for a nonminimal coupling between the phantom field and gravity, we present the slow-climb requirements, perform an analysis of the fluctuations, and finally we extract the overall conditions that are necessary in order to prevent eternality. Furthermore, we verify our results by solving explicitly the cosmological equations in a simple example of an exponential potential, formulating the classical motion plus the stochastic effect of the fluctuations through Langevin equations. Our analysis shows that phantom inflation can be finite without the need of additional exotic mechanisms.

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

Several Cardassian universe models including the original, modified polytropic and exponential Cardassian models are constrained by the latest Constitution Type Ia supernova data, the position of the first acoustic peak of CMB from the five years WMAP data and the size of baryonic acoustic oscillation peak from the SDSS data. Both the spatial flat and curved universe are studied, and we also take account of the possible bulk viscosity of the matter fluid in the flat universe case.

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

The possibility that the non-minimal coupling inflation could be eternal is investigated. We calculate the quantum fluctuation of the inflaton in a Hubble time and find that it has the same value as in the minimal case in the slow-roll limit. Armed with this result, we have studied some concrete non-minimal inflationary models including the chaotic inflation and the natural inflation while the inflaton is non-minimally coupled to the gravity and we find that these non-minimal inflations could be eternal in some parameter regions.

The class of Chaplygin gas models regarded as a candidate of dark energy can be realized by a scalar field, which could drive the variation of the fine structure constant $\alpha$ during the cosmic time. This phenomenon has been observed for almost ten years ago from the quasar absorption spectra and attracted many attentions. In this paper, we reconstruct the class of Chaplygin gas models to a kind of scalar fields and confront the resulting $\Delta\alpha/\alpha$ with the observational constraints. We found that if the present observational value of the equation of state of the dark energy was not exactly equal to -1, various parameters of the class of Chaplygin gas models are allowed to satisfy the observational constraints, as well as the equivalence principle is also respected.

In the holographic Ricci dark energy (RDE) model, the parameter $\alpha$ plays an important role in determining the evolutionary behavior of the dark energy. When $\alpha<1/2$, the RDE will exhibit a quintom feature, i.e., the equation of state of dark energy will evolve across the cosmological constant boundary $w=-1$. Observations show that the parameter $\alpha$ is indeed smaller than 1/2, so the late-time evolution of RDE will be really like a phantom energy. Therefore, it seems that the big rip is inevitable in this model. On the other hand, the big rip is actually inconsistent with the theoretical framework of the holographic model of dark energy. To avoid the big rip, we appeal to the extra dimension physics. In this Letter, we investigate the cosmological evolution of the RDE in the braneworld cosmology. It is of interest to find that for the far future evolution of RDE in a Randall-Sundrum braneworld, there is an attractor solution where the steady state (de Sitter) finale occurs, in stead of the big rip.

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

We investigate the viscous Ricci dark energy (RDE) model by assuming that there is bulk viscosity in the linear barotropic fluid and the RDE. In the RDE model without bulk viscosity, the universe is younger than some old objects at some redshifts. Since the age of the universe should be longer than any objects in the universe, the RDE model suffers the age problem, especially when we consider the object APM 08279+5255 at $z=3.91$, whose age is $t = 2.1$ Gyr. In this letter, we find that once the viscosity is taken into account, this age problem is alleviated.

We investigate the Ricci Dark Energy (RDE) in the braneworld models with a Gauss-Bonnet term in the Bulk. We solve the generalized Friedmann equation on the brane analytically and find that the universe will finally enter into a pure de Sitter spacetime in stead of the big rip that appears in the usual 4D Ricci dark energy model with parameter $\alpha<1/2$. We also consider the Hubble horizon as the IR cutoff in holographic dark energy model and find it can not accelerate the universe as in the usual case without interacting.

The Ricci dark energy (RDE) proposed to explain the accelerating expansion of the universe requires its parameter $\alpha < 1$, whose value will determine the behavior of RDE. In this Letter, we study the scalar perturbation of RDE with and without matter in the universe, and we find that in both cases, the perturbation is stable if $\alpha> 1/3$, which gives a lower bound for $\alpha$ theoretically.

Statefinder diagnostic is a useful method which can differ one dark energy model from each others. In this letter, we apply this method to a holographic dark energy model from Ricci scalar curvature, called the Ricci dark energy model(RDE). We plot the evolutionary trajectories of this model in the statefinder parameter-planes, and it is found that the parameter of this model plays a significant role from the statefinder viewpoint. In a very special case, the statefinder diagnostic fails to discriminate LCDM and RDE models, thus we apply a new diagnostic called the Om diagnostic proposed recently to this model in this case in Appendix A and it works well.

We examine different phenomenological interaction models for Dark Energy and Dark Matter by performing statistical joint analysis with observational data arising from the 182 Gold type Ia supernova samples, the shift parameter of the Cosmic Microwave Background given by the three-year Wilkinson Microwave Anisotropy Probe observations, the baryon acoustic oscillation measurement from the Sloan Digital Sky Survey and age estimates of 35 galaxies. Including the time-dependent observable, we add sensitivity of measurement and give complementary results for the fitting. The compatibility among three different data sets seem to imply that the coupling between dark energy and dark matter is a small positive value, which satisfies the requirement to solve the coincidence problem and the second law of thermodynamics, being compatible with previous estimates.

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

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

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

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We have a plan to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV energy region by discriminating them from cosmic-ray hadrons. We have deduced the attainable sensitivity of the muon detector array using our Monte Carlo simulation. We report here on the detailed procedure of our Monte Carlo simulation.

A general holographic relation between UV and IR cutoff of an effective field theory is proposed. Taking the IR cutoff relevant to the dark energy as the Hubble scale, we find that the cosmological constant is highly suppressed by a numerical factor and the fine tuning problem seems alleviative. We also use different IR cutoffs to study the case in which the universe is composed of matter and dark energy.

Using the data coming from the new 182 Gold type Ia supernova samples, the shift parameter of the Cosmic Microwave Background given by the three-year Wilkinson Microwave Anisotropy Probe observations, and the baryon acoustic oscillation measurement from the Sloan Digital Sky Survey, $H(z)$ and lookback time measurements, we have performed a statistical joint analysis of the interacting holographic dark energy model. Consistent parameter estimations show us that the interacting holographic dark energy model is a viable candidate to explain the observed acceleration of our universe.