We use three different data sets, specifically $H(z)$ measurements from cosmic chronometers, the HII-galaxy Hubble diagram, and reconstructed quasar-core angular-size measurements, to perform a joint analysis of three flat cosmological models: the $R_{\rm h}=ct$ Universe, $\Lambda$CDM, and $w$CDM. For $R_{\rm h}=ct$, the 1$\sigma$ best-fit value of the Hubble constant $H_0$ is $62.336\pm1.464$ $\mathrm{km \ s^{-1} \ Mpc^{-1}}$, which matches previous measurements ($\sim 63$ $\mathrm{km \ s^{-1} \ Mpc^{-1}}$) based on best fits to individual data sets. For $\Lambda$CDM, our inferred value of the Hubble constant, $H_0=67.013\pm2.578$ $\mathrm{km \ s^{-1} \ Mpc^{-1}}$, is more consistent with the ${\it Planck}$ optimization than the locally measured value using $\mbox{Cepheid}$ variables, and the matter density $\Omega_{\rm m}=0.347\pm0.049$ similarly coincides with its ${\it Planck}$ value to within 1$\sigma$. For $w$CDM, the optimized parameters are $H_0=64.718\pm3.088$ $\mathrm{km \ s^{-1} \ Mpc^{-1}}$, $\Omega_{\rm m}=0.247\pm0.108$ and $w=-0.693\pm0.276$, also consistent with ${\it Planck}$. A direct comparison of these three models using the Bayesian Information Criterion shows that the $R_{\rm h}=ct$ universe is favored by the joint analysis with a likelihood of $\sim 97\%$ versus $\lesssim 3\%$ for the other two cosmologies.

In this paper, we present the application of a new method measuring Hubble parameter $H(z)$ by using the anisotropy of luminosity distance($d_{L}$) of the gravitational wave(GW) standard sirens of neutron star(NS) binary system. The method has never been put into practice so far due to the lack of the ability of detecting GW. However, LIGO's success in detecting GW of black hole(BH) binary system merger announced the potential possibility of this new method. We apply this method to several GW detecting projects, including Advanced LIGO(aLIGO), Einstein Telescope(ET) and DECIGO, and evaluate its constraint ability on cosmological parameters of $H(z)$. It turns out that the $H(z)$ by aLIGO and ET is of bad accuracy, while the $H(z)$ by DECIGO shows a good one. We simulate $H(z)$ data at every 0.1 redshift span using the error information of $H(z)$ by DECIGO, and put the mock data into the forecasting of cosmological parameters. Compared with the previous data and method, we get an obviously tighter constraint on cosmological parameters by mock data, and a concomitantly higher value of Figure of Merit(FoM, the reciprocal of the area enclosed by the $2\sigma$ confidence region). For a 3-year-observation by standard sirens of DECIGO, the FoM value is as high as 170.82. If a 10-year-observation is launched, the FoM could reach 569.42. For comparison, the FoM of 38 actual observed $H(z)$ data(OHD) is 9.3. We also investigate the undulant universe, which shows a comparable improvement on the constraint of cosmological parameters. These improvement indicates that the new method has great potential in further cosmological constraints.

