We propose new massive gravity theories with 5 dynamical degrees of freedom. We evade uniqueness theorems regarding the form of the kinetic and potential terms by adopting the "generalized massive gravity" framework, where a global translation invariance is broken. By exploiting the rotation symmetry in the field space, we determine two novel classes of theories. The first one is an extension of generalized massive gravity with a non-minimal coupling. On the other hand, the second theory produces a mass term that is different from de Rham, Gabadadze, Tolley construction and trivially has 5 degrees of freedom. Both theories allows for stable cosmological solutions without infinite strong coupling, which are free of ghost and gradient instabilities.

We study thermal evolution of isolated neutron stars in scalar-tensor theories for the first time. Whether the rapid cooling due to the direct Urca process occurs or not is an interesting question in the viewpoint of the temperature observation of isolated neutron stars. The cooling effect is typically influenced by the proton fraction and the central density. If a fifth force is mediated due to modification of gravity, the relation between the central density and mass of neutron stars differs from one in general relativity, and the cooling curve is also naively expected to be varied. We find that an unscreened fifth force near the surface of neutron stars changes mass-central density relation, and the direct Urca process can be triggered even for neutron stars with smaller mass. This might enable us to test gravitational theories with neutron star cooling.

In this paper, we fully investigate cosmological scenarios in which dark matter is non-minimally coupled to an extra scalar degree of freedom. The interaction is realized by means of conformal and disformal terms in the transformed gravitational metric. Considering linear perturbation theory, we show that the growth rate of dark matter differs from the uncoupled case and that the well-known Kaiser formula undergoes modification. As a result, redshift space distortion measurements cease to be a direct probe of the linear growth rate of total matter, since the distortion factor has an extra, coupling-dependent term. We study the effect of the coupling in three cosmological models, two conformally and one disformally coupled, and forecast the constraints on the coupling, and other cosmological parameters, from future galaxy surveys.

We study cosmological applications of extended vector-tensor theories, whose Lagrangians contain up to two derivatives with respect to metric and vector field. We derive background equations under the assumption of homogeneous and isotropic universe and study the nature of cosmological perturbations on top of the background. We found an example of healthy cosmological solutions in non-trivial degenerate theory of gravity for the first time, where those perturbations do not suffer from any instabilities, that is, ghost and gradient instabilities.

We present an impact of coupling between dark matter and a scalar field, which might be responsible for dark energy, on measurements of redshift-space distortions. We point out that, in the presence of conformal and/or disformal coupling, linearized continuity and Euler equations for total matter fluid significantly deviate from the standard ones even in the sub-horizon scales. In such a case, a peculiar velocity of total matter field is determined not only by a logarithmic time derivative of its density perturbation but also by density perturbations for both dark matter and baryon, leading to a large modification of the physical interpretation of observed data obtained by measurements of redshift-space distortions. We reformulate galaxy two-point correlation function in the redshift space based on the modified continuity and Euler equations. We conclude from the resultant formula that the true value of the linear growth rate of large-scale structure cannot be necessarily constrained by single-redshift measurements of the redshift-space distortions, unless one observes the actual time-evolution of structure.

We study co-existence system of both bosonic and fermionic degrees of freedom. For such system with up to first derivatives in Lagrangian, we find Ostrogradsky-type ghost-free condition in Hamiltonian analysis, which is found to be the same with requiring that the equations of motion of fermions are first-order in Lagrangian formulation. When fermionic degrees of freedom are present, uniqueness of time evolution is not guaranteed a priori because of the Grassmann property. We confirm that the additional condition, which is introduced to close Hamiltonian analysis, also ensures the uniqueness of the time evolution of system.

Recently, several extensions of massive vector theory in curved space-time have been proposed in many literatures. In this paper, we consider the most general vector-tensor theories that contain up to two derivatives with respect to metric and vector field. By imposing a degeneracy condition of the Lagrangian in the context of ADM decomposition of space-time to eliminate an unwanted mode, we construct a new class of massive vector theories where five degrees of freedom can propagate, corresponding to three for massive vector modes and two for massless tensor modes. We find that the generalized Proca and the beyond generalized Proca theories up to the quartic Lagrangian, which should be included in this formulation, are degenerate theories even in curved space-time. Finally, introducing new metric and vector field transformations, we investigate the properties of thus obtained theories under such transformations.

