For a scalar field $\phi$ coupled to cold dark matter (CDM), we provide a general framework for studying the background and perturbation dynamics on the isotropic cosmological background. The dark energy sector is described by a Horndeski Lagrangian with the speed of gravitational waves equivalent to that of light, whereas CDM is dealt as a perfect fluid characterized by the number density $n_c$ and four-velocity $u_c^\mu$. For a very general interacting Lagrangian $f(n_c, \phi, X, Z)$, where $f$ depends on $n_c$, $\phi$, $X=-\partial^{\mu} \phi \partial_{\mu} \phi/2$, and $Z=u_c^{\mu} \partial_{\mu} \phi$, we derive the full linear perturbation equations of motion without fixing any gauge conditions. To realize a vanishing CDM sound speed for the successful structure formation, the interacting function needs to be of the form $f=-f_1(\phi, X, Z)n_c+f_2(\phi, X, Z)$. Employing a quasi-static approximation for the modes deep inside the sound horizon, we obtain analytic formulas for the effective gravitational couplings of CDM and baryon density perturbations as well as gravitational and weak lensing potentials. We apply our general formulas to several interacting theories and show that, in many cases, the CDM gravitational coupling around the quasi de-Sitter background can be smaller than the Newton constant $G$ due to a momentum transfer induced by the $Z$-dependence in $f_2$.

We investigate a dark energy scenario in which a canonical scalar field $\phi$ is coupled to the four velocity $u_{c}^{\mu}$ of cold dark matter (CDM) through a derivative interaction $u_{c}^{\mu} \partial_{\mu} \phi$. The coupling is described by an interacting Lagrangian $f(X, Z)$, where $f$ depends on $X=-\partial^{\mu} \phi \partial_{\mu} \phi/2$ and $Z=u_{c}^{\mu} \partial_{\mu} \phi$. We derive stability conditions of linear scalar perturbations for the wavelength deep inside the Hubble radius and show that the effective CDM sound speed is close to 0 as in the standard uncoupled case, while the scalar-field propagation speed is affected by the interacting term $f$. Under a quasi-static approximation, we also obtain a general expression of the effective gravitational coupling felt by the CDM perturbation. We study the late-time cosmological dynamics for the coupling $f \propto X^{(2-m)/2}Z^m$ and show that the gravitational coupling weaker than the Newton constant can be naturally realized for $m>0$ on scales relevant to the growth of large-scale structures. This allows the possibility for alleviating the tension of $\sigma_8$ between low- and high-redshift measurements.

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.

We provide a general framework for studying the dark energy cosmology in which a scalar field $\phi$ is nonminimally and kinetically coupled to Cold Dark Matter (CDM). The scalar-graviton sector is described by the action of Horndeski theories with the speed of gravitational waves equivalent to that of light, whereas CDM is treated as a perfect fluid given by a Schutz-Sorkin action. We consider two interacting Lagrangians of the forms $f_1(\phi,X)\rho_c (n_c)$ and $f_2 (n_c, \phi,X) J_c^{\mu} \partial_{\mu}\phi$, where $X=-\partial^{\mu} \phi \partial_{\mu} \phi/2$, $\rho_c$ and $n_c$ are the energy density and number density of CDM respectively, and $J_c^{\mu}$ is a vector field related to the CDM four velocity. We derive the scalar perturbation equations of motion without choosing any special gauges and identify conditions for the absence of ghosts and Laplacian instabilities on scales deep inside the sound horizon. Applying a quasi-static approximation in a gauge-invariant manner, we also obtain the effective gravitational couplings felt by CDM and baryons for the modes relevant to the linear growth of large-scale structures. In particular, the $n_c$ dependence in the coupling $f_2$ gives rise to an interesting possibility for realizing the gravitational coupling with CDM weaker than the Newton gravitational constant $G$.

In a class of generalized Proca theories, we study the existence of neutron star solutions with a nonvanishing temporal component of the vector field $A_\mu$ approaching 0 toward spatial infinity, as they may be the endpoints of tachyonic instabilities of neutron star solutions in general relativity with $A_{\mu}=0$. Such a phenomenon is called spontaneous vectorization, which is analogous to spontaneous scalarization in scalar-tensor theories with nonminimal couplings to the curvature or matter. For the nonminimal coupling $\beta X R$, where $\beta$ is a coupling constant and $X=-A_{\mu}A^{\mu}/2$, we show that there exist both 0-node and 1-node vector-field solutions, irrespective of the choice of the equations of state of nuclear matter. The 0-node solution, which is present only for $\beta=-{\cal O}(0.1)$, may be induced by some nonlinear effects such as the selected choice of initial conditions. The 1-node solution exists for $\beta=-{\cal O}(1)$, which suddenly emerges above a critical central density of star and approaches the general relativistic branch with the increasing central density. We compute the mass $M$ and radius $r_s$ of neutron stars for some realistic equations of state and show that the $M$-$r_s$ relations of 0-node and 1-node solutions exhibit notable difference from those of scalarized solutions in scalar-tensor theories. Finally, we discuss the possible endpoints of tachyonic instabilities.