In order to explore the generic properties of a backreaction model for explaining the accelerated expansion of the Universe, we exploit two metrics to describe the late time Universe. Since the standard FLRW metric cannot precisely describe the late time Universe on small scales, the template metric with an evolving curvature parameter $\kappa_{\mathcal{D}}(t)$ is employed. However, we doubt the validity of the prescription for $\kappa_{\mathcal{D}}$, which motivates us apply observational Hubble parameter data (OHD) to constrain parameters in dust cosmology. First, for FLRW metric, by getting best-fit constraints of $\Omega^{{\mathcal{D}}_0}_m = 0.25^{+0.03}_{-0.03}$, $n = 0.02^{+0.69}_{-0.66}$, and $H_{\mathcal{D}_0} = 70.54^{+4.24}_{-3.97}\ {\rm km \ s^{-1} \ Mpc^{-1}}$, the evolutions of parameters are explored. Second, in template metric context, by marginalizing over $H_{\mathcal{D}_0}$ as a prior of uniform distribution, we obtain the best-fit values of $n=-1.22^{+0.68}_{-0.41}$ and ${{\Omega}_{m}^{\mathcal{D}_{0}}}=0.12^{+0.04}_{-0.02}$. Moreover, we utilize three different Gaussian priors of $H_{\mathcal{D}_0}$, which result in different best-fits of $n$, but almost the same best-fit value of ${{\Omega}_{m}^{\mathcal{D}_{0}}}\sim0.12$. Also, the absolute constraints without marginalization of parameter are obtained: $n=-1.1^{+0.58}_{-0.50}$ and ${{\Omega}_{m}^{\mathcal{D}_{0}}}=0.13\pm0.03$. With these constraints, the evolutions of the effective deceleration parameter $q^{\mathcal{D}}$ indicate that the backreaction can account for the accelerated expansion of the Universe without involving extra dark energy component in the scaling solution context. Nevertheless, the results also verify that the prescription of $\kappa_{\mathcal{D}}$ is insufficient and should be improved.

In this paper, we present a new method of measuring Hubble parameter($H(z)$), making use of the anisotropy of luminosity distance($d_{L}$), and the analysis of gravitational wave(GW) of neutron star(NS) binary system. The method has never been put into practice before due to the lack of the ability of detecting GW. LIGO's success in detecting GW of black hole(BH) binary system merger announced the possibility of this new method. We apply this method to several GW detecting projects, including Advanced LIGO(Adv-LIGO), Einstein Telescope(ET) and DECIGO, finding that the $H(z)$ by Adv-LIGO and ET is of bad accuracy, while the $H(z)$ by DECIGO shows a good accuracy. We use the error information of $H(z)$ by DECIGO to simulate $H(z)$ data at every 0.1 redshift span, and put the mock data into the forecasting of cosmological parameters. Compared with the available 38 observed $H(z)$ data(OHD), mock data shows an obviously tighter constraint on cosmological parameters, and a concomitantly higher value of Figure of Merit(FoM). For a 3-year-observation by standard sirens of DECIGO, the FoM value is as high as 834.9. If a 10-year-observation is launched, the FoM could reach 2783.1. For comparison, the FoM of 38 actual observed $H(z)$ data is 9.3. These improvement indicates that the new method has great potential in further cosmological constraints.

The universe, with large-scale homogeneity, is locally inhomogeneous, clustering into stars, galaxies and larger structures. Such property is described by the smoothness parameter $\alpha$ which is defined as the proportion of matter in the form of intergalactic medium. If we take consideration of the inhomogeneities in small scale, there should be modifications of the cosmological distances compared to a homogenous model. Dyer and Roeder developed a second-order ordinary differential equation (D-R equation) that describes the angular diameter distance-redshift relation for inhomogeneous cosmological models. Furthermore, we may obtain the D-R equation for observational $H(z)$ data (OHD). The density-parameter $\Omega_{\rm M}$, the state of dark energy $\omega$, and the smoothness-parameter $\alpha$ are constrained by a set of OHD in a spatially flat $\Lambda$CDM universe as well as a spatially flat XCDM universe. By using of $\chi^2$ minimization method we get $\alpha=0.81^{+0.19}_{-0.20}$ and $\Omega_{\rm M}=0.32^{+0.12}_{-0.06}$ at $1\sigma$ confidence level. If we assume a Gaussian prior of $\Omega_{\rm M}=0.26\pm0.1$, we get $\alpha=0.93^{+0.07}_{-0.19}$ and $\Omega_{\rm M}=0.31^{+0.06}_{-0.05}$. For XCDM model, $\alpha$ is constrained to $\alpha\geq0.80$ but $\omega$ is weakly constrained around -1, where $\omega$ describes the equation of the state of the dark energy ($p_{\rm X}=\omega\rho_{\rm X}$). We conclude that OHD constrains the smoothness parameter more effectively than the data of SNe Ia and compact radio sources.

