We measure the expansion rate of the recent Universe and the calibration scale of the baryon acoustic oscillation (BAO) from low-redshift data. BAO relies on the calibration scale, i.e., the sound horizon at the end of drag epoch $r_d$, which often imposes a prior of the cosmic microwave background (CMB) measurement from the Planck satellite. In order to make really independent measurements of $H_0$, we leave $r_d$ completely free and use the BAO datasets combined with the 31 observational $H(z)$ data (OHD) and GW170817. For the two model-independent reconstructions of $H(z)$, we obtain $H_0=69.66^{+5.88}_{-6.63}$ km s$^{-1}$ Mpc$^{-1}$, $r_d=148.56^{+3.65}_{-4.08}$ Mpc in the cubic expansion, and $H_0=71.13\pm2.91$ km s$^{-1}$ Mpc$^{-1}$, $r_d=148.48^{+3.73}_{-3.74}$ Mpc in the polynomial expansion, and we find that the values of sound horizon $r_d$ are consistent with the estimate derived from the Planck CMB data assuming a flat $\Lambda$CDM model.

Cosmic Microwave Background (CMB) anisotropy encodes a lot of information about our Universe. In this paper we take the ground-based CMB observations (GCMB), including the South Pole Telescope (SPT), SPTpol and the Atacama Cosmology Telescope Polarimeter (ACTPol), as a new probe to the CMB anisotropy independent of two satellite observations, i.e. Wilkinson Microwave Anisotropy Probe (WMAP) and Planck. The combination of current GCMB data is consistent with WMAP and Planck. In the spatially flat $\Lambda$CDM model, the Hubble constant is $H_0=69.72\pm 1.63$ km/s/Mpc at $68\%$ confidence level (CL). Combining with baryon acoustic oscillation (BAO) and the Pantheon sample of Type Ia supernovae (SN), we find that $H_0=68.40\pm 0.58$ km/s/Mpc ($68\%$ CL) in the spatially flat $\Lambda$CDM cosmology which has a tension with local measurement given by Riess et al. in 2019 at $3.7\sigma$ level, and $\Omega_k=-0.0013\pm 0.0039$ and $N_{\rm{eff}}=2.90\pm 0.41$ ($68\%$ CL) in the extended cosmological models.

The Tianlai Cylinder Pathfinder is a radio interferometer array designed to test techniques for 21 cm intensity mapping in the post-reionization Universe, with the ultimate aim of mapping the large scale structure and measuring cosmological parameters such as the dark energy equation of state. Each of its three parallel cylinder reflectors is oriented in the north-south direction, and the array has a large field of view. As the Earth rotates, the northern sky is observed by drift scanning. The array is located in Hongliuxia, a radio-quiet site in Xinjiang, and saw its first light in September 2016. In this first data analysis paper for the Tianlai cylinder array, we discuss the sub-system qualification tests, and present basic system performance obtained from preliminary analysis of the commissioning observations during 2016-2018. We show typical interferometric visibility data, from which we derive the actual beam profile in the east-west direction and the frequency band-pass response. We describe also the calibration process to determine the complex gains for the array elements, either using bright astronomical point sources, or an artificial on site calibrator source, and discuss the instrument response stability, crucial for transit interferometry. Based on this analysis, we find a system temperature of about 90 K, and we also estimate the sensitivity of the array.

The detection of binary black hole coalescences by LIGO/Virgo has aroused the interest in primordial black holes (PBHs), because they could be both the progenitors of these black holes and a compelling candidate of dark matter (DM). PBHs are formed soon after the enhanced scalar perturbations re-enter horizon during radiation dominated era, which would inevitably induce gravitational waves as well. Searching for such scalar induced gravitational waves (SIGWs) provides an elegant way to probe PBHs. We perform the first direct search for the signals of SIGWs accompanying the formation of PBHs in North American Nanohertz Observatory for Gravitational waves (NANOGrav) 11-year data set. No statistically significant detection has been made, and hence we place a stringent upper limit on the abundance of PBHs at $95\%$ confidence level. In particular, less than one part in a million of the total DM mass could come from PBHs in the mass range of $[2 \times 10^{-3}, 7\times 10^{-1}] \Msun$.

Assuming that primordial black holes compose a fraction of dark matter, some of them may accumulate at the center of galaxy and perform a prograde or retrograde orbit against the gravity pointing towards the center exerted by the central massive black hole. If the mass of primordial black holes is of the order of stellar mass or smaller, such extreme mass ratio inspirals can emit gravitational waves and form a background due to incoherent superposition of all the contributions of the Universe. We investigate the stochastic gravitational-wave background energy density spectra from the directional source, the primordial black holes surrounding Sagittarius A$^\ast$ of the Milky Way, and the isotropic extragalactic total contribution, respectively. As will be shown, the resultant stochastic gravitational-wave background energy density shows different spectrum features such as the peak positions in the frequency domain for the above two kinds of sources. Detection of stochastic gravitational-wave background with such a feature may provide evidence for the existence of primordial black holes. Conversely, a null searching result can put constraints on the abundance of primordial black holes in dark matter.

The binary black holes (BBHs) formed near the supermassive black holes (SMBHs) in the galactic nuclei would undergo eccentricity excitation due to the gravitational perturbations from the SMBH and therefore merger more efficiently. In this paper, we study the coupling of the three body 1st post-Newtonian (pN) effects with the spin effects from the SMBH in the hierarchical triple system. The three body 1pN effects yielding the de-Sitter precession is usually decoupled in the secular dynamics, while it couples to the spin of SMBH through the Lense-Thirring precession of the outer orbital plane. This coupling includes both the precessions of the inner orbit angular momentum and the Runge-Lenz vector around the outer orbit angular momentum in a general reference frame. Our general argument on the coupling of the three body 1pN effects in three body systems could be extended to any other situation as long as the outer orbital plane evolves.

We analytically calculate the scalar induced gravitational waves (SIGWs) and find a log-dependent slope of SIGW in the infrared regions $(f<f_c)$, namely $n_{\mathrm{GW}}(f)=3-2/\ln(f_c/f)$, and $n_{\mathrm{GW}}(f)=2-2/\ln(f_c/f)$ near the peak if the power spectrum of scalar curvature perturbation is quite narrow, where $f_c$ is roughly the frequency at the peak of SIGW. Such a log-dependent slope can be taken as a new template for distinguishing SIGW from other sources.

PAON4 is an L-band (1250-1500 MHz) small interferometer operating in transit mode deployed at the Nan\c{c}ay observatory in France, designed as a prototype instrument for Intensity Mapping. It features four 5~meter diameter dishes in a compact triangular configuration, with a total geometric collecting area of $\sim75 \mathrm{m^2}$, and equipped with dual polarization receivers. A total of 36 visibilities are computed from the 8 independent RF signals by the software correlator over the full 250~MHz RF band. The array operates in transit mode, with the dishes pointed toward a fixed declination, while the sky drifts across the instrument. Sky maps for each frequency channel are then reconstructed by combining the time-dependent visibilities from the different baselines observed at different declinations. This paper presents an overview of the PAON4 instrument design and goals, as a prototype for dish arrays to map the Large Scale Structure in radio, using intensity mapping of the atomic hydrogen $21~\mathrm{cm}$ line. We operated PAON4 over several years and use data from observations in different periods to assess the array performance. We present preliminary analysis of a large fraction of this data and discuss crucial issues for this type of instrument, such as the calibration strategy, instrument response stability, and noise behaviour.

In this letter we calculate the scalar induced gravitational waves (SIGWs) accompanying the formation of primordial black hole during the radiation dominated era in three different gauges, i.e. synchronous gauge, Newton gauge and uniform curvature gauge, and we find that the energy density spectra of SIGWs, $\ogw(k)$, are identical in these three different gauges.

Adopting the low-redshift observational datasets, including the Pantheon sample of Type Ia supernovae, baryon acoustic oscillation measurements, and the tomographic Alcock-Paczynski method, we determine the Hubble constant to be $67.95^{+0.78}_{-1.03}$ , $69.81^{+2.22}_{-2.70}$ and $66.75^{+3.42}_{-4.23}$ km s$^{-1}$ Mpc$^{-1}$ at 68\% confidence level in the $\Lambda$CDM, $w$CDM and $w_0w_a$CDM models, respectively. Compared to the Hubble constant given by Riess et al. in 2019, we conclude that the new physics beyond the standard $\Lambda$CDM model is needed if all of these datasets are reliable.

Axion monodromy is a generic phenomenon in string compactifications. We investigate the production of gravitational-wave background from axion monodromy inflation during preheating. By performing lattice simulations using pseudo-spectral methods, we find that significant single-peak stochastic gravitational-wave backgrounds are generated during preheating for asymptotic linear $|\phi|$ and $\phi^{2/3}$ potentials in axion monodromy inflation and cuspy potentials. This may allow us to explore string compactifications through gravitational waves.

Recent theoretical studies suggest that stellar-mass binary black holes (BBHs) would merge more efficiently due to the Kozai-Lidov mechanism if these binaries form in the vicinity of supermassive black holes (SMBHs). Since SMBHs are likely rotating rapidly, we continue our earlier study on the generalization of the Kozai-Lidov formalism to include the spin of the SMBH and study the evolution of a nearby BBH. We find that the eccentricity and orbital inclination of the BBH is significantly affected, because the spin (i) forces the orbital plane of the center-of-mass of the BBH around the SMBH to precess (the Lense-Thirring effect) and (ii) imposes an additional gravitomagnetic force on the BBH. As a result, the merger time of the BBH could be significantly different. We calculate the waveform from the BBH in one representative example and study its detectability by a milli-Hertz GW detector, such as the Laser Interferometer Space Antenna (LISA). We find that the signal is distinguishable from that in the case without spin. Our results imply that the BBHs in the LISA band could potentially be used to probe the spin of the SMBHs in galaxy centers.

Primordial gravitational waves generated during inflation lead to the B-mode polarization in the cosmic microwave background and a stochastic gravitational wave background in the Universe. We will explore the current constraint on the tilt of primordial gravitational-wave spectrum, and forecast how the future observations can improve the current constraint.

