
A generalized BransDicke (GBD) theory is developed in this paper. The GBD
theory is obtained with replacing the Ricci scalar $R$ in the original
BransDicke (BD) action by an arbitrary function $f(R)$. In this theory, the
gravitational field equation and the motion equation of BD scalar field are
given. And then, we study the cosmological equations and their solutions. The
solutions show that they are consistent with the observational Hubble data. We
probe the effects of the BransDicke scalar field in the GBD, it is shown that
the modification to the cosmological quantities from the dynamical BD scalar
field is remarkable. In addition, the forms of the effective potential of
BransDicke field are investigated in the GBD model.

Tensions between several cosmic observations were found recently, such as the
inconsistent values of $H_{0}$ (or $\sigma_{8}$) were indicated by the
different cosmic observations. Introducing the massive neutrinos in
$\Lambda$CDM could potentially solve the tensions. Viable $f(R)$ gravity
producing $\Lambda$CDM background expansion with massive neutrinos is
investigated in this paper. We fit the current observational data: Planck2015
CMB, RSD, BAO and SNIa to constrain the mass of neutrinos in viable $f(R)$
theory. The constraint results at 95\% confidence level are: $\Sigma
m_\nu<0.202$ eV for the active neutrino case, $m_{\nu, sterile}^{eff}<0.757$ eV
with $N_{eff}<3.22$ for the sterile neutrino case. For the effects by the mass
of neutrinos, the constraint results on model parameter at 95\% confidence
level become $f_{R0}\times 10^{6}> 1.89$ and $f_{R0}\times 10^{6}> 2.02$
for two cases, respectively. It is also shown that the fitting values of
several parameters much depend on the neutrino properties, such as the cold
dark matter density, the cosmological quantities at matterradiation equality,
the neutrino density and the fraction of baryonic mass in helium. At last, the
constraint result shows that the tension between direct and CMB measurements of
$H_0$ gets slightly weaker in the viable $f(R)$ model than that in the base
$\Lambda$CDM model.

The coupled dark energy model provides a possible approach to mitigate the
coincidence problem of cosmological standard model. Here, the coupling term is
assumed as $\bar{Q}=3H\xi_x\bar{\rho}_x$, which is related to the interaction
rate and energy density of dark energy. We derive the background and
perturbation evolution equations for several coupled models. Then, we test
these models by currently available cosmic observations which include cosmic
microwave background radiation from Planck 2015, baryon acoustic oscillation,
type Ia supernovae, $f\sigma_8(z)$ data points from redshiftspace distortions,
and weak gravitational lensing. The constraint results tell us there is no
evidence of interaction at 2$\sigma$ level, it is very hard to distinguish
different coupled models from other ones.

Observations indicate that most of the universal matter are invisible and the
gravitational constant $G(t)$ maybe depends on the time. A theory of the
variational $G$ (VG) is explored in this paper, with naturally producing the
useful dark components in universe. We utilize the observational data: lookback
time data, modelindependent gamma ray bursts, growth function of matter linear
perturbations, type Ia supernovae data with systematic errors, CMB and BAO to
restrict the unified model (UM) of dark components in VG theory. Using the
bestfit values of parameters with the covariance matrix, constraints on the
variation of $G$ are $(\frac{G}{G_{0}})_{z=3.5}\simeq
1.0015^{+0.0071}_{0.0075}$ and $(\frac{\dot{G}}{G})_{today}\simeq
0.7252^{+2.3645}_{2.3645}\times 10^{13} yr^{1}$, the small uncertainties
around constants. Limit on the equation of state of dark matter is
$w_{0dm}=0.0072^{+0.0170}_{0.0170}$ with assuming $w_{0de}=1$ in unified
model, and dark energy is $w_{0de}=0.9986^{+0.0011}_{0.0011}$ with assuming
$w_{0dm}=0$ at prior. Restriction on UM parameters are
$B_{s}=0.7442^{+0.0137+0.0262}_{0.01320.0292}$ and
$\alpha=0.0002^{+0.0206+0.0441}_{0.02090.0422}$ with $1\sigma$ and $2\sigma$
confidence level. In addition, the effect of a cosmic string fluid on unified
model in VG theory are investigated. In this case it is found that the
$\Lambda$CDM ($\Omega_{s}=0$, $\beta=0$ and $\alpha=0$) is included in this
VGUM model at $1\sigma$ confidence level, and the larger errors are given:
$\Omega_{s}=0.0106^{+0.0312+0.0582}_{0.03050.0509}$ (dimensionless energy
density of cosmic string), $(\frac{G}{G_{0}})_{z=3.5}\simeq
1.0008^{+0.0620}_{0.0584}$ and $(\frac{\dot{G}}{G})_{today}\simeq
0.3496^{+26.3135}_{26.3135}\times 10^{13}yr^{1}$.

