Vainshtein Mechanism Research Articles

Scalar radiation with a quartic Galileon

The class of Galileon scalar fields theories encapsulate the Vainshtein screening mechanism, which is characteristic of a large range of infrared modified theories of gravity. Such theories can lead to testable departures from general relativity through fifth forces and new scalar modes of gravitational radiation. However, the inherent nonlinearity of the Vainshtein mechanism has limited analytic attempts to describe Galileon theories with both cubic and quartic interactions. To improve on this, we perform direct numerical simulations of the quartic Galileon model for a rotating binary source and infer the power spectrum of given multipoles. To tame numerical instabilities we utilize a low-pass filter, extending previous work on the cubic Galileon. Our findings show that the multipole expansion is well defined and under control. Moreover, our results confirm that despite being a nonlinear scalar, the dominant Galileon radiation is quadrupole, and we find a new scaling behavior deep inside the Vainshtein region. Published by the American Physical Society 2024

Revisiting Vainshtein screening for fast N-body simulations

We revisit a method to incorporate the Vainshtein screening mechanism in N-body simulations proposed by R. Scoccimarro in [1]. We further extend this method to cover a subset of Horndeski theories that evade the bound on the speed of gravitational waves set by the binary neutron star merger GW170817. The procedure consists of the computation of an effective gravitational coupling that is time and scale dependent, G eff (k,z), where the scale dependence will incorporate the screening of the fifth-force. This is a fast procedure that when contrasted to the alternative of solving the full equation of motion for the scalar field inside N-body codes, reduces considerably the computational time and complexity required to run simulations. To test the validity of this approach in the non-linear regime, we have implemented it in a COmoving Lagrangian Approximation (COLA) N-body code, and ran simulations for two gravity models that have full N-body simulation outputs available in the literature, nDGP and Cubic Galileon. We validate the combination of the COLA method with this implementation of the Vainshtein mechanism with full N-body simulations for predicting the boost function: the ratio between the modified gravity non-linear matter power spectrum and its General Relativity counterpart. This quantity is of great importance for building emulators in beyond-ΛCDM models, and we find that the method described in this work has an agreement of below 2% for scales down to k ≈ 3h/Mpc with respect to full N-body simulations.

Vainshtein screening in Horndeski theories nonminimally and kinetically coupled to ordinary matter

We study the Vainshtein screening mechanism in Horndeski theories in the presence of a scalar field $\phi$ nonminimally and kinetically coupled to ordinary matter field. A general interacting Lagrangian describing this coupling is characterized by energy transfer $f_1$ and momentum exchange $f_2$. For a spherically symmetric configurations on top of the cosmological background, we investigate the static perturbations in linear and nonlinear regimes with respect to the scalar field perturbation. In the former regime, the parametrized post-Newtonian parameter generally deviates from unity as long as the matter coupling or $G_{4,\phi}$ exists. On the other hand, in the latter regime, we show that the nonlinear self-interaction term of scalar field successfully activates the Vainshtein mechanism even in the presence of the couplings $f_1$ and $f_2$. The gravitational potentials recover the Newtonian behavior deep inside the Vainshtein radius. The bounds on coupling terms not to substantially change the Vainshtein radius are also given.

Fast full N-body simulations of generic modified gravity: derivative coupling models

We present mg-glam, a code developed for the very fast production of full N-body cosmological simulations in modified gravity (MG) models. We describe the implementation, numerical tests and first results of a large suite of cosmological simulations for two broad classes of MG models with derivative coupling terms — the Vainshtein- and Kmouflage-type models — which respectively features the Vainshtein and Kmouflage screening mechanism. Derived from the parallel particle-mesh code glam, mg-glam incorporates an efficient multigrid relaxation technique to solve the characteristic nonlinear partial differential equations of these models. For Kmouflage, we have proposed a new algorithm for the relaxation solver, and run the first simulations of the model to understand its cosmological behaviour. In a companion paper, we describe versions of this code developed for conformally-coupled MG models, including several variants of f(R) gravity, the symmetron model and coupled quintessence. Altogether, mg-glam has so far implemented the prototypes for most MG models of interest, and is broad and versatile. The code is highly optimised, with a tremendous (over two orders of magnitude) speedup when comparing its running time with earlier N-body codes, while still giving accurate predictions of the matter power spectrum and dark matter halo abundance. mg-glam is ideal for the generation of large numbers of MG simulations that can be used in the construction of mock galaxy catalogues and accurate emulators for ongoing and future galaxy surveys.

