We investigate a new mechanism to create large curvature perturbations on small scales due to parameter resonance in a single-field inflationary model with a small periodic structure upon the potential. After reentering the horizon, the amplified curvature perturbations can lead to observable primordial black holes as well as stochastic gravitational waves. The mass of primordial black holes and frequency of the induced gravitational waves depend on the model parameters. The resulted primordial black hole could constitute all dark matter or a fraction of dark matter in the universe, and corresponding stochastic gravitational waves fall in the frequency band measurable for the pulsar timing array and the space-based gravitational wave detectors.

The existing constraints from particle colliders reveal a suspicious but nonlethal meta-stability for our current electroweak vacuum of Higgs potential in the standard model of particle physics, which is, however, disfavored in the early Universe if the inflationary Hubble scale is larger than the instability scale when Higgs quartic self-coupling runs into negative value. Alternative to previous trials of acquiring a positive effective mass-squared from Higgs quadratic couplings to Ricci scalar or inflaton field, we propose a third approach to stabilize the Higgs potential in the early Universe by regarding Higgs as chameleon coupled to inflaton alone without conflicting to the present constraints on either Higgs or chameleon.

Strong gravitational lensing by galaxies provides us with a unique opportunity to understand the nature of gravity on galactic and extra-galactic scales. In this paper, we propose a new multimessenger approach using data from both gravitational wave (GW) and the corresponding electromagnetic (EM) counterpart to infer the constraint of the modified gravity (MG) theory denoted by the scale dependent phenomenological parameter. To demonstrate the robustness of this approach, we calculate the time-delay predictions by choosing various values of the phenomenological parameters and then compare them with that from the general relativity (GR). For the third generation ground-based GW observatory, with one typical strongly lensed GW+EM event, and assuming that the dominated error from the stellar velocity dispersions is 5\%, the GW time-delay data can distinguish an 18\% MG effect on a scale of tens of kiloparsecs with a $68\%$ confidence level. Assuming GR and a Singular Isothermal Sphere mass model, there exists a simplified consistency relationship between time-delay and imaging data. This relationship does not require for the velocity dispersion measurement, and hence can avoid major uncertainties. By using this relationship, the multimessenger approach is able to distinguish an $8\%$ MG effect. Our results show that the GW multimessenger approach can play an important role in revealing the nature of gravity on the galactic and extra-galactic scales.

The ultra-slow-roll (USR) inflationary models predict large-amplitude scalar perturbations at small scales which can lead to the primordial black hole production and scalar-induced gravitational waves. In general scalar perturbations in the USR models can only be obtained using numerical method because the usual slow-roll approximation breaks. In this work, we propose an analytical approach to estimate the scalar spectrum which is consistent with the numerical result. We find that the USR inflationary models predict a peak with power-law slopes in the scalar spectrum and energy spectrum of gravitational waves, and we derive the expression of the spectral indexes in terms of the inflationary potential. In turn, the inflationary potential near the USR regime can be reconstructed from the negative spectral index of the gravitational wave energy spectrum.

Both LISA and Taiji, planned space-based gravitational-wave detectors in orbit around the Sun, are expected to launch in 2030-2035. Assuming a one-year overlap, we explore a potential LISA-Taiji network to fast and accurately localize the gravitational-wave sources.

The recently proposed trans-Planckian censorship conjecture (TCC) imposes a strong constraint on the inflationary Hubble scale, of which the upper bound could be largely relaxed by considering a noninstantaneous reheating history. In this paper we will show that, if the primordial black holes (PBHs) are formed at reentry in the radiation-dominated era from the enhanced curvature perturbations at small scales, the TCC would impose a lower bound on the PBH mass $M_\mathrm{PBH}>\gamma(H_\mathrm{end}/10^9\,\mathrm{GeV})^2\,M_\odot$ regardless of the details for reheating history, where $\gamma$ is the collapse efficiency factor and $H_\mathrm{end}$ is the Hubble scale at the end of inflation. In particular, the current open window for PBHs to make up all the cold dark matter could be totally ruled out if the inflationary Hubble scale is larger than 10 TeV. For the case of PBHs formed in an early matter-dominated era, an upper mass bound is obtained.

We consider gravitational radiation and electromagnetic radiation from point mass binary with electric charges in a Keplerian orbit, and calculate the merger rate distribution of primordial black hole binaries with charges and a general mass function by taking into account gravitational torque and electromagnetic torque by the nearest primordial black hole. We apply the formalism to the extremal charged case and find that $\alpha=-(m_i+m_j)^2\partial^2 \ln {\cal R}(m_i,m_j)/\partial m_i \partial m_j=12/11$, which is independent of the mass function.

We propose a refined version of trans-Planckian censorship conjecture (TCC), which could be elaborated from the strong scalar weak gravity conjecture combined with some entropy bounds. In particular, no fine-tuning on the inflation model-building is required in the refined TCC, and it automatically passes the tests from those stringy examples that support the original TCC. Furthermore, our refined TCC could be consistent with hilltop eternal inflation.

If the black holes detected by LIGO/VIRGO are primordial black holes (PBHs) sourcing from a large primordial curvature perturbation on small scales, the corresponding induced gravitational waves (GWs) would peak at nanohertz that is detectable by the current and future observations of pulsar timing array (PTA). In this paper we show that with the mass function estimated from the merger rate of LIGO O1 and O2 events, the induced GWs from such a curvature perturbation with a Gaussian narrow peak at some small scale would be in a seemingly mild tension with current constraints from PTA. However, if the curvature perturbation is of local-type non-Gaussianity with a non-linear parameter $f_\text{NL}\gtrsim\mathcal{O}(10)$, the tension could be relieved. Nevertheless, such an induced GWs must be detectable by the Square Kilometer Array in a decade or less.

A space-based gravitational-wave detector, LISA, consists of a triangle of three spacecrafts with a separation distance of 2.5 million kilometers in a heliocentric orbit behind the Earth. Like LISA, Taiji is compose of a triangle of three spacecrafts with a separation distance of 3 million kilometers in a heliocentric orbit ahead of the Earth. They are expected to launch in 2030-2035. Assuming a one-year overlap, we propose the LISA-Taiji network in space to fast and accurately localize the gravitational-wave sources. We use the Fisher information matrix approach to analyze the sky localization for coalescing massive black hole binaries. For an equal-mass black hole binary located at redshift of 1 with a total intrinsic mass of $10^5 M_{\odot}$, the LISA-Taiji network may achieves about four orders of magnitude improvement on the event localization region compared to an individual detector. The precision measurement of sky location from the gravitational-wave signal may completely identify the host galaxy with low redshifts prior to the final black hole merger. Such the identification of the host galaxy is helpful for the follow-up change in electromagnetic emissions of the accretion disk when the massive black hole binary merges to a single massive black hole, and enables the coalescing massive black hole binaries to be used as a standard siren.

We study the general infrared behavior of the power spectrum of a stochastic gravitational wave background produced by stress tensor in the form bilinear in certain dynamical degrees of freedom. We find $\Omega_{\text{GW}}\propto k^3$ for a very wide class of the sources which satisfy a set of reasonable conditions. Namely, the $k^3$ scaling is universally valid when the source term is bounded in both frequency and time, is effective in a radiation-dominated stage, and for $k$ smaller than all the physical scales associated with the source, like the peak frequency, peak width, and time duration, etc.. We also discuss possible violations of these conditions and their physical implications.

We develop a formalism to calculate the merger rate density of primordial black hole binaries with a general mass function, by taking into account the merger history of primordial black holes. We apply the formalism to three specific mass functions, monochromatic, power-law and log-normal cases. In the former case, the merger rate density is dominated by the single-merger events, while in the latter two cases, the contribution of the multiple-merger events on the merger rate density can not be ignored. The effects of the merger history on the merger rate density depend on the mass function.

Recently it was shown that cosmic anisotropy can be well tested using either standard siren measurement of luminosity distance $d_\mathrm{L}(z)$ from gravitational-wave (GW) observation or dispersion measure ($\mathrm{DM}(z)$) from fast radio burst (FRB). It was also observed that the combined measurement of $d_\mathrm{L}(z)\cdot\mathrm{DM}(z)$ from the GW/FRB association system as suggested in some of FRB models is more effective to constrain cosmological parameters than $d_\mathrm{L}(z)$ or $\mathrm{DM}(z)$ separately due to its independence from Hubble constant. In this paper, we will show both theoretically and with simulation that, this upgraded sirens from combined GW/FRB observations could test cosmic anisotropy with a double relative sensitivity compared to the usual standard siren from GW observation alone.

We investigate the formation of primordial black holes (PBHs) from the collapse of spherically symmetric domain wall bubbles, which spontaneously nucleate via quantum tunneling during inflation. Since the tension of domain walls changes with time and so domain walls nucleate in a short time interval, the mass function of PBHs in general has a spike-like structure. In contrast to models in which PBHs produced from overdense regions, our model avoids the uncertainties of PBHs production mechanism. PBHs from domain walls with mass around $10^{20}\mathrm{g}$ may constitute all dark matter, those with mass around $10^{34}\mathrm{g}$ can explain the merger events of binary black holes detected by LIGO.

