We investigate observational constraints on the running vacuum model (RVM) of $\Lambda=3\nu (H^{2}+K/a^2)+c_0$ in the spatially curved universe, where $\nu$ is the model parameter, $K$ corresponds to the spatial curvature constant, and $c_{0}$ is a constant defined by the boundary conditions. As $\dot{\Lambda} \ne 0$, there are energy exchanges between vacuum, matter and radiation in RVM. We study the "geometrical degeneracy" of RVM on the CMB power spectra. By fitting the cosmological data, we find that the values of $\chi^2$ in RVM and $\Lambda$CDM are similar to each other for the non-flat universe. Explicitly, we obtain the constraints of $\nu\leq O(10^{-4})$ (68 $\%$ C.L.) and $|\Omega_K|\leq O(10^{-2})$ (95 $\%$ C.L.) in our study. In addition, we show that the cosmological constraints of $\Sigma m_{\nu}=0.416^{+0.311}_{-0.407}$ (RVM) and $\Sigma m_{\nu}=0.497^{+0.335}_{-0.387}$ ($\Lambda$CDM) at 95$\%$ C.L. for the neutrino mass sum are relaxed in both models in the spatially curved universe.

The scalar tensor theory contains a coupling function connecting the quantities in the Jordan and Einstein frames, which is constrained to guarantee a transformation rule between frames. We simulate the supernovae core collapse with different choices of coupling functions defined over the viable region of the parameter space and find that a generic inverse-chirp feature of the gravitational waves in the scalar tensor scenario.

The scalar-tensor theory can be formulated in both Jordan and Einstein frames, which are conformally related together with a redefinition of the scalar field. As the solution to the equation of the scalar field in the Jordan frame does not have the one-to-one correspondence with that in the Einstein frame, we give a criterion along with some specific models to check if the scalar field in the Einstein frame is viable or not by confirming whether this field is reversible back to the Jordan frame. We further show that the criterion in the first parameterized post-Newtonian approximation can be determined by the parameters of the osculating approximation of the coupling function in the Einstein frame and can be treated as a viable constraint on any numerical study in the scalar-tensor scenario. We also demonstrate that the Brans-Dicke theory with an infinite constant parameter $\omega_{\text{BD}}$ is a counterexample of the equivalence between two conformal frames due to the violation of the viable constraint.

In the light of the latest measurements on the total $e^+ + e^-$ flux by CALET and DAMPE experiments, we revisit the multicomponent leptonically decaying dark matter (DM) explanations to the cosmic-ray electron/positron excesses observed previously. Especially, we use the single and double-component DM models to explore the compatibility of the AMS-02 positron fraction with the new CALET or DAMPE data. It turns out that neither single nor double-component DM models are able to fit the AMS-02 positron fraction and DAMPE total $e^+ + e^-$ flux data simultaneously. On the other hand, for the combined AMS-02 and CALET dataset, both the single and double-component DM models can provide reasonable fits. If we further take into the diffuse $\gamma$-ray constraints from Fermi-LAT, only the double-component DM models are allowed.

We investigate the cosmological viability of the generalized proca theory. We first implement the background and linear perturbation equations of motion in the Boltzmann code and then study the constraints on the parameters of the generalized proca theory after running MCMC against the cosmological data set. With Planck + HST data, we obtain the constraint $h=0.7334_{-0.0269}^{+0.0246}$, which indicates that the tension between early universe and late time universe within this theory is removed. By adding other late-time data sets (BAO, RSD, etc.) we show that the tension is reduced, as the 2$\sigma$ allowed region for $h$ in Proca, $h=0.7041_{-0.0087}^{+0.0094}$, overlaps with the 2$\sigma$ region of the HST data.

We study a special running vacuum model (RVM) with $\Lambda = 3 \alpha H^2+3\beta H_0^4 H^{-2}+\Lambda_0$, where $\alpha$, $\beta$ and $\Lambda_0$ are the model parameters and $H$ is the Hubble one. This RVM has non-analytic background solutions for the energy densities of matter and radiation, which can only be evaluated numerically. From the analysis of the CMB power spectrum and baryon acoustic oscillation along with the prior of $\alpha>0$ to avoid having a negative dark energy density, we find that $\alpha<2.83\times 10^{-4}$ and $\beta=(-0.2^{+3.9}_{-4.5})\times 10^{-4}$ (95$\%$ C.L.). We show that the RVM fits the cosmological data comparably to the $\Lambda$CDM. In addition, we relate the fluctuation amplitude $\sigma_8$ to the neutrino mass sum $\Sigma m_\nu$.