We provide a general framework for studying the evolution of background and cosmological perturbations in the presence of a vector field $A_{\mu}$ coupled to cold dark matter (CDM). We consider an interacting Lagrangian of the form $Q f(X) T_c$, where $Q$ is a coupling constant, $f$ is an arbitrary function of $X=-A_{\mu}A^{\mu}/2$, and $T_c$ is a trace of the CDM energy-momentum tensor. The matter coupling affects the no-ghost condition and sound speed of linear scalar perturbations deep inside the sound horizon, while those of tensor and vector perturbations are not subject to modifications. The existence of interactions also modifies the no-ghost condition of CDM density perturbations. We propose a concrete model of coupled vector dark energy with the tensor propagation speed equivalent to that of light. In comparison to the $Q=0$ case, we show that the decay of CDM to the vector field leads to the phantom dark energy equation of state $w_{\rm DE}$ closer to $-1$. This alleviates the problem of observational incompatibility of uncoupled models in which $w_{\rm DE}$ significantly deviates from $-1$. The maximum values of $w_{\rm DE}$ reached during the matter era are bounded from the CDM no-ghost condition of future de Sitter solutions.

In $f(R)$ gravity and Brans-Dicke theory with scalar potentials, we study the structure of neutron stars on a spherically symmetric and static background for two equations of state: SLy and FPS. In massless BD theory, the presence of a scalar coupling $Q$ with matter works to change the star radius in comparison to General Relativity, while the maximum allowed mass of neutron stars is hardly modified for both SLy and FPS equations of state. In Brans-Dicke theory with the massive potential $V(\phi)=m^2 \phi^2/2$, where $m^2$ is a positive constant, we show the difficulty of realizing neutron star solutions with a stable field profile due to the existence of an exponentially growing mode outside the star. As in $f(R)$ gravity with the $R^2$ term, this property is related to the requirement of extra boundary conditions of the field at the surface of star. For the self-coupling potential $V(\phi)=\lambda \phi^4/4$, this problem can be circumvented by the fact that the second derivative $V_{,\phi \phi}=3\lambda\phi^2$ approaches 0 at spatial infinity. In this case, we numerically show the existence of neutron star solutions for both SLy and FPS equations of state and discuss how the mass-radius relation is modified as compared to General Relativity.

The gravitational-wave event GW170817 from a binary neutron star merger together with the electromagnetic counterpart showed that the speed of gravitational waves $c_t$ is very close to that of light for the redshift $z<0.009$. This places tight constraints on dark energy models constructed in the framework of modified gravitational theories. We review models of the late-time cosmic acceleration in scalar-tensor theories with second-order equations of motion (dubbed Horndeski theories) by paying particular attention to the evolution of dark energy equation of state and observables relevant to the cosmic growth history. We provide a gauge-ready formulation of scalar perturbations in full Horndeski theories and estimate observables associated with the evolution of large-scale structures, cosmic microwave background, and weak lensing by employing a so-called quasi-static approximation for the modes deep inside the sound horizon. In light of the recent observational bound of $c_t$, we also classify surviving dark energy models into four classes depending on different structure-formation patterns and discuss how they can be observationally distinguished from each other. In particular, the nonminimally coupled theories in which the scalar field $\phi$ has a coupling with the Ricci scalar $R$ of the form $G_4(\phi) R$, including $f(R)$ gravity, can be tightly constrained not only from the cosmic expansion and growth histories but also from the variation of screened gravitational couplings. The cross correlation of integrated Sachs-Wolfe signal with galaxy distributions can be a key observable for placing bounds on the relative ratio of cubic Galileon density to total dark energy density. The dawn of gravitational-wave astronomy will open up a new window to constrain nonminimally coupled theories further by the modified luminosity distance of tensor perturbations.

The gravitational-wave event GW170817, together with the electromagnetic counterpart, shows that the speed of tensor perturbations $c_T$ on the cosmological background is very close to that of light $c$ for the redshift $z<0.009$. In generalized Proca theories, the Lagrangians compatible with the condition $c_T=c$ are constrained to be derivative interactions up to cubic order, besides those corresponding to intrinsic vector modes. We place observational constraints on a dark energy model in cubic-order generalized Proca theories with intrinsic vector modes by running the Markov chain Monte Carlo (MCMC) code. We use the cross-correlation data of the integrated Sachs-Wolfe (ISW) signal and galaxy distributions in addition to the data sets of cosmic microwave background, baryon acoustic oscillations, type Ia supernovae, local measurements of the Hubble expansion rate, and redshift-space distortions. We show that, unlike cubic-order scalar-tensor theories, the existence of intrinsic vector modes allows the possibility for evading the ISW-galaxy anticorrelation incompatible with the current observational data. As a result, we find that the dark energy model in cubic-order generalized Proca theories exhibits a better fit to the data than the cosmological constant, even by including the ISW-galaxy correlation data in the MCMC analysis.

We study the cosmology in the presence of arbitrary couplings between $p$-forms in 4-dimensional space-time for a general action respecting gauge symmetry and parity invariance. The interaction between 0-form (scalar field $\phi$) and 3-form fields gives rise to an effective potential $V_{\rm eff}(\phi)$ for the former after integrating out the contribution of the latter. We explore the dynamics of inflation on an anisotropic cosmological background for a coupled system of 0-, 1-, and 2-forms. In the absence of interactions between 1- and 2-forms, we derive conditions under which the anisotropic shear endowed with nearly constant energy densities of 1- and 2-forms survives during slow-roll inflation for an arbitrary scalar potential $V_{\rm eff}(\phi)$. If 1- and 2-forms are coupled to each other, we show the existence of a new class of anisotropic inflationary solutions in which the energy density of 2-form is sustained by that of 1-form through their interactions. Our general analytic formulas for the anisotropic shear are also confirmed by the numerical analysis for a concrete inflaton potential.