We study gravitational waves (GWs) induced by non-Gaussian curvature perturbations. We calculate the density parameter per logarithmic frequency interval, $\Omega_\text{GW}(k)$, given that the power spectrum of the curvature perturbation $\mathcal{P}_\mathcal{R}(k)$ has a narrow peak at some small scale $k_*$, with a local-type non-Gaussianity, and constrain the nonlinear parameter $f_\text{NL}$ with the future LISA sensitivity curve as well as with constraints from the abundance of the primordial black holes (PBHs). We find that the non-Gaussian contribution to $\Omega_\text{GW}$ increases as $k^3$, peaks at $k/k_*=4/\sqrt{3}$, and has a sharp cutoff at $k=4k_*$. The non-Gaussian part can exceed the Gaussian part if $\mathcal{P}_\mathcal{R}(k)f_\text{NL}^2\gtrsim1$. If both a slope $\Omega_\text{GW}(k)\propto k^\beta$ with $\beta\sim3$ and the multiple-peak structure around a cutoff are observed, it can be recognized as a smoking gun of the primordial non-Gaussianity. We also find that if PBHs with masses of $10^{20}\text{g}$ to $10^{22}\text{g}$ are identified as cold dark matter of the Universe, the corresponding GWs must be detectable by LISA-like detectors, irrespective of the value of $\mathcal{P}_\mathcal{R}$ or $f_\text{NL}$.

We investigate the effect of the cuspiness of scalar potentials on the production of gravitational waves during oscillon formation after inflation. We consider a more general form of potentials with a mass parameter $M$, which repoduce cuspy potentials for fields much larger than $M$, and smooth potentials in the opposite limit. For cuspy potentials, nonsmooth oscillations of the inflaton induce an amplification of the inflaton fluctuations at the bottom of the potential, so that oscillons copiously form, which leads to a significant stochastic gravitational wave background with a double-peak spectrum. By varying the parameter $M$, we find that cuspy potentials yield stronger signals of gravitational waves and the generation of gravitational waves disappears for smooth potentials. Moreover, we calculate the equation of state after inflation and find the presence of a quasi matter-dominated stage right before the transition to the radiation-dominated stage.

We identify analytically a multiple-peak structure in the energy-density spectrum of induced gravitational waves (GWs) generated at second-order from a primordial scalar perturbations also with multiple($n$) peaks at small scales $k_{*i}$. The energy-density spectrum of induced GWs exhibits at most $C_{n+1}^2$ and at least $n$ peaks at wave-vectors $k_{ij}\equiv(k_{*i}+k_{*j})/\sqrt{3}$ due to resonant amplification, and, under the narrow-width approximation, it contains an universal factor that can be interpreted as a result of momentum conservation. We also extend these discussions to the case of non-Gaussian perturbations.

We develop an analytic formalism for computing the merger rate of primordial black hole binaries with a general mass function by taking into account the torques by the surrounding primordial black holes and linear density perturbations. We find that $\alpha=-(m_{i}+m_{j})^2\partial^{2} \ln {\mathcal R}(m_{i},m_{j})/\partial m_{i} \partial m_{j}=36/37$ is independent of the mass function. Moreover, the ratio of the merger rate density of primordial black hole binaries by taking into account the torques by the surrounding primordial black holes to by the nearest primordial black hole is independent of the masses of binaries.

We investigate the constraints on the reionization history of the Universe from a joint analysis of the cosmic microwave background and neutral hydrogen fraction data. The $\tanh$ parametrization and principal component analysis methods are applied to the reionization history respectively. The commonly used $\tanh$ parametrization is oversimplistic when the neutral hydrogen fraction data are taken into account. Using the principal component analysis method, the reconstructed reionization history is consistent with the neutral hydrogen fraction data. With the principal component analysis method, we reconstruct the neutral hydrogen fraction at $z=9.75$ as $x_{\text{HI}}=0.69^{+0.30}_{-0.32}$ for $6<z<20$ range reconstruction, and $x_{\text{HI}}=0.76^{+0.22}_{-0.27}$ for $6<z<30$ range reconstruction. These results suggest that the Universe began to reionize at redshift no later than $z=10$ at a $95\%$ confidence level.

We study a model of the emergent dark universe, which lives on the time-like hypersurface in a five-dimensional bulk spacetime. The holographic fluid on the hypersurface is assumed to play the role of the dark sector, mainly including the dark energy and apparent dark matter. Based on the modified Friedmann equations, we present a Markov-Chain-Monte-Carlo analysis with the observational data, including type Ia Supernova and the direct measurement of the Hubble constant. We obtain a good fitting result and the matter component turns out to be small enough, which matches well with our theoretical assumption that only the normal matter is required. After considering the fitting parameters, an effective potential of the model with a dynamical scalar field is reconstructed. The parameters in the swampland criteria are extracted, and they satisfy the criteria at the present epoch but are in tension with the criteria if the potential is extended to the future direction. The method to reconstruct the potential is helpful to study the swampland criteria of other models without an explicit scalar field.

We discuss a scenario that the dark matter in late time universe emerges as part of the holographic stress-energy tensor on the hypersurface in higher dimensional flat spacetime. Firstly we construct a toy model with a de Sitter hypersurface as the holographic screen in the flat bulk. After adding the baryonic matter on the screen, we assume that both of the dark matter and dark energy can be described by the Brown-York stress-energy tensor. From the Hamiltonian constraint equation in the flat bulk, we find an interesting relation between the dark matter and baryonic matter's energy density parameters, by comparing with the Lambda cold dark matter parameterization. We further compare this holographic embedding of emergent dark matter with traditional braneworld scenario and present an alternative interpretation as the holographic universe. It can be reduced to our toy constraint in the late time universe, with the new parameterization of the Friedmann equation. We also comment on the possible connection with Verlinde's emergent gravity, where the dark matter is treated as the elastic response of the baryonic matter on the de Sitter spacetime background. We show that from the holographic de Sitter model with elasticity, the Tully-Fisher relation and the dark matter distribution in the galaxy scale can be derived.

We investigate the sensitivity of reheating history on the abundance of primordial black holes (PBHs). Contrary to the monochromatic case of mass fraction of PBH, reheating era with different e-folding and equation-of-state can have a substantial impact on the abundance of PBH with an extended mass fraction. We demonstrate explicitly this reheating sensitivity in an illustrative model of single field inflation with a quasi-inflection point, and find that both the peak position and amplitude of the extended mass fraction as well as the abundance of PBH in DM can vary by many orders of magnitude, which adds another layer of uncertainty on the PBH scenarios as dark matter.

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.

The gravitational wave (GW) as a standard siren directly determines the luminosity distance from the gravitational waveform without reference to the specific cosmological model, of which the redshift can be obtained separately by means of the electromagnetic counterpart like GW events from binary neutron stars and massive black hole binaries (MBHBs). To see to what extent the standard siren can reproduce the presumed dipole anisotropy written in the simulated data of standard siren events from typical configurations of GW detectors, we find that (1) for the Laser Interferometer Space Antenna with different MBHB models during five-year observations, the cosmic isotropy can be ruled out at $3\sigma$ confidence level (C.L.) and the dipole direction can be constrained roughly around $20\%$ at $2\sigma$ C.L., as long as the dipole amplitude is larger than $0.03$, $0.06$ and $0.025$ for MBHB models Q3d, pop III and Q3nod with increasing constraining ability, respectively; (2) for Einstein Telescope with no less than $200$ standard siren events, the cosmic isotropy can be ruled out at $3\sigma$ C.L. if the dipole amplitude is larger than $0.06$, and the dipole direction can be constrained within $20\%$ at $3\sigma$ C.L. if the dipole amplitude is near $0.1$; (3) for the Deci-Hertz Interferometer Gravitational wave Observatory with no less than $100$ standard siren events, the cosmic isotropy can be ruled out at $3\sigma$ C.L. for dipole amplitude larger than $0.03$ , and the dipole direction can even be constrained within $10\%$ at $3\sigma$ C.L. if the dipole amplitude is larger than $0.07$. Our work manifests the promising perspective of the constraint ability on the cosmic anisotropy from the standard siren approach.

We study the hydrodynamics of bubble expansion in cosmological first-order phase transition in the Friedmann-Lema\^{\i}tre-Robertson-Walker (FLRW) background with probe limit. Different from previous studies for fast first-order phase transition in flat background, we find that, for slow first-order phase transition in FLRW background with a given peculiar velocity of the bubble wall, the efficiency factor of energy transfer into bulk motion of thermal fluid is significantly reduced, thus decreasing the previously-thought dominated contribution from sound wave to the stochastic gravitational-wave background.