We apply the first law of thermodynamics to the apparent horizon of the universe with the power-law corrected and non-extensive Tsallis entropies rather than the Bekenstein-Hawking one. We examine the cosmological properties in the two entropy models by using the CosmoMC package. In particular, the first numerical study for the cosmological observables with the power-law corrected entropy is performed. We also show that the neutrino mass sum has a non-zero central value with a relaxed upper bound in the Tsallis entropy model comparing with that in the $\Lambda$CDM one.

We study the matter power spectra in the viable $f(R)$ gravity models with the dynamical background evolution and linear perturbation theory by using the CosmoMC package. We show that these viable $f(R)$ models generally shorten the age of the universe and suppress the matter density fluctuation. We examine the allowed ranges of the model parameters and the constraints of the cosmological variables from the current observational data, and find that the dynamical evolution of $\rho_{DE}(z)$ plays an important role to constrain the neutrino masses.

We investigate the power spectra of the CMB temperature and matter density in the running vacuum model (RVM) with the time-dependent cosmological constant of $\Lambda = 3 \nu H^2 + \Lambda_0$, where $H$ is the Hubble parameter. In this model, dark energy decreases in time and decays to both matter and radiation. By using the Markov chain Monte Carlo method, we constrain the model parameter $\nu$ as well as the cosmological observables. Explicitly, we obtain $\nu \leq 1.54 \times 10^{-4}$ (68\% confidence level) in the RVM with the best-fit $\chi^2_{\mathrm{RVM}} = 13968.8$, which is slightly smaller than $\chi^2_{\Lambda \mathrm{CDM}} = 13969.8$ in the $\Lambda$CDM model of $\nu=0$.

We study the $\Lambda(\alpha)$CDM models with $\Lambda(\alpha)$ being a function of the time-varying fine structure constant $\alpha$. We give a close look at the constraints on two specific $\Lambda(\alpha)$CDM models with one and two model parameters, respectively, based on the cosmological observational measurements along with 313 data points for the time-varying $\alpha$. We find that the model parameters are constrained to be around $10^{-4}$, which are similar to the results discussed previously but more accurately.

We examine the running vacuum model with $\Lambda (H) = 3 \nu H^2 + \Lambda_0$, where $\nu$ is the model parameter and $\Lambda_0$ is the cosmological constant. From the data of the cosmic microwave background radiation, weak lensing and baryon acoustic oscillation along with the time dependent Hubble parameter $H(z)$ and weighted linear growth $f (z)\sigma_8(z)$ measurements, we find that $\nu=(1.37^{+0.72}_{-0.95})\times 10^{-4}$ with the best fitted $\chi^2$ value slightly smaller than that in the $\Lambda$CDM model.

We study quintessential inflation using a generalized exponential potential $V(\phi)\propto exp(-\lambda \phi^n/Mpl^n), n>1$, the model admits slow-roll inflation at early times and leads to close-to-scaling behaviour in the post inflationary era with an exit to dark energy at late times. We present detailed investigations of the inflationary stage in the light of the Planck 2015 results, study post-inflationary dynamics and analytically confirm the existence of an approximately scaling solution. Additionally, assuming that standard massive neutrinos are non-minimally coupled, makes the field $\phi$ dominant once again at late times giving rise to present accelerated expansion of the Universe. We derive observational constraints on the field and time-dependent neutrino masses. In particular, for $n=6 (8)$, the parameter $\lambda$ is constrained to be,$\log \lambda > -7.29 (-11.7)$; the model produces the spectral index of the power spectrum of primordial scalar (matter density) perturbations as $ n_s = 0.959 \pm 0.001 (0.961 \pm 0.001)$ and tiny tensor-to-scalar ratio, $r<1.72 \times 10^{-2} (2.32 \times 10^{-2})$ respectively. Consequently, the upper bound on possible values of the sum of neutrino masses $\Sigma m_{\nu} \lesssim 2.5$ eV significantly enhances compared to that in the standard $\Lambda$CDM model.

We investigate observational constraints on a specific one-parameter extension to the minimal quintessence model, where the quintessence field acquires a quadratic coupling to the scalar curvature through a coupling constant $\xi$. The value of $\xi$ is highly suppressed in typical tracker models if the late-time cosmic acceleration is driven at some field values near the Planck scale. We test $\xi$ in a second class of models in which the field value today becomes a free model parameter. We use the combined data from type-Ia supernovae, cosmic microwave background, baryon acoustic oscillations and matter power spectrum, to weak lensing measurements and find a best-fit value $\xi > 0.289$ where $\xi = 0$ is excluded outside the 95 per cent confidence region. The effective gravitational constant $G_{\rm eff}$ subject to the hint of a non-zero $\xi$ is constrained to $-0.003 < 1- G_{\rm eff}/G < 0.033$ at the same confidence level on cosmological scales, and can be narrowed down to $1- G_{\rm eff}/G < 2.2 \times 10^{-5}$ when combining with Solar System tests.

