We report pulsar timing observations carried out in L-band with NTSC's 40-meter Haoping Radio Telescope (HRT), which was constructed in 2014. The observations were carried out using the pulsar machine we developed. Timing observations toward millisecond pulsar J0437-4715 obtains a timing residual (r.m.s) of 397ns in the time span of 284 days. And our observations successfully detected Crab pulsar's glitch that happened on July 23rd, 2019.

We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties and both agree with TT(BIPM17), a post-corrected timescale produced by the Bureau International des Poids et Mesures (BIPM). We show that both methods could detect significant errors in TT(BIPM17) if they were present. We estimate the stability of the atomic clocks from which TT(BIPM17) is derived using observations of four rubidium fountain clocks at the US Naval Observatory. Comparing the power spectrum of TT(IPTA16) with that of these fountain clocks suggests that pulsar-based timescales are unlikely to contribute to the stability of the best timescales over the next decade, but they will remain a valuable independent check on atomic timescales. We also find that the stability of the pulsar-based timescale is likely to be limited by our knowledge of solar-system dynamics, and that errors in TT(BIPM17) will not be a limiting factor for the primary goal of the IPTA, which is to search for the signatures of nano-Hertz gravitational waves.

Employing multiple pulsars and using an appropriate algorithm to establish ensemble pulsar timescale can reduce the influences of various noises on the long-term stability of pulsar timescale, compared to a single pulsar. However, due to the low timing precision and the significant red noises of some pulsars, their participation in the construction of ensemble pulsar timescale is often limited. Inspired by the principle of solving non-stationary sequence modeling using co-integration theory, we puts forward an algorithm based on the co-integration theory to establish ensemble pulsar timescale. It is found that this algorithm can effectively suppress some noise sources if a co-integration relationship between different pulsar data exist. Different from the classical weighted average algorithm, the co-integration method provides the chances of the pulsar with significant red noises to attend the establishment of ensemble pulsar timescale. Based on the data from the North American Nanohertz Observatory for Gravitational Waves, we found that the co-integration algorithm can successfully reduce several timing noises and improve the long-term stability of the ensemble pulsar timescale.

Relic gravitational waves (RGWs) , a background originated during inflation, would give imprints on the pulsar timing residuals. This makes RGWs be one of important sources for detection using the method of pulsar timing. In this paper, we discuss the effects of RGWs on the single pulsar timing, and give quantitively the timing residuals caused by RGWs with different model parameters. In principle, if the RGWs are strong enough today, they can be detected by timing a single millisecond pulsar with high precision after the intrinsic red noise in pulsar timing residuals were understood, even though observing simultaneously multiple millisecond pulsars is a more powerful technique in extracting gravitational wave signals. We corrected the normalization of RGWs using observations of the cosmic microwave background (CMB), which leads to the amplitudes of RGWs being reduced by two orders of magnitude or so compared to our previous works. We made new constraints on RGWs using the recent observations from the Parkes Pulsar Timing Array, employing the tensor-to-scalar ratio $r=0.2$ due to the tensor-type polarization observations of CMB by BICEP2 as a referenced value even though it has been denied. Moreover, the constraints on RGWs from CMB and BBN (Big Bang nucleosynthesis) will also be discussed for comparison.

The pulsar timing residuals induced by gravitational waves from non-evolving single binary sources are affected by many parameters related to the relative positions of the pulsar and the gravitational wave sources. We will fully analyze the effects due to different parameters one by one. The standard deviations of the timing residuals will be calculated with a variable parameter fixing a set of other parameters. The orbits of the binary sources will be generally assumed to be elliptical. The influences of different eccentricities on the pulsar timing residuals will also studied in detail. We find that effects of the related parameters are quite different, and some of them present certain regularities.

In the non-standard model of relic gravitational waves (RGWs) generated in the early universe, the theoretical spectrum of is mainly described by an amplitude $r$ and a spectral index $\beta$, the latter usually being determined by the slope of the inflation potential. Pulsar timing arrays (PTAs) data have imposed constraints on the amplitude of strain spectrum for a power-law form as a phenomenological model. Applying these constraints to a generic, theoretical spectrum with $r$ and $\beta$ as independent parameters, we convert the PTAs constraint into an upper bound on the index $\beta$, which turns out to be less stringent than those upper bounds from BBN, CMB, and LIGO/VIRGO, respectively. Moreover, it is found that PTAs constrain the non-standard RGWs more stringent than the standard RGWs. If the condition of the quantum normalization is imposed upon a theoretical spectrum of RGWs, $r$ and $\beta$ become related. With this condition, a minimum requirement of the horizon size during inflation is greater than the Planck length results in an upper bound on $\beta$, which is comparable in magnitude to that by PTAs. When both PTAs and the quantum normalization are applied to a theoretical spectrum of RGWs, constraints can be obtained for other cosmic processes of the early universe, such as the reheating, a process less understood observationally so far. The resulting constraint is consistent with the slow-roll, massive scalar inflation model. The future SKA will be able to constrain RGWs further and might even detect RGWs, rendering an important probe to the very early universe.