Covariant effective action for generalized Proca theories

We investigate the quantum stability of generalised Proca theories in curved spacetime treating gravity as a dynamical field. To compute the quantum gravitational corrections, we evaluate the divergent part of the effective action at one-loop level using Vilkovisky-DeWitt formalism, which gives us a gauge invariant and gauge condition independent effective action. It is shown that the quantum corrections are suppressed by a UV cutoff parametrically higher than the Proca mass, if the coupling constants are restricted to lie in a certain range. Furthermore, it is shown that the quantum corrections remain suppressed even at scales where classical non-linearities dominate over kinetic terms, allowing Vainshtein mechanism to work.

Vainshtein mechanism in Generalised Massive Gravity

We present a non-linear analysis of perturbations around cosmological solutions in Generalised Massive gravity. This Lorentz invariant theory is an extension of de Rham, Gabadadze, Tolley massive gravity that propagates 5 degrees of freedom while allowing stable open FLRW cosmologies.For a minimal model that supports a self-accelerating background, we study the dynamics of non-linear perturbations. We find that the equation for the scalar graviton is distinct from the analogues in Horndeski and DHOST theories. We numerically solve the equation to find a new type of nonlinear solution for the scalar mode, and confirm the presence of a Vainshtein screening mechanism. We show that the PPN parameter approaches its GR value at solar system scales and satisfies the current bounds.

Fast generation of mock galaxy catalogues in modified gravity models with COLA

We investigate the viability of producing galaxy mock catalogues with COmoving Lagrangian Acceleration (COLA) simulations in Modified Gravity (MG) models employing the Halo Occupation Distribution (HOD) formalism. In this work, we focus on two theories of MG: f(R) gravity with the chameleon mechanism, and a braneworld model (nDGP) that incorporates the Vainshtein mechanism. We use a suite of full N-body simulations in MG as a benchmark to test the accuracy of COLA simulations. At the level of Dark Matter (DM), we show that COLA accurately reproduces the matter power spectrum up to k ∼ 1 h Mpc-1, while it is less accurate in reproducing the velocity field. To produce halo catalogues, we find that the rockstar halo-finder does not perform well with COLA simulations. On the other hand, using a simple Friends-of-Friends (FoF) finder and an empirical mass conversion from FoF to spherical over-density masses, we are able to produce halo catalogues in COLA that are in good agreement with those in N-body simulations. To consider the effects of the MG fifth force on the halo profile, we derive simple fitting formulae for the concentration-mass and the velocity dispersion-mass relations that we calibrate using rockstar halo catalogues in N-body simulations. We then use these results to extend the HOD formalism to modified gravity simulations in COLA. We use an HOD model with five parameters that we tune to obtain galaxy catalogues in redshift space. We find that despite the great freedom of the HOD model, MG leaves characteristic imprints in the redshift space power spectrum multipoles and these features are well captured by the COLA galaxy catalogues.

Vainshtein screening in bimetric cosmology

We demonstrate that the instabilities in linear cosmological perturbations in bimetric theory are the manifestation of the non-linear Vainshtein mechanism on an FRW background. The spin-2 mass serves as the cosmological Vainshtein scale in this case. This allows us to quantitatively address early universe cosmology. In particular, in a global analysis, we study data from the cosmic microwave background radiation and local measurements of the Hubble flow. We show that bimetric cosmology resolves the discrepancy in the local and early-time measurements of the Hubble scale via an effective phantom dark energy component.

