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.

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.

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.

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.

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.