Using the mathematical definitions of deceleration and jerk parameters we obtain a general differential equation for squared Hubble parameter. For a constant jerk, this differential equation leads to an exact function for Hubble parameter. By the aid of this exact Hubble function we can exactly reconstruct any other cosmographic parameters. We also obtained a general function for transition redshift as well as spacetime curvature. Our derived functions clearly impose a lower limit on the jerk parameter which is $j_{min}\geq-0.125$. Moreover, we found that the jerk parameter indicates the geometry of the spacetime i.e any deviation from $j=1$ imply to a non-flat spacetime. In other word $j\neq 1$ reefers to a dynamical, time varying, dark energy. From obtained Hubble function we recover the analogue of $\Lambda$CDM model. To constrain cosmographic parameters as well as transition redshift and spacetime curvature of the recovered $\Lambda$CDM model, we used Metropolis-Hasting algorithm to perform Monte Carlo Markov Chain analysis by using observational Hubble data obtained from cosmic chronometric (CC) technique, BAO data, Pantheon compilation of Supernovae type Ia, and their joint combination. The only free parameters are $H$, $A(\Omega_{m})$ and $j$. From joint analysis we obtained $H_{0}=69.9\pm 1.7$, $A(\sim\Omega_{0m})=0.279^{+0.013}_{-0.017}$, $B(\sim\Omega_{0X})=0.721^{+0.017}_{-0.013}$, $j_{0}=1.038^{+0.061}_{-0.023}$ and $z_{t}=0.706^{+0.031}_{-0.034}$.

We use recent 36 observational Hubble data (OHD) in the redshift range $0.07\leq z\leq 2.36$, latest \textgravedbl joint light curves\textacutedbl (JLA) sample, comprised of 740 type Ia supernovae (SNIa) in the redshift range $0.01\leq z \leq 1.30$, and their joint combination datasets to constrain anisotropic Bianchi type I (BI) dark energy (DE) model. To estimate model parameters, we apply Hamiltonian Monte Carlo technique. We also compute the covariance matrix for BI dark energy model by considering different datasets to compare the correlation between model parameters. To check the acceptability of our fittings, all results are compared with those obtained from 9 year WMAP as well as Planck (2015) collaboration. Our estimations show that at 68\% confidence level the dark energy equation of state (EOS) parameter for OHD or JLA datasets alone varies between quintessence and phantom regions whereas for OHD+JLA dataset this parameter only varies in phantom region. It is also found that the current cosmic anisotropy is of order $\sim10^{-3}$ which imply that the OHD and JLA datasets do not put tight constraint on this parameter. Therefore, to constraint anisotropy parameter, it is necessary to use high redshif dataset namely cosmic microwave background (CMB). Moreover, from the calculation of $p$-value associated with $\chi^{2}$ statistic we observed that non of the $\omega \mbox{BI}$ and flat $\omega\mbox{CDM}$ models rule out by OHD or JLA datasets. The deceleration parameter is obtained as $q=-0.46^{+0.89 +0.36}_{-0.41 -0.37}$, $q=-0.619^{+0.12 +0.20}_{-0.095 -0.24}$, and $q=-0.52^{+0.080 +0.014}_{-0.046 -0.15}$ for OHD, SNIa, and OHD+SNIa data respectively.

The present study reveals observational constraints on the coupling between dark components of anisotropic Bianchi type I universe. We assume interaction between dark matter and dark energy and split the continuity equation with inclusion of interaction term $\Gamma$. Two scenarios have been considered (i) when coupling between dark components is constant and (ii) when it is a function of redshift ($z$). Metropolis-Hasting algorithm has been used to perform Monte Carlo Markov Chain (MCMC) analysis by using observational Hubble data obtained from cosmic chronometric (CC) technique, cosmic microwave background (CMB) baryon acoustic oscillation (BAO), Pantheon compilation of Supernovae type Ia (SNIa), their joint combination and a gaussian prior on the Hubble parameter $H_{0}$. It is obtained that the combination of all databases plus $H_{0}$ prior marginalized over a present dark energy density gives stringent constraints on the current value of coupling as $-0.001<\delta<0.041$ in constant coupling model and $-0.042<\delta<0.053$ in varying coupling model at 68\% confident level. In general, for both models, we found $\omega^{X}\approx -1$ and $\delta(\delta_{0})\approx 0$ which indicate that still recent data favor uncoupled $\Lambda$CDM model. Our estimations show that in constant coupling model $H_{0}=73.9^{+1.5}_{-0.95}$ which naturally solves the Hubble tension problem. In other word, in constant coupling model, we did not find any disagreement between our estimated $H_{0}$ and those reported by Hubble space telescope (HST) and large scale structure (LSS) experiments.