Vainshtein screening for slowly rotating stars

We study the Vainshtein mechanism in the context of slowly rotating stars in scalar-tensor theories. While the Vainshtein screening is well established for spherically symmetric spacetimes, we examine its validity in the axisymmetric case for slowly rotating sources. We show that the deviations from the general relativity solution are small in the weak-field approximation outside the star: the solution for the frame-dragging function is the same as in general relativity at leading order. Moreover, in most cases the corrections are suppressed by powers of the Vainshtein radius provided that the screening operates in spherical symmetry. Outside the Vainshtein radius, the frame dragging function receives corrections that are not suppressed by the Vainshtein radius, but which are still subleading. This suggests that the Vainshtein mechanism in general can be extended to slowly rotating stars and that it works analogously to the static case inside the Vainshtein radius. We also study relativistic stars and show that for some theories the frame-dragging function in vacuum does not receive corrections at all, meaning that the screening is perfect outside the star.

New graviton mass bound from binary pulsars

In Einstein's general relativity, gravity is mediated by a massless metric field. The extension of general relativity to consistently include a mass for the graviton has profound implications for gravitation and cosmology. Salient features of various massive gravity theories can be captured by Galileon models, the simplest of which is the cubic Galileon. The presence of the Galileon field leads to additional gravitational radiation in binary pulsars where the Vainshtein mechanism is less suppressed than its fifth-force counterpart, which deserves a detailed confrontation with observations. We prudently choose fourteen well-timed binary pulsars, and from their intrinsic orbital decay rates we put a new bound on the graviton mass, $m_g \lesssim 2 \times 10^{-28}\,{\rm eV}/c^2$ at the 95% confidence level, assuming a flat prior on $\ln m_g$. It is equivalent to a bound on the graviton Compton wavelength $\lambda_g \gtrsim 7 \times 10^{21}\,{\rm m}$. Furthermore, we extensively simulate times of arrival for pulsars in orbit around stellar-mass black holes and the supermassive black hole at the Galactic center, and investigate their prospects in probing the cubic Galileon theory in the near future.

Testing Brans-Dicke gravity with screening by scalar gravitational wave memory

The Brans-Dicke theory of gravity is one of the oldest ideas to extend general relativity by introducing a nonminimal coupling between the scalar field and gravity. The Solar System tests put tight constraints on the theory. In order to evade these constraints, various screening mechanisms have been proposed. These screening mechanisms allow the scalar field to couple to matter as strongly as gravity in low density environments while suppressing it in the Solar System. The Vainshtein mechanism, which is found in various modified gravity models such as massive gravity, braneworld models and scalar tensor theories, suppresses the scalar field efficiently in the vicinity of a massive object. This makes it difficult to test these theories from gravitational wave observations. We point out that the recently found scalar gravitational wave memory effect, which is caused by a permanent change in spacetime geometry due to the collapse of a star to a back hole can be significantly enhanced in the Brans-Dicke theory of gravity with the Vainshtein mechanism. This provides a possibility to detect scalar gravitational waves by a network of three or more gravitational wave detectors.

Relativistic stars in a cubic Galileon universe

We study relativistic stars in Hordenski theories that evade the gravitational wave constraints and exhibit the Vainshtein mechanism, focusing on a model based on the cubic Galileon Lagrangian. We derive the scalar field profile for static spherically symmetric objects in asymptotically de Sitter space-time with a linear time dependence. The exterior solution matches to the black hole solution found in the literature. Due to the Vainshtein mechanism, the stellar structure is indistinguishable from that of General Relativity with the same central density as long as the radius of the star is shorter than the Vainshtein radius. On the other hand, the scalar field is not suppressed beyond the Vainshtein radius. These solutions have an additional integration constant in addition to the mass of the star.