We study the dynamics of dark energy in the presence of a 2-form field coupled to a canonical scalar field $\phi$. We consider the coupling proportional to $e^{-\mu \phi/M_{\rm pl}} H_{\alpha \beta \gamma}H^{\alpha \beta \gamma}$ and the scalar potential $V(\phi) \propto e^{-\lambda \phi/M_{\rm pl}}$, where $H_{\alpha \beta \gamma}$ is the 2-form field strength, $\mu, \lambda$ are constants, and $M_{\rm pl}$ is the reduced Planck mass. We show the existence of an anisotropic matter-dominated scaling solution followed by a stable accelerated fixed point with a non-vanishing shear. Even if $\lambda \geq {\cal O}(1)$, it is possible to realize the dark energy equation of state $w_{\rm DE}$ close to $-1$ at low redshifts for $\mu \gg \lambda$. The existence of anisotropic hair and the oscillating behavior of $w_{\rm DE}$ are key features for distinguishing our scenario from other dark energy models like quintessence.

In quadratic-order degenerate higher-order scalar-tensor (DHOST) theories compatible with gravitational-wave constraints, we derive the most general Lagrangian allowing for tracker solutions characterized by $\dot{\phi}/H^p={\rm constant}$, where $\dot{\phi}$ is the time derivative of a scalar field $\phi$, $H$ is the Hubble expansion rate, and $p$ is a constant. While the tracker is present up to the cubic-order Horndeski Lagrangian $L=c_2X-c_3X^{(p-1)/(2p)} \square \phi$, where $c_2, c_3$ are constants and $X$ is the kinetic energy of $\phi$, the DHOST interaction breaks this structure for $p \neq 1$. Even in the latter case, however, there exists an approximate tracker solution in the early cosmological epoch with the nearly constant field equation of state $w_{\phi}=-1-2p\dot{H}/(3H^2)$. The scaling solution, which corresponds to $p=1$, is the unique case in which all the terms in the field density $\rho_{\phi}$ and the pressure $P_{\phi}$ obey the scaling relation $\rho_{\phi} \propto P_{\phi} \propto H^2$. Extending the analysis to the coupled DHOST theories with the field-dependent coupling $Q(\phi)$ between the scalar field and matter, we show that the scaling solution exists for $Q(\phi)=1/(\mu_1 \phi+\mu_2)$, where $\mu_1$ and $\mu_2$ are constants. For the constant $Q$, i.e., $\mu_1=0$, we derive fixed points of the dynamical system by using the general Lagrangian with scaling solutions. This result can be applied to the model construction of late-time cosmic acceleration preceded by the scaling $\phi$-matter-dominated epoch.

In cubic-order Horndeski theories where a scalar field $\phi$ is coupled to nonrelativistic matter with a field-dependent coupling $Q(\phi)$, we derive the most general Lagrangian having scaling solutions on the isotropic and homogenous cosmological background. For constant $Q$ including the case of vanishing coupling, the corresponding Lagrangian reduces to the form $L=Xg_2(Y)-g_3(Y)\square \phi$, where $X=-\partial_{\mu}\phi\partial^{\mu}\phi/2$ and $g_2, g_3$ are arbitrary functions of $Y=Xe^{\lambda \phi}$ with constant $\lambda$. We obtain the fixed points of the scaling Lagrangian for constant $Q$ and show that the $\phi$-matter-dominated-epoch ($\phi$MDE) is present for the cubic coupling $g_3(Y)$ containing inverse power-law functions of $Y$. The stability analysis around the fixed points indicates that the $\phi$MDE can be followed by a stable critical point responsible for the cosmic acceleration. We propose a concrete dark energy model allowing for such a cosmological sequence and show that the ghost and Laplacian instabilities can be avoided even in the presence of the cubic coupling.

In scalar-vector-tensor (SVT) theories with parity invariance, we perform a gauge-ready formulation of cosmological perturbations on the flat Friedmann-Lema\^{i}tre-Robertson-Walker (FLRW) background by taking into account a matter perfect fluid. We derive the second-order action of scalar perturbations and resulting linear perturbation equations of motion without fixing any gauge conditions. Depending on physical problems at hand, most convenient gauges can be chosen to study the development of inhomogeneities in the presence of scalar and vector fields coupled to gravity. This versatile framework, which encompasses Horndeski and generalized Proca theories as special cases, is applicable to a wide variety of cosmological phenomena including nonsingular cosmology, inflation, and dark energy. By deriving conditions for the absence of ghost and Laplacian instabilities in several different gauges, we show that, unlike Horndeski theories, it is possible to evade no-go arguments for the absence of stable nonsingular bouncing/genesis solutions in both generalized Proca and SVT theories. We also apply our framework to the case in which scalar and vector fields are responsible for dark energy and find that the separation of observables relevant to the evolution of matter perturbations into tensor, vector, and scalar sectors is transparent in the unitary gauge. Unlike the flat gauge chosen in the literature, this result is convenient to confront SVT theories with observations associated with the cosmic growth history.