We use Gaussian processes in combination with MCMC method to place constraints on cosmological parameters of three dark energy models including flat and curved FRW and Bianchi type I spacetimes. To do so, we use recently compiled 36 measurements of the Hubble parameter $H(z)$ in the redshifts intermediate $0.07\leqslant z \leqslant 2.36$. Moreover, we use these models to estimate the redshift of the deceleration-acceleration transition. We consider two Gaussian priors for current value of the Hubble constant i.e $H_{0}=73\pm1.74 (68\pm 2.8)$ km/s/Mpc to investigate the effect of the assumed $H_{0}$ on our parameters estimations. For statistical analysis we use NUTS sampler which is an extension of Hamiltonian Monte Carlo algorithm to generate MCMC chains for parameters of dark energy models. To compare the considered cosmologies, we perform Akaike information criterion (AIC) and Bayes factor ($\Psi$). In general, when we compared our results with 9 years WMAP as well as Planck 2015 Collaboration, we found that Bianchi type I model is slightly fits better to the observational Hubble data with respect to the non-flat FRW model.

It is well known that flat FRW metric is a special case of Bianchi type I spacetime. In this paper, we use 38 Hubble parameter, $H(z)$, measurements at intermediate redshifts $0.07\leq z\leq 2.36$ and its joint combination with the latest \textgravedbl joint light curves\textacutedbl (JLA) sample, comprised of 740 type Ia supernovae in the redshift range $z \epsilon [0.01, 1.30]$ to constrain the parameters of Bianchi type I dark energy model. We also use the same datasets to constrain flat $\Lambda$CDM Model. In both cases, we specifically address the expansion rate $H_{0}$ as well as the transition redshift $z_{t}$ determinations out of these measurements. In both Models we found that using joint combination of datasets gives rise to lower values for model parameters. Also to compare the considered cosmologies, we have made Akaike information criterion (AIC) and Bayes factor ($\Psi$) tests.

We study the constraints on dark energy equation of state $\omega^{X}$ and the deceleration parameter $q$ from the recent observational data including Hubble data and the cosmic microwave background (CMB) radiation by using a model-independent deceleration parameter $q(z)=1/2-a/(1+z)^b$ and dark energy equation of state $\omega^{X}=\omega_{0}+\omega_{1}z/(1+z)$ in the scope of anisotropic Bianchi type I space-time. For the cases of Hubble dataset, CMB data, and their combination, our results indicate that the constraints on transition redshift $z_{\ast}$ are $0.62^{+1.45}_{-0.56}$, $0.34^{+0.13}_{-0.06}$, and $0.60^{+0.20}_{-0.10}$ respectively.

We study the interaction between dark energy (DE) and dark matter (DM) in the scope of anisotropic bianchi type I space-time. First we derive the general form of the dark energy equation of state parameter (EoS) in both non-interacting and interacting cases and then we examine it's future by applying a hyperbolic scale factor. It is shown that in non-interacting case, depending on the value of the anisotropy parameter $K$, the dark energy EoS parameter is varying from phantom to quintessence whereas in interacting case EoS parameter vary in quintessence region. However, in both cases the dark energy EoS parameter $\omega^{de}$, ultimately (i. e at $z=-1$) tends to the cosmological constant ($\omega^{de}=-1$). Moreover, we fixed the cosmological bound on the anisotropy parameter $K$ by using the recent observational data of Hubble parameter.

In this paper we study the evolution of the equation of state of viscous dark energy in the scope of Bianchi type III space-time. We consider the case when the dark energy is minimally coupled to the perfect fluid as well as direct interaction with it. The viscosity and the interaction between the two fluids are parameterized by constants $\zeta_{0}$ and $\sigma$ respectively. We have made a detailed investigation on the cosmological implications of this parametrization. To differentiate between different dark energy models, we have performed a geometrical diagnostic by using the statefinder pair $\{s, r\}$.

In this paper we have investigated the general form of viscous and non-viscous dark energy equation of state (EoS) parameter in the scope of anisotropic Bianchi type I space-time. We show that the presence of bulk viscosity causes transition of $\omega^{de}$ from quintessence to phantom but the phantom state is an unstable state (as expected) and EoS of DE tends to $-1$ at late time. Then we show this phantomic description of the viscous dark energy and reconstruct the potential of the phantom scalar field. It is found that bulk viscosity pushes the universe to a darker region. We have also shown that at late time $q\sim-\Omega^{de}$.