We investigate gravitational Cherenkov radiation in a healthy branch of background solutions in the ghost-free bigravity model. In this model, because of the modification of dispersion relations, each polarization mode can possess subluminal phase velocities, and the gravitational Cherenkov radiation could be potentially emitted from a relativistic particle. In the present paper, we derive conditions for the process of the gravitational Cherenkov radiation to occur and estimate the energy emission rate for each polarization mode. We found that the gravitational Cherenkov radiation emitted even from an ultrahigh energy cosmic ray is sufficiently suppressed for the graviton's effective mass less than $100\,{\rm eV}$, and the bigravity model with dark matter coupled to the hidden metric is therefore consistent with observations of high energy cosmic rays.

In this paper, we scrutinize very closely the cosmology in the proxy theory to massive gravity obtained in Phys. Rev. D84 (2011) 043503. This proxy theory was constructed by covariantizing the decoupling limit Lagrangian of massive gravity and represents a subclass of Horndeski scalar-tensor theory. Thus, this covariantization unifies two important classes of modified gravity theories, namely massive gravity and Horndeski theories. We go beyond the regime which was studied in Phys. Rev. D84 (2011) 043503 and show that the theory does not admit any homogeneous and isotropic self-accelerated solutions. We illustrate that the only attractor solution is flat Minkowski solution, hence this theory is less appealing as a dark energy model. We also show that the absence of de Sitter solutions is tightly related to the presence of shift symmetry breaking interactions.

Quasidilaton massive gravity is an extension of massive General Relativity to a theory with additional scale invariance and approximate internal Galilean symmetry. The theory has a novel self-accelerated solution with the metric indistinguishable (in the decoupling limit) from the de Sitter space, its curvature set by the graviton mass. The spectra of tensor, vector, and scalar perturbations on this solution contain neither ghosts, nor gradient instabilities or superluminal modes, for a range of the parameter space. This represents an example of a self-accelerated solution with viable perturbations, attainable within a low energy effective field theory.

We demonstrate that the general second-order scalar-tensor theories, which have attracted attention as possible modified gravity models to explain the late time cosmic acceleration, could be strongly constrained from the argument of the gravitational Cherenkov radiation. To this end, we consider the purely kinetic coupled gravity and the extended galileon model. In these models, the propagation speed of tensor mode could be less than the speed of light, which puts very strong constraints from the gravitational Cherenkov radiation.

Detections of molecular lines, mainly from H2$ and CO, reveal molecular material in planetary nebulae. Observations of a variety of molecules suggest that the molecular composition in these objects differs from that found in interstellar clouds or in circumstellar envelopes. The success of the models, which are mostly devoted to explain molecular densities in specific planetary nebulae, is still partial, however. The present study aims at identifying the influence of stellar and nebular properties on the molecular composition of planetary nebulae by means of chemical models. A comparison of theoretical results with those derived from the observations may provide clues to the conditions that favor the presence of a particular molecule. A self-consistent photoionization numerical code was adapted to simulate cold molecular regions beyond the ionized zone. The code was used to obtain a grid of models and the resulting column densities are compared with those inferred from observations. Our models show that the inclusion of an incident flux of X-rays is required to explain the molecular composition derived for planetary nebulae. We also obtain a more accurate relation for the N(CO)/N(H2) ratio in these objects. Molecular masses obtained by previous works in the literature were then recalculated, showing that these masses can be underestimated by up to three orders of magnitude. We conclude that the problem of the missing mass in planetary nebulae can be solved by a more accurate calculation of the molecular mass.

The cross-correlation between the integrated Sachs-Wolfe (ISW) effect and the large scale structure (LSS) is a powerful tool to constrain dark energy and alternative theories of gravity. In this paper, we obtain observational constraints on kinetic gravity braiding from the ISW-LSS cross-correlation. We find that the late-time ISW effect in the kinetic gravity braiding model anti-correlates with large scale structures in a wide range of parameters, which clearly demonstrates how one can distinguish modified gravity theories from the LCDM model using the ISW effect. In addition to the analysis based on a concrete model, we investigate a future prospect of the ISW-LSS cross-correlation by using a phenomenological parameterization of modified gravity models.