The scalar-vector-tensor theories with second-order equations of motion can accommodate both Horndeski and generalized Proca theories as specific cases. In the presence of a perfect fluid, we study the cosmology in such a most general scheme of scalar-vector-tensor theories with parity invariance by paying particular attention to the application to dark energy. We obtain a closed-form expression of the background equations of motion by using coefficients appearing in the second-order action of scalar perturbations. We also derive conditions for the absence of ghost and Laplacian instabilities of tensor and vector perturbations and show that the existence of matter does not substantially modify the stabilities of dynamical degrees of freedom in the small-scale limit. On the other hand, the sound speed of scalar perturbations is affected by the presence of matter. Employing the quasi-static approximation for scalar perturbations deep inside the sound horizon, we derive analytic expressions of Newtonian and weak lensing gravitational potentials as well as two scalar perturbations arising from the scalar and vector fields. We apply our general framework to dark energy theories with the tensor propagation speed equivalent to the speed of light and show that the observables associated with the growth of matter perturbations and weak lensing potentials are generally affected by intrinsic vector modes and by interactions between scalar and vector fields.

We study the cosmology on the Friedmann-Lemaitre-Robertson-Walker background in scalar-vector-tensor theories with a broken $U(1)$ gauge symmetry. For parity-invariant interactions arising in scalar-vector-tensor theories with second-order equations of motion, we derive conditions for the absence of ghosts and Laplacian instabilities associated with tensor, vector, and scalar perturbations at linear order. This general result is applied to the computation of the primordial tensor power spectrum generated during inflation as well as to the speed of gravity relevant to dark energy. We also construct a concrete inflationary model in which a temporal vector component $A_0$ contributes to the dynamics of cosmic acceleration besides a scalar field $\phi$ through their kinetic mixings. In this model, we show that all the stability conditions of perturbations can be consistently satisfied during inflation and subsequent reheating.

In scalar-vector-tensor theories with $U(1)$ gauge invariance, it was recently shown that there exists a new type of hairy black hole (BH) solutions induced by a cubic-order scalar-vector interaction. In this paper, we derive conditions for the absence of ghosts and Laplacian instabilities against odd-parity perturbations on a static and spherically symmetric background for most general $U(1)$ gauge-invariant scalar-vector-tensor theories with second-order equations of motion. We apply those conditions to hairy BH solutions arising from the cubic-order coupling and show that the odd-parity stability in the gravity sector is always ensured outside the event horizon with the speed of gravity equivalent to that of light. We also study the case in which quartic-order interactions are present in addition to the cubic coupling and obtain conditions under which black holes are stable against odd-parity perturbations.

The generalized Proca theories with second-order equations of motion can be healthily extended to a more general framework in which the number of propagating degrees of freedom remains unchanged. In the presence of a quartic-order nonminimal coupling to gravity arising in beyond-generalized Proca theories, the speed of gravitational waves $c_t$ on the Friedmann-Lemaitre-Robertson-Walker (FLRW) cosmological background can be equal to that of light $c$ under a certain condition. By using this condition alone, we show that the speed of gravitational waves in the vicinity of static and spherically symmetric black holes is also equivalent to $c$ for the propagation of odd-parity perturbations along both radial and angular directions. As a by-product, the black holes arising in our beyond-generalized Proca theories are plagued by neither ghost nor Laplacian instabilities against odd-parity perturbations. We show the existence of both exact and numerical black hole solutions endowed with vector hairs induced by the quartic-order coupling.

The Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theories up to quartic order are the general scheme of scalar-tensor theories allowing the possibility for realizing the tensor propagation speed $c_t$ equivalent to 1 on the isotropic cosmological background. We propose a dark energy model in which the late-time cosmic acceleration occurs by a simple k-essence Lagrangian analogous to the ghost condensate with cubic and quartic Galileons in the framework of GLPV theories. We show that a wide variety of the variation of the dark energy equation of state $w_{\rm DE}$ including the entry to the region $w_{\rm DE}<-1$ can be realized without violating conditions for the absence of ghosts and Laplacian instabilities. The approach to the tracker equation of state $w_{\rm DE}=-2$ during the matter era, which is disfavored by observational data, can be avoided by the existence of a quadratic k-essence Lagrangian $X^2$. We study the evolution of nonrelativistic matter perturbations for the model $c_t^2=1$ and show that the two quantities $\mu$ and $\Sigma$, which are related to the Newtonian and weak lensing gravitational potentials respectively, are practically equivalent to each other, such that $\mu \simeq \Sigma>1$. For the case in which the deviation of $w_{\rm DE}$ from $-1$ is significant at a later cosmological epoch, the values of $\mu$ and $\Sigma$ tend to be larger at low redshifts. We also find that our dark energy model can be consistent with the bounds on the deviation parameter $\alpha_{\rm H}$ from Horndeski theories arising from the modification of gravitational law inside massive objects.

We study relativistic star solutions in second-order generalized Proca theories characterized by a $U(1)$-breaking vector field with derivative couplings. In the models with cubic and quartic derivative coupling, the mass and radius of stars become larger than those in general relativity for negative derivative coupling constants. This phenomenon is mostly attributed to the increase of star radius induced by a slower decrease of the matter pressure compared to general relativity. There is a tendency that the relativistic star with a smaller mass is not gravitationally bound for a low central density and hence dynamically unstable, but that with a larger mass is gravitationally bound. On the other hand, we show that the intrinsic vector-mode couplings give rise to general relativistic solutions with a trivial field profile, so the mass and radius are not modified from those in general relativity.

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.