The possibility that primordial black holes (PBHs) represent all of the dark matter (DM) in the Universe and explain the coalescences of binary black holes detected by LIGO/Virgo has attracted a lot of attention. PBHs are generated by the enhancement of scalar perturbations which inevitably produce the induced gravitational waves (GWs). We calculate the induced GWs up to the third-order correction which not only enhances the amplitude of induced GWs, but also extends the cutoff frequency from $2k_*$ to $3k_*$. Such effects of the third-order correction lead to an around $10\%$ increase of the signal-to-noise ratio (SNR) for both LISA and pulsar timing array (PTA) observations, and significantly widen the mass range of PBHs in the stellar mass window accompanying detectable induced GWs for PTA observations including IPTA, FAST and SKA. On the other hand, the null detections of the induced GWs by LISA and PTA experiments will exclude the possibility that all of the DM is comprised of PBHs and the GW events detected by LIGO/Virgo are generated by PBHs.

This paper employs a simple model, considering just geometry and linear or quadratic limb darkening, to fit Kepler transit data via a Markov Chain Monte Carlo (MCMC) methodology for Kepler-1b, 5b, 8b, 12b, 77b, 428b, 491b, 699b, 706b, and 730b. Additional fits were made of the systems using the more sophisticated modeller Winfitter, which gives results in general agreement with the simpler model. Analysis of data with longer integration times showed biasing of the derived parameters, as expected from the literature,} leading to larger estimates for radii and reducing estimates of the system inclination.

We report a ground X-ray calibration of two X-ray telescope prototypes at the PANTER X-ray Test Facility, of the Max-Planck-Institute for Extraterrestrial Physics, in Neuried, Germany. The X-ray telescope prototypes were developed by the Institute of Precision Optical Engineering (IPOE) of Tongji University, in a conical Wolter-I configuration, using thermal glass slumping technology. Prototype #1 with 3 layers and Prototype #2 with 21 layers were tested to assess the prototypes' on-axis imaging performance. The measurement of Prototype #1 indicates a Half Power Diameter (HPD) of 82" at 1.49 keV. As for Prototype #2, we performed more comprehensive measurements of on-axis angular resolution and effective area at several energies ranging from 0.5-10 keV. The HPD and effective area are 111" and 39 cm^2 at 1.49 keV, respectively, at which energy the on-axis performance of the prototypes is our greatest concern.

We update the constraints on the power spectra of primordial curvature perturbation and tensor perturbation including Planck data 2015 (P15) and recently released BICEP2/Keck data (BK15), Baryon Acoustic Oscillation data and the Type Ia supernovae data. We find that the upper limits of tensor-to-scalar ratio are $0.061$, $0.064$ and $0.072$ at $95\%$ confidence level (CL) in the $\Lambda$CDM+$r$, $\Lambda$CDM+$r$+$\alpha_s$ and $\Lambda$CDM+$r$+$\alpha_s$+$\beta_s$ models respectively, where $\alpha_s$ and $\beta_s$ are the running of scalar spectral index and running of running. The inflation model with a concave potential is favored at more than $95\%$ CL. In addition, parametrizing the slow-roll parameter $\epsilon\sim 1/N^p$, where $N$ is the e-folding number before the end of inflation and taken in the range of $[50,60]$ and $[14,75]$ respectively, we conclude that the inflation model with a monomial potential $V(\phi)\sim \phi^n$ is disfavored at more than $95\%$ CL, and both the Starobinsky inflation model and brane inflation model are still consistent with the data.

We investigate how the next generation gravitational-wave (GW) detectors, such as Einstein Telescope (ET) and Cosmic Explorer (CE), can be used to distinguish primordial black holes (PBHs) from astrophysical black holes (ABHs). Since a direct detection of sub-solar mass black holes can be taken as the smoking gun for PBHs, we figure out the detectable limits of the abundance of sub-solar mass PBHs in cold dark matter by the targeted search for sub-solar mass PBH binaries and binaries containing a sub-solar mass PBH and a super-solar mass PBH, respectively. On the other hand, according to the different redshift evolutions of merger rate for PBH binaries and ABH binaries, we forecast the detectable event rate distributions for the PBH binaries and ABH binaries by ET and CE respectively, which can serve as a method to distinguish super-solar mass PBHs from ABHs.

The advent of gravitational wave (GW) and multi-messenger astronomy has stimulated the research on the formation mechanisms of binary black holes (BBHs) observed by LIGO/Virgo. In literature, the progenitors of these BBHs could be stellar-origin black holes (sBHs) or primordial black holes (PBHs). In this paper we calculate the Stochastic Gravitational-Wave Background (SGWB) from BBHs, covering the astrophysical and primordial scenarios separately, together with the one from binary neutron stars (BNSs). Our results indicate that PBHs contribute a stronger SGWB than that from sBHs, and the total SGWB from both BBHs and BNSs has a high possibility to be detected by the future observing runs of LIGO/Virgo and LISA. On the other hand, the SGWB from BBHs and BNSs also contributes an additional source of confusion noise to LISA's total noise curve, and then weakens LISA's detection abilities. For instance, the detection of massive black hole binary (MBHB) coalescences is one of the key missions of LISA, and the largest detectable redshift of MBHB mergers can be significantly reduced.

The tomographic Alcock-Paczynski (AP) test is a robust large-scale structure (LSS) measurement that receives little contamination from the redshift space distortion (RSD). It has placed tight cosmological constraints by using small and intermediate clustering scales of the LSS data. However, previous works have neglected the cross-correlation among different redshift bins, which could cause the statistical uncertainty being underestimated by $\sim$20\%. In this work, we further improve this method by including this multi-redshifts full correlation. We apply it to the SDSS DR12 galaxies sample and find out that, for $\Lambda$CDM, the combination of AP with the Planck+BAO dataset slightly reduces (within 1-$\sigma$) $\Omega_m$ to $0.304\pm0.007$ (68.3\% CL). This then leads to a larger $H_0$ and also mildly affects $\Omega_b h^2$, $n_s$ and the derived parameters $z_*$, $r_*$, $z_{re}$ but not $\tau$, $A_s$ and $\sigma_8$. For the flat $w$CDM model, our measurement gives $\Omega_m=0.301\pm 0.010$ and $w=-1.090\pm 0.047$, where the additional AP measurement reduces the error budget by $\sim 25\%$. When including more parameters into the analysis, the AP method also improves the constraints on $\Omega_k$, $\sum m_\mu$, $N_{\rm eff}$ by $20-30\%$. Early universe parameters such as $dn_s/d{\rm ln}k$ and $r$, however, are unaffected. Assuming the dark energy equation of state $w=w_0+w_a \frac{z}{1+z}$, the Planck+BAO+SNIa+$H_0$+AP datasets prefer a dynamical dark energy at $\approx1.5 \sigma$ CL. Finally, we forecast the cosmological constraints expected from the DESI galaxy survey and find that combining AP with CMB+BAO method would improve the $w_0$-$w_a$ constraint by a factor of $\sim 10$.

The leading order effect of the spin of a distant supermassive black hole (SMBH) on the orbit of compact binary revolving around it appears at the 1.5 post Newtonian (PN) expansion through coupling with the inner and outer orbital angular momenta, and makes an oscillating characteristic contribution to the secular evolution of the compact binary eccentricity over a long time scale compared to the Kozai-Lidov dynamics caused by a non-spinning SMBH.

The neutral hydrogen 21cm line is potentially a very powerful probe of the observable universe, and a number of on-going experiments are trying to detect it at cosmological distances. However, the presence of strong foreground radiations such as the galactic synchrotron radiation, galactic free-free emission and extragalactic radio sources make it a very challenging task. For the design of 21cm experiments and analysis of their data, simulation is an essential tool, and good sky foreground model is needed. With existing data the whole sky maps are available only in low angular resolutions or for limited patches of sky, which is inadequate in the simulation of these new 21cm experiments. In this paper, we present the method of constructing a high resolution self-consistent sky model at low frequencies, which incorporates both diffuse foreground and point sources. Our diffuse map is constructed by generating physical foreground components including the galactic synchrotron emission and galactic free-free emission. The point source sample is generated using the actual data from the NRAO VLA Sky Survey (NVSS) and the Sydney University Molonglo Sky Survey (SUMSS) where they are available and complete in flux limit, and mock point sources according to statistical distributions. The entire model is made self-consistent by removing the integrated flux of the point sources from the diffuse map so that this part of radiation is not double counted. We show that with the point sources added, a significant angular power is introduced in the mock sky map, which may be important for foreground subtraction simulations. Our sky maps and point source catalogues are available to download.

In this paper we present the enhanced X-ray Timing and Polarimetry mission - eXTP. eXTP is a space science mission designed to study fundamental physics under extreme conditions of density, gravity and magnetism. The mission aims at determining the equation of state of matter at supra-nuclear density, measuring effects of QED, and understanding the dynamics of matter in strong-field gravity. In addition to investigating fundamental physics, eXTP will be a very powerful observatory for astrophysics that will provide observations of unprecedented quality on a variety of galactic and extragalactic objects. In particular, its wide field monitoring capabilities will be highly instrumental to detect the electro-magnetic counterparts of gravitational wave sources. The paper provides a detailed description of: (1) the technological and technical aspects, and the expected performance of the instruments of the scientific payload; (2) the elements and functions of the mission, from the spacecraft to the ground segment.