We use the redshift Hubble parameter $H(z)$ data derived from relative galaxy ages, distant type Ia supernovae (SNe Ia), the Baryonic Acoustic Oscillation (BAO) peak, and the Cosmic Microwave Background (CMB) shift parameter data, to constrain cosmological parameters in the Undulant Universe. We marginalize the likelihood functions over $h$ by integrating the probability density $P\propto e^{-\chi^2/2}$. By using the Markov Chain Monte Carlo (MCMC) technique, we obtain the best fitting results and give the confidence regions on the $b-\Omega_{\rm m0}$ plane. Then we compare their constraints. Our results show that the $H(z)$ data play a similar role with the SNe Ia data in cosmological study. By presenting the independent and joint constraints, we find that the BAO and CMB data play very important roles in breaking the degeneracy compared with the $H(z)$ and SNe Ia data alone. Combined with the BAO or CMB data, one can improve the constraints remarkably. The SNe Ia data sets constrain $\Omega_{\rm m0}$ much tighter than the $H(z)$ data sets, but the $H(z)$ data sets constrain $b$ much tighter than the SNe Ia data sets. All these results show that the Undulant Universe approaches the $\Lambda \rm$CDM model. We expect more $H(z)$ data to constrain cosmological parameters in future.

We use the Radial Baryon Acoustic Oscillation (RBAO) measurements, distant type Ia supernovae (SNe Ia), the observational $H(z)$ data (OHD) and the Cosmic Microwave Background (CMB) shift parameter data to constrain cosmological parameters of $\Lambda$CDM and XCDM cosmologies and further examine the role of OHD and SNe Ia data in cosmological constraints. We marginalize the likelihood function over $h$ by integrating the probability density $P\propto e^{-\chi^{2}/2}$ to obtain the best fitting results and the confidence regions in the $\Omega_{m}-\Omega_{\Lambda}$ plane.With the combination analysis for both of the {\rm $\Lambda$}CDM and XCDM models, we find that the confidence regions of 68.3%, 95.4% and 99.7% levels using OHD+RBAO+CMB data are in good agreement with that of SNe Ia+RBAO+CMB data which is consistent with the result of Lin et al's work. With more data of OHD, we can probably constrain the cosmological parameters using OHD data instead of SNe Ia data in the future.

It has been revealed recently that, in the scale free range, i.e. from the scale of the onset of nonlinear evolution to the scale of dissipation, the velocity and mass density fields of cosmic baryon fluid are extremely well described by the self-similar log-Poisson hierarchy. As a consequence of this evolution, the relations among various physical quantities of cosmic baryon fluid should be scale invariant, if the physical quantities are measured in cells on scales larger than the dissipation scale, regardless the baryon fluid is in virialized dark halo, or in pre-virialized state. We examine this property with the relation between the Compton parameter of the thermal Sunyaev-Zel'dovich effect, $y(r)$, and X-ray luminosity, $L_{\rm x}(r)$, where $r$ being the scale of regions in which $y$ and $L_{\rm x}$ are measured. According to the self-similar hierarchical scenario of nonlinear evolution, one should expect that 1.) in the $y(r)$-$L_x(r)$ relation, $y(r)=10^{A(r)}[L_{\rm x}(r)]^{\alpha(r)}$, the coefficients $A(r)$ and $\alpha(r)$ are scale-invariant; 2.) The relation $y(r)=10^{A(r)}[L_{\rm x}(r)]^{\alpha(r)}$ given by cells containing collapsed objects is also available for cells without collapsed objects, only if $r$ is larger than the dissipation scale. These two predictions are well established with a scale decomposition analysis of observed data, and a comparison of observed $y(r)$-$L_x(r)$ relation with hydrodynamic simulation samples. The implication of this result on the characteristic scales of non-gravitational heating is also addressed.

In this work, we use observations of the Hubble parameter from the differential ages of passively evolving galaxies and the recent detection of the Baryon Acoustic Oscillations (BAO) at $z_1=0.35$ to constrain the Dvali-Gabadadze-Porrati (DGP) universe. For the case with a curvature term, we set a prior $h=0.73\pm0.03$ and the best-fit values suggest a spatially closed Universe. For a flat Universe, we set $h$ free and we get consistent results with other recent analyses.