In generalized Proca theories with vector-field derivative couplings, a bunch of hairy black hole solutions have been derived on a static and spherically symmetric background. In this paper, we formulate the odd-parity black hole perturbations in generalized Proca theories by expanding the corresponding action up to second order and investigate whether or not black holes with vector hair suffer ghost or Laplacian instabilities. We show that the models with cubic couplings $G_3(X)$, where $X=-A_{\mu}A^{\mu}/2$ with a vector field $A_{\mu}$, do not provide any additional stability condition as in General Relativity. On the other hand, the exact charged stealth Schwarzschild solution with a nonvanishing longitudinal vector component $A_1$, which originates from the coupling to the Einstein tensor $G^{\mu\nu}A_\mu A_\nu$ equivalent to the quartic coupling $G_4(X)$ containing a linear function of $X$, is unstable in the vicinity of the event horizon. The same instability problem also persists for hairy black holes arising from general quartic power-law couplings $G_4(X) \supset \beta_4 X^n$ with the nonvanishing $A_1$, while the other branch with $A_1=0$ can be consistent with conditions for the absence of ghost and Laplacian instabilities. We also discuss the case of other exact and numerical black hole solutions associated with intrinsic vector-field derivative couplings and show that there exists a wide range of parameter spaces in which the solutions suffer neither ghost nor Laplacian instabilities against odd-parity perturbations.

We present a family of exact black-hole solutions on a static spherically symmetric background in second-order generalized Proca theories with derivative vector-field interactions coupled to gravity. We also derive non-exact solutions in power-law coupling models including vector Galileons and numerically show the existence of regular black holes with a primary hair associated with the longitudinal propagation. The intrinsic vector-field derivative interactions generally give rise to a secondary hair induced by non-trivial field profiles. The deviation from General Relativity is most significant around the horizon and hence there is a golden opportunity for probing the Proca hair by the measurements of gravitational waves (GWs) in the regime of strong gravity.

We derive the profile of a vector field coupled to matter on a static and spherically symmetric background in the context of generalized Proca theories. The cubic Galileon self-interaction leads to the suppression of a longitudinal vector component due to the operation of the Vainshtein mechanism. For quartic and sixth-order derivative interactions, the solutions consistent with those in the continuous limit of small derivative couplings correspond to the branch with the vanishing longitudinal mode. We compute the corrections to gravitational potentials outside a compact body induced by the vector field in the presence of cubic, quartic, and sixth-order derivative couplings, and show that the models can be consistent with local gravity constraints under mild bounds on the temporal vector component. The quintic Galileon interaction does not allow regular solutions of the longitudinal mode for a rapidly decreasing matter density outside the body.

We study static and spherically symmetric black hole (BH) solutions in second-order generalized Proca theories with nonminimal vector field derivative couplings to the Ricci scalar, the Einstein tensor, and the double dual Riemann tensor. We find concrete Lagrangians which give rise to exact BH solutions by imposing two conditions of the two identical metric components and the constant norm of the vector field. These exact solutions are described by either Reissner-Nordstr\"{o}m (RN), stealth Schwarzschild, or extremal RN solutions with a non-trivial longitudinal mode of the vector field. We then numerically construct BH solutions without imposing these conditions. For cubic and quartic Lagrangians with power-law couplings which encompass vector Galileons as the specific cases, we show the existence of BH solutions with the difference between two non-trivial metric components. The quintic-order power-law couplings do not give rise to non-trivial BH solutions regular throughout the horizon exterior. The sixth-order and intrinsic vector-mode couplings can lead to BH solutions with a secondary hair. For all the solutions, the vector field is regular at least at the future or past horizon. The deviation from General Relativity induced by the Proca hair can be potentially tested by future measurements of gravitational waves in the nonlinear regime of gravity.

The beyond-generalized Proca theories are the extension of second-order massive vector-tensor theories (dubbed generalized Proca theories) with two transverse vector modes and one longitudinal scalar besides two tensor polarizations. Even with this extension, the propagating degrees of freedom remain unchanged on the isotropic cosmological background without an Ostrogradski instability. We study the cosmology in beyond-generalized Proca theories by paying particular attention to the dynamics of late-time cosmic acceleration and resulting observational consequences. We derive conditions for avoiding ghosts and instabilities of tensor, vector, and scalar perturbations and discuss viable parameter spaces in concrete models allowing the dark energy equation of state smaller than $-1$. The propagation speeds of those perturbations are subject to modifications beyond the domain of generalized Proca theories. There is a mixing between scalar and matter sound speeds, but such a mixing is suppressed during most of the cosmic expansion history without causing a new instability. On the other hand, we find that derivative interactions arising in beyond-generalized Proca theories give rise to important modifications to the cosmic growth history. The growth rate of matter perturbations can be compatible with the redshift-space distortion data due to the realization of gravitational interaction weaker than that in generalized Proca theories. Thus, it is possible to distinguish the dark energy model in beyond-generalized Proca theories from the counterpart in generalized Proca theories as well as from the $\Lambda$CDM model.