We report the constraints of $H_0$ obtained from Wilkinson Microwave Anisotropy Probe (WMAP) combined with the latest baryonic acoustic oscillations (BAO) measurements. We use the BAO measurements from the 6dF Galaxy Survey (6dFGS), the SDSS DR7 main galaxies sample (MGS), the BOSS DR12 galaxies and the eBOSS DR14 quasars. Adding the recent BAO measurements to the cosmic microwave background (CMB) data from WMAP, we constrain cosmological parameters $\Omega_m=0.298\pm0.005$, $H_0=68.36^{+0.53}_{-0.52}$ km s$^{-1}$ Mpc$^{-1}$, $\sigma_8=0.8170^{+0.0159}_{-0.0175}$ in a spatially flat $\Lambda$ cold dark matter ($\Lambda$CDM) model, and $\Omega_m=0.302\pm0.008$, $H_0=67.63\pm1.30$ km s$^{-1}$ Mpc$^{-1}$, $\sigma_8=0.7988^{+0.0345}_{-0.0338}$ in a spatially flat $w$CDM model, respectively. The combined constraint on $w$ from CMB and BAO in a spatially flat $w$CDM model is $w=-0.96\pm0.07$. Our measured $H_0$ results prefer a value lower than 70 km s$^{-1}$ Mpc$^{-1}$, consistent with the recent data on CMB constraints from Planck (2018), but in $3.1\sim 3.5\sigma$ tension with local measurements of Riess et al. (2018) in $\Lambda$CDM and $w$CDM framework, respectively. Compared with the WMAP alone analysis, the WMPA+BAO analysis reduces the error bar by 75.4% in $\Lambda$CDM model and 95.3% in $w$CDM model.

In the measurements of cosmic microwave background polarizations, three frequency channels are necessary for discriminating the primordial B-modes from the polarized dust and the synchrotron emission. We carry out an optimistic estimate on the sensitivity to the detection of primordial gravitational waves using the cosmic microwave background B-modes only, and explore how to reach the thresholds for the tensor-to-scalar ratio $r$ in the theoretically well-motivated inflation models. For example, Lyth bound implies $r \simeq 2\times10^{-3}$, a corrected Lyth bound shows $r \simeq 7\times10^{-4}$, and some typical inflation models gives $r\simeq4\times10^{-5}$. Taking into account the up-to-date constraints on $r$, i.e. $r_{0.05}<0.07$ at $95\%$ confidence, we find that the consistency relation $n_t=-r/8$ in the canonical single-field slow-roll inflation cannot be distinguished from the scale invariance, due to the cosmic variance in the power spectrum of cosmic microwave background B-modes. The cosmic variance places an inevitable limit on the measurements of the tensor spectral index, i.e. $\sigma_{n_t}\simeq0.01$ for $2\leqslant\ell\leqslant \ell_\text{max}=300$.

We present the distance priors from the finally released $Planck~ \textrm{TT,TE,EE}+\textrm{lowE}$ data in 2018. The uncertainties are around $40\%$ smaller than those from $Planck$ 2015 TT$+$lowP. In order to check the validity of these new distance priors, we adopt the distance priors to constrain the cosmological parameters in different dark energy models, including the $\Lambda$CDM model, the $w$CDM model and the CPL model, and conclude that the distance priors provide consistent constraints on the relevant cosmological parameters compared to those from the full $Planck$ 2018 data release.

Up to now several gravitational-wave events from the coalescences of black hole binaries have been reported by LIGO/VIRGO, and imply that black holes should have an extended mass function. We work out the merger rate distribution of primordial-black-hole binaries with a general mass function by taking into account the torques by all primordial black holes and linear density perturbations. In the future, many more coalescences of black hole binaries are expected to be detected, and the one-dimensional and two-dimensional merger rate distributions will be crucial for reconstructing the mass function of primordial black holes.

The Tianlai Pathfinder is designed to demonstrate the feasibility of using a wide field of view radio interferometers to map the density of neutral hydrogen in the Universe after the Epoch of Reionizaton. This approach, called 21~cm intensity-mapping, promises an inexpensive means for surveying the large-scale structure of the cosmos. The Tianlai Pathfinder presently consists of an array of three, 15~m $\times$ 40~m cylinder telescopes and an array of sixteen, 6~m diameter dish antennas located in a radio-quiet part of western China. The two types of arrays were chosen to determine the advantages and disadvantages of each approach. The primary goal of the Pathfinder is to make 3D maps by surveying neutral hydrogen over large areas of the sky %$20,000 {\rm deg}^2$ in two different redshift ranges: first at $1.03 > z > 0.78$ ($700 - 800$~MHz) and later at $0.21 > z > 0.12$ ($1170 - 1270$~MHz). The most significant challenge to $21$~cm intensity-mapping is the removal of strong foreground radiation that dwarfs the cosmological signal. It requires exquisite knowledge of the instrumental response, i.e. calibration. In this paper, we provide an overview of the status of the Pathfinder and discuss the details of some of the analysis that we have carried out to measure the beam function of both arrays. We compare electromagnetic simulations of the arrays to measurements, discuss measurements of the gain and phase stability of the instrument, and provide a brief overview of the data processing pipeline.

In this paper, we simulate Super-Eddington accreting massive black holes (SEAMBHs) as the candles to probe cosmology for the first time. SEAMBHs have been demonstrated to be able to provide a new tool for estimating cosmological distance. Thus, we create a series of mock data sets of SEAMBHs, especially in the high redshift region, to check their abilities to probe the cosmology. To fulfill the potential of the SEAMBHs on the cosmology, we apply the simulated data to three projects. The first is the exploration of their abilities to constrain the cosmological parameters, in which we combine different data sets of current observations such as the cosmic microwave background from {\it Planck} and type Ia supernovae from Joint Light-curve Analysis (JLA). We find that the high redshift SEAMBHs can help to break the degeneracies of the background cosmological parameters constrained by {\it Planck} and JLA, thus giving much tighter constraints of the cosmological parameters. The second uses the high redshift SEAMBHs as the complements of the low redshift JLA to constrain the early expansion rate and the dark energy density evolution in the cold dark matter frame. Our results show that these high redshift SEAMBHs are very powerful on constraining the early Hubble rate and the evolution of the dark energy density; thus they can give us more information about the expansion history of our Universe, which is also crucial for testing the $\Lambda$CDM model in the high redshift region. Finally, we check the SEAMBH candles' abilities to reconstruct the equation of state of dark energy at high redshift. In summary, our results show that the SEAMBHs, as the rare candles in the high redshift region, can provide us a new and independent observation to probe cosmology in the future.

We constrain the spectral running by combining cosmic microwave background (CMB) data, baryon acoustic oscillation (BAO) data and the constraint from primordial black holes (PBHs). We find that the constraint from PBHs has a significant impact on the running of running of scalar spectral index, and a power-law scalar power spectrum without running is consistent with observational data once the constraint from PBHs is taken into account. In addition, from the constraints on the slow-roll parameters, the derived tensor spectral index in the single-field slow-roll inflation model is quite small, namely $|n_t|\lesssim 9.3\times 10^{-3}$ which will be very difficult to be measured by CMB data only in the future, and the absolute value of derived running of tensor spectral index is not larger than $2.1\times 10^{-4}$ at 95% confidence level.

Extracting the neutral hydrogen (HI) signal is a great challenge for cosmological 21cm experiments, both the astrophysical foregrounds and the receiver noise are typically several orders of magnitude greater than the 21cm signal. However, the different properties of the 21cm signal, foreground, and noise can be exploited to separate these components. The foregrounds are generally smooth or correlated over the frequency space along a line of sight (l.o.s.), while both the 21cm signal and the noise varies stochastically along the same l.o.s. The foreground can be removed by filtering out the smooth component in the frequency space. The receiver noise is basically uncorrelated for observations at different times, hence for surveys they are also uncorrelated in the different directions, while the 21cm signal which traces the large scale structure are correlated up to certain scales. In this exercise, we apply Wiener filters in frequency and angular space, to extract the 21cm signals. We found that the method works well. The inaccurate knowledge about beam could degrade the reconstruction, but the overall result is still good, showing that the method is fairly robust.

Radio astronomical observation below 30 MHz is hampered by the refraction and absorption of the ionosphere, and the radio frequency interference (RFI), so far high angular resolution sky intensity map is not available. An interferometer array on lunar orbit provides a perfect observatory in this frequency band: it is out of the ionosphere and the Moon helps to block the RFIs from the Earth. The satellites can make observations on the far side of the Moon and then send back the data on the near side part of the orbit. However, for such array the traditional imaging algorithm is not applicable: the field of view is very wide (almost whole sky), and for baselines distributed on a plane, there is a mirror symmetry between the two sides of the plane. A further complication is that for each baseline, the Moon blocks part of the sky, but as the satellites orbit the Moon, both the direction of the baseline and the blocked sky change, so even imaging algorithms which can deal with non-coplanar baseline may not work in this case. Here we present an imaging algorithm based on solving the linear mapping equations relating the sky intensity to the visibilities. We show that the mirror symmetry can be broken by the three dimensional baseline distribution generated naturally by the precession of the orbital plane of the satellites. The algorithm is applicable and good maps could be reconstructed, even though for each baseline the sky blocking by the Moon is different. We also investigate how the map-making is affected by inhomogeneous baseline distributions.

We study a model of inflation in which a scalar field $\chi$ is non-minimally coupled to Starobinsky's $R^2$ gravity. After transforming it to the Einstein frame, a new scalar field, the scalaron $\phi$, will appear and couple to $\chi$ with a nontrivial field metric, while $\chi$ acquires a positive mass via the non-minimal coupling. Initially inflation occurs along the $\phi$ direction with $\chi$ trapped near its origin by this induced mass. After $\phi$ crosses a critical value, it starts rolling down rapidly and proceeds todamped oscillations around an effective local minimum determined by the value of $\chi$, while inflation still continues, driven by the $\chi$ field at this second stage where the effect of the non-minimal coupling becomes negligible. The presence of the damped oscillations during the transition from the first to second stage of inflation causes enhancement and oscillation features in the power spectrum of the curvature perturbation. Assuming that the oscillations may be treated perturbatively, we calculate these features by using the $\delta N$ formalism, and discuss its observational implications to large scale CMB anomalies or primordial black hole formation, depending on the scale of the features.

Advanced LIGO's discovery of gravitational-wave events is stimulating extensive studies on the origin of binary black holes. Assuming that the gravitational-wave events can be explained by binary primordial black hole mergers, we utilize the upper limits on the stochastic gravitational-wave background given by Advanced LIGO as a new observational window to independently constrain the abundance of primordial black holes in dark matter. We show that Advanced LIGO's first observation run gives the best constraint on the primordial black hole abundance in the mass range $1 M_\odot \lesssim M_\text{PBH}\lesssim 100 M_\odot$, pushing the previous microlensing and dwarf galaxy dynamics constraints tighter by 1 order of magnitude. Moreover, we discuss the possibility to detect the stochastic gravitational-wave background from primordial black holes, in particular from subsolar mass primordial black holes, by Advanced LIGO in the near future.