In recent years, superfluid dark matter (SfDM) has become a competitive model of emergent modified Newtonian dynamics (MOND) scenario: MOND phenomenons naturally emerge as a derived concept due to an extra force mediated between baryons by phonons as a result of axionlike particles condensed as superfluid at galactic scales; Beyond galactic scales, these axionlike particles behave as normal fluid without phonon-mediated MOND-like force between baryons, therefore SfDM also maintains the usual success of $\Lambda$CDM at cosmological scales. In this paper, we use gravitational waves (GWs) to probe the relevant parameter space of SfDM. GWs through Bose-Einstein condensate (BEC) could propagate with a speed slightly deviation from the speed-of-light due to the change in the effective refractive index, which depends on the SfDM parameters and GW-source properties. We find that Five hundred meter Aperture Spherical Telescope (FAST), Square Kilometre Array (SKA) and International Pulsar Timing Array (IPTA) are the most promising means as GW probe of relevant parameter space of SfDM. Future space-based GW detectors are also capable of probing SfDM if a multimessenger approach is adopted.

We study the production of gravitational waves during oscillations of the inflaton around the minimum of a cuspy potential after inflation. We find that a cusp in the potential can trigger copious oscillon formation, which sources a characteristic energy spectrum of gravitational waves with double peaks. The discovery of such a double-peak spectrum could test the underlying inflationary physics.

We present a fast numerical screened halo model algorithm (CHAM) for modeling non-linear power spectrum for the alternative models to LCDM. This method has three obvious advantages. First of all, it is not being restricted to a specific dark energy/modified gravity model. In principle, all of the screened scalar-tensor theories can be applied. Second, the least assumptions are made in the calculation. Hence, the physical picture is very easily understandable. Third, it is very predictable and does not rely on the calibration from N-body simulation. As an example, we show the case of Hu-Sawicki f(R) gravity. In this case, the typical CPU time with the current parallel Python script (8 threads) is roughly within $10$ minutes. The resulting spectra are in a good agreement with N-body data within a few percentage accuracy up to k~1 h/Mpc.

We investigate the impact of different redshift distributions of random samples on the baryon acoustic oscillations (BAO) measurements of $D_V(z)r_\mathrm{d}^\mathrm{fid}/r_\mathrm{d}$ from the two-point correlation functions (2PCF) of galaxies in the Data Release 12 (DR12) of the Baryon Oscillation Spectroscopic Survey (BOSS). Big surveys, such as BOSS, usually assign redshifts to the random samples by randomly drawing values from the measured redshift distributions of the data, which would necessarily introduce fiducial signals of fluctuations into the random samples, weakening the signals of BAO, if the cosmic variance cannot be ignored. We propose a smooth function of redshift distribution that fits the data well to populate the random galaxy samples. The resulting cosmological parameters match the input parameters of the mock catalogue very well. The significance of BAO signals has been improved by $0.33\sigma$ for a low-redshift (LOWZ) sample and by $0.03\sigma$ for a constant-stellar-mass (CMASS) sample, though the absolute values do not change significantly. Given the precision of the measurements of current cosmological parameters, it would be appreciated for the future improvements on the measurements of galaxy clustering.

The gravitational waves from compact binary systems are viewed as a standard siren to probe the evolution of the universe. This paper summarizes the potential and ability to use the gravitational waves to constrain the cosmological parameters and the dark sector interaction in the Gaussian process methodology. After briefly introducing the method to reconstruct the dark sector interaction by the Gaussian process, the concept of standard sirens and the analysis of reconstructing the dark sector interaction with LISA are outlined. Furthermore, we estimate the constraint ability of the gravitational waves on cosmological parameters with ET. The numerical methods we use are Gaussian process and the Markov-Chain Monte-Carlo. Finally, we also forecast the improvements of the abilities to constrain the cosmological parameters with ET and LISA combined with the \it Planck.

We investigate in details the gravitational wave (GW) from the first-order phase transition (PT) in the extended standard model of particle physics with a dimension-six operator, which is capable of exhibiting the recently discovered slow first-order PT in addition to the usually studied fast first-order PT. To simplify the discussion, it is sufficient to work with an example of a toy model with the sextic term, and we propose an unified description for both slow and fast first-order PTs. We next study the full one-loop effective potential of the model with fixed/running renormalization-group (RG) scales. Compared to the prediction of GW energy density spectrum from the fixed RG scale, we find that the presence of running RG scale could amplify the peak amplitude by amount of one order of magnitude while shift the peak frequency to the lower frequency regime, and the promising regime of detection within the sensitivity ranges of various space-based GW detectors shrinks down to a lower cut-off value of the sextic term rather than the previous expectation.

We investigate constraints on Lorentz invariance violation in the neutrino sector from a joint analysis of big bang nucleosynthesis and the cosmic microwave background. The effect of Lorentz invariance violation during the epoch of big bang nucleosynthesis changes the predicted helium-4 abundance, which influences the power spectrum of the cosmic microwave background at the recombination epoch. In combination with the latest measurement of the primordial helium-4 abundance, the Planck 2015 data of the cosmic microwave background anisotropies give a strong constraint on the deformation parameter since adding the primordial helium measurement breaks the degeneracy between the deformation parameter and the physical dark matter density.

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

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

We perform a forecast analysis of the ability of the LISA space-based interferometer to reconstruct the dark sector interaction using gravitational wave (GW) standard sirens at high redshift. We employ Gaussian process methods to reconstruct the distance-redshift relation in a model independent way. We adopt simulated catalogues of standard sirens given by merging massive black hole binaries (MBHBs) visible by LISA, with an electromagnetic (EM) counterpart detectable by future telescopes. The catalogues are constructed considering three different astrophysical scenarios for the evolution of MBHB mergers based on the semi-analytic model of E. Barausse (2012). We first use these standard siren datasets to assess the potential of LISA in reconstructing a possible interaction between vacuum dark energy and dark matter. Then we combine the LISA cosmological data with supernovae data simulated for the Dark Energy Survey (DES). We consider two scenarios distinguished by the time duration of the LISA mission: 5 and 10 years. Using only LISA standard siren data, the dark sector interaction can be well reconstructed from redshift $z\sim1$ to $z\sim3$ (5 yr) and $z\sim1$ up to $z\sim5$ (10 yr), though the reconstruction is inefficient at lower redshift. When combined with the DES datasets, the interaction is well reconstructed in the whole redshift region from $z\sim0$ to $z\sim3$ (5 yr) and $z\sim0$ to $z\sim5$ (10 yr). MBHB standard sirens can thus be used to constrain the dark sector interaction at redshift ranges not reachable by usual supernovae datasets which probe only the $z\lesssim 1.5$ range. GW standard sirens will not only constitute a complementary and alternative way, with respect to familiar EM observations, to probe the cosmic expansion, but will also provide new tests to constrain possible deviations from the standard $\Lambda$CDM dynamics, especially at high redshift.

We investigate the constraint ability of the gravitational wave (GW) as the standard siren on the cosmological parameters by using the third-generation gravitational wave detector: the Einstein Telescope. We simulate the luminosity distances and redshift measurements from 100 to 1000 GW events. We use two different algorithms to constrain the cosmological parameters. For the Hubble constant $H_0$ and dark matter density parameter $\Omega_m$, we adopt the Markov chain Monte Carlo approach. We find that with about 500-600 GW events we can constrain the Hubble constant with an accuracy comparable to \textit{Planck} temperature data and \textit{Planck} lensing combined results, while for the dark matter density, GWs alone seem not able to provide the constraints as good as for the Hubble constant; the sensitivity of 1000 GW events is a little lower than that of \textit{Planck} data. It should require more than 1000 events to match the \textit{Planck} sensitivity. Yet, for analyzing the more complex dynamical property of dark energy, i.e., the equation of state $w$, we adopt a new powerful nonparametric method: the Gaussian process. We can reconstruct $w$ directly from the observational luminosity distance at every redshift. In the low redshift region, we find that about 700 GW events can give the constraints of $w(z)$ comparable to the constraints of a constant $w$ by \textit{Planck} data with type Ia supernovae. Those results show that GWs as the standard sirens to probe the cosmological parameters can provide an independent and complementary alternative to current experiments.

We develop a new model-independent method to probe the constancy of the speed of light $c$. In our method, the degeneracy between the cosmic curvature and the speed of light can be eliminated, which makes the test more natural and general. Combining the independent observations of Hubble parameter $H(z)$ and luminosity distance $d_L(z)$, we use the model-independent smoothing technique, Gaussian processes, to reconstruct them and then detect variation of the speed of light. We find no signal of deviation from the present value of the speed of light $c_0$. Moreover, to demonstrate the improvement in probing the constancy of the speed of light from future experiments, we produce a series of simulated data. The Dark Energy Survey will be able to detect $\Delta c /c_0 \sim 1\%$ at $\sim 1.5\sigma$ confidence level and $\Delta c /c_0 \sim 2\%$ at $\sim 3\sigma$ confidence level. If the errors are reduced to one-tenth of the expected DES ones, it can detect a $\Delta c /c_0 \sim 0.1\%$ variation at $\sim 2\sigma$ confidence level.