A non-abelian $SU(2)$ gauge field with a non-minimal Horndeski coupling to gravity gives rise to a de Sitter solution followed by a graceful exit to a radiation-dominated epoch. In this Horndeski Yang-Mills (HYM) theory we derive the second-order action for tensor perturbations on the homogeneous and isotropic quasi de Sitter background. We find that the presence of the Horndeski non-minimal coupling to the gauge field inevitably introduces ghost instabilities in the tensor sector during inflation. Moreover, we also find Laplacian instabilities for the tensor perturbations deep inside the Hubble radius during inflation. Thus, we conclude that the HYM theory does not provide a consistent inflationary framework due to the presence of ghosts and Laplacian instabilities.

In Gleyzes-Langlois-Piazza-Vernizzi (GLPV) scalar-tensor theories, which are outside the domain of second-order Horndeski theories, it is known that there exists a solid angle deficit singularity in the case where the parameter $\alpha_{\rm H}$ characterizing the deviation from Horndeski theories approaches a non-vanishing constant at the center of a spherically symmetric body. Meanwhile, it was recently shown that second-order generalized Proca theories with a massive vector field $A^{\mu}$ can be consistently extended to beyond-generalized Proca theories, which recover shift-symmetric GLPV theories in the scalar limit $A^{\mu} \to \nabla^{\mu} \chi$. In beyond-generalized Proca theories up to quartic-order Lagrangians, we show that solid angle deficit singularities are generally absent due to the existence of a temporal vector component. We also derive the vector-field profiles around a compact object and show that the success of the Vainshtein mechanism operated by vector Galileons is not prevented by new interactions in beyond generalized Proca theories.

We consider higher-order derivative interactions beyond second-order generalized Proca theories that propagate only the three desired polarizations of a massive vector field besides the two tensor polarizations from gravity. These new interactions follow the similar construction criteria to those arising in the extension of scalar-tensor Horndeski theories to Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theories. On the isotropic cosmological background, we show the existence of a constraint with a vanishing Hamiltonian that removes the would-be Ostrogradski ghost. We study the behavior of linear perturbations on top of the isotropic cosmological background in the presence of a matter perfect fluid and find the same number of propagating degrees of freedom as in generalized Proca theories (two tensor polarizations, two transverse vector modes, and two scalar modes). Moreover, we obtain the conditions for the avoidance of ghosts and Laplacian instabilities of tensor, vector, and scalar perturbations. We observe key differences in the scalar sound speed, which is mixed with the matter sound speed outside the domain of generalized Proca theories.

We consider the finite interactions of the generalized Proca theory including the sixth-order Lagrangian and derive the full linear perturbation equations of motion on the flat Friedmann-Lema\^{i}tre-Robertson-Walker background in the presence of a matter perfect fluid. By construction, the propagating degrees of freedom (besides the matter perfect fluid) are two transverse vector perturbations, one longitudinal scalar, and two tensor polarizations. The Lagrangians associated with intrinsic vector modes neither affect the background equations of motion nor the second-order action of tensor perturbations, but they do give rise to non-trivial modifications to the no-ghost condition of vector perturbations and to the propagation speeds of vector and scalar perturbations. We derive the effective gravitational coupling $G_{\rm eff}$ with matter density perturbations under a quasi-static approximation on scales deep inside the sound horizon. We find that the existence of intrinsic vector modes allows a possibility for reducing $G_{\rm eff}$. In fact, within the parameter space, $G_{\rm eff}$ can be even smaller than the Newton gravitational constant $G$ at the late cosmological epoch, with a peculiar phantom dark energy equation of state (without ghosts). The modifications to the slip parameter $\eta$ and the evolution of growth rate $f\sigma_8$ are discussed as well. Thus, dark energy models in the framework of generalized Proca theories can be observationally distinguished from the $\Lambda$CDM model according to both cosmic growth and expansion history. Furthermore, we study the evolution of vector perturbations and show that outside the vector sound horizon the perturbations are nearly frozen and start to decay with oscillations after the horizon entry.

In beyond-generalized Proca theories including the extension to theories higher than second order, we study the role of a spatial component $v$ of a massive vector field on the anisotropic cosmological background. We show that, as in the case of the isotropic cosmological background, there is no additional ghostly degrees of freedom associated with the Ostrogradski instability. In second-order generalized Proca theories we find the existence of anisotropic solutions on which the ratio between the anisotropic expansion rate $\Sigma$ and the isotropic expansion rate $H$ remains nearly constant in the radiation-dominated epoch. In the regime where $\Sigma/H$ is constant, the spatial vector component $v$ works as a dark radiation with the equation of state close to $1/3$. During the matter era, the ratio $\Sigma/H$ decreases with the decrease of $v$. As long as the conditions $|\Sigma| \ll H$ and $v^2 \ll \phi^2$ are satisfied around the onset of late-time cosmic acceleration, where $\phi$ is the temporal vector component, we find that the solutions approach the isotropic de Sitter fixed point ($\Sigma=0=v$) in accordance with the cosmic no-hair conjecture. In the presence of $v$ and $\Sigma$ the early evolution of the dark energy equation of state $w_{\rm DE}$ in the radiation era is different from that in the isotropic case, but the approach to the isotropic value $w_{\rm DE}^{{\rm (iso)}}$ typically occurs at redshifts $z$ much larger than 1. Thus, apart from the existence of dark radiation, the anisotropic cosmological dynamics at low redshifts is similar to that in isotropic generalized Proca theories. In beyond-generalized Proca theories the only consistent solution to avoid the divergence of a determinant of the dynamical system corresponds to $v=0$, so $\Sigma$ always decreases in time.