We investigate the matter density perturbation $\delta_m$ and power spectrum $P(k)$ in the running vacuum model (RVM) with the cosmological constant being a function of the Hubble parameter, given by $\Lambda = \Lambda_0 + 6 \sigma H H_0+ 3\nu H^2$, in which the linear and quadratic terms of $H$ would originate from the QCD vacuum condensation and cosmological renormalization group, respectively. Taking the dark energy perturbation into consideration, we derive the evolution equation for $\delta_m$ and find a specific scale $d_{cr}=2 \pi/k_{cr}$, which divides the evolution of the universe into the sub and super-interaction regimes, corresponding to $k \ll k_{cr}$ and $k \gg k_{cr}$, respectively. For the former, the evolution of $\delta_m$ has the same behavior as that in the $\Lambda$CDM model, while for the latter, the growth of $\delta_m$ is frozen (greatly enhanced) when $\nu + \sigma >(<)0$ due to the couplings between radiation, matter and dark energy. It is clear that the observational data rule out the cases with $\nu<0$ and $\nu + \sigma <0$, while the allowed window for the model parameters is extremely narrow with $\nu, |\sigma| \lesssim \mathcal{O}(10^{-7})$.

We investigate the running cosmological constant model with dark energy linearly proportional to the Hubble parameter, $\Lambda = \sigma H + \Lambda_0$, in which the $\Lambda$CDM limit is recovered by taking $\sigma=0$. We derive the linear perturbation equations of gravity under the Friedmann-Lema\"itre-Robertson-Walker cosmology, and show the power spectra of the CMB temperature and matter density distribution. By using the Markov chain Monte Carlo method, we fit the model to the current observational data and find that $\sigma H_0/ \Lambda_0 \lesssim 2.63 \times 10^{-2}$ and $6.74 \times 10^{-2}$ for $\Lambda(t)$ coupled to matter and radiation-matter, respectively, along with constraints on other cosmological parameters.

In order to solve the fine-tuning problem of the cosmological constant, we propose a simple model with the vacuum energy non-minimally coupled to the inflaton field. In this model, the vacuum energy decays to the inflaton during pre-inflation and inflation eras, so that the cosmological constant effectively deflates from the Planck mass scale to a much smaller one after inflation and plays the role of dark energy in the late-time of the universe. We show that our deflationary scenario is applicable to arbitrary slow-roll inflation models. We also take two specific inflation potentials to illustrate our results.

We consider generic models of quintessence and we investigate the influence of massive neutrino matter with field-dependent masses on the matter power spectrum. In case of minimally coupled neutrino matter, we examine the effect in tracker models with inverse power-law and double exponential potentials. We present detailed investigations for the scaling field with a steep exponential potential, non-minimally coupled to massive neutrino matter, and we derive constraints on field-dependent neutrino masses from the observational data.

A long-lived decaying dark matter as a resolution to Fermi, PAMELA and ATIC anomalies is investigated in the framework of split supersymmetry (SUSY) without R-parity, where the neutralino is regarded as the dark matter and the extreme fine-tuned couplings for the long-lived neutralino are naturally evaded in the usual approach.

We investigate the observational constraints on the modified gravity, which combines the $R^{2-q}$ inflation with the power-law (exponential) type of the viable $f(R)$ dark energy models. We discuss the difference between the combined model and $R^{2-q}$ gravity in the inflationary epoch and obtain the constraints on the deviation power $q$ as well as the parameters in $f(R)$ by using the CosmoMC package. The allowed ranges of the spectral index and tensor-to-scalar ratio from the Planck data are highly restricted, resulting in $q < 2.66 \times 10^{-2}$ and $2.17 \times 10^{-2}$ for the power-law and exponential types of $f(R)$ gravity, respectively.

We investigate two classes of models of quintessential inflation, based upon canonical as well as noncanonical scalar fields. In particular, introducing potentials steeper than the standard exponential, we construct models that can give rise to a successful inflationary phase, with signatures consistent with Planck 2015 results. Additionally, using nonminimal coupling of the scalar field with massive neutrino matter, we obtain the standard thermal history of the Universe, with late-time cosmic acceleration as the last stage of evolution. In both cases, inflation and late-time acceleration are connected by a tracker solution.

We revisit the multi-component leptonically decaying dark matter (DM) scenario to explain the possible electron/positron excesses with the recently updated AMS-02 data. We find that both the single- and two-component DM models can fit the positron fraction and $e^+/e^-$ respective fluxes, in which the two-component ones provide better fits. However, for the single-component models, the recent AMS-02 data on the positron fraction limit the DM cutoff to be smaller than 1 TeV, which conflicts with the high-energy behavior of the AMS-02 total $e^++e^-$ flux spectrum, while the two-component DM models do not possess such a problem. We also discuss the constraints from the Fermi-LAT measurement of the diffuse $\gamma$-ray spectrum. We show that the two-component DM models are consistent with the current DM lifetime bounds. In contrast, the best-fit DM lifetimes in the single-component models are actually excluded.