Resolving the van Dam-Veltman-Zakharov and strong coupling problems in massive gravity and bigravity

As is well known, both massive gravity and bigravity exhibit the linear van Dam-Veltman-Zakharov (vDVZ) discontinuity that is cured classically by the nonlinear Vainshtein mechanism due to certain low scale strongly coupled interactions. Here we show how both the vDVZ and strong coupling problems can be removed by embedding 4D covariant massive gravity into a certain 5D warped geometry. The 4D theory is a nonlinear strongly coupled massive gravity, that is being coupled to a 5D bulk theory that generates a bulk graviton mass via a one loop diagram. This induced mass leads to an additional 4D kinetic term for the 4D longitudinal mode, even on flat space. Due to this kinetic term the 4D massive theory becomes weakly coupled all the way up to a high energy scale set by the bulk cosmological constant. The same effect leads to a suppression of the interactions of the 4D longitudinal mode with a 4D matter stress-tensor, thus removing the vDVZ discontinuity. The proposed mechanism has a pure 4D holographic interpretation: a 4D nonlinear massive gravity mixes to a non-conserved symmetric tensor of a 4D CFT that has a cutoff; the latter mixing generates a large kinetic term for the longitudinal mode, and this makes the longitudinal mode be weakly coupled to a matter stress-tensor, and weakly self-coupled, all the way up to the scale of the 4D CFT cutoff.

Screened and unscreened solutions for relativistic star in de Rham-Gabadadze-Tolley massive gravity

We study the static and spherical symmetric (SSS) configurations in the non-minimal model of the de Rham-Gabadadze-Tolley (dRGT) massive gravity with a flat reference metric. Considering the modified Tolman-Oppenheimer-Volkoff (TOV) equation, the Bianchi identity, and energy-momentum conservation, we find a new algebraic equation for the radial coordinate of the reference metric. We demonstrate that this equation suggests an absence of the Vainshtein mechanism in the minimal model of the dRGT massive gravity, while it has two branches of solutions where one connects with the Schwarzschild space-time and another implies the significant deviation from the asymptotically flat space-time in the non-minimal model. We also briefly discuss the boundary conditions for the relativistic stars in the dRGT massive gravity and a potential relation with the mass-radius relation of the stars.

Cosmological dynamics and double screening of DBI-Galileon gravity

We investigate cosmological dynamics and screening mechanism of the Dirac-Born-Infeld (DBI) Galileon model. The model has been divided into two regimes, one has positive signs in front of scalar field kinetic terms so-called the DBI galileon, another one has negative signs and it is dubbed as the DBIonic galileon. We find de Sitter solution and evolution of the Universe starting from radiation dominated era to late-time accelerated expansion in the DBI galileon model without the presence of potential term. In one of the attractors, the ghost and Laplacian instabilities vanishes for the whole evolution. We find mixing of screening mechanisms between the Vainshtein mechanism and the DBIonic screening mechanism in the DBIonic galileon model, in which a scale changing between these two mechanisms depends on a mass of a source.

Lunar laser ranging constraints on nonminimally coupled dark energy and standard sirens

In dark energy models where a scalar field $\phi$ is coupled to the Ricci scalar $R$ of the form $e^{-2Q (\phi-\phi_0)/M_{\rm pl}}R$, where $Q$ is a coupling constant, $\phi_0$ is today's value of $\phi$, and $M_{\rm pl}$ is the reduced Planck mass, we study how the recent Lunar Laser Ranging (LLR) experiment places constraints on the nonminimal coupling from the time variation of gravitational coupling. Besides a potential of the light scalar responsible for cosmic acceleration, we take a cubic Galileon term into account to suppress fifth forces in over-density regions of the Universe. Even if the scalar-matter interaction is screened by the Vainshtein mechanism, the time variation of gravitational coupling induced by the cosmological background field $\phi$ survives in the solar system. For a small Galileon coupling constant $\beta_3$, there exists a kinetically driven $\phi$-matter-dominated-epoch ($\phi$MDE) prior to cosmic acceleration. In this case, we obtain the stringent upper limit $Q \le 3.4 \times 10^{-3}$ from the LLR constraint. For a large $\beta_3$ without the $\phi$MDE, the coupling $Q$ is not particularly bounded from above, but the cosmological Vainshtein screening strongly suppresses the time variation of $\phi$ such that the dark energy equation of state $w_{\rm DE}$ reaches the value close to $-1$ at high redshifts. We study the modified gravitational wave propagation induced by the nonminimal coupling to gravity and show that, under the LLR bound, the difference between the gravitational wave and luminosity distances does not exceed the order $10^{-5}$ over the redshift range $0<z<100$. In dark energy models where the Vainshtein mechanism is at work through scalar derivative self-interactions, it is difficult to probe the signature of nonminimal couplings from the observations of standard sirens.