We introduce a model-independent approach to the null test of the cosmic curvature which is geometrically related to the Hubble parameter $H(z)$ and luminosity distance $d_L(z)$. Combining the independent observations of $H(z)$ and $d_L(z)$, we use the model-independent smoothing technique, Gaussian processes, to reconstruct them and determine the cosmic curvature $\Omega_K^{(0)}$ in the null test relation. The null test is totally geometrical and without assuming any cosmological model. We show that the cosmic curvature $\Omega_K^{(0)}=0$ is consistent with current observational data sets, falling within the $1\sigma$ limit. To demonstrate the effect on the precision of the null test, we produce a series of simulated data of the models with different $\Omega_K^{(0)}$. Future observations in better quality can provide a greater improvement to constrain or refute the flat universe with $\Omega_K^{(0)}=0$.

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

It was shown recently that, without jeopardizing the success of the $\Lambda$ cold dark matter model on cosmic scales, the modified Newtonian dynamics (MOND) can be derived as an emergent phenomenon when axionlike dark matter particles condense into superfluid on the galactic scales. We propose in this paper a Dirac-Born-Infeld (DBI) scalar field conformally coupled to the matter components. To maintain the success of MOND phenomenon of dark matter superfluid on the galactic scales, the fifth force introduced by the DBI scalar should be screened on the galactic scales. It turns out that the screening effect naturally leads to a simple explanation for a longstanding puzzle that the MOND critical acceleration coincides with present Hubble scale. This galactic coincidence problem is solved, provided that the screened DBI scalar also plays the role of dark energy on the cosmic scales.

The measured masses of the Higgs boson and top quark indicate that the effective potential of the standard model either develops an unstable electroweak vacuum or stands stable all the way up to the Planck scale. In the latter case in which the top quark mass is about $2\sigma$ below its present central value, the Higgs boson can be the inflaton with the help of a large nonminimal coupling to curvature in four dimensions. We propose a scenario in which the Higgs boson can be the inflaton in a five-dimensional Gauss-Bonnet braneworld model to solve both the unitarity and stability problems which usually plague Higgs inflation. We find that in order for Higgs inflation to happen successfully in the Gauss-Bonnet regime, the extra dimension scale must appear roughly in the range between the TeV scale and the instability scale of standard model. At the tree level, our model can give rise to a naturally small nonminimal coupling $\xi\sim\mathcal{O}(1)$ for the Higgs quartic coupling $\lambda\sim\mathcal{O}(0.1)$ if the extra dimension scale lies at the TeV scale. At the loop level, the inflationary predictions at the tree level are preserved. Our model can be confronted with future experiments and observations from both particle physics and cosmology.

We investigate the influence on the inflationary predictions from the reheating processes characterized by the $e$-folding number $N_{\mathrm{reh}}$ and the effective equation-of-state parameter $w_{\mathrm{reh}}$ during the reheating phase. For the first time, reheating processes can be constrained in the $N_{\mathrm{reh}}\!-\!w_{\mathrm{reh}}$ plane from Planck 2015. We find that for Higgs inflation with a nonminimal coupling to gravity, the predictions are insensitive to the reheating phase for current CMB measurements. We also find that the spontaneously broken SUSY inflation and axion monodromy inflation with $\phi^{2/3}$ potential, which with instantaneous reheating lie outside or at the edge of the $95\%$ confidence region in the $n_s\!-\!r$ plane from Planck 2015 TT,TE,EE$+$lowP, can well fit the data with the help of reheating processes. Future CMB experiments would put strong constraints on reheating processes.

The simple assumption of an instantaneous reionization of the Universe may bias estimates of cosmological parameters. In this paper a model-independent principal component method for the reionization history is applied to give constraints on the cosmological parameters from recent Planck 2015 data. We find that the Universe are not completely reionized at redshifts $z \ge 8.5$ at 95% CL. Both the reionization optical depth and the matter fluctuation amplitude are higher than but consistent with those obtained in the standard instantaneous reionization scheme. The high estimated value of the matter fluctuation amplitude strengthens the tension between Planck CMB observations and some astrophysical data, such as cluster counts and weak lensing. The tension can significantly be relieved if the neutrino masses are allowed to vary. Thanks to a high scalar spectral index, the low-scale spontaneously broken SUSY inflationary model can fit the data well, which is marginally disfavored at 95% CL in the Planck analysis.

We obtain the reduced CMB data $\{l_A, R, z_*\}$ from WMAP9, WMAP9+BICEP2, Planck+WP and Planck+WP+BICEP2 for the $\Lambda$CDM and $w$CDM models with or without spatial curvature. We then use these reduced CMB data in combination with low-redshift observations to put constraints on cosmological parameters. We find that including BICEP2 results in a higher value of the Hubble constant especially when the equation of state of dark energy and curvature are allowed to vary. For the $\Lambda$CDM model with curvature, the estimate of the Hubble constant with Planck+WP+Lensing is inconsistent with the one derived from Planck+WP+BICEP at about 1.3 $\sigma$ confidence level.

We present a nonparametric approach to reconstruct the interaction between dark energy and dark matter directly from SNIa Union 2.1 data using Gaussian processes, which is a fully Bayesian approach for smoothing data. In this method, once the equation of state ($w$) of dark energy is specified, the interaction can be reconstructed as a function of redshift. For the decaying vacuum energy case with $w=-1$, the reconstructed interaction is consistent with the standard $\Lambda$CDM model, namely, there is no evidence for the interaction. This also holds for the constant $w$ cases from $-0.9$ to $-1.1$ and for the Chevallier-Polarski-Linder (CPL) parametrization case. If the equation of state deviates obviously from $-1$, the reconstructed interaction exists at $95\%$ confidence level. This shows the degeneracy between the interaction and the equation of state of dark energy when they get constraints from the observational data.

Constrains of dark energy (DE) at high redshift from current and mock future observational data are obtained. It is found that present data give poor constraints of DE even beyond redshift z=0.4, and mock future 2298 type Ia supernove data only give a little improvement of the constraints. We analyze in detail why constraints of DE decrease rapidly with the increasing of redshift. Then we try to improve the constraints of DE at high redshift. It is shown that the most efficient way is to improve the error of observations.

Recently released Planck data favor a lower value of the Hubble constant and a higher value of the fraction matter density in the standard $\Lambda$CDM model, which are discrepant with some of the low-redshift measurements. Within the context of this cosmology, we examine the consistency of the estimated values for the Hubble constant and fraction matter density with redshift tomography. Using the SNe Ia, Hubble parameter, BAO and CMB data, which are divided into three bins, we find no statistical evidence for any tension in the three redshift bins, although there exists a 1.4$\sigma$ deviation of the Hubble constant in the middle redshift from the one in the high redshift bin.

By using the cubic spline interpolation method, we reconstruct the shape of the primordial scalar and tensor power spectra from the recently released {\it Planck} temperature and BICEP2 polarization cosmic microwave background data. We find that the vanishing scalar index running ($\dd n_s/\dd\ln k$) model is strongly disfavored at more than $3\sigma$ confidence level on the $k=0.0002$ Mpc$^{-1}$ scale. Furthermore, the power-law parameterization gives a blue-tilt tensor spectrum, no matter using only the first 5 bandpowers $n_t = 1.20^{+0.56}_{-0.64} (95% {\rm CL})$ or the full 9 bandpowers $n_t = 1.24^{+0.51}_{-0.58} (95% {\rm CL})$ of BICEP2 data sets. Unlike the large tensor-to-scalar ratio value ($r\sim0.20$) under the scale-invariant tensor spectrum assumption, our interpolation approach gives $r_{0.002} < 0.060 (95% {\rm CL})$ by using the first 5 bandpowers of BICEP2 data. After comparing the results with/without BICEP2 data, we find that {\it Planck} temperature with small tensor amplitude signals and BICEP2 polarization data with large tensor amplitude signals dominate the tensor spectrum reconstruction on the large and small scales, respectively. Hence, the resulting blue tensor tilt actually reflects the tension between {\it Planck} and BICEP2 data.