We consider a massive vector field with derivative interactions that propagates only the 3 desired polarizations (besides two tensor polarizations from gravity) with second-order equations of motion in curved space-time. The cosmological implications of such generalized Proca theories are investigated for both the background and the linear perturbation by taking into account the Lagrangian up to quintic order. In the presence of a matter fluid with a temporal component of the vector field, we derive the background equations of motion and show the existence of de Sitter solutions relevant to the late-time cosmic acceleration. We also obtain conditions for the absence of ghosts and Laplacian instabilities of tensor, vector, and scalar perturbations in the small-scale limit. Our results are applied to concrete examples of the general functions in the theory, which encompass vector Galileons as a specific case. In such examples, we show that the de Sitter fixed point is always a stable attractor and study viable parameter spaces in which the no-ghost and stability conditions are satisfied during the cosmic expansion history.

For a massive vector field with derivative self-interactions, the breaking of the gauge invariance allows the propagation of a longitudinal mode in addition to the two transverse modes. We consider generalized Proca theories with second-order equations of motion in a curved space-time and study how the longitudinal scalar mode of the vector field gravitates on a spherically symmetric background. We show explicitly that cubic-order self-interactions lead to the suppression of the longitudinal mode through the Vainshtein mechanism. Provided that the dimensionless coupling of the interaction is not negligible, this screening mechanism is sufficiently efficient to give rise to tiny corrections to gravitational potentials consistent with solar-system tests of gravity. We also study the quartic interactions with the presence of non-minimal derivative coupling with the Ricci scalar and find the existence of solutions where the longitudinal mode completely vanishes. Finally, we discuss the case in which the effect of the quartic interactions dominates over the cubic one and show that local gravity constraints can be satisfied under a mild bound on the parameters of the theory.

In Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theories, it is known that the conical singularity arises at the center of a spherically symmetric body ($r=0$) in the case where the parameter $\alpha_{{\rm H}4}$ characterizing the deviation from the Horndeski Lagrangian $L_4$ approaches a non-zero constant as $r \to 0$. We derive spherically symmetric solutions around the center in full GLPV theories and show that the GLPV Lagrangian $L_5$ does not modify the divergent property of the Ricci scalar $R$ induced by the non-zero $\alpha_{{\rm H}4}$. Provided that $\alpha_{{\rm H}4}=0$, curvature scalar quantities can remain finite at $r=0$ even in the presence of $L_5$ beyond the Horndeski domain. For the theories in which the scalar field $\phi$ is directly coupled to $R$, we also obtain spherically symmetric solutions inside/outside the body to study whether the fifth force mediated by $\phi$ can be screened by non-linear field self-interactions. We find that there is one specific model of GLPV theories in which the effect of $L_5$ vanishes in the equations of motion. We also show that, depending on the sign of a $L_5$-dependent term in the field equation, the model can be compatible with solar-system constraints under the Vainshtein mechanism or it is plagued by the problem of a divergence of the field derivative in high-density regions.

We propose a dark energy model where a scalar field $\phi$ has nonlinear self-interactions in the presence of a dilatonic coupling with the Ricci scalar. This belongs to a sub-class of theories beyond Horndeski, which accommodates covariant Galileons and Brans-Dicke theories as specific cases. We derive conditions under which the scalar sound speed squared $c_{\rm s}^2$ is positive from the radiation era to today. Since $c_{\rm s}^2$ remains to be smaller than the order of 1, the deviation from Horndeski theories does not cause heavy oscillations of gauge-invariant gravitational potentials. In this case, the evolution of matter perturbations at low redshifts is similar to that in the coupled dark energy scenario with an enhanced gravitational interaction. On the spherically symmetric background with a matter source, the existence of field self-interactions suppresses the propagation of fifth force inside a Vainshtein radius. We estimate an allowed parameter space in which the model can be compatible with solar system constraints while satisfying conditions for the cosmological viability of background and perturbations.

In a class of Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theories, we derive both vacuum and interior Schwarzschild solutions under the condition that the derivatives of a scalar field $\phi$ with respect to the radius $r$ vanish. If the parameter $\alpha_{\rm H}$ characterizing the deviation from Horndeski theories approaches a non-zero constant at the center of a spherically symmetric body, we find that the conical singularity arises at $r=0$ with the Ricci scalar given by $R=-2\alpha_{\rm H}/r^2$. This originates from violation of the geometrical structure of four-dimensional curvature quantities. The conical singularity can disappear for the models in which the parameter $\alpha_{\rm H}$ vanishes in the limit that $r \to 0$. We propose explicit models without the conical singularity by properly designing the classical Lagrangian in such a way that the main contribution to $\alpha_{\rm H}$ comes from the field derivative $\phi'(r)$ around $r=0$. We show that the extension of covariant Galileons with a diatonic coupling allows for the recovery of general relativistic behavior inside a so-called Vainshtein radius. In this case, both the propagation of a fifth force and the deviation from Horndeski theories are suppressed outside a compact body in such a way that the model is compatible with local gravity experiments inside the solar system.