We constrain the scalar field dark energy model with an inverse power-law potential, i.e., $V(\phi)\propto {\phi}^{-\alpha}$ ($\alpha>0$), from a set of recent cosmological observations by compiling an updated sample of Hubble parameter measurements including 30 independent data points. Our results show that the constraining power of the updated sample of $H(z)$ data with the HST prior on $H_0$ is stronger than those of the SCP Union2 and Union2.1 compilations. A recent sample of strong gravitational lensing systems is also adopted to confine the model even though the results are not significant. A joint analysis of the strong gravitational lensing data with the more restrictive updated Hubble parameter measurements and the Type Ia supernovae data from SCP Union2 indicates that the recent observations still can not distinguish whether dark energy is a time-independent cosmological constant or a time-varying dynamical component.

We investigate the viable exponential $f(R)$ gravity in the metric formalism with $f(R)=-\beta R_s (1-e^{-R/R_s})$. The latest sample of the Hubble parameter measurements with 23 data points is used to place bounds on this $f(R)$ model. A joint analysis is also performed with the luminosity distances of Type Ia supernovae and baryon acoustic oscillations in the clustering of galaxies, and the shift parameters from the cosmic microwave background measurements, which leads to $0.240<\Omega_m^0<0.296$ and $\beta>1.47$ at 1$\sigma$ confidence level. The evolutions of the deceleration parameter $q(z)$ and the effective equations of state $\omega_{de}^{eff}(z)$ and $\omega_{tot}^{eff}(z)$ are displayed. By taking the best-fit parameters as prior values, we work out the transition redshift (deceleration/acceleration) $z_T$ to be about 0.77. It turns out that the recent observations are still unable to distinguish the background dynamics in the $\Lambda$CDM and exponential $f(R)$ models.

We investigate the matter power spectra in the power law and exponential types of viable $f(R)$ theories along with massive neutrinos. The enhancement of the matter power spectrum is found to be a generic feature in these models. In particular, we show that in the former type, such as the Starobinsky model, the spectrum is magnified much larger than the latter one, such as the exponential model. A greater scale of the total neutrino mass, $\Sigma m_{\nu}$, is allowed in the viable $f(R)$ models than that in the $\Lambda$CDM one. We obtain the constraints on the neutrino masses by using the CosmoMC package with the modified MGCAMB. Explicitly, we get $\Sigma m_{\nu} < 0.451 \ (0.214)\ \mathrm{eV}$ at 95% C.L. in the Starobinsky (exponential) model, while the corresponding one for the $\Lambda$CDM model is $\Sigma m_{\nu} < 0.200\ \mathrm{eV}$. Furthermore, by treating the effective number of neutrino species $N_{\mathrm{eff}}$ as a free parameter along with $\Sigma m_{\nu}$, we find that $N_{\mathrm{eff}} = 3.78^{+0.64}_{-0.84} (3.47^{+0.74}_{-0.60})$ and $\Sigma m_{\nu} = 0.533^{+0.254}_{-0.411}$ ($< 0.386) \ \mathrm{eV}$ at 95% C.L. in the Starobinsky (exponential) model.

We study teleparallel gravity in the \emph{original} Kaluza-Klein (KK) scenario. Our calculation of the KK reduction of teleparallel gravity indicates that the 5-dimensional torsion scalar $^{(5)}T$ generates the non-Brans-Dicke type effective Lagrangian in 4-dimension due to an additional coupling between the derivative of the scalar field and torsion, but the result is equivalent to that in general relativity. We also discuss the cosmological behavior in the FLRW universe based on the effective teleparallel gravity.

We study a two-component dark matter (DM) model in which the two Majorana fermionic DM components with nearly degenerate masses are stabilized by an $Z_2$ symmetry and interact with the right-handed muon and tau only via real Yukawa couplings, together with an additional $Z_2$-odd singly-charged scalar. In this setup, the decay from the heavy DM to the lighter one via the transition electric dipole yields the 3.55 keV X-ray signal observed recently. The Yukawa couplings in the dark sector are assumed to be hierarchical, so that the observed DM relic abundance can be achieved with the leading s-wave amplitudes without a fine-tuning. We also consider the constraints from flavor physics, DM direct detections and collider searches, respectively.