A generic second-order scalar-tensor theory contains a nonlinear derivative self-interaction of the scalar degree of freedom $\phi$ \`{a} la Galileon models, which allows for the Vainshtein screening mechanism. We investigate this effect on subhorizon scales in a cosmological background, based on the most general second-order scalar-tensor theory. Our analysis takes into account all the relevant nonlinear terms and the effect of metric perturbations consistently. We derive an explicit form of Newton's constant, which in general is time-dependent and hence is constrained from observations, as suggested earlier. It is argued that in the most general case the inverse-square law cannot be reproduced on the smallest scales. Some applications of our results are also presented.

We investigate the structure of halos in the sDGP (self-accelerating branch of the Dvali-Gavadadze-Porrati braneworld gravity) model and the galileon modified gravity model on the basis of the static and spherically symmetric solutions of the collisionless Boltzmann equation, which reduce to the singular isothermal sphere model and the King model in the limit of Newtonian gravity. The common feature of these halos is that the density of a halo in the outer region is larger (smaller) in the sDGP (galileon) model, respectively, in comparison with Newtonian gravity. This comes from the suppression (enhancement) of the effective gravity at large distance in the sDGP (galileon) model, respectively. However, the difference between these modified gravity models and Newtonian gravity only appears outside the halo due to the Vainshtein mechanism, which makes it difficult to distinguish between them. We also discuss the case in which the halo density profile is fixed independently of the gravity model for comparison between our results and previous work.

Molecular hydrogen emission is commonly observed in planetary nebulae. Images taken in infrared H2 emission lines show that at least part of the molecular emission is produced inside the ionised region. In the best-studied case, the Helix nebula, the H2 emission is produced inside cometary knots (CKs), comet-shaped structures believed to be clumps of dense neutral gas embedded within the ionised gas. Most of the H2 emission of the CKs seems to be produced in a thin layer between the ionised diffuse gas and the neutral material of the knot, in a mini photodissociation region (PDR). However, PDR models published so far cannot fully explain all the characteristics of the H2 emission of the CKs. In this work, we use the photoionisation code \textsc{Aangaba} to study the H2 emission of the CKs, particularly that produced in the interface H^+/H^0 of the knot, where a significant fraction of the H2 1-0S(1) emission seems to be produced. Our results show that the production of molecular hydrogen in such a region may explain several characteristics of the observed emission, particularly the high excitation temperature of the H2 infrared lines. We find that the temperature derived from H2 observations even of a single knot, will depend very strongly on the observed transitions, with much higher temperatures derived from excited levels. We also proposed that the separation between the H_alpha and NII peak emission observed in the images of CKs may be an effect of the distance of the knot from the star, since for knots farther from the central star the NII line is produced closer to the border of the CK than H_alpha.

We study cosmological consequences of a kinetic gravity braiding model, which is proposed as an alternative to the dark energy model. The kinetic braiding model we study is characterized by a parameter n, which corresponds to the original galileon cosmological model for n=1. We find that the background expansion of the universe of the kinetic braiding model is the same as the Dvali-Turner's model, which reduces to that of the standard cold dark matter model with a cosmological constant (LCDM model) for n equal to infinity. We also find that the evolution of the linear cosmological perturbation in the kinetic braiding model reduces to that of the LCDM model for n=\infty. Then, we focus our study on the growth history of the linear density perturbation as well as the spherical collapse in the nonlinear regime of the density perturbations, which might be important in order to distinguish between the kinetic braiding model and the LCDM model when n is finite. The theoretical prediction for the large scale structure is confronted with the multipole power spectrum of the luminous red galaxy sample of the Sloan Digital Sky survey. We also discuss future prospects of constraining the kinetic braiding model using a future redshift survey like the WFMOS/SuMIRe PFS survey as well as the cluster redshift distribution in the South Pole Telescope survey.

In the present work, we use a photoionisation code to study the H2 emission of the Helix nebula (NGC 7293) cometary knots, particularly that produced in the interface H+/H0 of the knot, where a significant fraction of the H2 1-0 S(1) emission seems to be produced. Our results show that the production of molecular hydrogen in such region may explain several characteristics of the observed emission, particularly the high excitation temperature of the H2 infrared lines.