By using a Markov Chain Monte Carlo simulation, we investigate cosmological constraints on the ghost dark energy (GDE) model in the framework of the Brans-Dicke (BD) theory. A combination of the latest observational data of the cosmic microwave background radiation data from seven-year WMAP, the baryon acoustic oscillation data form the SDSS, the supernovae type Ia data from the Union2 and the X-ray gas mass fraction data from the Chandra X-ray observations of the largest relaxed galaxy clusters are used to perform constraints on GDE in the BD cosmology. In this paper, we consider both flat and non-flat universes together with interaction between dark matter and dark energy. The main cosmological parameters are obtained as: $\Omega_{\rm b}h^2= 0.0223^{+0.0016}_{-0.0013}$, $\Omega_{\rm c}h^2=0.1149^{+0.0088}_{-0.0104}$ and $\Omega_{\rm k}=0.0005^{+0.0025}_{-0.0073}$. In addition, the Brans-Dicke parameter $\omega$ is estimated as $1/\omega\simeq 0.002$.

In this paper we consider quintessence reconstruction of interacting holographic dark energy in a non-flat background. As system's IR cutoff we choose the radius of the event horizon measured on the sphere of the horizon, defined as $L=ar(t)$. To this end we construct a quintessence model by a real, single scalar field. Evolution of the potential, $V(\phi)$, as well as the dynamics of the scalar field, $\phi$, are obtained according to the respective holographic dark energy. The reconstructed potentials show a cosmological constant behavior for the present time. We constrain the model parameters in a flat universe by using the observational data, and applying the Monte Carlo Markov chain simulation. We obtain the best fit values of the holographic dark energy model and the interacting parameters as $c=1.0576^{+0.3010+0.3052}_{-0.6632-0.6632}$ and $\zeta=0.2433^{+0.6373+0.6373}_{-0.2251-0.2251}$, respectively. From the data fitting results we also find that the model can cross the phantom line in the present universe where the best fit value of of the dark energy equation of state is $w_D=-1.2429$.

We constrain holographic dark energy (HDE) with time varying gravitational coupling constant in the framework of the modified Friedmann equations using cosmological data from type Ia supernovae, baryon acoustic oscillations, cosmic microwave background radiation and X-ray gas mass fraction. Applying a Markov Chain Monte Carlo (MCMC) simulation, we obtain the best fit values of the model and cosmological parameters within $1\sigma$ confidence level (CL) in a flat universe as: $\Omega_{\rm b}h^2=0.0222^{+0.0018}_{-0.0013}$, $\Omega_{\rm c}h^2 =0.1121^{+0.0110}_{-0.0079}$, $\alpha_{\rm G}\equiv \dot{G}/(HG) =0.1647^{+0.3547}_{-0.2971}$ and the HDE constant $c=0.9322^{+0.4569}_{-0.5447}$. Using the best fit values, the equation of state of the dark component at the present time $w_{\rm d0}$ at $1\sigma$ CL can cross the phantom boundary $w=-1$.

We consider the effect of a logarithmic f(R) theory, motivated by the form of the one-loop effective action arising from gluons in curved spacetime, on the structure of relativistic stars. In addition to analysing the consistency constraints on the potential of the scalar degree of freedom, we discuss the possibility of observational features arising from a fifth force in the vicinity of the neutron star surface. We find that the model exhibits a chameleon effect that completely suppresses the effect of the modification on scales exceeding a few radii, but close to the surface of the neutron star, the deviation from General Relativity can significantly affect the surface redshift that determines the shift in absorption (or emission) lines. We also use the method of perturbative constraints to solve the modified Tolman-Oppenheimer-Volkov equations for normal and self-bound neutron stars (quark stars).