We present a unified dark fluid model to describe the possible evolutionary behavior of $\Delta N_\mathrm{eff}$ in dark radiation. This model can be viewed as an interacting model for the dark sectors, in which dark matter interacts with dark radiation. We show that the evolution of $\Delta N_\mathrm{eff}$ can be nicely explained without some drawbacks, such as the blowup of $\Delta N_\mathrm{eff}$ and the non-vanishing interaction at the late time.

The multi-component decaying dark matter (DM) scenario is investigated to explain the possible excesses in the positron fraction by PAMELA and recently confirmed by AMS-02, and in the total $e^+ +e^-$ flux observed by Fermi-LAT. By performing the $\chi^2$ fits, we find that two DM components are already enough to give a reasonable fit of both AMS-02 and Fermi-LAT data. The best-fitted results show that the heavier DM component with its mass 1.5 TeV dominantly decays through the $\mu$-channel, while the lighter one of 100 GeV mainly through the $\tau$-channel. As a byproduct, the fine structure around 100 GeV observed by AMS-02 and Fermi-LAT can be naturally explained by the dropping due to the lighter DM component. With the obtained model parameters by the fitting, we calculate the diffuse $\gamma$-ray emission spectrum in this two-component DM scenario, and find that it is consistent with the data measured by Fermi-LAT. We also construct a microscopic particle DM model to naturally realize the two-component DM scenario, and point out an interesting neutrino signal which is possible to be measured in the near future by IceCube.

We study the post-Newtonian limit in the teleparallel equivalent of General Relativity with a scalar field which non-minimally couples to gravity. The metric perturbation is obtained from the vierbein field expansion with respect to the Minkowski background. Due to the structure of the teleparallel gravity Lagrangian, the potential of the scalar field shows no effect to the parametrized post-Newtonian parameters, and compatible results with Solar System observations are found.

Generation of large-scale magnetic fields in inflationary cosmology is studied in teleparallelism, where instead of the scalar curvature in general relativity, the torsion scalar describes the gravity theory. In particular, we investigate a coupling of the electromagnetic field to the torsion scalar during inflation, which leads to the breaking of conformal invariance of the electromagnetic field. We demonstrate that for a power-law type coupling, the current magnetic field strength of $\sim 10^{-9}$ G on 1 Mpc scale can be generated, if the backreaction effects and strong coupling problem are not taken into consideration.

We study future singularity in teleparallel dark energy models, particularly its behavior and its (non)occurrence in the observationally viable models. For the models with a general self-potential of the scalar field, we point out that both at early times and in the future near the singularity the behavior of dark energy can be described by the analytic solutions of the scalar field we obtained for the model with no self-potential. As to the (non)occurrence in the viable models, we consider a natural binding-type self-potential, the quadratic potential, when fitting observational data, and illustrate the constraining region up to the $3\sigma$ confidence level as well as the region where a singularity will occur. As a result, the singularity region is outside the $3\sigma$ constraint. Thus, although the future singularity problem potentially exists in teleparallel dark energy models, the observationally viable models may not suffer this problem.

We study the cosmological perturbations in teleparallel dark energy models in which there is a dynamical scalar field with a non-minimal coupling to gravity. We find that the propagating degrees of freedom are the same as in quintessence cosmology despite that variables of the perturbed vierbein field are greater than those in metric theories. We numerically show some evident discrepancy from general relativity in the evolutions of the perturbations on all scales of the universe. We also demonstrate that the gravitational interactions are enhanced during the unique tracker evolutions in these models.

We analyze the stability of the Einstein static closed and open universe in two types of exponential $f(T)$ gravity theories. We show that the stable solutions exist in these two models. In particular, we find that large regions of parameter space in equation of state $w=p/\rho$ for the stable universe are allowed in the $f(T)$ theories.

We study the cosmological effect of the simple scalar-torsion ($0^+$) mode in Poincar\'{e} gauge theory of gravity. We find that for the non-constant (affine) curvature case, the early evolution of the torsion density $\rho_T$ has a radiation-like asymptotic behavior of $a^{-4}$ with $a$ representing the scale factor, along with the stable point of the torsion pressure ($P_T$) and density ratio $P_T/\rho_T\rightarrow 1/3$ in the high redshift regime $(z \gg 0)$, which is different from the previous result in the literature. We use the Laurent expansion to resolve the solution. We also illustrate our result by the execution of numerical computations.

We propose the simplest model of teleparallel dark energy with purely a non-minimal coupling to gravity but no self-potential, a single model possessing various interesting features: simplicity, self-potential-free, the guaranteed late-time cosmic acceleration driven by the non-minimal coupling to gravity, tracker behavior of the dark energy equation of state at earlier times, a crossing of the phantom divide at a late time, and the existence of a finite-time future singularity. We find the analytic solutions of the dark-energy scalar field respectively in the radiation, matter, and dark energy dominated eras, thereby revealing the above features. We further illustrate possible cosmic evolution patterns and present the observational constraint of this model obtained by numerical analysis and data fitting.