We report a tight Hubble constant constraint $67.78^{+1.21}_{-1.86}$ km s$^{-1}$ Mpc$^{-1}$ (2.26\% precision) derived from galaxy redshift surveys. We combine the BAO measurements from 6dFGS, the SDSS DR7 main galaxies, the BOSS DR12 galaxies, and eBOSS DR14 quasars, and also apply the tomographic Alcock-Paczynski (AP) method to the BOSS DR12 galaxies, to place constraints on $H_0$ in the spatially flat $\Lambda$CDM framework. Our result is fully consistent with the CMB constraints from Planck, but in $2.58\sigma$ tension with local measurements of Riess et al. 2016. Compared with the BAO alone constraint, the BAO+AP combined result reduces the error bar by 32%. This shows the strong power of the tomographic AP method in extracting cosmological information from galaxy redshift surveys.

The dispersion relation of graviton is a fundamental issue for fundamental physics about gravity. In this paper we investigate how the modified dispersion relation of graviton affects the cosmic microwave background (CMB) power spectra, in particular the B-mode polarization. Our results will be useful to test the dispersion relation of graviton at the energy scale around $10^{-29}$ eV.

We update the constraints on the cosmological parameters by adopting the Planck data released in 2015 and Baryon Acoustic Oscillation (BAO) measurements including the new DR14 quasar sample measurement at redshift $z=1.52$, and we conclude that the based six-parameter $\Lambda$CDM model is preferred. Exploring some extensions to the $\Lambda$CDM models, we find that the equation of state of dark energy reads $w=-1.036\pm 0.056$ in the $w$CDM model, the effective relativistic degrees of freedom in the Universe is $N_\text{eff}=3.09_{-0.20}^{+0.18}$ in the $N_\text{eff}+\Lambda$CDM model and the spatial curvature parameter is $\Omega_k=(1.8\pm1.9)\times 10^{-3}$ in the $\Omega_k+\Lambda$CDM model at $68\%$ confidence level (C.L.), and the $95\%$ C.L. upper bounds on the sum of three active neutrinos masses are $\sum m_\nu<0.16$ eV for the normal hierarchy (NH) and $\sum m_\nu<0.19$ eV for the inverted hierarchy (IH) with $\Delta\chi^2\equiv \chi^2_\text{NH}-\chi^2_\text{IH}=-1.25$.

The early reionization (ERE) is supposed to be a physical process which happens after recombination, but before the instantaneous reionization caused by the first generation of stars. We investigate the effect of the ERE on the temperature and polarization power spectra of cosmic microwave background (CMB), and adopt principal components analysis (PCA) to model-independently reconstruct the ionization history during the ERE. In addition, we also discuss how the ERE affects the cosmological parameter estimates, and find that the ERE does not impose any significant influences on the tensor-to-scalar ratio $r$ and the neutrino mass at the sensitivities of current experiments. The better CMB polarization data can be used to give a tighter constraint on the ERE and might be important for more precisely constraining cosmological parameters in the future.

We use Planck data released in 2015 to constrain the abundance of primordial black holes (PBHs) in dark matter in two different reionization models (one is the instantaneous reionization and the other is the asymmetric reionization), and significantly improve the upper limits on the abundance of PBHs from WMAP 3-year data by around two orders of magnitude. Furthermore, these new limits imply that the event rates of mergers of PBH binaries (Gpc$^{-3}$ yr$^{-1}$) are less than $0.002$ for $M_\text{pbh}=30M_\odot$, $5$ for $M_\text{pbh}=10M_\odot$ and $2000$ for $M_\text{pbh}=2M_\odot$ at $95\%$ confidence level (C.L.), and thus GW150914 seems very unlikely produced by the merger of a PBH binary.

We propose a new parameterization to measure the neutrino mass hierarchy, namely $\Delta=(m_3-m_1)/(m_1+m_3)$ which is dimensionless and varies in the range $[-1,1]$. Taking into account the results of neutrino oscillation experiments, $\Delta$ is the unique parameter for determining all the masses of neutrinos, and a positive (negative) sign of $\Delta$ denotes the normal (inverted) mass hierarchy. Adopting the currently available cosmic observations, we find that the normal mass hierarchy is slightly favored, and the mass of lightest neutrino is less than $0.030$ eV for the normal mass hierarchy and $0.024$ eV for the inverted mass hierarchy at $95\%$ confidence level.

The effects of heavy fields modulate the scalar power spectrum during inflation. We analytically calculate the modulations of the scalar power spectrum from a heavy field with a separable monomial potential, i.e. V(\phi)~\phi^n. In general the modulation is characterized by a power-law oscillation which is reduced to the logarithmic oscillation in the case of n=2.

The LAMP (Lightweight Asymmetry and Magnetism Probe) X-ray telescope is a mission concept to measure the polarization of X-ray astronomical sources at 250 eV via imaging mirrors that reflect at incidence angles near the polarization angle, i.e., 45 deg. Hence, it will require the adoption of multilayer coatings with a few nanometers d-spacing in order to enhance the reflectivity. The nickel electroforming technology has already been successfully used to fabricate the high angular resolution imaging mirrors of the X-ray telescopes SAX, XMM-Newton, and Swift/XRT. We are investigating this consolidated technology as a possible technique to manufacture focusing mirrors for LAMP. Although the very good reflectivity performances of this kind of mirrors were already demonstrated in grazing incidence, the reflectivity and the scattering properties have not been tested directly at the unusually large angle of 45 deg. Other possible substrates are represented by thin glass foils or silicon wafers. In this paper we present the results of the X-ray reflectivity campaign performed at the BEAR beamline of Elettra - Sincrotrone Trieste on multilayer coatings of various composition (Cr/C, Co/C), deposited with different sputtering parameters on nickel, silicon, and glass substrates, using polarized X-rays in the spectral range 240 - 290 eV.

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

We try to reconcile the tension between the local 2.4\% determination of Hubble constant and its global determination by $\textit{Planck}$ CMB data and BAO data through modeling the dark energy variously. We find that the chi-square is significantly reduced by $\Delta\chi^2_\text{all}=-6.76$ in the redshift-binned dark energy model where the $68\%$ limits of the equation of state of dark energy read $w(0\leq z\leq 0.1)=-1.958_{-0.508}^{+0.509}$, $w(0.1< z\leq 1.5)=-1.006_{-0.082}^{+0.092}$, and here $w(z>1.5)$ is fixed to $-1$.

We provide the latest constraints on the power spectra of both scalar and tensor perturbations from the CMB data (including \textit{Planck}~2015, BICEP2 \& \textit{Keck Array} experiments) and the new BAO scales from SDSS-III BOSS observation. We find that the inflation model with a convex potential is not favored and both the inflation model with a monomial potential and the natural inflation model are marginally disfavored at around $95\%$ confidence level. But both the Brane inflation model and the Starobinsky inflation model fit the data quite well.

Both WMAP and Planck data show a significant odd-multipole preference in the large scales of the cosmic microwave background (CMB) temperature anisotropies. If this pattern originates from cosmological effects, then it can be considered a crucial clue for a violation in the cosmological principle. By defining various direction dependent statistics in the full-sky Planck 2015 maps (see, for instance, Naselsky et al. (2012); W. Zhao (2014)), we found that the CMB parity asymmetry has a preferred direction, which is independent of the choices of the statistics. In particular, this preferred axis is strongly aligned with those in the CMB quadrupole and octopole, as well as that in the CMB kinematic dipole, which hints to their non-cosmological origin. In realistic observations, the foreground residuals are inevitable, and should be properly masked out in order to avoid possible misinterpretation of the results. In this paper, we extend our previous analyses to the masked Planck 2015 data. By defining a similar direction dependent statistic in the masked map, we find a preferred direction of the CMB parity asymmetry, in which the axis also coincides with that found in the full-sky analysis. Therefore, our conclusions on the CMB parity violation and its directional properties are confirmed.

Considering the mass splitting between three active neutrinos, we represent the new constraints on the sum of neutrino mass $\sum m_\nu$ by updating the anisotropic analysis of Baryon Acoustic Oscillation (BAO) scale in the CMASS and LOWZ galaxy samples from Data Release 12 of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS DR12). Combining the BAO data of 6dFGS, MGS, LOWZ and CMASS with $\textit{Planck}$~2015 data of temperature anisotropy and polarizations of Cosmic Microwave Background (CMB), we find that the $95\%$ C.L. upper bounds on $\sum m_\nu$ refer to $\sum m_{\nu,\rm NH}<0.18$ eV for normal hierarchy (NH), $\sum m_{\nu,\rm IH}<0.20$ eV for inverted hierarchy (IH) and $\sum m_{\nu,\rm DH}<0.15$ eV for degenerate hierarchy (DH) respectively, and the normal hierarchy is slightly preferred than the inverted one ($\Delta \chi^2\equiv \chi^2_{\rm NH}-\chi^2_{\rm IH} \simeq -3.4$). In addition, the additional relativistic degrees of freedom and massive sterile neutrinos are neither favored at present.

We update the distance priors by adopting $Planck~ \textrm{TT,TE,EE}+\textrm{lowP}$ data released in 2015, and our results impose at least $30\%$ tighter constraints than those from $Planck~ \textrm{TT}+\textrm{lowP}$. Combining the distance priors with the combination of supernova Union~2.1 compilation of 580 SNe (Union~2.1) and low redshift Baryon Acoustic Oscillation (BAO) data, we constrain the cosmological parameters in the freely binned dark energy (FBDE) and FBDE$+\Omega_k$ models respectively, and find that the equations of state of dark energy in both models are consistent with $w=-1$. Furthermore, we show that the tension with the BAO data at $z=2.34$ from Ly$\alpha$ forest (Ly$\alpha$F) auto-correlation and Combined Ly$\alpha$F cannot be relaxed in the FBDE and FBDE$+\Omega_k$ models.