Screenings in modified gravity: a perturbative approach

We present a formalism to study screening mechanisms in modified theories of gravity through perturbative methods in different cosmological scenarios. We consider Einstein-frame posed theories that are recast as Jordan-frame theories, where a known formalism is employed, although the resulting nonlinearities of the Klein–Gordon equation acquire an explicit coupling between matter and the scalar field, which is absent in Jordan-frame theories. The obtained growth functions are then separated into screening and non-screened contributions to facilitate their analysis. This allows us to compare several theoretical models and to recognize patterns that can be used to distinguish models and their screening mechanisms. In particular, we find anti-screening features in the symmetron model. In contrast, chameleon-type theories in both the Jordan and Einstein frames always present a screening behaviour. Up to third order in perturbation, we find no anti-screening behaviour in theories with a Vainshtein mechanism, such as the Dvali Gabadadze Porrati braneworld model and the cubic Galileon.

Anti-screening of the Galileon force around a disk center hole

The Vainshtein mechanism is known as an efficient way of screening the fifth force around a matter source in modified gravity. This has been verified mainly in highly symmetric matter configurations. To study how the Vainshtein mechanism works in a less symmetric setup, we numerically solve the scalar field equation around a disk with a hole at its center in the cubic Galileon theory. We find, surprisingly, that the Galileon force is enhanced, rather than suppressed, in the vicinity of the hole. This anti-screening effect is larger for a thinner, less massive disk with a smaller hole. At this stage, our setup is only of academic interest and its astrophysical consequences are unclear, but this result implies that the Vainshtein screening mechanism around less symmetric matter configurations are quite nontrivial.

Scalar gravitational radiation from binaries: Vainshtein mechanism in time-dependent systems

We develop a full four-dimensional numerical code to study scalar gravitational radiation emitted from binary systems and probe the Vainshtein mechanism in situations that break the static and spherical symmetry, relevant for binary pulsars as well as black holes and neutron stars binaries. The present study focuses on the cubic Galileon which arises as the decoupling limit of massive theories of gravity. Limitations associated with the numerical methods prevent us from reaching a physically realistic hierarchy of scales; nevertheless, within this context we observe the same power law scaling of the radiated power as previous analytic estimates, and confirm a strong suppression of the power emitted in the monopole and dipole as compared with quadrupole radiation. Following the trend to more physically realistic parameters, we confirm the suppression of the power emitted in scalar gravitational radiation and the recovery of general relativity with good accuracy. This paves the way for future numerical work, probing more generic, physically relevant situations and sets of interactions that may exhibit the Vainshtein mechanism.

Horndeski gravity without screening in binary pulsars

We test the subclasses of Horndeski gravity without Vainshtein mechanism in the strong field regime of binary pulsars. We find the rate of energy losses via the gravitational radiation predicted by such theories and compare our results with observational data from quasi-circular binaries PSR J1738+0333, PSR J0737-3039, PSR J1012+5307. In addition, we consider few specific cases: the hybrid metric-Palatini f(R)-gravity and massive Brans-Dicke theory.