The pulsar timing residuals induced by gravitational waves from non-evolving single binary sources with general elliptical orbits will be analyzed. For different orbital eccentricities, the timing residuals present different properties. The standard deviations of the timing residuals induced by a fixed gravitational wave source will be calculated for different values of the eccentricity. We will also analyze the timing residuals of PSR J0437-4715 induced by one of the best known single gravitational wave sources, the supermassive black hole binary in the blazar OJ287.

We obtained the analytic solutions of relic gravitational waves (RGWs) for the slow-roll inflation with a power-law form potential of the scalar field, $V=\lambda\phi^n$. Based on a reasonable range of $n$ constrained by cosmic microwave background (CMB) observations, we give tight constraints of the tensor-to-scalar ratio $r$ and the inflation expansion index $\beta$ for the fixed scalar spectral index $n_s$. Even though, the spectrum of RGWs in low frequencies is hardly depends on any parameters, the high frequency parts will be affected by several parameters, such as $n_s$, the reheating temperature $T_{RH}$ and the index $\beta_s$ describing the expansion from the end of inflation to the reheating process. We analyzed in detail all the factors which would affect the spectrum of RGWs in high frequencies including the quantum normalization. We found that the future GW detectors SKA, eLISA, BBO and DECIGO are promising to catch the signals of RGWs. Furthermore, BBO and DECIGO have the potential not only to distinguish the spectra with different parameters but also to examine the validity of the quantum normalization.

According to the CMB observations, Mielczarek (\cite{Mielczarek}) evaluated the reheating temperature, which could help to determine the history of the Universe. In this paper, we recalculate the reheating temperature using the new data from WMAP 7 observations. Based on that, we list the approximate solutions of relic gravitational waves (RGWs) for various frequency bands. With the combination of the quantum normalization of RGWs when they are produced and the CMB observations, we obtain the relation between the tensor-to-scalar ratio $r$ and the inflation index $\beta$ for a given scalar spectral index $n_s$. As a comparison, the diagram $r-\beta$ in the slow-roll inflation model is also given. Thus, the observational limits of $r$ from CMB lead to the constraints on the value of $\beta$. Then, we illustrate the energy density spectrum of RGWs with the quantum normalization for different values of $r$ and the corresponding $\beta$. For comparison, the energy density spectra of RGWs with parameters based on slow-roll inflation are also discussed. We find that the values of $n_s$ affect the spectra of RGWs sensitively in the very high frequencies. Based on the current and planed gravitational wave detectors, we discuss the detectabilities of RGWs.

The decaying vacuum model (DV), treating dark energy as a varying vacuum, has been studied well recently. The vacuum energy decays linearly with the Hubble parameter in the late-times, $\rho_\Lambda(t) \propto H(t)$, and produces the additional matter component. We constrain the parameters of the DV model using the recent data-sets from supernovae, gamma-ray bursts, baryon acoustic oscillations, CMB, the Hubble rate and x-rays in galaxy clusters. It is found that the best fit of matter density contrast $\Omega_m$ in the DV model is much lager than that in $\Lambda$CDM model. We give the confidence contours in the $\Omega_m-h$ plane up to $3\sigma$ confidence level. Besides, the normalized likelihoods of $\Omega_m$ and $h$ are presented, respectively. %

Motivated by Yang-Mills dark energy model, we propose a new model by introducing a logarithmic correction. we find that this model can avoid the coincidence problem naturally and gives an equation of state $w$ smoothly crossing -1 if an interaction between dark energy and dark matter exists. It has a stable tracker solution as well. To confront with observations based on the combined data of SNIa, BAO, CMB and Hubble parameter, we obtain the best fit values of the parameters with $1\sigma, 2\sigma, 3\sigma$ errors for the noncoupled model: $\Omega_m=0.276\pm0.008^{+0.016+0.024}_{-0.015-0.022}$, $h=0.699\pm0.003\pm0.006\pm0.008$, and for the coupled model with a decaying rate $\gamma=0.2$: $\Omega_m=0.291\pm0.004^{+0.008+0.012}_{-0.007-0.011}$, $h=0.701\pm0.002\pm0.005\pm0.007$. In particular, it is found that the non-coupled model has a dynamic evolution almost undistinguishable to $\Lambda$CDM at the late-time Universe.