We update the cosmological parameter estimation for three non-vanilla models by a joint analysis of \CCCP\ X-ray cluster, the newly released \Planck\ CMB data as well as some external data sets, such as baryon acoustic oscillation measurements from the 6dFGS, SDSS DR7 and BOSS DR9 surveys, and Hubble Space Telescope $H_0$ measurement. First of all, we find that X-ray cluster data sets strongly favor a non-zero summed neutrino mass at more than 3$\sigma$ confidence level in these non-vanilla models. And then, we reveal some tensions between X-ray cluster and {\it Planck} data in some cosmological parameters. For the matter power spectrum amplitude $\sigma_8$, X-ray cluster data favor a lower value compared with {\it Planck}. Because of the strong $\sigma_8-\sum m_{\nu}$ degeneracy, this tension could beyond 2$\sigma$ confidence level when the summed neutrino mass $\sum m_{\nu}$ is allowed to vary. For the CMB lensing amplitude $A_L$, the addition of X-ray cluster data results in a 3$\sigma$ deviation from the vanilla model. Furthermore, {\it Planck}+X-ray data prefer a large Hubble constant and phantom-like dark energy equation of state, which are in $2\sigma$ tension with those from WMAP7+X-ray data. Finally, we find that these tensions/descrepencies could be relaxed in some sense by adding a $9\%$ systematic shift in the cluster mass functions.

We study the possibly existing anisotropy in the accelerating expansion universe with the Union2 Type Ia supernovae data and Gamma-ray burst data. We construct a direction-dependent dark energy model and constrain the anisotropy direction and strength of modulation. We find that the maximum anisotropic deviation direction is $(l,\,b)=(126^{\circ},\,13^{\circ})$ (or equivalently $(l,\,b)=(306^{\circ},\,-13^{\circ})$), and the current anisotropy level is $g_0=0.030_{+0.010}^{-0.030}$ ($1\sigma$ confidence level with Union2 data). Our results do not show strong evidence for the anisotropic dark energy model. We also discuss potential methods that may distinguish the peculiar velocity field from the anisotropic dark energy model.

The difference between vacuum energy of quantum fields in Minkowski space and in Friedmann-Robterson-Walker universe might be related to the observed dark energy. The vacuum energy of the Veneziano ghost field introduced to solve the $U(1)_A$ problem in QCD is of the form, $ H+ {\cal O}(H^2)$. Based on this, we study the dynamical evolution of a phenomenological dark energy model whose energy density is of the form $\alpha H+\beta H^2$. In this model, the universe approaches to a de Sitter phase at late times. We fit the model with current observational data including SnIa, BAO, CMB, BBN, Hubble parameter and growth rate of matter perturbation. It shows that the universe begins to accelerate at redshift $z\sim 0.75$ and this model is consistent with current data. In particular, this model fits the data of growth factor well as the $\Lambda CDM$ model.

The deceleration parameter q as the diagnostic of the cosmological accelerating expansion is investigated. By expanding the luminosity distance to the fourth order of redshift and the so-called y-redshift in two redshift bins and fitting the SNIa data (Union2), the marginalized likelihood distribution of the current deceleration parameter shows that the cosmic acceleration is still increasing, but there might be a tendency that the cosmic acceleration will slow down in the near future. We also fit the Hubble evolution data together with SNIa data by expanding the Hubble parameter to the third order, showing that the present decelerating expansion is excluded within $2\sigma$ error. Further exploration on this problem is also approached in a non-parametrization method by directly reconstructing the deceleration parameter from the distance modulus of SNIa, which depends neither on the validity of general relativity nor on the content of the universe or any assumption regarding cosmological parameters. More accurate observation datasets and more effective methods are still in need to make a clear answer on whether the cosmic acceleration will keep increasing or not.

In this paper we study the possibly existing anisotropy in the acceleration expansion by use of the full sample of Union2 data. Using the hemisphere comparison method to search for a preferred direction, we take the deceleration parameter q0 as the diagnostic to quantify the anisotropy level in the wCDM model. We find that the maximum accelerating expansion direction is $(l,b)=(314_{+20^{\circ}}^{\circ-13^{\circ}},28_{+11^{\circ}}^{\circ-33^{\circ}})$,with the maximum anisotropy level of $\delta q_{0,max}/\bar{q}_0=0.79_{+0.27}^{-0.28}$, and that the anisotropy is more prominent when only low redshift data ($z\leq0.2$) are used. We also discuss this issue in the $CPL$ parameterized model, showing a similar result.

In this paper we study the non-Gaussian features of the primordial fluctuations in loop quantum cosmology with the inverse volume corrections. The detailed analysis is performed in the single field slow-roll inflationary models. However, our results reflect the universal characteristics of bispectrum in loop quantum cosmology. The main corrections to the scalar bispectrum come from two aspects: one is the modifications to the standard Bunch-Davies vacuum, the other is the corrections to the background dependent variables, such as slow-roll parameters. Our calculations show that the loop quantum corrections make $f_{{\rm NL}}$ of the inflationary models increase 0.1%. Moreover, we find that two new shapes arise, namely $\mathcal F_{1}$ and $\mathcal F_{2}$. The former gives a unique loop quantum feature which is less correlated with the local, equilateral and single types, while the latter is highly correlated with the local one.

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.

The model-independent piecewise parametrizations (0-spline, linear-spline and cubic-spline) are used to estimate constraints of equation of state of dark energy ($w_{de}$) from current observational data (including SNIa, BAO and Hubble parameter) and the simulated future data. A combination of fitting results of $w_{de}$ from these three spline methods reveal essential properties of real equation of state $w_{de}$. It is shown that $w_{de}$ beyond redshift $z\sim0.5$ is poorly constrained from current data, and the mock future $\sim2300$ supernovae data give poor constraints of $w_{de}$ beyond $z\sim1$. The fitting results also indicate that there might exist a rapid transition of $w_{de}$ around $z\sim0.5$. The difference between three spline methods in reconstructing and constraining $w_{de}$ has also been discussed.

To constrain cosmological parameters, one often makes a joint analysis with different combinations of observational data sets. In this paper we take the figure of merit (FoM) for Dark Energy Task Force fiducial model (CPL model) to estimate goodness of different combinations of data sets, which include 11 widely-used observational data sets (Type Ia Supernovae, Observational Hubble Parameter, Baryon Acoustic Oscillation, Cosmic Microwave Background, X-ray Cluster Baryon Mass Fraction, and Gamma-Ray Bursts). We analyze different combinations and make a comparison for two types of combination based on two types of basic combinations, which are often adopted in the literatures. We find two sets of combinations, which have strong ability to constrain the dark energy parameters, one has the largest FoM, the other contains less observational data with a relative large FoM and a simple fitting procedure.

In this note we study the linear dynamics of scalar graviton in a de Sitter background in the infrared limit of the healthy extension of Ho\v{r}ava-Lifshitz gravity with the dynamical critical exponent $z=3$. Both our analytical and numerical results show that the non-zero Fourier modes of scalar graviton oscillate with an exponentially damping amplitude on the sub-horizon scale, while on the super-horizon scale, the phases are frozen and they approach to some asymptotic values. In addition, as the case of the non-zero modes on super-horizon scale, the zero mode also initially decays exponentially and then approaches to an asymptotic constant value.

We study a phenomenological dark energy model which is rooted in the Veneziano ghost of QCD. In this dark energy model, the energy density of dark energy is proportional to Hubble parameter and the proportional coefficient is of the order $\Lambda^3_{QCD}$, where $\Lambda_{QCD}$ is the mass scale of QCD. The universe has a de Sitter phase at late time and begins to accelerate at redshift around $z_{acc}\sim0.6$. We also fit this model and give the constraints on model parameters, with current observational data including SnIa, BAO, CMB, BBN and Hubble parameter data. We find that the squared sound speed of the dark energy is negative, which may cause an instability. We also study the cosmological evolution of the dark energy with interaction with cold dark matter.

Using the full radiation transfer function, we numerically calculate the CMB angular bispectrum seeded by the compensated magnetic scalar density mode. We find that, for the string inspired primordial magnetic fields characterized by index $n_B=-2.9$ and mean-field amplitude $B_{\lam}=9{\rm nG}$, the angular bispectrum is dominated by two primordial magnetic shapes. The first magnetic shape looks similar to the one from local-type primordial curvature perturbations, so both the amplitude and profile of the Komatsu-Spergel estimator (reduced bispectrum) seeded by this shape are almost the same as those of the primary CMB anisotropies. However, for different parameter sets ($l_1,l_2$), this "local-type" reduced bispectrum oscillates around different asymptotic values in the high-$l_3$ regime because of the effect of the Lorentz force, which is exerted by the primordial magnetic fields on the charged baryons. This feature is different from the standard case where all modes approach to zero asymptotically in the high-$l$ limit. On the other hand, the second magnetic shape appears only in the primordial magnetic field model. The amplitude of the Komatsu-Spergel estimator sourced by the second shape diverges in the low-$l$ regime because of the negative slope of shape. In the high-$l$ regime, this amplitude is approximately equal to that of the first estimator, but with a reversal phase.