By taking into account the contamination of foreground radiations, we employ the Fisher matrix to forecast the future sensitivity on the tilt of power spectrum of primordial tensor perturbations for several ground-based (AdvACT, CLASS, Keck/BICEP3, Simons Array, SPT-3G), balloon-borne (EBEX, Spider) and satellite (CMBPol, COrE, LiteBIRD) experiments of B-mode polarizations. For the fiducial model $n_t=0$, our results show that the satellite experiments give good sensitivity on the tensor tilt $n_t$ to the level $\sigma_{n_t}\lesssim0.1$ for $r\gtrsim2\times10^{-3}$, while the ground-based and balloon-borne experiments give worse sensitivity. By considering the BICEP2/Keck Array and Planck (BKP) constraint on the tensor-to-scalar ratio $r$, we see that it is impossible for these experiments to test the consistency relation $n_t=-r/8$ in the canonical single-field slow-roll inflation models.

In de Rham-Gabadadze-Tolley (dRGT) massive gravity and bigravity, a non-minimal matter coupling involving both metrics generically re-introduces the Boulware--Deser (BD) ghost. A non-minimal matter coupling via a simple, yet specific composite metric has been proposed, which eliminates the BD ghost below the strong coupling scale. Working explicitly in the metric formulation and for arbitrary spacetime dimensions, we show that this composite metric is the unique consistent non-minimal matter coupling below the strong coupling scale, which emerges out of two diagnostics, namely, absence of Ostrogradski ghosts in the decoupling limit and absence of the BD ghost from matter quantum loop corrections.

We expand the Einstein-Hilbert action with a positive cosmological constant up to the fourth order in terms of gravity fluctuations, and then use the in-in formalism to calculate the four-point correlation function for gravitational waves, including both contact and exchange diagrams, generated during a period of exactly de Sitter expansion. In addition, we also present the general properties of the $n$-point function of graviton in terms of both circularly and linearly polarized states.

Imposing that the excursion distance of inflaton in field space during inflation be less than the Planck scale, we derive an upper bound on the tensor-to-scalar ratio at the CMB scales, i.e. $r_{*,max}$, in the general canonical single-field slow-roll inflation model, in particular the model with non-negligible running of the spectral index $\alpha_s$ and/or the running of running $\beta_s$. We find that $r_{*,max}\simeq 7\times 10^{-4}$ for $n_s=0.9645$ without running and running of running, and $r_{*,max}$ is significantly relaxed to the order of ${\cal O}(10^{-2}\sim 10^{-1})$ in the inflation model with $\alpha_s$ and/or $\beta_s\sim +{\cal O}(10^{-2})$ which are marginally preferred by the Planck 2015 data.

In this paper, we constrain the tilt of the power spectrum of relic gravitational waves by combining the data from BICEP2/Keck array and Planck (BKP) and the Laser Interferometer Gravitational-Waves Observatory (LIGO). From the data of BKP B-modes, the constraint on the tensor tilt is $n_t=0.66^{+1.83}_{-1.44}$ at the $68%$ confidence level. By further adding the LIGO upper limit on the energy density of gravitational waves, the constraint becomes $n_t=-0.76^{+1.37}_{-0.52}$ at the $68%$ confidence level. We conclude that there is no evidence for a blue-tilted power spectrum of relic gravitational waves and either sign of the index of tensor power spectrum is compatible with the data.

We make a joint analysis of BICEP2 and recently released Planck HFI 353 GHz dust polarization data, and find that there is no evidence for the primordial gravitational waves and the bound on the tensor-to-scalar ratio becomes $r<0.083$ at $95%$ confidence level in the base $\Lambda$CDM + tensor model. Extending to the model with running of scalar spectral index, the bound is a little bit relaxed to $r<0.116$ at $95%$ confidence level. Our results imply that the inflation model with a single monomial potential is marginally disfavored at around $95%$ confidence level. Especially, the $m^2\phi^2/2$ inflation model is disfavored at more than $2\sigma$ level. However, the Starobinsky inflation model gives a nice fit.

In this paper we constrain the cosmological parameters, in particular the tilt of tensor power spectrum, by adopting Background Imaging of Cosmic Extragalactic Polarization (B2), Planck released in 2013 (P13) and Wilkinson Microwaves Anisotropy Probe 9-year Polarization (WP) data. We find that a blue tilted tensor power spectrum is preferred at more than $3\sigma$ confidence level if the data from B2 are assumed to be totally interpreted as the relic gravitational waves, but a scale invariant tensor power spectrum is consistent with the data once the polarized dust is taken into account. The recent Planck 353 GHz HFI dust polarization data imply that the B2 data are perfectly consistent with there being no gravitational wave signal.

Even though the Hubble constant cannot be significantly determined by the low-redshift Baryon Acoustic Oscillation (BAO) data alone, it can be tightly constrained once the high-redshift BAO data are combined. Combining BAO data from 6dFGS, BOSS DR11 clustering of galaxies, WiggleZ and $z=2.34$ from BOSS DR11 quasar Lyman-$\alpha$ forest lines, we get $H_0=68.17^{+1.55}_{-1.56}$ km s$^{-1}$ Mpc$^{-1}$. In addition, adopting the the simultaneous measurements of $H(z)$ and $D_A(z)$ from the two-dimensional two-point correlation function from BOSS DR9 CMASS sample and two-dimensional matter power spectrum from SDSS DR7 sample, we obtain $H_0=68.11\pm1.69$ km s$^{-1}$ Mpc$^{-1}$. Finally, combining all of the BAO datasets, we conclude $H_0=68.11\pm 0.86$ km s$^{-1}$ Mpc$^{-1}$, a 1.3% determination.

In this paper we propose a new inflation model named (p,q) inflation model in which the inflaton potential contains both positive and negative powers of inflaton field in the polynomial form. We will show that such a simple model can easily generate a large amplitude of tensor perturbation and a negative running of spectral index with large absolute value.

We derive the analytic formula of the growth index for $f(R)$ dark energy model where the effect on the growth of matter density perturbation $\delta_m$ from modified gravity (MG) is encoded in the effective Newton coupling constant $G_{\rm eff}$ in MG (or equivalently $g\equiv {G_{\rm eff}/ G}$). Based on the analytic formula, we propose that the parameter $g$ can be directly figured out by comparing the observed growth rate $f_g\equiv d\ln\delta_m/d\ln a$ to the prediction of $f_g$ in general relativity.

Since the temperature fluctuations in cosmic microwave background (CMB) on large-angular scales probe length scales that were super-horizon sized at photon decoupling and hence insensitive to microphysical processes, the low-multipole CMB data are supposed to be a good probe to the physics of the primordial Universe. In this letter we will constrain the cosmological parameters in the base $\Lambda$CDM model with tensor perturbations by only using low-multipole CMB data, including Background Imaging of Cosmic Extragalactic Polarization (B2), Planck released in 2013 (P13) and Wilkinson Microwaves Anisotropy Probe 9-year data (W9). We find that a red tilted power spectrum of relic gravitational waves is compatible with the data, but a blue tilted power spectrum of scalar perturbations on the large scales is preferred at around $2\sigma$ confidence level.

Recently Background Imaging of Cosmic Extragalactic Polarization (B2) discovered the relic gravitational waves at $7.0\sigma$ confidence level. However, the other cosmic microwave background (CMB) data, for example Planck data released in 2013 (P13), prefer a much smaller amplitude of the primordial gravitational waves spectrum if a power-law spectrum of adiabatic scalar perturbations is assumed in the six-parameter $\Lambda$CDM cosmology. In this paper, we explore whether the $w$CDM model and the running spectral index can relax the tension between B2 and other CMB data. In particular, we find that a positive running of running of spectral index is preferred at $1.7\sigma$ level from the combination of B2, P13 and WMAP Polarization data.

Even though Planck data released in 2013 (P13) is not compatible with Background Imaging of Cosmic Extragalactic Polarization (B2) and some local cosmological observations, including Supernova Legacy Survey (SNLS) samples and $H_0$ prior from Hubble Space Telescope (HST) etc, Wilkinson Microwaves Anisotropy Probe 9-year data (W9) is consistent with all of them in the base six-parameter $\Lambda$CDM+tensor cosmology quite well. In this letter, we adopt the combinations of B2+W9 and B2+W9+SNLS+BAO+HST to constrain the cosmological parameters in the base six-parameter $\Lambda$CDM+tensor model with $n_t=-r/8$, where r and $n_t$ are the tensor-to-scalar ratio and the tilt of relic gravitational wave spectrum, and BAO denotes Baryon Acoustic Oscillation. We find that the Harrison-Zel'dovich (HZ) scale invariant scalar power spectrum is consistent with both data combinations, chaotic inflation is marginally disfavored by the data at around $2\sigma$ level, but the power-law inflation model and the inflation model with inverse power-law potential can fit the data nicely.

In this paper we constrain the tilt of the spectra of primordial gravitational waves from Background Imaging of Cosmic Extragalactic Polarization (BICEP2) data only. We find $r=0.21_{-0.10}^{+0.04}$ and $n_t=-0.06_{-0.23}^{+0.25}$ (at $68\%$ C.L.) which implies that a scale-invariant primordial gravitational waves spectra is consistent with BICEP2 nicely. Our results provide strong evidence for supporting inflation model, and the alternative models, for example the ekpyrotic model which predicts $n_t=2$, are ruled out at more than $5\sigma$ significance.

Relic gravitational waves (RGWs) leave well-understood imprints on the anisotropies in the temperature and polarization of cosmic microwave background (CMB) radiation. In the TT and TE information channels, which have been well observed by WMAP and Planck missions, RGWs compete with density perturbations mainly at low multipoles. It is dangerous to include high-multipole CMB data in the search for gravitational waves, as the spectral indices may not be constants. In this paper, we repeat our previous work [W.Zhao & L.P.Grishchuk, Phys.Rev.D {\bf 82}, 123008 (2010)] by utilizing the Planck TT and WMAP TE data in the low-multipole range $\ell\le100$. We find that our previous result is confirmed {with higher confidence}. The constraint on the tensor-to-scalar ratio from Planck TT and WMAP TE data is $r\in [0.06,~0.60]$ (68% C.L.) with the maximum likelihood at around $r\sim 0.2$. Correspondingly, the spectral index at the pivot wavenumber $k_*=0.002$Mpc$^{-1}$ is $n_s=1.13^{+0.07}_{-0.08}$, which is larger than 1 at more than $1\sigma$ level. So, we conclude that the new released CMB data indicate a stronger hint for the RGWs with the amplitude $r\sim 0.2$, which is hopeful to be confirmed by the imminent BICEP and Planck polarization data.