We study the cosmological perturbations for the possible inflation scenario in the teleparallel equivalence of general relativity specified with parallelizable topological conditions. By acquiring the identical physical observables to general relativity under the teleparallel formalism, we perform a 3+1 decomposition of the vierbein field, which can be interpreted as the time gauge fixing between coordinate and tangent frames. We also extend our discussion to the higher-order action, $f(T)$ gravity.

We investigate the equation of state (EoS) of the scalar-torsion mode in Poincar\'{e} gauge theory of gravity. We concentrate on two cases with the constant curvature solution and positive kinetic energy, respectively. In the former, we find that the torsion EoS has different values in the various stages of the universe. In particular, it behaves like the radiation (matter) EoS of $w_r=1/3$ ($w_m=0$) in the radiation (matter) dominant epoch, while in the late time the torsion density is supportive for the accelerating universe. In the latter, our numerical analysis shows that in general the EoS has an asymptotic behavior in the high redshift regime, while it could cross the phantom divide line in the low redshift regime.

We explore the generation of large-scale magnetic fields from inflation in teleparallelism, in which the gravitational theory is described by the torsion scalar instead of the scalar curvature in general relativity. In particular, we examine the case that the conformal invariance of the electromagnetic field during inflation is broken by a non-minimal gravitational coupling between the torsion scalar and the electromagnetic field. It is shown that for a power-law type coupling, the magnetic field on 1 Mpc scale with its strength of $\sim 10^{-9}$ G at the present time can be generated.

The curvature singularity in viable f(R) gravity models is examined when the background density is dense. This singularity could be eliminated by adding the $R^{2}$ term in the Lagrangian. Some of cosmological consequences, in particular the source for the scalar mode of gravitational waves, are discussed.

We investigate the generation of large-scale magnetic fields due to the breaking of the conformal invariance in the electromagnetic field through the $CPT$-even dimension-six Chern-Simons-like effective interaction with a fermion current by taking account of the dynamical Kalb-Ramond and scalar fields in inflationary cosmology. It is explicitly demonstrated that the magnetic fields on 1Mpc scale with the field strength of $\sim 10^{-9}$G at the present time can be induced.

We study magnetic dipole moments of right-handed neutrinos in a keV neutrino dark matter model. This model is a simple extension of the standard model with only right-handed neutrinos and a pair of charged particles added. One of the right-handed neutrinos is the candidate of dark matter with a keV mass. Some bounds on the dark matter magnetic dipole moment and model parameters are obtained from cosmological observations.

We use data from Type Ia Supernovae (SNIa), Baryon Acoustic Oscillations (BAO), and Cosmic Microwave Background (CMB) observations to constrain the recently proposed teleparallel dark energy scenario based on the teleparallel equivalent of General Relativity, in which one adds a canonical scalar field, allowing also for a nonminimal coupling with gravity. Using the power-law, the exponential and the inverse hyperbolic cosine potential ansatzes, we show that the scenario is compatible with observations. In particular, the data favor a nonminimal coupling, and although the scalar field is canonical the model can describe both the quintessence and phantom regimes.

We explore thermodynamics of the apparent horizon in $f(T)$ gravity with both equilibrium and non-equilibrium descriptions. We find the same dual equilibrium/non-equilibrium formulation for $f(T)$ as for $f(R)$ gravity. In particular, we show that the second law of thermodynamics can be satisfied for the universe with the same temperature of the outside and inside the apparent horizon.

Using the "teleparallel" equivalent of General Relativity as the gravitational sector, which is based on torsion instead of curvature, we add a canonical scalar field, allowing for a nonminimal coupling with gravity. Although the minimal case is completely equivalent to standard quintessence, the nonminimal scenario has a richer structure, exhibiting quintessence-like or phantom-like behavior, or experiencing the phantom-divide crossing. The richer structure is manifested in the absence of a conformal transformation to an equivalent minimally-coupled model.

We study gravitational waves in viable $f(R)$ theories under a non-zero background curvature. In general, an $f(R)$ theory contains an extra scalar degree of freedom corresponding to a massive scalar mode of gravitational wave. For viable $f(R)$ models, since there always exits a de-Sitter point where the background curvature in vacuum is non-zero, the mass squared of the scalar mode of gravitational wave is about the de-Sitter point curvature $R_{d}\sim10^{-66}eV^{2}$. We illustrate our results in two types of viable $f(R)$ models: the exponential gravity and Starobinsky models. In both cases, the mass will be in the order of $10^{-33}eV$ when it propagates in vacuum. However, in the presence of matter density in galaxy, the scalar mode can be heavy. Explicitly, in the exponential gravity model, the mass becomes almost infinity, implying the disappearance of the scalar mode of gravitational wave, while the Starobinsky model gives the lowest mass around $10^{-24}eV$, corresponding to the lowest frequency of $10^{-9}$ Hz, which may be detected by the current and future gravitational wave probes, such as LISA and ASTROD-GW.