Five-dimensional spacetimes of two orbifold 3-branes are studied, by assuming that {\em the two 3-branes are spatially homogeneous, isotropic, and independent of time}, following the so-called "bulk-based" approach. The most general form of the metric is obtained, and the corresponding field equations are divided into three groups, one is valid on each of the two 3-branes, and the third is valid in the bulk. The Einstein tensor on the 3-branes is expressed in terms of the discontinuities of the first-order derivatives of the metric coefficients. Thus, once the metric is known in the bulk, the distribution of the Einstein tensor on the two 3-branes is uniquely determined. As applications, we consider two different cases, one is in which the bulk is locally $AdS_{5}$, and the other is where it is vacuum. In some cases, it is shown that the universe is first decelerating and then accelerating. The global structure of the bulk as well as the 3-branes is also studied, and found that in some cases the solutions may represent the collision of two orbifold 3-branes. The applications of the formulas to the studies of the cyclic universe and the cosmological constant problem are also pointed out.

We examine in some detail the influence of the systematics in different data sets including type Ia supernova sample, baryon acoustic oscillation data and the cosmic microwave background information on the fitting results of the Chevallier-Polarski-Linder parametrization. We find that the systematics in the data sets does influence the fitting results and leads to different evolutional behavior of dark energy. To check the versatility of Chevallier-Polarski-Linder parametrization, we also perform the analysis on the Wetterich parametrization of dark energy. The results show that both the parametrization of dark energy and the systematics in data sets influence the evolutional behavior of dark energy.

We investigate dynamical behavior of the equation of state of dark energy $w_{de}$ by employing the linear-spline method in the region of low redshifts from observational data (SnIa, BAO, CMB and 12 $H(z)$ data). The redshift is binned and $w_{de}$ is approximated by a linear expansion of redshift in each bin. We leave the divided points of redshift bins as free parameters of the model, the best-fitted values of divided points will represent the turning positions of $w_{de}$ where $w_{de}$ changes its evolving direction significantly (if there exist such turnings in our considered region). These turning points are natural divided points of redshift bins, and $w_{de}$ between two nearby divided points can be well approximated by a linear expansion of redshift. We find two turning points of $w_{de}$ in $z\in(0,1.8)$ and one turning point in $z\in (0,0.9)$, and $w_{de}(z)$ could be oscillating around $w=-1$. Moreover, we find that there is a $2\sigma$ deviation of $w_{de}$ from -1 around $z=0.9$ in both correlated and uncorrelated estimates.

It is possible that there exist some interactions between dark energy (DE) and dark matter (DM), and a suitable interaction can alleviate the coincidence problem. Several phenomenological interacting forms are proposed and are fitted with observations in the literature. In this paper we investigate the possible interaction in a way independent of specific interacting forms by use of observational data (SNe, BAO, CMB and Hubble parameter). We divide the whole range of redshift into a few bins and set the interacting term $\delta(z)$ to be a constant in each redshift bin. We consider four parameterizations of the equation of state $w_{de}$ for DE and find that $\delta(z)$ is likely to cross the non-interacting ($\delta=0$) and have an oscillation form. It suggests that to study the interaction between DE and DM, more general phenomenological forms of the interacting term should be considered.

Singularities in $(3+1)$-dimensional Horava-Lifshitz (HL) theory of gravity are studied. These singularities can be divided into scalar, non-scalar curvature, and coordinate singularities. Because of the foliation-preserving diffeomorphisms of the theory, the number of scalars that can be constructed from the extrinsic curvature tensor $K_{ij}$, the 3-dimensional Riemann tensor and their derivatives is much large than that constructed from the 4-dimesnional Riemann tensor and its derivatives in general relativity (GR). As a result, even for the same spacetime, it may be singular in the HL theory but not in GR. Two representative families of solutions with projectability condition are studied, one is the (anti-) de Sitter Schwarzschild solutions, and the other is the Lu-Mei-Pope (LMP) solutions written in a form satisfying the projectability condition - the generalized LMP solutions. The (anti-) de Sitter Schwarzschild solutions are vacuum solutions of both HL theory and GR, while the LMP solutions with projectability condition satisfy the HL equations coupled with an anisotropic fluid with heat flow. It is found that the scalars $K$ and $K_{ij}K^{ij}$ are singular only at the center for the de Sitter Schwarzschild solution, but singular at both the center and $ r = (3M/|\Lambda|)^{1/3}$ for the anti-de Sitter Schwarzschild solution. The singularity at $ r = (3M/|\Lambda|)^{1/3}$ is absent in GR. In addition, all the generalized LMP solutions have two scalar curvature singularities, located at either $r = 0$ and $r=r_{s} > 0$, or $r=r_{1}$ and $r= r_{2}$ with $r_{2} > r_{1} > 0$, or $r=r_{s} > 0$ and $r = \infty$, depending on the choice of the free parameter $\lambda$.

We use the Constitution supernova, the baryon acoustic oscillation, the cosmic microwave background, and the Hubble parameter data to analyze the evolution property of dark energy. We obtain different results when we fit different baryon acoustic oscillation data combined with the Constitution supernova data to the Chevallier-Polarski-Linder model. We find that the difference stems from the different values of $\Omega_{m0}$. We also fit the observational data to the model independent piecewise constant parametrization. Four redshift bins with boundaries at $z=0.22$, 0.53, 0.85 and 1.8 were chosen for the piecewise constant parametrization of the equation of state parameter $w(z)$ of dark energy. We find no significant evidence for evolving $w(z)$. With the addition of the Hubble parameter, the constraint on the equation of state parameter at high redshift isimproved by 70%. The marginalization of the nuisance parameter connected to the supernova distance modulus is discussed.

In this paper we study the entropy perturbations in N-flation by using the $\d\ma{N}$ formalism. We calculate the entropy corrections to the power spectrum of the overall curvature perturbation $P_{\z}$. We obtain an analytic form of the transfer coefficient $T^2_{\ma{R}\ma{S}}$, which describes the correlation between the curvature and entropy perturbations, and investigate its behavior numerically. It turns out that the entropy perturbations cannot be neglected in N-flation, because the amplitude of entropy components is approximately in the same order as the adiabatic one at the end of inflation $T^2_{\ma{R}\ma{S}}\sim\ma{O}(1)$. The spectral index $n_S$ is calculated and it becomes smaller after the entropy modes are taken into account, i.e., the spectrum becomes redder, compared to the pure adiabatic case. Finally we study the modified consistency relation of N-flation, and find that the tensor-to-scalar ratio ($r\simeq0.006$) is greatly suppressed by the entropy modes, compared to the pure adiabatic one ($r\simeq0.017$) at the end of inflation.

We investigate the linear cosmological perturbations of Ho\v{r}ava-Lifshitz gravity in a FRW universe without any matter. Our results show that a new gauge invariant dynamical scalar mode emerges, due to the gauge transformation under the "foliation-preserving" diffeomorphism and "projectability condition", and it can produce a scale invariant power spectrum. In the infrared regime with $\lambda=1$, the dynamical scalar degree of freedom turns to be a non-dynamical one at the linear order level.

A new dark energy model, named "agegraphic dark energy", has been proposed recently, based on the so-called K\'{a}rolyh\'{a}zy uncertainty relation, which arises from quantum mechanics together with general relativity. In this note, we extend the original agegraphic dark energy model by including the interaction between agegraphic dark energy and pressureless (dark) matter. In the interacting agegraphic dark energy model, there are many interesting features different from the original agegraphic dark energy model and holographic dark energy model. The similarity and difference between agegraphic dark energy and holographic dark energy are also discussed.

The constraint on the total energy in a given spatial region is given from holography by the mass of a black hole which just fits in that region, which leads to an UV/IR relation: the maximal energy density in that region is proportional to $M_p^2/L^2$, where $M_p$ is the Planck mass and $L$ is the spatial scale of that region under consideration. Assuming the maximal black hole in the universe is formed through gravitational collapse of perturbations in the universe, then the "Jeans" scale of the perturbations gives a causal connection scale $R_{\rm CC}$. For gravitational perturbations, $R^{-2}_{\rm CC}={\rm Max}(\dot H+2H^2, -\dot H)$ for a flat universe. We study the cosmological dynamics of the corresponding vacuum energy density by choosing the causal connection scale as the IR cutoff in the UV/IR relation, in the cases of the vacuum energy density as an independently conserved energy component and an effective dynamical cosmological constant, respectively. It turns out that only the case with the choice $R_{\rm CC}^{-2}= \dot H +2H^2$, could be consistent with the current cosmological observations when the vacuum density appears as an independently conserved energy component. In this case, the model is called holographic Ricci scalar dark energy model in the literature.

In this work, we consider the cosmological constraints on the new agegraphic dark energy (NADE) proposed in arXiv:0708.0884, by using the observational data of type Ia supernovae (SNIa), the shift parameter from cosmic microwave background (CMB) and the baryon acoustic oscillation (BAO) peak from large scale structures (LSS). Thanks to its special analytic features in the radiation-dominated and matter-dominated epochs, NADE is a single-parameter model in practice because once the single model parameter $n$ is given, all other physical quantities of NADE can be determined correspondingly. The joint analysis gives the best-fit value (with $1\sigma$ uncertainty) $n=2.716^{+0.111}_{-0.109}$, and the derived $\Omega_{m0}$, $\Omega_{q0}$ and $w_{q0}$ (with $1\sigma$ uncertainties) are $0.295^{+0.020}_{-0.020}$, $0.705^{+0.020}_{-0.020}$ and $-0.794^{+0.006}_{-0.005}$, respectively. In addition, we find that the coincidence problem could be solved naturally in the NADE model provided that $n$ is of order unity.