In this paper, we study cosmologies containing the reduced modified Chaplygin
gas (RMCG) fluid which is reduced from the modified Chaplygin gas
$p=A\rhoB\rho^{\alpha}$ for the value of $\alpha=1/2$. In this special case,
dark cosmological models can be realized for different values of model
parameter $A$. We investigate the viabilities of these dark cosmological models
by discussing the evolutions of cosmological quantities and using the currently
available cosmic observations. It is shown that the special RMCG model ($A=0$
or $A=1$) which unifies the dark matter and dark energy should be abandoned.
For $A=1/3$, RMCG which unifies the dark energy and dark radiation is the
favorite model according to the objective Akaike information criteria. In the
case of $A<0$, RMCG can achieve the features of the dynamical quintessence and
phantom models, where the evolution of the universe is not sensitive to the
variation of model parameters.

In this paper, a unified dark fluid with constant adiabatic sound speed is
decomposed into cold dark matter interacting with vacuum energy. Based on
Markov chain Monte Carlo method, we constrain this model by jointing the
geometry and dynamical measurement. The geometry test includes cosmic microwave
background radiation from \textit{Planck}, baryon acoustic oscillation, and
type Ia supernovae; the dynamic measurement is $f\sigma_8(z)$ data points which
is obtained from the growth rate via redshiftspace distortion, and
$\sigma_8(z)$ is the rootmeansquare amplitude of the density contrast
$\delta$ at the comoving $8h^{1}$ Mpc scale. The jointed constraint shows that
$\alpha$ = $ 0.000662_{ 0.000662}^{+ 0.000173}$ and $\sigma_8$ = $ 0.824_{
0.0166}^{+ 0.0128}$. The CMB and matter power spectra are both similar for the
case of $\alpha=$ mean value and that of $\alpha=0$. However, the evolutionary
curves of $f\sigma_8(z)$ are different. This means that, to some extent, the
data points of the growth rate could break the degeneracy of the dark energy
models.

Motivated by the work [K. Karami, J. Fehri, {{\it Phys. Lett. B}} {\bf 684},
61 (2010)] and [A. Sheykhi, {{\it Phys. Lett. B}} {\bf 681}, 205 (2009)], we
generalize their work to the new holographic dark energy model with
$\rho_D=\frac{3\phi^2}{4\omega}(\mu H^2+\nu\dot{H})$ in the framework of
BransDicke cosmology. Concretely, we study the correspondence between the
quintessence, tachyon, Kessence, dilaton scalar field and Chaplygin gas model
with the new holographic dark energy model in the nonflat BransDicke
universe. Furthermore, we reconstruct the potentials and dynamics for these
models. By analysis we can show that for new holographic quintessence and
Chaplygin gas models, if the related parameters to the potentials satisfy some
constraints, the accelerated expansion can be achieved in BransDicke
cosmology. Especially the counterparts of fields and potentials in general
relativity can describe accelerated expansion of the universe. It is worth
stressing that not only can we give some new results in the framework of
BransDicke cosmology, but also the previous results of the new holographic
dark energy in Einstein gravity can be included as special cases given by us.

Motivated by recent works [1,2], the influences of selfgravitational
interaction on the Hawking radiation are studied both for (3 + 1) and (4 + 1)
dimensional black holes in $z = 4$ Ho\u{r}avaLifshitz gravity. It is found
that the tunneling entropies $S_{B(3 + 1)}$ and $S_{B(4 + 1)}$ independent on
particle's mass are consistent with the thermodynamical entropies $S_{BH(3 +
1)}$ and $S_{BH(4 + 1)}$, respectively. There is a very visible degree of
uniformity between thermodynamics and quantum tunneling in $z = 4$
Ho\u{r}avaLifshitz gravity. It suggests that the entropies contained corrected
terms could be explained well by the selfgravitational interaction of Hawking
radiation. The study of tunneling process may shed light on understanding the
Ho\u{r}avaLifshitz gravity.