Recently released Planck data implies a smaller Hubble constant $H_0$ than that from Hubble Space Telescope project (HST) and a larger percentage of the matter components $\Omega_m$ compared to Supernova Legacy Survey (SNLS) in $\Lambda$CDM model. In this paper we found that even though the tension on $H_0$ between Planck and HST can be relaxed if the dark radiation is introduced ($\Delta N_{\rm eff}=0.536_{-0.224}^{+0.229}$ at $68\%$ CL from the datasets of Planck+WMAP Polarization (WP)+baryon acoustic oscillation (BAO)+the combination of supernova Union2.1 compilation of 580 SNe (Union2.1)+HST), $\Omega_m$ from Planck is still not nicely compatible with that from SNLS. The tensions between Planck and other astrophysical datasets can be significantly relaxed in $w$CDM model, and the combination of these datasets prefers a phantom-like dark energy at more than $95\%$ CL: $w=-1.15\pm 0.07$ and $w=-1.16\pm 0.06$ at $68\%$ CL from Planck+WP+BAO+Union2.1+HST and Planck+WP+BAO+SNLS+HST respectively. From the statistical point of view, there is no evidence for a time-evolving equation of state ($\Delta \chi^2=-0.3$ compared to a constant equation of state for the combination of Planck+WP+BAO+SNLS+HST).

Motivated by the ultraviolet complete theory of quantum gravity, for example the string theory, we investigate a polynomial $f(R )$ inflation model in detail. We calculate the spectral index and tensor-to-scalar ratio in the $f(R )$ inflation model with the form of $f(R )=R+{R^2\over 6M^2}+{\lambda_n\over 2n}{R^n\over (3M^2)^{n-1}}$. Compared to Planck 2013, we find that $R^n$ term should be exponentially suppressed, i.e. $|\lambda_n|\lesssim 10^{-2n+2.6}$.

It has recently been shown that the graviton can consistently gain a constant mass without introducing the Boulware-Deser ghost. We propose a gravity model where the graviton mass is set by a scalar field and prove that this model is free of the Boulware-Deser ghost by analyzing its constraint system and showing that two constraints arise. We also initiate the study of the model's cosmic background evolution and tentatively discuss possible cosmological implications of this model. In particular, we consider a simple scenario where the scalar field setting the graviton mass is identified with the inflaton and the graviton mass evolves from a high to a low energy scale, giving rise to the current cosmic acceleration.

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

We extend the four-dimensional de Rham-Gabadadze-Tolley (dRGT) massive gravity model to a general scalar massive-tensor theory in arbitrary dimensions, coupling a dRGT massive graviton to multiple scalars and allowing for generic kinetic and mass matrix mixing between the massive graviton and the scalars, and derive its Hamiltonian formulation and associated constraint system. When passing to the Hamiltonian formulation, two different sectors arise: a general sector and a special sector. Although obtained via different ways, there are two second class constraints in either of the two sectors, eliminating the BD ghost. However, for the special sector, there are still ghost instabilities except for the case of two dimensions. In particular, for the special sector with one scalar, there is a "second BD ghost".

We constrain models of single field inflation with the pre-\textit{Planck} CMB data. The data used here is the 9-year Wilkinson Microwave Anisotropy Probe (\textit{WMAP}) data, South Pole Telescope (SPT) data and Atacama Cosmology Telescope (ACT) data. By adding in running of spectral index parameter, we find that the $\chi^2$ is improved by a factor of $\Delta \chi^2=8.44$, which strongly indicates the preference of this parameter from current data. In addition, we find that the running of spectral index $\alpha_s$ does not change very much even if we switch to different pivot scales, which suggests that the power law expansion of power spectrum is accurate enough till the 1st order term. Furthermore, we find that the joint constraints on $r-n_{s}$ give very tight constraints on single-field inflation models, and the models with power law potential $\phi^{p}$ can only survive if $0.9 \lesssim p \lesssim 2.1$, so a large class of inflation models have already been ruled out before \textit{Planck} data. Finally, we use the $f_{NL}$ data to constrain the non-trivial sound speed $c_s$. We find that the current constraint is dominated by the power spectrum constraints which have some inconsistency with the constraints from $f_{NL}$. This poses important questions of consistency between power spectrum and bispectrum of \textit{WMAP} data.

In this paper, we compare brane inflation models with the Planck data and the pre-Planck data (which combines WMAP, ACT, SPT, BAO and H0 data). The Planck data prefer a spectral index less than unity at more than 5\sigma confidence level, and a running of the spectral index at around 2\sigma confidence level. We find that the KKLMMT model can survive at the level of 2\sigma only if the parameter $\beta$ (the conformal coupling between the Hubble parameter and the inflaton) is less than $\mathcal{O}(10^{-3})$, which indicates a certain level of fine-tuning. The IR DBI model can provide a slightly larger negative running of spectral index and red tilt, but in order to be consistent with the non-Gaussianity constraints from Planck, its parameter also needs fine-tuning at some level.

As a simplest extension to the mass-term curvaton model, the curvaton model with a polynomial potential can relax the restricted constraint from PLANCK due to the non-linear dynamics of curvaton field before it decays. We find that there is still a big room for producing a large negative $g_{\rm NL}$, but not positive $g_{\rm NL}$. For example, we only need around 10% "tuning" for $-g_{\rm NL}>10^4$.

We investigate the non-Gaussian signatures of ultra slow-roll inflation. The bispectrum and the trispectrum are calculated with general initial conditions. The trispectrum is of local shape, as in the case of the bispectrum. We show that the prediction of local non-Gaussianity is robust again generalizing the Bunch-Davies vacuum. The Suyama-Yamaguchi relation is saturated in this scenario.

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

We study the Standard Model Higgs field as a source for the primordial curvature perturbation, particularly in the curvaton and modulated reheating scenario. We conclude that the Higgs cannot play as a curvaton due to the small energy density when it decays, however the modulated reheating by Higgs can be a viable scenario. In the latter case, the non-Gaussianity is inevitably generated and strongly constrains the type of potential of inflaton field and Higgs-dependent interaction term. For the quadratic potential of the inflaton field with decay rate which non-linearly depends on the Higgs vacuum expectation value, the contribution of Higgs field to the primordial curvature perturbation must be less than 8%.

In this paper we investigate the perturbation theory of the asymptotically safe inflation and we find that all modes of gravitational waves perturbation become ghosts in order to achieve a large enough number of e-folds. Formally we can calculate the power spectrum of gravitational waves perturbation, but we find that it is negative. It indicates that there is serious trouble with the asymptotically safe inflation.

We make a comprehensive investigation of the observational effect of the inflation consistency relation. We focus on the general single-field inflation model with the consistency relation $r=-8c_s n_t$, and investigate the observational constraints of sound speed $c_s$ by using the Seven-Year WMAP data, the BICEP tensor power spectrum data, and the constraints on $f_{\rm NL}^{\rm equil.}$ and $f_{\rm NL}^{\rm orth.}$ from the Five-Year WMAP observations. We find that the constraints on the tensor-to-scalar ratio $r$ is much tighter if $c_s$ is small, since a large tilt $n_t$ is strongly constrained by the observations. We obtain $r<0.37, 0.27$ and 0.09 ($dn_s/d\ln k=0$) for $c_s$=1, 0.1 and 0.01 models at 95.4% confidence level. When taking smaller values of $c_s$, the positive correlation between $r$ and $n_s$ also leads to slightly tighter constraint on the upper bound of $n_s$, while the running of scalar spectral index $dn_s/d\ln k$ is generally unaffected. For the sound speed $c_s$, it is not well constrained if only the CMB power spectrum data is used, while the constraints are obtainable by taking $f_{\rm NL}^{\rm equil.}$ and $f_{\rm NL}^{\rm orth.}$ priors into account. With the constraining data of $f_{\rm NL}^{\rm equil.}$ and $f_{\rm NL}^{\rm orth.}$, we find that, $c_s\lesssim 0.01$ region is excluded at 99.7% CL, and the $c_s=1$ case (the single-field slow-roll inflation) is slightly disfavored at 68.3% CL. In addition, the inclusion of $f_{\rm NL}^{\rm equil.}$ and $f_{\rm NL}^{\rm orth.}$ into the analysis can improve the constraints on $r$ and $n_s$. We further discuss the implications of our constraints on the test of inflation models.

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

We derive the Boltzmann equation in the synchronous gauge for massive neutrinos with a deformed dispersion relation. Combining the 7-year WMAP data with lower-redshift measurements of the expansion rate, we give constraints on the deformation parameter and find that the deformation parameter is strong degenerate with the physical dark matter density instead of the neutrino mass. Our results show that there is no evidence for Lorentz invariant violation in the neutrino sector. The ongoing Planck experiment could provide improved constraints on the deformation parameter.

We explore the cosmological implications of five modified gravity (MG) models by using the recent cosmological observational data, including the recently released SNLS3 type Ia supernovae sample, the cosmic microwave background anisotropy data from the Wilkinson Microwave Anisotropy Probe 7-yr observations, the baryon acoustic oscillation results from the Sloan Digital Sky Survey data release 7, and the latest Hubble constant measurement utilizing the Wide Field Camera 3 on the Hubble Space Telescope. The MG models considered include the Dvali-Gabadadze-Porrati(DGP) model, two $f(R)$ models, and two $f(T)$ models. We find that compared with the $\Lambda$CDM model, MG models can not lead to a appreciable reduction of the $\chi^2_{min}$. The analysis of AIC and BIC shows that the simplest cosmological constant model($\Lambda$CDM) is still most preferred by the current data, and the DGP model is strongly disfavored. In addition, from the observational constraints, we also reconstruct the evolutions of the growth factor in these models. We find that the current available growth factor data are not enough to distinguish these MG models from the $\Lambda$CDM model.