In this note, we propose a new model of agegraphic dark energy based on the K\'{a}rolyh\'{a}zy relation, where the time scale is chosen to be the conformal time $\eta$ of the Friedmann-Robertson-Walker (FRW) universe. We find that in the radiation-dominated epoch, the equation-of-state parameter of the new agegraphic dark energy $w_q=-1/3$ whereas $\Omega_q=n^2a^2$; in the matter-dominated epoch, $w_q=-2/3$ whereas $\Omega_q=n^2a^2/4$; eventually, the new agegraphic dark energy dominates; in the late time $w_q\to -1$ when $a\to\infty$, and the new agegraphic dark energy mimics a cosmological constant. In every stage, all things are consistent. The confusion in the original agegraphic dark energy model proposed in arXiv:0707.4049 disappears in this new model. Furthermore, $\Omega_q\ll 1$ is naturally satisfied in both radiation-dominated and matter-dominated epochs where $a\ll 1$. In addition, we further extend the new agegraphic dark energy model by including the interaction between the new agegraphic dark energy and background matter. In this case, we find that $w_q$ can cross the phantom divide.

We consider the evolution of the vacuum energy in the DGP model according to the holographic principle under the assumption that the relation linking the IR and UV cut-offs still holds in this scenario. The model is studied when the IR cut-off is chosen to be the Hubble scale $H^{-1}$, the particle horizon $R_{\rm ph}$ and the future event horizon $R_{\rm eh}$, respectively. And the two branches of the DGP model are also taken into account. Through numerical analysis, we find that in the cases of $H^{-1}$ in the (+) branch and $R_{\rm eh}$ in both branches, the vacuum energy can play the role of dark energy. Moreover, when considering the combination of the vacuum energy and the 5D gravity effect in both branches, the equation of state of the effective dark energy may cross -1, which may lead to the Big Rip singularity. Besides, we constrain the model with the Type Ia supernovae and baryon oscillation data and find that our model is consistent with current data within $1\sigma$, and that the observations prefer either a pure holographic dark energy or a pure DGP model

A new dark energy model, named as ``agegraphic dark energy'', has been proposed by one of us (R. G. Cai) in arXiv:0707.4049, based on the K\'{a}rolyh\'{a}zy uncertainty relation, which arises from the quantum mechanics together with general relativity. Then, in arXiv:0707.4052, it has been extended by including the interaction between the agegraphic dark energy and the pressureless (dark) matter. In this note, we investigate the agegraphic dark energy models without and with interaction by means of statefinder diagnostic and w-w' analysis.

In this paper, we use the Type Ia supernova data as well as the CMB and LSS data to constrain the agegraphic dark energy model recently proposed by Cai. Due to its peculiar nature, the parameter $n$ of this model cannot be well constrained by the SNIa data, while the other parameter $\Omega_{m0}$ can be constrained to be $0.34\pm0.04$. When combined with CMB and LSS data, the range of $1\sigma$ confidence level for $n$ is greatly narrowed, albeit still very large. The best fit result is $\Omega_{m0}=0.28\pm0.02$, which is consistent with most observations like WMAP and SDSS, and $n=3.4$, of which a meaningful range of confidence level can not be obtained due to the fact that the contours are not closed. Despite of this result, we conclude that for $n>1$ this model is consistent with SNIa, CMB and LSS observations. Furthermore, the fitting results indicate a generalized definition for the agegraphic dark energy.

Quantum mechanics together with general relativity leads to the K\'arolyh\'azy relation and a corresponding energy density of quantum fluctuations of space-time. Based on the energy density we propose a dark energy model, in which the age of the universe is introduced as the length measure. This dark energy is consistent with astronomical data if the unique numerical parameter in the dark energy model is taken to be a number of order one. The dark energy behaves like a cosmological constant at early time and drives the universe to an eternally accelerated expansion with power-law form at late time. In addition, we point out a subtlety in this kind of dark energy model.

Weyl's idea on scale invariance was resurrected by Cheng in 1988. The requirement of local scale invariance leads to a completely new vector field, which we call the ``Cheng-Weyl vector field''. The Cheng-Weyl vector field couples only to a scalar field and the gravitational field naturally. It does not interact with other known matters in the standard model of particle physics. In the present work, the (generalized) Cheng-Weyl vector field coupled with the scalar field and its cosmological application are investigated. A mixture of the scalar field and a so-called ``cosmic triad'' of three mutually orthogonal Cheng-Weyl vector fields is regarded as the dark energy in the universe. The cosmological evolution of this ``mixed'' dark energy model is studied. We find that the effective equation-of-state parameter of the dark energy can cross the phantom divide $w_{de}=-1$ in some cases; the first and second cosmological coincidence problems can be alleviated at the same time in this model.

Motivated by the AMS project, we assume that after the Big Bang or inflation epoch, antimatter was repelled onto one brane which is separated from our brane where all the observational matter resides. It is suggested that CP may be spontaneously broken, the two branes would correspond to ground states for matter and antimatter respectively. Generally a complex scalar field which is responsible for the spontaneous CP violation, exists in the space between the branes and causes a repulsive force against the gravitation. A possible potential barrier prevents the mater(antimatter) particles to enter the space between two branes. However, by the quantum tunnelling, a sizable anti-matter flux may come to our brane. In this work by considering two possible models, i.e. the naive flat space-time and Randall-Sundrum models and using the observational data on the visible matter in our universe as inputs, we derive the antimatter flux which would be observed by the AMS detector.

We investigate the age constraint on the agegraphic dark energy model by using two old galaxies (LBDS 53W091 and LBDS 53W069) and the old high redshift quasar APM 08279+5255. We find that the agegraphic dark energy model can easily accommodate LBDS 53W091 and LBDS 53W069. To accommodate APM 08279+5255, one can take the reduced Hubble parameter as large as h=0.64, when the fraction matter energy density $\Omega_{m0}\approx 0.22$.

In this paper we first obtain Friedmann equations for the $(n-1)$-dimensional brane embedded in the $(n+1)$-dimensional bulk, with intrinsic curvature term of the brane included in the action (DGP model). Then, we show that one can always rewrite the Friedmann equations in the form of the first law of thermodynamics, $dE=TdS+WdV$, at apparent horizon on the brane, regardless of whether there is the intrinsic curvature term on the brane or a cosmological constant in the bulk. Using the first law, we extract the entropy expression of the apparent horizon on the brane. We also show that in the case without the intrinsic curvature term, the entropy expressions are the same by using the apparent horizon on the brane and by using the bulk geometry. When the intrinsic curvature appears, the entropy of apparent horizon on the brane has two parts, one part follows the $n$-dimensional area formula on the brane, and the other part is the same as the entropy in the case without the intrinsic curvature term. As an interesting result, in the warped DGP model, the entropy expression in the bulk and on the brane are not the same. This is reasonable, since in this model gravity on the brane has two parts, one induced from the $(n+1)$-dimensional bulk gravity and the other due to the intrinsic curvature term on the brane.

One of the puzzles of the dark energy problem is the (first) cosmological coincidence problem, namely, {\em why does our universe begin the accelerated expansion recently? why are we living in an epoch in which the dark energy density and the dust matter energy density are comparable?} On the other hand, cosmological observations hint that the equation-of-state parameter (EoS) of dark energy crossed the phantom divide $w_{de}=-1$ in the near past. Many dark energy models whose EoS can cross the phantom divide have been proposed. However, to our knowledge, these models with crossing the phantom divide only provide the possibility that $w_{de}$ can cross -1. They do not answer another question, namely, {\em why crossing the phantom divide occurs recently?} Since in many existing models whose EoS can cross the phantom divide, $w_{de}$ undulates around -1 randomly, {\em why are we living in an epoch $w_{de}<-1$?} This can be regarded as the second cosmological coincidence problem. In this work, the cosmological evolution of the vector-like dark energy interacting with background perfect fluid is investigated. We find that the first and second cosmological coincidence problems can be alleviated at the same time in this scenario.

The gravitational collapse of a spherically symmetric star, made of a dust fluid, $\rho_{DM}$, in a background of dark energy, $p = w\rho,\; (w < -1/3)$, is studied. It is found that when only dark energy is present, black holes are {\em never} formed. When both of them are present, black holes can be formed, due to the condensation of the dust fluid. Initially the dust fluid may not play an important role, but, as time increases, it will dominate the collapse and finally leads to formation of black holes. This result remains true even when the interaction between the dust fluid and dark energy does not vanish. When $w < -1$ (phantoms), some models can also be interpreted as representing the death of a white hole that ejects both dust and phantoms. The ejected matter re-collapses to form a black hole.