A solution to the coincidence and Big Rip problems on the bases of an anisotropic space-time is proposed. To do so, we study the interaction between viscous dark energy and dark matter in the scope of the Bianchi type-I Universe. We parameterize the viscosity and the interaction between the two fluids by constants $\zeta_{0}$ and $\sigma$ respectively. A detailed investigation on the cosmological implications of this parametrization has been made. We have also performed a geometrical diagnostic by using the statefinder pairs $\{s, r\}$ and $\{q, r\}$ in order to differentiate between different dark energy models. Moreover, we fit the coupling parameter $\sigma$ as well as the Hubble's parameter $H_{0}$ of our model by minimizing the $\chi^{2}$ through the age differential method, involving a direct measurement of $H$.

Phantom dark energy is a proposal that explains the current observations that mildly favor the equation of state of dark energy $\omega^{de}$ crossing -1 at 68% confidence level. However, phantom fields are generally ruled out by ultraviolet quantum instabilities. To overcome this discrepancy, in this paper we propose a mechanism to show that how the presence of bulk viscosity in the cosmic fluid can temporarily drive the fluid into the phantom region ($\omega<-1$). As time is going on, phantom decays and ultimately $\omega^{de}$ approaches to -1. Then we show these quintessence and phantom descriptions of non-viscous and viscous dark energy and reconstruct the potential of these two scalar fields. Also a diagnostic for these models are performed by using the statefinder pairs $\{s,r\}$. All results are obtained in an anisotropic space-time which is a generalization of FLRW universe.

Some new exact solutions of Einstein's field equations have come forth within the scope of a spatially homogeneous and anisotropic Bianchi type-III space-time filled with barotropic fluid and dark energy by considering a variable deceleration parameter. We consider the case when the dark energy is minimally coupled to the perfect fluid as well as direct interaction with it. Under the suitable condition, the anisotropic models approach to isotropic scenario. We also find that during the evolution of the universe, the equation of state (EoS) for dark energy $\omega^{(de)}$, in both cases, tends to -1 (cosmological constant, $\omega^{(de)} = -1$), by displaying various patterns as time increases, which is consistent with recent observations. The cosmic jerk parameter in our derived models are in good agreement with the recent data of astrophysical observations under appropriate condition. It is observed that the universe starts from an asymptotic Einstein static era and reaches to the $\Lambda$CDM model. So from recently developed Statefinder parameters, the behavior of different stages of the universe has been studied. The physical and geometric properties of cosmological models are also discussed.

A model of a cloud formed by massive strings is used as a source of Bianchi type II. We assumed that the expansion $(\theta)$ in the model is proportional to the shear $(\sigma)$. To get exact solution, we have considered the equation of state of the fluid to be in the stiff form. It is found that the bulk viscosity plays a very important rule in the history of the universe. In presence of bulk viscosity the particles dominate over strings whereas in absence of it, strings dominate over the particles which is not in consistence with the recent observations. Also we observe that the viscosity caused the expansion of the universe to be accelerating. Our models are evolving from an early decelerating phase to a late time accelerating phase. The physical and geometrical behavior of these models are discussed.

We study the evolution of the dark energy parameter within the scope of a spatially non-flat and isotropic Friedmann-Robertson-Walker (FRW) model filled with barotropic fluid and bulk viscous stresses. We have obtained cosmological solutions which exhibit without a big rip singularity. It is concluded that in both non-interacting and interacting cases non-flat open universe crosses the phantom region. We find that during the evolution of the universe, the equation of state (EoS) for dark energy $\omega_{D}$ changes from $\omega^{eff}_{D} < -1$ to $\omega^{eff}_{D} > -1$, which is consistent with recent observations.

A model of a cloud formed by massive strings is used as a source of LRS Bianchi type II with time decaying vacuum energy density $\Lambda$. To construct string cosmological models we have used the energy-momentum tensor for such string as formulated by Letelier (1983). The high nonlinear field equations have been solved for two types of strings, (i) Massive string and (ii) Nambu string. The expansion $\theta$ in the model is assumed to be proportional to the shear $\sigma$. This condition leads to $A=\beta B^{m}$, where $A$ and $B$ are the metric coefficients, $m$ is a constant and $\beta$ is an integrating constant. Our models are in accelerating phase which is consistent to the recent observations of supernovae type Ia. The physical and geometrical behavior of these models are also discussed.