From the viewpoint of no-cloning theorem we postulate a relation between the current accelerated expansion of our universe and the inflationary expansion in the very early universe. It implies that the fate of our universe should be in a state with accelerated expansion. Quantitatively we find that the no-cloning theorem leads to a lower bound on the cosmological constant which is compatible with observations.

We explore the ultimate fate of the Universe by using a divergence-free parametrization for dark energy $w(z)=w_0+w_a({\ln (2+z)\over 1+z}-\ln2)$. Unlike the CPL parametrization, this parametrization has well behaved, bounded behavior for both high redshifts and negative redshifts, and thus can genuinely cover many theoretical dark energy models. After constraining the parameter space of this parametrization by using the current cosmological observations, we find that, at the 95.4% confidence level, our Universe can still exist at least 16.7 Gyr before it ends in a big rip. Moreover, for the phantom energy dominated Universe, we find that a gravitationally bound system will be destroyed at a time $t \simeq P\sqrt{2|1+3w(-1)|}/[6\pi |1+w(-1)|]$, where $P$ is the period of a circular orbit around this system, before the big rip.

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.

Testing the so-called consistency relations plays an important role for distinguishing the different classes of inflation models. In this paper, we investigate the possible testing based on the potential observations of the cosmic microwave background (CMB) radiation, including the planned CMBPol mission and the ideal CMB experiment where only the reduced cosmic weak lensing contamination for the B-mode polarization is considered. We find that for the canonical single-field inflation and the phantom inflation, the consistency relations are quite hard to be tested: the testing is possible only if $r>0.14$ for CMBPol mission, and $r>0.06$ for the ideal experiment. However, the situation could become much better for the general Lorentz invariant single-field inflation with large non-gaussian signal and the two-field inflation with strong correlation between the adiabatic and the isocurvature perturbations. We find that for these two classes of inflation the testing is possible if $r\gtrsim 10^{-2}$ or even smaller for both CMB experiments.

We derive the spectral indices and their runnings of local form $f_{NL}$ and $g_{NL}$ from a DBI isocurvature field and we find that the indices are suppressed by the sound speed $c_s$. This effect can be interpreted by the Lorentz boost from the viewpoint in the frame where brane is moving.

We probe the cosmic acceleration by using the recently released SNLS3 sample of 472 type Ia supernovae. Combining this type Ia supernovae dataset with the cosmic microwave background anisotropy data from the Wilkinson Microwave Anisotropy Probe 7-yr observations, the baryon acoustic oscillation results from the Sloan Digital Sky Survey data release 7, and the Hubble constant measurement from the Wide Field Camera 3 on the Hubble Space Telescope, we measure the dark energy equation of state $w$ and the deceleration parameter $q$ as functions of redshift by using the Chevallier-Polarski-Linder parametrization. Our result is consistent with a cosmological constant at 1$\sigma$ confidence level, without evidence for the recent slowing down of the cosmic acceleration. Furthermore, we consider three binned parametrizations ($w$ is piecewise constant in redshift $z$) based on different binning methods. The similar results are obtained, i.e., the $\Lambda$CDM model is still nicely compatible with current observations.

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.

It is possible that the primordial non-Gaussianity is dominated by the higher order terms, such as that set by $g_{NL}$, not $f_{NL}$. In this paper we derive the spectral index and running of $g_{NL}$ from a single isocurvature scalar field. The scale dependences of non-Gaussianity parameters are detectable if the mass of isocurvature field is not too small compared to the Hubble parameter during inflation. In addition, we also apply our results to investigate the curvaton model with near quadratic potential in detail.

We derive the evolution equation of growth factor for the matter over-dense perturbation in $f(T)$ gravity. For instance, we investigate its behavior in power law model at small redshift and compare it to the prediction of $\Lambda$CDM and dark energy with the same equation of state in the framework of Einstein general relativity. We find that the perturbation in $f(T)$ gravity grows slower than that in Einstein general relativity if $\p f/\p T>0$ due to the effectively weakened gravity.

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.

Adopting the horizon-crossing approximation, we derive the spectral index of $f_{NL}$ in general N-flation model. Axion N-flation model is taken as a typical model for generating a large $f_{NL}$ which characterizes the size of local form bispectrum. We find that its tilt $n_{f_{NL}}$ is negligibly small when all inflatons have the same potential, but a negative detectable $n_{f_{NL}}$ can be achieved in the axion N-flation with different decay constants for different inflatons. The measurement of $n_{f_{NL}}$ can be used to support or falsify the axion N-flation in the near future.

We derive the spectral index of $f_{NL}$ and its running from isocurvature single field and investigate the curvaton models with a negative spectral index of $f_{NL}$ in detail. In particular, a numerical study of the axion-type curvaton model is illustrated, and we find that the spectral index of $f_{NL}$ is negative and its absolute value is maximized around $\sigma_*=\pi f/2$ for the potential $V(\sigma)=m^2f^2(1-\cos{\sigma\over f})$. The spectral index of $f_{NL}$ can be ${\cal O}(-0.1)$ for the axion-type curvaton model. A convincing detection of a positive $n_{f_{NL}}$ will rule out the axion-type curvaton model. In addition, we also give a general discussion about the detectable parameter space for the curvaton model with a polynomial potential.

We calculate the trispectrum in ghost inflation where both the contact diagram and scale-exchange diagram are taken into account. The shape of trispectrum is discussed carefully and we find that the local form is absent in ghost inflation. In general, for the non-local shape trispectrum there are not analogous parameters to $\tau_{NL}^{loc.}$ and $g_{NL}^{loc.}$ which can completely characterize the size of local form trispectrum.

We demonstrate novel features in the behavior of the second and third order non-linearity parameters of the curvature perturbation, namely, $f_{NL}$ and $g_{NL}$, arising from non-linear motion of curvaton field. We investigate two classes of potentials for the curvaton - the first has tiny oscillations super-imposed upon the quadratic potential. The second is characterized by a single 'feature' separating two quadratic regimes with different mass scales. The feature may either be a bump or a flattening of the potential. In the case of the oscillatory potential we find that as the width and height of superimposed oscillations increase, both $f_{NL}$ and $g_{NL}$ deviate strongly from their expected values from a quadratic potential. $f_{NL}$ changes sign from positive to negative as the oscillations in the potential become more prominent. Hence, this model can be severely constrained by convincing evidence from observations that $f_{NL}$ is positive. $g_{NL}$, on the other hand, acquires very large negative values. For the the single feature potential, we find that $f_{NL}$ and $g_{NL}$ exhibit oscillatory behavior as a function of the parameter that controls the feature.

In this paper we compute the sizes of equilateral and orthogonal shape bispectrum for the general Lorentz invariant single-field inflation. The stability of field theory implies a non-negative square of sound speed which leads to a consistency relation between the sizes of orthogonal and equilateral shape bispectrum, namely $f_{NL}^{orth.}\lesssim -0.054 f_{NL}^{equil.}$. In particular, for the single-field Dirac-Born-Infeld (DBI) inflation, the consistency relation becomes $f_{NL}^{orth.}\simeq 0.070 f_{NL}^{equil.}\lesssim 0$. These consistency relations are also valid in the mixed scenario where the quantum fluctuations of some other light scalar fields contribute to a part of total curvature perturbation on the super-horizon scale and may generate a local form bispectrum. A distinguishing prediction of the mixed scenario is $\tau_{NL}^{loc.}>({6\over 5}f_{NL}^{loc.})^2$. Comparing these consistency relations to WMAP 7yr data, there is still a big room for the Lorentz invariant inflation, but DBI inflation has been disfavored at more than 68% CL.

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.

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.

In this work, we explore the cosmological consequences of the recently released Constitution sample of 397 Type Ia supernovae (SNIa). By revisiting the Chevallier-Polarski-Linder (CPL) parameterization, we find that, for fitting the Constitution set alone, the behavior of dark energy (DE) significantly deviate from the cosmological constant $\Lambda$, where the equation of state (EOS) $w$ and the energy density $\rho_{\Lambda}$ of DE will rapidly decrease along with the increase of redshift $z$. Inspired by this clue, we separate the redshifts into different bins, and discuss the models of a constant $w$ or a constant $\rho_{\Lambda}$ in each bin, respectively. It is found that for fitting the Constitution set alone, $w$ and $\rho_{\Lambda}$ will also rapidly decrease along with the increase of $z$, which is consistent with the result of CPL model. Moreover, a step function model in which $\rho_{\Lambda}$ rapidly decreases at redshift $z\sim0.331$ presents a significant improvement ($\Delta \chi^{2}=-4.361$) over the CPL parameterization, and performs better than other DE models. We also plot the error bars of DE density of this model, and find that this model deviates from the cosmological constant $\Lambda$ at 68.3% confidence level (CL); this may arise from some biasing systematic errors in the handling of SNIa data, or more interestingly from the nature of DE itself. In addition, for models with same number of redshift bins, a piecewise constant $\rho_{\Lambda}$ model always performs better than a piecewise constant $w$ model; this shows the advantage of using $\rho_{\Lambda}$, instead of $w$, to probe the variation of DE.

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 study the axion-type curvaton model, with emphasis on the large field regime where analytic results are very difficult to obtain. We evaluate the tensor-scalar ratio $r$ using WMAP normalization and the non-linearity parameters $f_{NL}$ and $g_{NL}$ by solving the equations numerically using the $\delta N$ formalism. We compare them with results for the curvaton with quadratic potential. We find that $r$ is much smaller for the axion-type case compared to the result from the quadratic potential, with the difference increasingly more pronounced at larger field values. $g_{NL}$ is found to be positive, unlike the quadratic case where it is negative, and the amplitude of $g_{NL}$ is much larger. Moreover, there is a nearly linear scaling between $g_{NL}$ and $f_{NL}$, with small deviation from linearity at large field values. The slope between $g_{NL}$ and $f_{NL}$ depends on the parameters characterizing the axion-type curvaton model. We further consider a mixed scenario where both the inflaton and the curvaton contribute to the primordial power spectrum and the non-Gaussianity parameters are found to be much larger than those in the case with quadratic potential.