We consider a brane cosmological model with energy exchange between brane and bulk. Parameterizing the energy exchange term by the scale factor and Hubble parameter, we are able to exactly solve the modified Friedmann equation on the brane. In this model, the equation of state for the effective dark energy has a transition behavior changing from $w_{de}^{eff}>-1$ to $w_{de}^{eff}<-1$, while the equation of state for the dark energy on the brane has $w>-1$. Fitting data from type Ia supernova, Sloan Digital Sky Survey and Wilkinson Microwave Anisotropy Probe, our universe is predicted now in the state of super-acceleration with $w_{de0}^{eff}=-1.21$.

Recently a lot of attention has been given to building dark energy models in which the equation-of-state parameter $w$ can cross the phantom divide $w=-1$. However, to our knowledge, these models with crossing the phantom divide only provide the possibility that $w$ can cross -1. They do not answer another question: {\em why crossing phantom divide occurs recently?} Since in many existing models whose equation-of-state parameter can cross the phantom divide, $w$ undulates around -1 randomly, {\em why are we living in an epoch $w<-1$?} This can be regarded as the second cosmological coincidence problem. In this note, we propose a possible approach to alleviate this problem within a hybrid dark energy model.

Recently, many dark energy models whose equation-of-state parameter can cross the phantom divide $w_{de}=-1$ have been proposed. In a previous paper [Class. Quant. Grav. {\bf 22}, 3189 (2005); hep-th/0501160], we suggest such a model named hessence, in which a non-canonical complex scalar field plays the role of dark energy. In this work, the cosmological evolution of the hessence dark energy is investigated. We consider two cases: one is the hessnece field with an exponential potential, and the other is with a (inverse) power law potential. We separately investigate the dynamical system with four different interaction forms between hessence and background perfect fluid. It is found that the big rip never appears in the hessence model, even in the most general case, beyond particular potentials and interaction forms.

Using the evolution history of the universe, one can make constraint on the parameter space of dynamic dark energy models. We discuss two different parameterized dark energy models. Our results further restrict the combined constraints obtained from supernova and WMAP observations. From the allowed parameter space, it is found that our universe will experience an eternal acceleration. We also estimate the bound on the physically relevant regions both in the re-inflationary and inflationary phases.

Recently the 1/R gravity has been proposed in order to explain the accelerated expansion of the universe. However, it was argued that the 1/R gravity conflicts with solar system tests. While this statement is true if one views the 1/R gravity as an effective theory, we find that this difficulty might be avoided if one treats the 1/R term as a correction to the scalar curvature term in the high curvature limit $R \gg \mu^2$.

Recently a lot of attention has been drawn to build dark energy model in which the equation-of-state parameter $w$ can cross the phantom divide $w=-1$. One of models to realize crossing the phantom divide is called quintom model, in which two real scalar fields appears, one is a normal scalar field and the other is a phantom-type scalar field. In this paper we propose a non-canonical complex scalar field as the dark energy, which we dub ``hessence'', to implement crossing the phantom divide, in a similar sense as the quintom dark energy model. In the hessence model, the dark energy is described by a single field with an internal degree of freedom rather than two independent real scalar fields. However, the hessence is different from an ordinary complex scalar field, we show that the hessence can avoid the difficulty of the Q-balls formation which gives trouble to the spintessence model (An ordinary complex scalar field acts as the dark energy). Furthermore, we find that, by choosing a proper potential, the hessence could correspond to a Chaplygin gas at late times.

Recent type Ia supernovas data seemingly favor a dark energy model whose equation of state $w(z)$ crosses -1 very recently, which is a much more amazing problem than the acceleration of the universe. In this paper we show that it is possible to realize such a crossing without introducing any phantom component in a Gauss-Bonnet brane world with induced gravity, where a four dimensional curvature scalar on the brane and a five dimensional Gauss-Bonnet term in the bulk are present. In this realization, the Gauss-Bonnet term and the mass parameter in the bulk play a crucial role.

A solution of codimension 2 brane is found for which 4 dimensional Friedmann cosmology is recovered on the brane with time dependent tension, in the Einstein frame. The effective parameter $p/\rho$ of equation of state on the brane can be quintessence like, de Sitter like or phantom like, depending on integration constants of the solution.

In this paper we present a hybrid model of k-essence and chameleon, named as k-chameleon. In this model, due to the chameleon mechanism, the directly strong coupling between the k-chameleon field and matters (cold dark matters and baryons) is allowed. In the radiation dominated epoch, the interaction between the k-chameleon field and background matters can be neglected, the behavior of the k-chameleon therefore is the same as that of the ordinary k-essence. After the onset of matter domination, the strong coupling between the k-chameleon and matters dramatically changes the result of the ordinary k-essence. We find that during the matter-dominated epoch, only two kinds of attractors may exist: one is the familiar {\bf K} attractor and the other is a completely {\em new}, dubbed {\bf C} attractor. Once the universe is attracted into the {\bf C} attractor, the fraction energy densities of the k-chameleon $\Omega_{\phi}$ and dust matter $\Omega_m$ are fixed and comparable, and the universe will undergo a power-law accelerated expansion. One can adjust the model so that the {\bf K} attractor do not appear. Thus, the k-chameleon model provides a natural solution to the cosmological coincidence problem.

We study a cosmological model in which phantom dark energy is coupled to dark matter by phenomenologically introducing a coupled term to the equations of motion of dark energy and dark matter. This term is parameterized by a dimensionless coupling function $\delta$, Hubble parameter and the energy density of dark matter, and it describes an energy flow between the dark energy and dark matter. We discuss two cases: one is the case where the equation-of-state $\omega_e$ of the dark energy is a constant; the other is that the dimensionless coupling function $\delta$ is a constant. We investigate the effect of the interaction on the evolution of the universe, the total lifetime of the universe, and the ratio of the period when the universe is in the coincidence state to its total lifetime. It turns out that the interaction will produce significant deviation from the case without the interaction.

We investigate the cosmological evolution of an interacting phantom energy model in which the phantom field has interaction with the dark matter. We discuss the existence and stability of scaling solutions for two types of specific interactions. One is motivated by the conformal transformation in string theory and the other is motivated by analogy with dissipation. In the former case, there exist no scaling solutions. In the latter case, there exist stable scaling solutions, which may give a phenomenological solution of the coincidence problem. Furthermore, the universe either accelerates forever or ends with a singularity, which is determined by not only the model parameters but also the initial velocity of the phantom field.

We extend Green's function method developed by Stewart and Gong to calculate the power spectrum of density perturbation in the case with a time-dependent sound speed, and explicitly give the power spectrum and spectral index up to second-order corrections in the slow-roll expansion. The case of tachyon inflation is included as a special case.

We discuss an inflation model, in which the inflation is driven by a single scalar field with exponential potential on a warped DGP brane. In contrast to the power law inflation in standard model, we find that the inflationary phase can exit spontaneously without any mechanism. The running of the index of scalar perturbation spectrum can take an enough large value to match the observation data, while other parameters are in a reasonable region.

An elegant approach, which incorporates the effect of the stringy spacetime uncertainty relation, to calculate power spectra of fluctuations during inflation has been suggested by Brandenberger and Ho. In this approach, one of important features is the appearance of an upper bound on the comoving momentum $k$, at which the stringy spacetime uncertainty relation is saturated. As a result, the time-dependent upper bound leads us to choose naturally a set of initial vacua for each mode, in which the stringy uncertainty relation is saturated. In this note, with that set of vacua we calculate power spectrum of curvature fluctuation for a power law inflation, up to the leading order of a parameter describing the spacetime noncommutativity. It turns out that this choice of initial vacuum has a significant effect on the power spectrum of fluctuations.

If the source of the current accelerating expansion of the universe is a positive cosmological constant, Banks and Fischler argued that there exists an upper limit of the total number of e-foldings of inflation. We further elaborate on the upper limit in the senses of viewing the cosmological horizon as the boundary of a cavity and of the holographic D-bound in a de Sitter space. Assuming a simple evolution model of inflation, we obtain an expression of the upper limit in terms of the cosmological constant, the energy density of inflaton when the inflation starts, the energy density as the inflation ends, and reheating temperature. We discuss how the upper limit is modified in the different evolution models of the universe. The holographic D-bound gives more high upper limit than the entropy threshold in the cavity. For the most extremal case where the initial energy density of inflation is as high as the Planck energy, and the reheating temperature is as low as the energy scale of nucleosynthesis, the D-bound gives the upper limit as 146 and the entropy threshold as 122. For reasonable assumption in the simplest cosmology, the holographic D-bound leads to a value about 85, while the cavity model gives a value around 65 for the upper limit, which is close to the value in order to solve the flatness problem and horizon problem in the hot big bang cosmology.