In this paper we investigate the interaction between dark matter and dark
energy by considering two different interacting scenarios, i.e. the cases of
constant interaction function and variable interaction function. By fitting the
current observational data to constrain the interacting models, it is found
that the interacting strength is nonvanishing, but weak for the case of
constant interaction function, and the interaction is not obvious for the case
of variable interaction function. In addition, for seeing the influence from
interaction we also investigate the evolutions of interaction function,
effective state parameter for dark energy and energy density of dark matter. At
last some geometrical quantities in the interacting scenarios are discussed.

Time variable cosmological constant (TVCC) of holographic origin with
interaction in BransDicke theory is discussed in this paper. We investigate
some characters for this model, and show the evolutions of deceleration
parameter and equation of state (EOS) for dark energy. It is shown that in this
scenario an accelerating universe can be obtained and the evolution of EOS for
dark energy can cross over the boundary of phantom divide. In addition, a
geometrical diagnostic method, jerk parameter is applied to this model to
distinguish it with cosmological constant.

We use the Markov ChainMonte Carlo method to investigate global constraints
on the generalized holographic (GH) dark energy with flat and nonflat universe
from the current observed data: the Union2 dataset of type supernovae Ia
(SNIa), highredshift GammaRay Bursts (GRBs), the observational Hubble data
(OHD), the cluster Xray gas mass fraction, the baryon acoustic oscillation
(BAO), and the cosmic microwave background (CMB) data. The most stringent
constraints on the GH model parameter are obtained. In addition, it is found
that the equation of state for this generalized holographic dark energy can
cross over the phantom boundary wde =1.

The evolution of universe in BransDicke (BD) theory is discussed in this
paper.
Considering a parameterized scenario for BD scalar field
$\phi=\phi_{0}a^{\alpha}$ which plays the role of gravitational "constant" $G$,
we apply the Markov Chain Monte Carlo method to investigate a global
constraints on BD theory with a selfinteracting potential according to the
current observational data: Union2 dataset of type supernovae Ia (SNIa),
highredshift GammaRay Bursts (GRBs) data, observational Hubble data (OHD),
the cluster Xray gas mass fraction, the baryon acoustic oscillation (BAO), and
the cosmic microwave background (CMB) data. It is shown that an expanded
universe from deceleration to acceleration is given in this theory, and the
constraint results of dimensionless matter density $\Omega_{0m}$ and parameter
$\alpha$ are, $\Omega_{0m}=0.286^{+0.037+0.050}_{0.0390.047}$ and
$\alpha=0.0046^{+0.0149+0.0171}_{0.01710.0206}$ which is consistent with the
result of current experiment exploration, $\mid\alpha\mid \leq 0.132124$. In
addition, we use the geometrical diagnostic method, jerk parameter $j$, to
distinguish the BD theory and cosmological constant model in Einstein's theory
of general relativity.

We use the Markov Chain Monte Carlo method to investigate a global
constraints on the modified Chaplygin gas (MCG) model as the unification of
dark matter and dark energy from the latest observational data: the Union2
dataset of type supernovae Ia (SNIa), the observational Hubble data (OHD), the
cluster Xray gas mass fraction, the baryon acoustic oscillation (BAO), and the
cosmic microwave background (CMB) data. In a flat universe, the constraint
results for MCG model are,
$\Omega_{b}h^{2}=0.02263^{+0.00184}_{0.00162}$ ($1\sigma$)
$^{+0.00213}_{0.00195}$ $(2\sigma)$,
$B_{s}=0.7788^{+0.0736}_{0.0723}$ ($1\sigma$) $^{+0.0918}_{0.0904}$
$(2\sigma)$,
$\alpha=0.1079^{+0.3397}_{0.2539}$ ($1\sigma$) $^{+0.4678}_{0.2911}$
$(2\sigma)$,
$B=0.00189^{+0.00583}_{0.00756}$ ($1\sigma$) $^{+0.00660}_{0.00915}$
$(2\sigma)$, and $H_{0}=70.711^{+4.188}_{3.142}$ ($1\sigma$)
$^{+5.281}_{4.149}$ $(2\sigma)$.