In general a weakly self-interacting curvaton field is expected and the curvaton potential takes the polynomial form. The curvaton potential can be dominated by the self-interaction term during the period of inflation if the curvaton field stays at a large vacuum expectation value. We use the $\delta {\cal N}$ formalism to calculate the primordial curvature perturbation in the various possible scenarios which make the curvaton model much richer.

We argue that the typical energy density of a light scalar field should not be less than $H^4$ in the inflationary Universe. This requirement implies that the non-Gaussianity parameter $f_{NL}$ is typically bounded by the tensor-scalar ratio $r$ from above, namely $f_{NL}\lesssim 518\cdot r^{1\over 4}$. If $f_{NL}=10^2$, inflation occurred around the GUT scale.

In general there are a large number of light scalar fields in the theories going beyond standard model, such as string theory, and some of them can be taken as the candidates of curvatons. For simplicity, we assume all of curvatons have the same decay rate and suddenly decay into radiation at the same time. In order to distinguish this scenario from the more general case, we call it "N-vaton". We use $\delta {\cal N}$ formalism to calculate the primordial power spectrum and bispectrum in N-vaton model and investigate various bounds on the non-Gaussianity parameter $f_{NL}$. A red tilted primordial power spectrum and a large value of $f_{NL}$ can be naturally obtained if the curvature perturbation generated by inflaton also makes a significant contribution to the primordial power spectrum. As a realistic N-vaton model, we suppose that the axions in the KKLT compactifications of Type IIB string theory are taken as curvatons and a rich phenomenology is obtained.

We investigate the curvaton dynamics and the non-linearity parameters in curvaton model with potential slightly deviating from the quadratic form in detail. The non-linearity parameter $g_{NL}$ will show up due to the curvaton self-interaction. We also point out that the leading order of non-quadratic term in the curvaton potential can be negative, for example in the axion-type curvaton model. If a large positive $g_{NL}$ is detected, the axion-type curvaton model will be preferred.

A red tilted primordial power spectrum is preferred by WMAP five-year data and a large positive local-type non-Gaussianity $f_{NL}$ might be observed as well. In this short note we find that a red tilted and large non-Gaussian primordial power spectrum cannot be naturally obtained in curvaton model, because $f_{NL}$ is related to the initial condition of inflation.

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.

We conjecture an intrinsic UV cutoff for the validity of the effective field theory with a large number of species coupled to gravity. In four dimensions such a UV cutoff takes the form $\Lambda=\sqrt{\lambda/ N}M_p$ for $N$ scalar fields with the same potential $\lambda \phi_i^4$, $i=1,...,N$. This conjecture implies that the assisted chaotic inflation or N-flation might be in the swampland, not in the landscape. Similarly a UV cutoff $\Lambda=gM_p/\sqrt{N}$ is conjectured for the U(1) gauge theory with $N$ species.

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

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 value of scalar field coupled to gravity should be less than the Planck scale in the consistent theory of quantum gravity. It provides a theoretic constraint on the equation of state parameter for the quintessence. In some cases our theoretic constraints are more stringent than the constraints from the present experiments.

We conjecture that the inflation models with trans-Planckian excursions in the field space should be in the swampland. We check this conjecture in a few examples and investigate the constraints on the spectral index for the slow-roll inflation model in string landscape where the variation of inflaton during the period of inflation is less than the Planck scale $M_p$. A red primordial power spectrum with a lower bound on the spectral index is preferred. Both the tensor-scalar ratio and the running can be ignored.

We investigate whether the eternal chaotic inflation can be achieved when the weak gravity conjecture is taken into account. We show that even the assisted chaotic inflation with potential $\lambda\phi^4$ or $m^2\phi^2$ can not be eternal. The effective field theory description for the inflaton field breaks down before inflation reaches the eternal regime. We also find that the total number of e-folds is still bounded by the inflationary entropy for the assisted inflation.

We not only reconstruct the slow-roll parameters for fit to the running spectral index from WMAP three-year data in the usual slow-roll inflation model and noncommutative inflation model, but also investigate the evolution of these slow-roll parameters. Requiring slow-roll inflation lasts more than 20 e-folds after CMB scales leave the horizon suggests that the potential at the last stage of inflation takes the form $V(\phi)=V_0(1+{\eta_c\over 2}{(\phi-\phi_c)^2\over M_p^2})$, where $\eta_c$ is a constant. A natural mechanism to end inflation at $\phi=\phi_c$ is hybrid type inflation.

Our universe is born of a tunnelling from nothing in quantum cosmology. Nothing here can be interpreted as a state with zero entropy. As a reliable modification of the Hartle-Hawking wave function of the universe, the improved Hartle-Hawking wave function proposed by Firouzjahi, Sarangi and Tye gives many interesting observational consequences which we explore in this paper. Fruitful observations are obtained for chaotic inflation, including a detectable spatial curvature and a negligible tunnelling probability for eternal chaotic inflation. And we find that the tensor-scalar ratio and the spatial curvature for brane inflation type models should be neglected.

WMAP three-year data favors a red power spectrum at the level of 2 standard deviations, which provides a stringent constraint on the inflation models. In this note we use this data to constrain brane inflation models and find that KKLMMT model can not fit WMAP+SDSS data at the level of 1 standard deviation and a fine-tuning, eight parts in thousand at least, is needed at the level of 2 standard deviation.

Three-year data of WMAP implies not only a negative running of the spectral index with large absolute value, but also a large positive running of running of the spectral index with order of the magnitude $10^{-2}$. We calculate the running of running in usual inflation model and noncommutative inflation model. A large tensor-scalar ratio $r\geq 1.23$ is needed in order to fit the WMAP data in the noncommutative inflation model, which roughly saturates the observational upper bound on it.

A model independent analysis shows that the running of the spectral index of the three year WMAP results can be nicely realized in noncommutative inflation. We also re-examine some concrete noncommutative inflation models. We find that a large tensor-scalar ratio is required, corresponding to a low number of e-folds before the end of inflation in some simple models.

In the noncommutative space-time the fine tuning in KKLMMT model can be significantly released and a nice running of spectral index fitting the WMAP three year data can be achieved. The fitting results show that the noncommutative mass scale is roughly $5\times 10^{14}$ Gev. The string mass scale is higher than the noncommutative scale unless the string coupling is smaller than $10^{-10}$.

We investigate the non-Gaussianity of the brane inflation which happens in the same throat in the framework of the generalized KKLMMT model. When we take the constraints from non-Gaussianity into account, various consequences are discussed including the bound on the string coupling, such as the string coupling is larger than 0.08 and the effective string scale on the brane is larger than $1.3 \times 10^{-4} M_p$ in KKLMMT model.

String/M theory inspires an uncertainty relation between space and time which deviates from general relativity. It is possible to explore this deviation from cosmological observations, in particular from the CMB fluctuation spectrum. This paper extends some previous observations to more general inflation schemes, we find that the noncommutative spacetime effects always suppress the power spectrum, of both the scalar and tensor perturbations, and may provide a large enough running of the spectral index to fit the WMAP data in the inflation model.

We discuss the anthropic principle when applied to the holographic dark energy. We find that if the amplitude of the density fluctution is variable, the holographic dark energy fares better than the cosmological constant. More generally, the anthropic predictions agree better with observation for dark energy with $w_\d=p_\d/\rho_\d$ decreasing over time.

In this paper, we use the type Ia supernova data to constrain the model of holographic dark energy. For $d=1$, the best fit result is $\Omega^0_m = 0.25$, the equation of the state of the holographic dark energy $w^0_\d = -0.91$ and the transition between the decelerating expansion and the accelerating expansion happened when the cosmological red-shift was $z_T=0.72$. If we set $d$ as a free parameter, the best fit results are $d=0.21$, $\Omega_m^0=0.46$, $w_\d^0=-2.67$, which sounds like a phantom today, and the transition redshift is $z_T=0.28$.

We study the model for holographic dark energy in a spatially closed universe, generalizing the proposal in hep-th/0403127 for a flat universe. We provide independent arguments for the choice of the parameter $c=1$ in the holographic dark energy model. On the one hand, $c$ can not be less than 1, to avoid violating the second law of thermodynamics. On the other hand, observation suggests $c$ be very close to 1, it is hard to justify a small deviation of $c$ from 1, if $c>1$.

The first year results of WMAP tentatively indicate running of the spectral index as well as a deficit of power in the low multipoles in the CMB spectrum. The former can be rather easily understood in the noncommutative inflation model, and the latter, as we shall show in this paper, still appears to be an anomaly, even though the noncommutative inflation model already suppresses the low multipoles to a certain degree. By fitting the power spectrum, we determine the string scale to be $l_s\sim 4\times 10^{-29}$cm.

Very recent CMB data of WMAP offers an opportunity to test inflation models, in particular, the running of spectral index is quite new and can be used to rule out some models. We show that an noncommutative spacetime inflation model gives a good explanation of these new results. In fitting the data, we also obtain a relationship between the noncommutative parameter (string scale) and the ending time of inflation.

There has been some evidence that the fine structure "constant" $\alpha$ may vary with time. We point out that this variation can be described by a scalar field in some supergravity theory in our toy model, for instance, the N=8 extended supergravity in four dimensions which can be accommodated in M-theory.

We present a compilation of radio, infrared, optical and hard X-ray (2-10kev) data for a sample of 90 Seyfert 2 galaxies(Sy2s) with spectropolarimetric observations (41 Sy2s with detection of polarized broad lines (PBL) and 49 without PBL). Compared to Sy2s without PBL, Sy2s with PBL tend to be earlier-type spirals, and show warmer mid-infrared color and significant excess of emissions (including the hard X-ray(2-10kev), [O {\sc iii}]$\lambda$5007, infrared (25 $\mu$m) and radio). Our analyses indicate that the majority of Sy2s without PBL are those sources having less powerful AGN activities, most likely caused by low accretion rate. It implies that the detectability of the polarized broad emission lines in Sy2s may depend on their central AGN activities in most cases. Based on the available data, we find no compelling evidence for the presence of two types of Sy2s, one type of them has been proposed to be intrinsically different from Sy2s claimed in Unification Model.