We demonstrate that there does exist an equilibrium description of thermodynamics on the apparent horizon in the expanding cosmological background for a wide class of modified gravity theories with the Lagrangian density $f(R, \phi, X)$, where $R$ is the Ricci scalar and $X$ is the kinetic energy of a scalar field $\phi$. This comes from a suitable definition of an energy momentum tensor of the "dark" component obeying the local energy conservation law in the Jordan frame. It is shown that the equilibrium description in terms of the horizon entropy $S$ is convenient because it takes into account the contribution of the horizon entropy $\hat{S}$ in non-equilibrium thermodynamics as well as an entropy production term.

We study the cosmological evolutions of the equation of state for dark energy $w_{\mathrm{DE}}$ in the exponential and logarithmic as well as their combination $f(T)$ theories. We show that the crossing of the phantom divide line of $w_{\mathrm{DE}} = -1$ can be realized in the combined $f(T)$ theory even though it cannot be in the exponential or logarithmic $f(T)$ theory. In particular, the crossing is from $w_{\mathrm{DE}} > -1$ to $w_{\mathrm{DE}} < -1$, in the opposite manner from $f(R)$ gravity models. We also demonstrate that this feature is favored by the recent observational data.

We study the equation of state for dark energy and explicitly demonstrate that the future crossings of the phantom divide line $w_{\mathrm{DE}}=-1$ are the generic feature in the existing viable $f(R)$ gravity models. We also explore the future evolution of the cosmological horizon entropy and illustrate that the cosmological horizon entropy oscillates with time due to the oscillatory behavior of the Hubble parameter. The important cosmological consequence is that in the future, the sign of the time derivative of the Hubble parameter changes from negative to positive in these viable $f(R)$ gravity models.

We study the observational constraints on the exponential gravity model of f(R)=-beta*Rs(1-e^(-R/Rs)). We use the latest observational data including Supernova Cosmology Project (SCP) Union2 compilation, Two-Degree Field Galaxy Redshift Survey (2dFGRS), Sloan Digital Sky Survey Data Release 7 (SDSS DR7) and Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP7) in our analysis. From these observations, we obtain a lower bound on the model parameter beta at 1.27 (95% CL) but no appreciable upper bound. The constraint on the present matter density parameter is 0.245< Omega_m^0<0.311 (95% CL). We also find out the best-fit value of model parameters on several cases.

We explore the cosmological evolution in the exponential gravity $f(R)=R +c_1 (1-e^{- c_2 R})$ ($c_{1, 2} = \mathrm{constant}$). We summarize various viability conditions and explicitly demonstrate that the late-time cosmic acceleration following the matter-dominated stage can be realized. We also study the equation of state for dark energy and confirm that the crossing of the phantom divide from the phantom phase to the non-phantom (quintessence) one can occur. Furthermore, we illustrate that the cosmological horizon entropy globally increases with time.

We show that in the exponential $f(T)$ model studied by E. Linder [arXiv:1005.3039, Phys.Rev.D 81, 127301 (2010)], it is impossible to have the crossing of the phantom divide line $w_{\mathrm{DE}}=-1$.

We study the CPT-even dimension-six Chern-Simons-like term by including dynamical Kalb-Ramond and scalar fields to examine the cosmological birefringence. We show that the combined effect of neutrino current and Kalb-Ramond field could induce a sizable rotation polarization angle in the cosmic microwave background radiation polarization.

We investigate thermodynamics of the apparent horizon in $f(R)$ gravity in the Palatini formalism with non-equilibrium and equilibrium descriptions. We demonstrate that it is more transparent to understand the horizon entropy in the equilibrium framework than that in the non-equilibrium one. Furthermore, we show that the second law of thermodynamics can be explicitly verified in both phantom and non-phantom phases for the same temperature of the universe outside and inside the apparent horizon.

We show that it is possible to obtain a picture of equilibrium thermodynamics on the apparent horizon in the expanding cosmological background for a wide class of modified gravity theories with the Lagrangian density $f(R, \phi, X)$, where $R$ is the Ricci scalar and $X$ is the kinetic energy of a scalar field $\phi$. This comes from a suitable definition of an energy momentum tensor of the "dark" component that respects to a local energy conservation in the Jordan frame. In this framework the horizon entropy $S$ corresponding to equilibrium thermodynamics is equal to a quarter of the horizon area $A$ in units of gravitational constant $G$, as in Einstein gravity. For a flat cosmological background with a decreasing Hubble parameter, $S$ globally increases with time, as it happens for viable $f(R)$ inflation and dark energy models. We also show that the equilibrium description in terms of the horizon entropy $S$ is convenient because it takes into account the contribution of both the horizon entropy $\hat{S}$ in non-equilibrium thermodynamics and an entropy production term.