We review the effective field theory of modified gravity in which the Lagrangian involves three dimensional geometric quantities appearing in the 3+1 decomposition of space-time. On the flat isotropic cosmological background we expand a general action up to second order in the perturbations of geometric scalars, by taking into account spatial derivatives higher than two. Our analysis covers a wide range of gravitational theories-- including Horndeski theory/its recent generalizations and the projectable/non-projectable versions of Ho\v{r}ava-Lifshitz gravity. We derive the equations of motion for linear cosmological perturbations and apply them to the calculations of inflationary power spectra as well as the dark energy dynamics in Galileon theories. We also show that our general results conveniently recover stability conditions of Ho\v{r}ava-Lifshitz gravity already derived in the literature.

We consider perturbations of a static and spherically symmetric background endowed with a metric tensor and a scalar field in the framework of the effective field theory of modified gravity. We employ the previously developed 2+1+1 canonical formalism of a double Arnowitt-Deser-Misner (ADM) decomposition of space-time, which singles out both time and radial directions. Our building block is a general gravitational action that depends on scalar quantities constructed from the 2+1+1 canonical variables and the lapse. Variation of the action up to first-order in perturbations gives rise to three independent background equations of motion, as expected from spherical symmetry. The dynamical equations of linear perturbations follow from the second-order Lagrangian after a suitable gauge fixing. We derive conditions for the avoidance of ghosts and Laplacian instabilities for the odd-type perturbations. We show that our results not only incorporates those derived in the most general scalar-tensor theories with second-order equations of motion (the Horndeski theories) but they can be applied to more generic theories beyond Horndeski.

We study the cosmology of an extended version of Horndeski theories with second-order equations of motion on the flat Friedmann-Lema\^{i}tre-Robertson-Walker (FLRW) background. In addition to a dark energy field $\chi$ associated with the gravitational sector, we take into account multiple scalar fields $\phi_I$ ($I=1,2\cdots,N-1$) characterized by the Lagrangians $P^{(I)}(X_I)$ with $X_I=\partial_{\mu}\phi_I\partial^{\mu}\phi_I$. These additional scalar fields can model the perfect fluids of radiation and non-relativistic matter. We derive propagation speeds of scalar and tensor perturbations as well as conditions for the absence of ghosts. The theories beyond Horndeski induce non-trivial modifications to all the propagation speeds of $N$ scalar fields, but the modifications to those for the matter fields $\phi_I$ are generally suppressed relative to that for the dark energy field $\chi$. We apply our results to the covariantized Galileon with an Einstein-Hilbert term in which partial derivatives of the Minkowski Galileon are replaced by covariant derivatives. Unlike the covariant Galileon with second-order equations of motion in general space-time, the scalar propagation speed square $c_{s1}^2$ associated with the field $\chi$ becomes negative during the matter era for late-time tracking solutions, so the two Galileon theories can be clearly distinguished at the level of linear cosmological perturbations.

We study how the Vainshtein mechanism operates in the most general scalar-tensor theories with second-order equations of motion. The field equations of motion, which can be also applicable to most of other screening scenarios proposed in literature, are generally derived in a spherically symmetric space-time with a matter source. In the presence of a field coupling to the Ricci scalar, we clarify conditions under which the Vainshtein mechanism is at work in a weak gravitational background. We also obtain the solutions of the field equation inside a spherically symmetric body and show how they can be connected to exterior solutions that accommodate the Vainshtein mechanism. We apply our general results to a number of concrete models such as the covariant/extended Galileons and the DBI Galileons with Gauss-Bonnet and other terms. In these models the fifth force can be suppressed to be compatible with solar-system constraints, provided that non-linear field kinetic terms coupled to the Einstein tensor do not dominate over other non-linear field self-interactions.

In second-order scalar-tensor theories we study how the Vainshtein mechanism works in a spherically symmetric background with a matter source. In the presence of the field coupling $F(\phi)=e^{-2Q\phi}$ with the Ricci scalar $R$ we generally derive the Vainshtein radius within which the General Relativistic behavior is recovered even for the coupling $Q$ of the order of unity. Our analysis covers the models such as the extended Galileon and Brans-Dicke theories with a dilatonic field self-interaction. We show that, if these models are responsible for the cosmic acceleration today, the corrections to gravitational potentials are generally small enough to be compatible with local gravity constraints.

We study the evolution of matter density perturbations in Galileon cosmology where the late-time cosmic acceleration can be realized by a field kinetic energy. We obtain full perturbation equations at linear order in the presence of five covariant Lagrangians ${cal L}_i$ ($i=1,...,5$) satisfying a Galilean symmetry in the flat space-time. The equations for a matter perturbation as well as an effective gravitational potential are derived under a quasi-static approximation on sub-horizon scales. This approximation can reproduce full numerical solutions with high accuracy for the wavelengths relevant to large-scale structures. For the model parameters constrained by the background expansion history of the Universe the growth rate of matter perturbations is larger than that in the LCDM model, with the growth index $gamma$ today typically smaller than 0.4. We also find that, even on very large scales associated with the Integrated-Sachs-Wolfe (ISW) effect in Cosmic Microwave Background (CMB) temperature anisotropies, the effective gravitational potential exhibits a temporal growth during the transition from the matter era to the epoch of cosmic acceleration. These properties are useful to distinguish the Galileon model from the LCDM in future high-precision observations.