With LIGO having achieved its design sensitivity and the LIGO S5 strain data being available, constraints on the relic gravitational waves (RGWs) becomes realistic. The analytical spectrum of RGWs generated during inflation depends sensitively on the initial condition, which is generically described by the index $\beta$, the running index $\alpha_t$, and the tensor-to-scalar ratio $r$. By the LIGO S5 data of the cross-correlated two detectors, we obtain constraints on the parameters $(\beta, \alpha_t,r)$. As a main result, we have computed the theoretical signal-to noise ratio (SNR) of RGWs for various values of $(\beta, \alpha_t, r)$, using the cross-correlation for the given pair of LIGO detectors. The constraints by the indirect bound on the energy density of RGWs by BBN and CMB have been obtained, which turn out to be still more stringent than LIGO S5.

We review the quantum Yang-Mills condensate (YMC) dark energy models. As the effective Yang-Mills Lagrangian is completely determined by the quantum field theory, there is no adjustable parameter in the model except the energy scale. In this model, the equation-of-state (EOS) of the YMC dark energy, $w_y > -1$ and $w_y < -1$, can both be naturally realized. By studying the evolution of various components in the model, we find that, in the early stage of the universe, dark energy tracked the evolution of the radiation, i.e. $w_y \to 1/3$. However, in the late stage, $w_y$ naturally runs to the critical state with $w_y = -1$, and the universe transits from matter-dominated into dark energy dominated stage only at recently $z \sim 0.3$. These characters are independent of the choice of the initial condition, and the cosmic coincidence problem is avoided in the models. We also find that, if the possible interaction between YMC and dust matter is considered, the late time attractor solution may exist. In this case, the EOS of YMC must evolve from $w_y>0$ into $w_y < -1$, which is slightly suggested by the observations. At the same time, the total EOS in the attractor solution is $w_{tot} = -1$, the universe being the de Sitter expansion in the late stage, and the cosmic big rip is naturally avoided. These features are all independent of the interacting forms.

We study the impact of a running index $\alpha_t$ on the spectrum of relic gravitational waves (RGWs) over the whole range of frequency $(10^{-18}\sim 10^{10})$ Hz and reveal its implications in RGWs detections and in cosmology. Analytical calculations show that, although the spectrum of RGWs on low frequencies is less affected by $\alpha_t\ne 0$, but, on high frequencies, the spectrum is modified substantially. Investigations are made toward potential detections of the $\alpha_t$-modified RGWs for several kinds of current and planned detectors. The Advanced LIGO will likely be able to detect RGWs with $\alpha_t\ge 0$ for inflationary models with the inflation index $\beta=-1.956$ and the tensor-scalar ratio $r= 0.55$. The future LISA can detect RGWs for a much broader range of ($\alpha_t$, $\beta$, $r$), and will have a better chance to break a degeneracy between them. Constraints on $\alpha_t$ are estimated from several detections and cosmological observations. Among them, the most stringent one is from the bound of the Big Bang nucleosynthesis (BBN), and requires $\alpha_t < 0.008$ rather conservatively for any reasonable ($\beta$, $r$), preferring a nearly power-law spectrum of RGWs. In light of this result, one would expect the scalar running index $\alpha_s$ to be of the same magnitude as $\alpha_t$, if both RGWs and scalar perturbations are generated by the same scalar inflation.

As an extension of $\Lambda$CDM, the decaying vacuum model (DV) describes the dark energy as a varying vacuum whose energy density decays linearly with the Hubble parameter in the late-times, $\rho_\Lambda(t) \propto H(t)$, and produces the matter component. We examine the high-$z$ cosmic age problem in the DV model, and compare it with $\Lambda$CDM and the Yang-Mills condensate (YMC) dark energy model. Without employing a dynamical scalar field for dark energy, these three models share a similar behavior of late-time evolution. It is found that the DV model, like YMC, can accommodate the high-$z$ quasar APM 08279+5255, thus greatly alleviates the high-$z$ cosmic age problem. We also calculate the Statefinder $(r,s)$ and the {\it Om} diagnostics in the model. It is found that the evolutionary trajectories of $r(z)$ and $s(z)$ in the DV model are similar to those in the kinessence model, but are distinguished from those in $\Lambda$CDM and YMC. The ${\it Om}(z)$ in DV has a negative slope and its height depends on the matter fraction, while YMC has a rather flat ${\it Om}(z)$, whose magnitude depends sensitively on the coupling.

The quantum effective Yang-Mills condensate (YMC) dark energy model has some distinguished features that it naturally solves the coincidence problem and, at the same time, is able to give an equation of state $w$ crossing -1. In this work we further employ the Statefinder pair $(r,s)$ introduced by Sahni et al to diagnose the YMC model for three cases: the non-coupling, the YMC decaying into matter only, and the YMC decaying into both matter and radiation. The trajectories $(r,s)$ and $(r,q)$, and the evolutions $r(z)$, $s(z)$ are explicitly presented. It is found that, the YMC model in all three cases has $r\simeq 1$ for $ z < 10$ and $s\simeq 0$ for $z<5$ with only small deviations $\simeq 0.02$, quite close to the cosmological constant model (LCDM), but is obviously differentiated from other dark energy models, such as quiesence, kinessence etc.