An explicit model of $F(R)$ gravity with realizing a crossing of the phantom divide is reconstructed. In particular, it is shown that the Big Rip singularity may appear in the reconstructed model of $F(R)$ gravity. Such a Big Rip singularity could be avoided by adding $R^2$ term or non-singular viable $F(R)$ theory to the model because phantom behavior becomes transient.

We explicitly show that the equations of motion for modified gravity theories of $F(R)$-gravity, the scalar-Gauss-Bonnet gravity, $F(\mathcal{G})$-gravity and the non-local gravity are equivalent to the Clausius relation in thermodynamics. In addition, we discuss the relation between the expression of the entropy and the contribution from the modified gravity as well as the matter to the definition of the energy flux (heat).

We investigate the oscillating effective equation of state (EoS) of the universe around the phantom divide in the framework of $F(R)$ gravity. We illustrate the behavior of $F(R)$ with realizing multiple crossings of the phantom divide.

We study thermodynamics of the apparent horizon in $F(R)$ gravity. In particular, we demonstrate that a $F(R)$ gravity model with realizing a crossing of the phantom divide can satisfy the second law of thermodynamics in the effective phantom phase as well as non-phantom one.

We reconstruct an explicit model of modified gravity in which a crossing of the phantom divide can be realized. It is shown that the (finite-time) Big Rip singularity appears in the model of modified gravity (i.e., in the so-called Jordan frame), whereas that in the corresponding scalar field theory obtained through the conformal transformation (i.e., in the so-called Einstein frame) the singularity becomes the infinite-time one. Furthermore, we investigate the relations between the scalar field theories with realizing a crossing of the phantom divide and the corresponding modified gravitational theories by using the inverse conformal transformation. It is demonstrated that the scalar field theories describing the non-phantom phase (phantom one with the Big Rip) can be represented as the theories of real (complex) $F(R)$ gravity through the inverse (complex) conformal transformation. We also study a viable model of modified gravity in which the transition from the de Sitter universe to the phantom phase can occur. In addition, we explore the stability for the obtained solutions of the crossing of the phantom divide under a quantum correction coming from conformal anomaly.

We discuss the generation of large-scale magnetic fields due to the breaking of the conformal invariance in the electromagnetic field through the $CPT$-even dimension-six Chern-Simons-like effective interaction with a fermion current in inflationary cosmology. It is shown that the magnetic fields on 1Mpc scale with the field strength of $\sim 10^{-9}$G at the present time can be generated even for the scale of the effective interaction being the Planck scale.

We study the models with radiative neutrino mass generation and explore the relation between the neutrino masses and dark energy. In these models, the pseudo-Nambu-Goldston bosons (pNGBs) arise at two-loop level via the Majorana neutrino masses. In particular, we demonstrate that the potential energy of the pNGB can be the dark energy potential and the observed value of the equation of state (EoS) parameter of the universe, $i.e.$, $w\simeq -1$, can be realized.

We study a new scenario for baryogenesis due to the spontaneous breaking of the $CPT$ invariance through the interaction between a baryon current and a hypermagnetic helicity. The hypermagnetic helicity (Chern-Simons number) of $U(1)_Y$ provides a $CPT$ violation background for the generation of baryons via sphaleron processes, which protects these baryons from the sphaleron wash-out effect in thermal equilibrium. It is shown that if the present amplitude of the resultant magnetic fields are sufficiently large, for a wide range mass scale (from TeV to the Planck scale), the observational magnitude of the baryon asymmetry of the Universe can be realized.

We demonstrate that the cosmological birefringence can arise from CPT conserving effect, originated from the CPT-even dimension-six Chern-Simons-like term. We show that a sizable rotation polarization angle in the data of the cosmic microwave background radiation polarization can be induced.

We review our recent work on the cosmological birefringence. We propose a new type of effective interactions in terms of the $CPT$-even dimension-six Chern-Simons-like term to generate the cosmological birefringence. We use the neutrino number asymmetry to induce a non-zero rotation polarization angle in the data of the cosmic microwave background radiation polarization.

We study a new type of effective interactions in terms of the $CPT$-even dimension-six Chern-Simons-like term, which could originate from superstring theory, to generate the cosmological birefringence. We use the neutrino number asymmetry to induce a sizable rotation polarization angle in the data of the cosmic microwave background radiation polarization. The combined effect of the new term and the neutrino asymmetry provides an alternative way to understand the birefringence.