Next Article in Journal
Simulation of Closed Timelike Curves in a Darwinian Approach to Quantum Mechanics
Next Article in Special Issue
Unimodular Approaches to the Cosmological Constant Problem
Previous Article in Journal
Tests of Charge–Parity Symmetry and Lepton Flavor Conservation in the Top Quark Sector
Previous Article in Special Issue
Particles of a de Sitter Universe
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Modified Gravity Approaches to the Cosmological Constant Problem

by
Foundational Aspects of Dark Energy (FADE) Collaboration
,
Heliudson Bernardo
1,
Benjamin Bose
2,3,4,
Guilherme Franzmann
4,5,
Steffen Hagstotz
6,7,
Yutong He
5,8,
Aliki Litsa
8 and
Florian Niedermann
5,*
1
Department of Physics, Ernest Rutherford Physics Building, McGill University, 3600 Rue Université, Montréal, QC H3A 2T8, Canada
2
Département de Physique Théorique, Université de Genève, 24 quai Ernest Ansermet, 1211 Genève, Switzerland
3
Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK
4
Basic Research Community for Physics e.V., Veddeler Brückenstrasse 122, 20539 Hamburg, Germany
5
Nordita, KTH Royal Institute of Technology and Stockholm University, Hannes Alfvéns väg 12, SE-106 91 Stockholm, Sweden
6
Universitäts-Sternwarte, Fakultät für Physik, Ludwig-Maximilians Universität München, Scheinerstraße 1, D-81679 München, Germany
7
Excellence Cluster ORIGINS, Boltzmannstraße 2, D-85748 Garching, Germany
8
The Oskar Klein Centre for Cosmoparticle Physics, Stockholm University Roslagstullsbacken 21A, SE-106 91 Stockholm, Sweden
*
Author to whom correspondence should be addressed.
Universe 2023, 9(2), 63; https://doi.org/10.3390/universe9020063
Submission received: 31 October 2022 / Revised: 4 January 2023 / Accepted: 13 January 2023 / Published: 20 January 2023
(This article belongs to the Special Issue Cosmological Constant)

Abstract

:
The cosmological constant and its phenomenology remain among the greatest puzzles in theoretical physics. We review how modifications of Einstein’s general relativity could alleviate the different problems associated with it that result from the interplay of classical gravity and quantum field theory. We introduce a modern and concise language to describe the problems associated with its phenomenology, and inspect no-go theorems and their loopholes to motivate the approaches discussed here. Constrained gravity approaches exploit minimal departures from general relativity; massive gravity introduces mass to the graviton; Horndeski theories lead to the breaking of translational invariance of the vacuum; and models with extra dimensions change the symmetries of the vacuum. We also review screening mechanisms that have to be present in some of these theories if they aim to recover the success of general relativity on small scales as well. Finally, we summarize the statuses of these models in their attempts to solve the different cosmological constant problems while being able to account for current astrophysical and cosmological observations.

 

Contents
 
1.      Introduction...................................................................................................................................................................................... 2
 
2.      What Is the Problem after All?............................................................................................................................................................ 4
 
3.      How to Modify Gravity.......................................................................................................................................................... 5
3.1.      Self-Tuning................................................................................................................................................................. 5
            3.1.1.      A Failed Start......................................................................................................................................................6
            3.1.2.      Weinberg’s Argument......................................................................................................................................... 8
            3.1.3.      Beyond Weinberg’s Argument.................................................................................................................................. 9
            3.1.4.      Loopholes..................................................................................................................................................... 11
3.2.      Screening Mechanisms........................................................................................................................................................ 12
            3.2.1.      Potential Screening........................................................................................................................................... 13
            3.2.2.      Derivative Screening.......................................................................................................................................... 15
 
4.      Modified Gravity Approaches............................................................................................................................................................ 17
4.1.      What Do We Want?............................................................................................................................................................ 17
4.2.      Constraining Gravity........................................................................................................................................................ 17
            4.2.1.      The Global Vacuum Energy Sequestering............................................................................................................18
            4.2.2.      Local Sequestering............................................................................................................................................ 19
            4.2.3.      Nonlocal Approach............................................................................................................................................. 20
            4.2.4.      Unimodular Gravity............................................................................................................................................ 22
4.3.      Massive Gravity............................................................................................................................................................. 23
            4.3.1.      Degravitation................................................................................................................................................. 23
            4.3.2.      Linear Massive Gravity........................................................................................................................................ 24
            4.3.3.      Nonlinear Theories............................................................................................................................................ 25
4.4.      Self-Tuning with Horndeski Theories.......................................................................................................................... 25
            4.4.1.      The Fab-4..................................................................................................................................................... 27
            4.4.2.      Well-Tempered Self-Tuning..................................................................................................................................... 27
            4.4.3.      Outlook.................................................................................................................................................................... 28
4.5.      Braneworld Models............................................................................................................................................................... 28
            4.5.1.      Our Universe as a Cosmic String in Six Dimensions....................................................................................... 29
            4.5.2.      Supersymmetric Large Extra Dimensions............................................................................................................. 31
            4.5.3.      Outlook.................................................................................................................................................................. 32
 
5.      Conclusions—You Can’t Always Get What You Want, or Can You?........................................................................................ 33
 
References................................................................................................................................................................................................... 37

1. Introduction

The cosmological constant problem (CCP) is one of the most persistent puzzles in theoretical physics. It appears at the interface between quantum mechanics and gravity, and it seemingly contradicts one of the major building blocks of modern physics, which is that scales in nature decouple. We hope that it is possible to understand physics at low energies without a detailed knowledge of the ultraviolet (UV)-complete theory. Formally, this decoupling of scales is expressed through the enormously successful framework of effective field theories (EFT), which is challenged by the CCP.
Viewed solely as a metric theory for spacetime, Einstein’s general relativity (GR) allows for a cosmological constant Λ that is left unfixed by the principles and symmetries of the theory. Such a parameter may compete with the energy-momentum tensor (EMT) of local sources in the equations of motion to produce certain solutions, as was the case of Einstein’s cosmological static solution [1,2]. In fact, the weak strength of the gravitational interaction is such that all local tests of GR are compatible with Λ = 0 , singling out cosmological observations to constrain Λ .
The current cosmological concordance model, Λ CDM, is fully based on classical GR and employs a cosmological constant (CC) to describe the present period of accelerated expansion necessary to fit both early- [3] and late-time (e.g., [4,5,6,7,8,9,10], and see [11] for a recent review) observational data. Although there are other cosmological models in which the driver of the current epoch of accelerated expansion, herein referred to as dark energy (DE), is not described by a cosmological constant, standard GR with a CC is still the preferred model once many datasets are analyzed together [11].
However, problems appear once we consider the standard model of particle physics (SM) together with GR. The SM is based on quantum field theory (QFT) in Minkowski spacetime. Within the QFT paradigm, the local energy-momentum tensor of field systems might be nonzero even in the vacuum state, in which case they have the same structure as the EMT of a CC term [12]. Hence, the existence of the SM fields, inferred from local experiments, implies that there are quantum-matter-induced contributions to the cosmological constant. This means that if DE is a CC, then the Λ -fitted value (also known as effective cosmological constant) also includes contributions from the SM quantum fields, generically referred to as the vacuum energy density, ρ vac .
Typically, any vacuum energy contribution is enormous and would completely dominate the gravitational dynamics of the Universe. This was first noted by Nernst in 1916 [13], and Pauli famously quipped that “the Universe would not even reach to the moon” [14]1 given the large vacuum energy from quantum fields. The problem was formally analyzed in a cosmological context by Zel’dovich in the 1970s [19,20]. By the late 1980s, due to the perturbative UV incompleteness of GR [21,22], it was expected that a quantum version of the theory would fix Λ . Prior to the discovery of the accelerated cosmological expansion, popular attempts to achieve this were based on supergravity theories, which typically do not allow a positive cosmological constant. The seminal review [23] summarizes the status of the problem before the early 1990s.
Perhaps the most straightforward procedure for solving the mismatch is to postulate a bare, purely classical contribution inherent to GR to the cosmological constant, Λ b, that would cancel part of ρ vac / M Pl 2 . However, the modern view of QFT within the effective field theory approach is that such a proposal is radiatively unstable or extremely fine-tuned, as discussed later. There are many reviews on other approaches (see, for instance, [24,25]), and the present work aims at reviewing modified gravity theories to alleviate or solve the cosmological constant problem. As we will discuss in this review, the problem actually has different facets, and it is a priori unclear whether all of them can be addressed at once or not.
Since the problem occurs when trying to calculate the gravitational effect of QFT vacuum energies, one approach to tackle the issue is to modify gravity itself. As shown in Section 2 and Section 3, there is little space to solve the CCP within classical GR. Setting aside anthropic arguments and discussions about fine-tuning [26,27,28,29,30,31,32,33], modifying GR not only offers a new view on the CCP but also provides new phenomenological signatures. Generally speaking, modified gravity theories introduce extra fields and/or constraints on GR such that the coupling of matter to the metric is modified and/or there are extra universal couplings with the extra fields in a way that is consistent with observations. It is natural then to ask if and how the CCP manifests itself in such theories, and this is one of the motivations for developing them.
Another motivation for studying modified gravity approaches is the fact that some of these theories are typical examples of how to evade Weinberg’s no-go theorem, as discussed in Section 3.1.2, and its complement, reviewed in Section 3.1.3. We highlight possible loopholes in the assumptions of the no-go theorems in Section 3.1.4 and use them as a systematic guiding principle for the modified gravity theories considered throughout the review. The idea of this review is not to cover all the modified gravity literature. Instead, we make a narrow choice of models that we deem promising and exemplify the discussed loopholes with an emphasis on how the CCP is addressed in each one of them. For other reviews on modified gravity and the cosmological constant issues, see [23,24,34,35,36,37,38,39,40,41,42].
We start by carefully introducing the different contributions and aspects of the cosmological constant problem in Section 2. Any possible solution to the cosmological constant problem is severely constrained by powerful no-go theorems, which we recap in Section 3 together with potential ways to circumvent the theorems while taking into consideration various observational constraints. In Section 4.1, we summarize the task at hand from both a theoretical and phenomenological point of view before diving into various approaches in Section 4. We discuss the extent to which the various modified gravity theories can solve parts of the problem in Section 5.
Conventions: Unless otherwise stated, we set c = = 1 and use the metric signature ( + + + ) . The gravitational coupling is denoted by M Pl 2 in all sections, apart from Section 4.2.1 and Section 4.2.2 on sequestering, where κ 2 and κ 2 ( x ) are used to indicate the promotion of the Planck mass to a variable. We use L m or S m for a generic matter component, and only specify its content as L ( . . . ) or S ( . . . ) when necessary.
Disclaimer: This review does not intend, in any way, to present an exhaustive collection of modified gravity approaches. Instead, the models were chosen to exemplify the different ways one can bypass no-go theorems surrounding the cosmological constant problem, and they are subjected to the authors’ personal biases and preferences on the subject.

2. What Is the Problem after All?

Much discussion has been focused on the CCP in the past. There is some confusion in the literature about what the problem is, as there are different distinct issues that need to be tackled. To make it clear what the challenges are, and how modifications of gravity can help sort them out, we briefly review them below:
  • The weight of vacuum. The gravitating vacuum energy at the level of the classical Einstein equations receives contributions from the vacuum energy of the fields in the SM, ρ vac . QFT calculations of the latter for a given energy scale μ , defined by the renormalization scale at which we trust our theories, seem to result in a very high value for the vacuum energy that scales as m 4 ln ( m 2 / μ 2 ) , where m is the mass of the particle for which the vacuum energy is being computed [24]. This typically differs by many orders of magnitude from the actual value associated with the accelerated expansion of the Universe at cosmological scales, which is roughly ( 10 3 eV ) 4 (the mass of the top quark, for instance, is roughly 10 11 eV , leading to a 56-orders-of-magnitude gap). One of the problems is to reconcile this discrepancy. Following standard practice, we call this the old cosmological constant problem (old-CCP).
  • Phase transitions. The potential energy for a constant field configuration contributes as vacuum energy, but, for instance, the Higgs field’s effective potential depends on the background temperature. As the Universe cools down, the potential changes its shape, shifting its global minimum. Therefore, its contribution to the vacuum energy also changes. This problem is a classical instability due to the change in the contribution coming from the potential energy of the Higgs field before and after its phase transition. In fact, this is just one example of what happens generally for any phase transition, including the one due to quantum chromodynamics (QCD) and possibly others from unification theories. This is called the classical cosmological constant problem [24] (class-CCP).
  • Dark energy. The Universe is undergoing a period of accelerated expansion [43,44] that can be explained by a nonzero vacuum energy in the form of a positive CC. Thus, we need to explain where this positive-cosmological-constant-like term is coming from at the cosmological level when using Einstein equations. We will call this the dark energy problem (DEP).2
  • UV sensitivity. Finally, another problem comes from the fact that the vacuum energy computed in QFT is UV-sensitive, despite being possibly the most infrared (IR) quantity one could conceive (as it is a constant throughout spacetime). In particular, there are two ways in which the UV sensitivity of ρ vac manifests:
    It is directly connected with the Higgs’ mass UV sensitivity. As we saw above, the vacuum energy scales with the mass of particles, and the Higgs’ mass squared is itself highly UV-sensitive (quadratic in cutoff). This is the hierarchy problem (see [47,48] for a recent review), and it manifests as an even worse sensitivity in the computation of the vacuum energy (since it is quartic in the cutoff).
    As we change the QFT cutoff by increasing it to higher-energy scales, we might be able to excite new fields with higher masses, again disturbing the fixing of the CC performed at lower scales. Thus, even if all the masses were not fine-tuned in the SM, would still be sensitive to new fields showing up at new higher-energy EFTs.
    Thus, in short, once we change the energy scale in which we are computing the vacuum energy, the radiative corrections from higher-order loop corrections in QFT will shift the value of ρ vac . We will refer to this issue as the new cosmological constant problem (new-CCP).
Note that these problems do not have the same nature. The weight of vacuum can be easily fixed with the presence of a bare CC in Einstein equations, which is then fixed by measurements (as any other fundamental constant of nature is fixed by applying the renormalization program). This cancellation between the bare CC and the vacuum energy could be seen as a coincidence, but poses no problem to the theory. Thus, we take the stance that this has never been an issue, despite it still being widely referred to as such, both in the physics literature and also in popularized accounts about the CCP.
However, fixing the value of the bare CC remains a potential issue as phase transitions unfold and shift its value, and the final value of the effective vacuum energy can be adjusted only once by a bare CC. On the other hand, if an approach can guarantee that phase transitions do not spoil this adjustment, then, naturally, the final value of the effective CC could correspond to the one providing the correct phenomenology associated with DE, sorting out the DE problem. Thus, the old, the classical, and the dark energy problems can all be understood at the classical level of Einstein equations and are directly related to phenomenology.
Note, however, that the dark energy problem may or may not be solved in these approaches. That depends on whether these models allow to have some residual-like CC term that would induce an accelerated expansion. Nonetheless, in principle, this problem can be tackled independently by combining these approaches with quintessential models that would drive the dark energy dynamics (see [49,50,51] for a review).
Finally, the UV sensitivity is a deeper issue. Within GR, even if one manages to classically solve the other problems, this sensitivity would force the solution to be fine-tuned (see also [52]). As we consider models that modify gravity for which there is a self-tuning mechanism that effectively prevents the vacuum energy from gravitating, these models only directly solve the old- and the classical-CCPs. By this, they tacitly accept that the UV sensitivity is not necessarily an issue, in the sense that they do not attempt to fix it, leaving QFT calculations untouched. Instead, gravity is modified such that the breakdown of the UV–IR decoupling does not lead to dramatic observable effects, thereby disposing of the need for fine-tuning.3

3. How to Modify Gravity

There are two features shared by almost all modified gravity approaches to the CCPs. First, all models discussed in this review rely on the idea of making gravity less sensitive to the presence of vacuum energy. We will refer to this general idea as self-tuning and introduce it in Section 3.1. Second, self-tuning is typically implemented in terms of additional degrees of freedom, which, in general, mediate a gravitational force. This requires the presence of a mechanism that suppresses this additional force’s phenomenology at observationally accessible scales. We provide a review of common screening mechanisms in Section 3.2.

3.1. Self-Tuning

The key observation is that, due to the equivalence principle, gravity is sourced by vacuum energy, raising the need to cancel it against a highly fine-tuned cosmological constant (see Section 2). Therefore, if we modify our gravitational theory in such a way that vacuum energy either drops out of Einstein equations or is strongly suppressed, this avoids the need to fine-tune the CC to make its value compatible with observations. Instead, we assume that the CC takes on a generic value set by the cutoff of the theory. In other words, we accept that the cosmological constant is highly sensitive to unknown UV physics but have our gravitational theory turn a blind eye4.
This rather general idea comes in many flavors and under different names. Common notions are “self-tuning”, “shielding of the CC”, or “self-adjustment mechanism”. A model-independent implementation is provided by the “degravitation proposal” which promotes the Newtonian constant to a high-pass filter (see Section 4.3.1). For simplicity, we will refer to the general idea as self-tuning but resort to the notion commonly used in the literature when reviewing individual models.
To explain the concept of self-tuning and gain an idea of common obstacles, we will first consider a somewhat naive and insufficient implementation of self-tuning in Section 3.1.1. In a more general context, the reason for its failure is formalized by Weinberg’s no-go theorem reviewed in Section 3.1.2. A complementary argument based on a spectral decomposition will be provided in Section 3.1.3. This will guide us towards a small handful of loopholes that can tell us how to modify gravity in order to implement self-tuning. They are later used to motivate the different models reviewed in Section 4.

3.1.1. A Failed Start

As a warm-up, we consider a simple scalar field toy model that illustrates the idea of self-tuning and highlights the essential difficulties we face when implementing it. We will then generalize the discussion, using model-independent arguments, when reviewing two complementary no-gos in the literature in Section 3.1.2 and Section 3.1.3.
Consider an action containing GR, a (dimensionless) scalar field ϕ , and a matter sector S m [ e β ϕ g μ ν , Ψ ] with field content collectively denoted by Ψ that exhibits a conformal coupling to ϕ (controlled by β )5
S = M Pl 2 2 d 4 x g R + d 4 x g M 2 2 μ ϕ μ ϕ V ( ϕ ) + S m [ e β ϕ g μ ν , Ψ ] ,
where M Pl is the Planck mass and M sets the coupling strength of the scalar ϕ . In other words, we can think of ϕ as mediating an additional gravitational force. Note that we define the theory in the “Einstein frame” where the action for the metric is the Einstein–Hilbert action. The gravitational equations of motion are
M 2 ϕ = V ( ϕ ) β 2 T ,
M Pl 2 G μ ν = T μ ν ( ϕ ) + T μ ν ,
where we defined the EMT of the scalar field as
T μ ν ( ϕ ) = M 2 μ ϕ ν ϕ + g μ ν 1 2 M 2 g α β α ϕ β ϕ V ( ϕ ) ,
and that of the matter sector as
T μ ν = 2 g δ S m δ g μ ν .
Taking the trace of Equation (2b) returns
M Pl 2 R + T ( ϕ ) + T = 0 .
As usual, we expect the vacuum energy ρ vac to contribute to the matter energy density. To make this more explicit, we can decompose
T μ ν = τ μ ν ρ vac e 2 β ϕ g μ ν ,
where τ μ ν corresponds to a localized source with asymptotic fall-off. In particular, τ μ ν = 0 in the vacuum.
The old and new CCP in this language boil down to the statement that ρ vac , due to Equation (5), needs to be finely tuned to a small value (irrespective of the cutoff of the theory or the loop-order) to avoid making R unacceptably large. The classical CCP, on the other hand, is the observation that this tuning would be spoiled in a phase transition (unless the phase transition itself is tuned not to make a contribution to ρ vac ). Of course, this reasoning only holds if there is no cancellation between T ( ϕ ) and T at play. Now, self-tuning is based on the the idea that such a cancellation does take place. As we will explain, it can be enforced through the ansatz (at this stage, other definitions differing by gradient terms are possible too):
T + T ( ϕ ) = V ( ϕ ) + β 2 T .
For the mechanism to work, a sufficient assumption is that the scalar field vacuum is translationally invariant,
μ ϕ vac = 0 .
The reason is that when the system settles to its vacuum configuration, the LHS of Equation (2a) vanishes, enforcing ( V ( ϕ ) β / 2 T ) 0 . The self-tuning Equation (7) along with Equation (6) then imposes the condition 4 ρ vac e 2 β ϕ T ( ϕ ) | ϕ = ϕ vac = 0 , which, in turn, cancels vacuum energy from Equation (5) (and similarly from Equation (2b)). In particular, if the vacuum energy changes during a phase transition or by integrating in or out new particles, leading to ρ vac ρ ^ vac = ρ vac + Δ ρ vac , the field ϕ is pushed away from its stationary point and settles to a new configuration ϕ ^ vac , which now imposes 4 ρ ^ vac e β ϕ ^ vac T ( ϕ ) | ϕ = ϕ ^ vac = 0 . As a consequence, in the absence of localized matter sources, i.e., for τ μ ν = 0 , the Minkowski metric remains a solution despite the presence of huge vacuum energy. This is indeed the definition of self-tuning.
In short, we have—somewhat naively—imposed self-tuning by supplementing the system (2) with the ad hoc Equation (7). Can this be performed consistently without overconstraining the system? To find out, we try to derive a potential V ( ϕ ) such that all equations can be satisfied simultaneously. Translational invariance in Equation (8) implies that close to the vacuum configuration we can neglect gradients and time derivatives; specifically, μ ϕ ν ϕ V ( ϕ ) , and thus T ( ϕ ) 4 V ( ϕ ) . Substituting this back into the self-tuning Equation (7) and using Equation (6) yields
4 V ( ϕ ) V ( ϕ ) 4 1 β 2 ρ vac e 2 β ϕ ,
which is a differential equation for V ( ϕ ) . Its general solution takes the form
V ( ϕ ) ρ vac e 2 β ϕ + V 0 e 4 ϕ for β 2 V 0 e 4 ϕ for β = 2 ,
where V 0 is an integration constant. There are two interesting cases:
  • For β 2 , the potential directly depends on ρ vac in such a way that it exactly cancels the vacuum energy in both equations in Equation (2) (or Equation (5) equivalently). This is fine-tuning again and it would not be robust under a change of the loop-order or a shift in the EFT cutoff, and it also means the cancellation is no longer dynamical. Therefore, this model does not address any of the CCPs discussed in Section 2.
  • For β = 2 , the potential does not depend on ρ vac . Unfortunately, we have still arrived at an impasse. Our self-tuning mechanism assumed that ϕ settles down to a constant vacuum configuration ϕ vac ; however, the potential in Equation (10) has no minimum at finite field values. Instead, it is a runaway potential, which approaches its minimum only asymptotically as ϕ . Setting V 0 = 0 to avoid this conclusion would constitute a fine-tuning again. The reason is that V 0 receives huge (UV-sensitive) radiative corrections because only the combination V 0 + ρ vac contributes to Equation (2).
We could hope that the running itself is not a problem. After all, it seems to exponentially suppress the vacuum energy contribution to T ρ vac e 4 ϕ g μ ν as ϕ . The problem with that is that it is not only the vacuum energy that is being suppressed.For example, the conformal coupling with ϕ implies that the Higgs mass term is replaced through m H 2 H 2 m H 2 H 2 e 4 ϕ . Now, as ϕ runs down its potential, the Higgs mass goes to zero, which—unless the running can be made very slow—is not compatible with observations. Thus, neither setting V 0 = 0 nor accepting the running seems to be a valid option. In the next section, we will find that this failure can be explained by a very general symmetry argument.

3.1.2. Weinberg’s Argument

Weinberg in his seminal paper on the cosmological constant problem investigates self-tuning and finds that under a few key assumptions, self-tuning is not possible without fine-tuning [23]. Here, we review his argument (see, also, [41,52]).
We consider the Lagrangian L ( g μ ν , ϕ i ) of a metric field g μ ν and a set of self-tuning scalars ϕ i . The presence of additional matter fields Ψ that give rise to vacuum energy is understood. This approach covers the toy model from the previous section but is far more general. In particular, it does not assume the gravitational sector to be described by GR or canonical kinetic terms for ϕ i . What we do assume is that the vacuum of the theory is translationally invariant,
( ϕ vac ) i = const and ( g vac ) μ ν = const .
While this is in agreement with our expectation from typical field theories,6 we will discuss an example of self-tuning in Section 4.3 and Section 4.4, where this assumption is violated. In any event, general coordinate invariance tells us that L has to transform as a density (to make d 4 x L invariant). This is guaranteed if L takes on the form [41]7
L vac L ( ( g vac ) μ ν , ( ϕ vac ) i ) = g vac ρ vac ( ( ϕ vac ) i , ) ,
where the ellipsis indicates the dependence on further SM vacuum expectation values. Note that ρ vac , as defined here, also includes the contribution from a bare cosmological constant and the self-tuning field’s potential.
In this general context, self-tuning can be realized through the equation
2 g μ ν L g μ ν = i f i ( ϕ ) L ϕ i ,
where f i are a set of general functions that depend on the self-tuning fields ϕ i . As a quick sanity check, with these definitions, we recover Equation (9) from Equation (1) for f 1 = 4 . The idea is the same as before. The ϕ i are subject to their equations of motion,
L ϕ i μ L ( μ ϕ i ) = 0 .
Therefore, as the ϕ i approach their constant vacuum value, L / ϕ i 0 , which, due to Equation (13), implies g μ ν L vac / g μ ν = 0 . From Equation (12), we see that this requires the vanishing of ρ vac . It indeed looks as if the ϕ i have “self-tuned” to make the vacuum energy vanish.
To learn more about the viability of the mechanism, we further constrain the functional form of L vac . The key observation is that Equation (13) implies the scaling symmetry (with ϵ 1 )
δ ϵ ( g vac ) μ ν = 2 ϵ ( g vac ) μ ν and δ ϵ ( ϕ vac ) i = ϵ f i ( ϕ vac ) ,
under which the vacuum theory is invariant. Indeed, it is easy to check that δ ϵ L vac = 0 if we use Equation (13) along with Equation (11). Applying an appropriate field transformation ϕ i ϕ ˜ i , we can write Equation (15) as (suppressing the subscript “vac”)
δ ϵ g μ ν = 2 ϵ g μ ν , δ ϵ ϕ ˜ 0 = ϵ , and δ ϵ ϕ ˜ i 0 = 0 .
As δ ϵ e 2 ϕ ˜ 0 g μ ν = 0 , we conclude that L vac = L ( e 2 ϕ ˜ 0 g μ ν ) . This can be made compatible with Equation (12), if we factorize ρ vac ( ϕ i ) e 4 ϕ ˜ 0 ρ ˜ vac ( ϕ ˜ i 0 ) . We therefore find,
L vac = det ( g μ ν ) e 4 ϕ ˜ 0 ρ ˜ vac ( ϕ ˜ i 0 ) ,
which, indeed, takes the right form if we use det ( g μ ν ) e 4 ϕ ˜ 0 = det ( e 2 ϕ ˜ 0 g μ ν ) . Moreover, this recovers our previous result in Section 3.1.1 (for β = 2 ) if we identify ρ ˜ vac = V 0 + ρ vac . As before, setting ρ ˜ vac = 0 (or setting it to its observed value) is fine-tuning. For ρ ˜ vac 0 , on the other hand, we again find that Equation (17) describes a runaway potential, which contradicts the initial assumption in Equation (11) and therefore concludes the no-go argument.
We summarize the main assumptions used for the argument:
(i)
The scalar field vacuum is constant and thus invariant under translations (see Equation (11));
(ii)
The vacuum geometry is constant (this assumption is stronger than maximal symmetry);
(iii)
The theory can be formulated in terms of a Lagrange density L . The covariant integration measure is provided by d 4 x g as in (12).
We will discuss loopholes in Section 3.1.4.

3.1.3. Beyond Weinberg’s Argument

Weinberg’s no-go theorem can be complemented by an independent argument that does not assume the translational invariance of the scalar field vacuum, can be generalized to maximally symmetric vacua with positive curvature, and allows contact with phenomenology [57]. It builds on a spectral decomposition of the gravitational exchange amplitude,
A = s = 0 , 2 0 d μ ρ s ( μ ) x , y T ¯ ( x ) α β G α β γ δ ( s ) ( x , y ; μ ) T ( y ) γ δ .
Here, T ¯ μ ν is assumed to be a source that falls off at infinity. Physically, it plays the role of a local mass distribution that is accessible to observations such as a finite set of point masses. Its counterpart T μ ν represents a generic source, which, in particular, is not assumed to fall off at infinity and will later be identified with vacuum energy. The gravitational exchange between both sources is described in terms of a propagator G α β γ δ ( s ) ( μ ) . The spectral density ρ s ( μ ) then gives a particular weight to contributions with different mass-squared μ = m 2 and bosonic spin s { 0 , 2 } 8. It can be decomposed into massive and massless contributions as ρ s ( μ ) = ρ s , m ( μ ) + ρ ¯ s δ ( μ ) . For example, GR only contains a massless spin-2 field, which is realized through ρ ¯ 2 = 1 / M Pl 2 , ρ ¯ 0 = 0 , and ρ s , m ( μ ) = 0 , and gives rise to the familiar exchange amplitude between two localized sources after substituting the propagator for a massless graviton (for a detailed derivation of G α β γ δ ( 2 ) ( μ ) for μ = 0 and μ 0 , see, for example, Section 2.2 in [37]):
A GR = 1 M Pl 2 T ¯ α β 1 T α β 1 2 T ¯ 1 T ,
where we use 1 f as a shorthand for the convolution of f with the Green function of the d’Alembertian. On the other hand, a single massive spin-2 particle with mass m and Planckian coupling is introduced through ρ 2 , m = δ ( m 2 μ ) / M Pl 2 . It contributes to Equation (18) as
A FP = 1 M Pl 2 T ¯ α β 1 + m 2 T α β 1 3 T ¯ 1 + m 2 T ,
where the factor 1 / 3 corresponds to the notorious Fierz–Pauli tensor structure [58] (see also Section 4.3.2). Finally, a scalar mediator will lead to a contribution with trivial tensor structure T ¯ 1 ¯ + m 2 T .
This general formulation covers a wide range of gravitational models that rely on introducing additional gravitational degrees of freedom. It builds on a minimal set of assumptions:
(i)
We assume classical and quantum stability, which requires μ 0 (absence of tachyons) and ρ s ( μ ) > 0 (unitarity);
(ii)
The vacuum geometry is either Minkowski or (quasi) de Sitter;
(iii)
A weakly coupling assumption underlies the explicit derivation of G α β γ δ ( s ) ( μ ) and limits the approach to linear source couplings;
(iv)
The spin-0 and spin-2 propagators take on their canonical form.
Now, self-tuning describes the absence of a gravitational response in the presence of vacuum energy, explicitly
A | T μ ν = ρ vac g μ ν = ! 0 .
This is a constraint that applies in the deep infrared as vacuum energy is a spacetime constant and therefore only couples to the zero (four-)momentum mode. On a Minkowski background, this follows from the fact that the momentum space representation of vacuum energy is T μ ν η μ ν δ ( 4 ) ( p ) . This infrared condition has to be supplemented with a phenomenological constraint that enforces the approximate recovery of Equation (19) for localized sources and hence applies in the ultraviolet where 9. To be specific, we demand10
A A GR ϵ = 1 M Pl 2 T ¯ α β 1 T α β 1 2 ( 1 ϵ ) T ¯ 1 T ,
where the parameter ϵ < 10 5 parametrizes a small deviation from the GR tensor structure in Equation (19) allowed by solar system observations [59,60]. This approach is complementary to Weinberg’s argument because it does not require the scalar field vacuum to be invariant under translations. In addition, it does not rely on a Lagrangian formulation, nor is it restricted to the classical limit.
Both conditions in Equation (21) can be readily evaluated on a Minkowski background. The self-tuning constraint yields11  ρ ¯ 2 = 6 ρ ¯ 0 . We note that it does not depend on the massive particle spectrum, which, in our language, shows that Fierz–Pauli massive gravity has a built-in self-tuning mechanism (see Section 4.3). In any event, the phenomenological constraint yields
ρ ¯ 2 p 2 + 0 + d μ ρ 2 , m ( μ ) p 2 + μ p 1 p 2 M Pl ,
1 3 0 d μ ρ 2 , m ( μ ) p 2 + μ + 2 ρ ¯ 0 p 2 + 2 0 + d μ ρ 0 ( μ ) p 2 + μ p ϵ p 2 M Pl .
Using the positivity of each term in (22b), we find ρ ¯ 0 ϵ / M Pl 2 and p 2 μ ρ 2 , m p 2 + μ ϵ / M Pl 2 . With this, we conclude from Equation (22a) that ρ ¯ 2 1 / M Pl 2 . The conditions on ρ ¯ 0 and ρ ¯ 2 are, however, in contradiction with the self-tuning relation ρ ¯ 2 = 6 ρ ¯ 0 . In other words, self-tuning under the above assumptions is incompatible with phenomenology. This concludes the argument. As shown in [57], the same conclusion can be reached for de Sitter vacua.

3.1.4. Loopholes

Is this generalized approach, alongside Weinberg’s argument, heralding the end of self-tuning models? It certainly rules out many direct attempts at implementing the mechanism, but, as is true with every no-go theorem, these, too, are limited by their assumptions. They allow us to formulate loopholes that themselves can serve as the seeds of future research. Here, we single out three possibilities.
  • Models that rely on nonlinear couplings to restore GR (violating assumption (iii) of the amplitude argument). In particular, this includes models that build on the Vainshtein mechanism presented in Section 3.2.2. In this review, we discuss two examples: massive gravity in Section 4.3 and scalar–tensor theories such as Fab-4 in Section 4.4. Whether the no-go arguments can be generalized to this case remains to be seen. It should also be noted that Vainshtein screening comes with its own set of problems pertaining to its UV sensitivity in the nonlinear regime [61,62,63]. Both models also circumnavigate Weinberg’s complementary no-go by relying on a scalar field vacuum that breaks translational invariance (violating assumption (i) in Section 3.1.2)a. In the case of massive gravity, this happens in the Stuckelberg sector of the theory.
  • Another key assumption in both theorems is that vacuum energy, and, thus, also the corresponding geometry, is maximally symmetric. In particular, the vacuum does not break local Lorentz invariance. As we will argue in Section 4.5, in braneworld models, the vacuum energy that arises from SM matter loops (spontaneously) breaks Lorentz invariance in the directions orthogonal to the brane (violating assumption (ii) of both arguments).
  • For the generalized no-go, we assumed that propagators take on their canonical form as could be derived in canonical local and ghost-free field theories of a given spin. Accordingly, this loophole relies on introducing nonstandard propagators (violating assumption (iv) of the amplitude argument). This is, for example, how sequestering in Section 4.2 avoids the no-go. It also avoids Weinberg’s argument by introducing either global variables or employing a four-form contribution to the volume measure (violating assumption (iii) in Section 3.1.2).
To provide an explicit example of how the last loophole can be exploited, we follow [57] and discuss the decapitation idea first introduced in a string theory context [64]. To that end, we consider an exchange amplitude of the form
A = A GR + 1 6 M Pl 2 T ¯ 1 T A ghost .
The second term is the contribution of a canonical scalar field. Its numerical coefficient is chosen such that it cancels with A GR for a vacuum energy source. As we have seen before, this is not possible without violating the observational constraint in Equation (21b) when considering localized sources. This is where the third term A ghost comes to the rescue. For localized sources, it exactly cancels the second term while it vanishes for vacuum energy:
A ghost = 1 6 M Pl 2 T ¯ 1 T for T μ ν ( x ) local 0 for T μ ν = ρ vac g μ ν .
As a result, we recover exact GR for every source that is localized, but realize self-tuning for vacuum energy sources. We call it a ghost term because it leads to a negative contribution in Equation (23). However, it does not lead to an instability, as could be expected for generic ghost fields, because it is either vanishing or exactly compensated by the healthy scalar mode (this is somewhat analogous to Faddeev–Popov ghosts in gauge quantum field theory, which also do not lead to instabilities due to an cancellation with healthy but unphysical modes). Writing A ghost = 1 6 M Pl 2 x , y T ¯ ( x ) G decap ( x , y ) T ( y ) , the nontrivial question is what type of propagator G decap will lead to this exotic behavior. As it turns out, a nonlocal construction can achieve it, explicitly [64]
G decap ( x , y ) = G 0 ( x , y ) 1 V z G 0 ( z , y ) 1 V z G 0 ( x , z ) + 1 V 2 z 1 , z 2 G 0 ( z 1 , z 2 ) ,
where G 0 denotes the standard Green function of the d’Alembertian and V is the (appropriately regularized) spacetime volume. It is referred to as a decapitated propagator because it does not couple to vacuum loops. To show that it fulfills Equation (24), the property x G decap ( x , y ) = y G decap ( x , y ) = 0 can be used. It was demonstrated in [57] that the sequestering model (discussed in Section 4.2.1) provides a possible field theory realization of the decapitation idea.
Of course, the list of loopholes is not complete, and other models are conceivable that, for example, incorporate a mild tachyonic instability, consider an anti-de Sitter vacuum, or find different ways of deviating from the canonical field theory framework.

3.2. Screening Mechanisms

As we have seen, by challenging the assumption of GR we can hope to address aspects of the CCPs. A crucial constraint coming from observations on such avenues is that the proposed modification to gravity be undetectable in our local environment. Specifically, the numerous precision tests of gravity in the Solar System only allow extremely mild deviations from GR at these scales (see, for example, [65,66,67,68,69]). This poses a problem for scalar tensor theories which generically predict a supplementary force on top of the standard GR prediction, usually referred to as the “fifth force”. This is of particular relevance for scalar tensor theories such as the Fab-4 (see Section 4.4) or the Dvali–Gabadadze–Porrati braneworld model (DGP) [70,71] (see Section 4.5), which can be written as a scalar tensor theory in the limit that the five-dimensional degrees of freedom decouple (see [72] for an example).
The scalar tensor theories we will consider here and in this review will fall under the broad class of Horndeski theories [73], which can be written as
S H = d 4 x g [ G 2 ( ϕ , X ) G 3 ( ϕ , X ) ϕ + G 4 ( ϕ , X ) R + G 4 , X ( ϕ , X ) [ ( ϕ ) 2 ( μ ν ϕ ) 2 ] + G 5 ( ϕ , X ) G μ ν μ ν ϕ 1 6 G 5 , X ( ϕ , X ) [ ( ϕ ) 3 3 ϕ ( μ ν ϕ ) 2 + 2 ( μ ν ϕ ) 3 ] ] ,
where each G i ( ϕ , X ) , i = 2 , 3 , 4 , 5 is a free function of the scalar field ϕ and its canonical kinetic term X = ( ϕ ) 2 / 2 , and G i , X ( ϕ , X ) = G i / X , with □ being the D’Alembertian operator. R and G μ ν are the Ricci scalar and Einstein tensor components, respectively. This class ensures the resulting equations of motion are at most second-order and nondegenerate which makes them Ostrogradski-ghost-free (see [74] for a discussion) and theoretically viable. We note that one can have healthy theories with higher-than-second-order derivatives by considering degenerate theories or certain field symmetries [75,76,77,78,79].
As an illustrative example of the problem, we can consider the archetypal Jordan–Brans–Dicke theory [80,81] which is contained in the Horndeski class. We can write this theory in the Einstein frame, with a minimal coupling of the scalar field to the gravitational sector, i.e., it does not multiply the Ricci scalar, but a modified matter sector as
L JBD = M Pl 2 R 2 1 2 g μ ν μ ϕ ν ϕ V ( ϕ ) + L m ( A 2 ( ϕ ) g μ ν , Ψ M ) ,
where A ( ϕ ) is the conformal transformation between minimally and nonminimally coupled frames, ϕ is the scalar degree of freedom, and Ψ M are the collective matter fields.
If we consider now a perfect fluid
T μ ν = ρ m + p c 2 u μ u ν + p g μ ν ,
with rest frame mass density ρ m and we assume a small pressure p ρ c 2 , in a perturbed Minkowski spacetime,
d s 2 = ( 1 + 2 Φ ) d t 2 + ( 1 2 Ψ ) ( d x 2 + d y 2 + d z 2 ) ,
we can write the 00-th Einstein field equation in the Newtonian limit as
2 Φ = 1 2 M Pl 2 ρ m 1 2 2 ϕ ,
where 2 is the Euclidean Laplacian operator. Note that we ignored pressure contributions and only considered first-order contributions in the metric perturbation Φ . We can also write down the geodesic equation in the nonrelativistic limit ( d x i / d t 1 ), which gives the acceleration a test particle experience:
a = Φ β ϕ a N + a 5 ,
where we defined β ( ϕ ) [ ln A ( ϕ ) ] , a prime denoting a scalar field derivative, and defined the Newtonian acceleration and an additional acceleration associated with a “5th force” sourced by ϕ . One basic criterion for any viable gravitational theory is that the scalar field’s contribution to the acceleration is negligible, or screened, in environments such as the Solar System12. We will discuss two generic ways of ensuring this holds true: screening through a well-chosen potential for the field and screening that employs derivatives of the field. We refer the interested reader to [36,82,83] for more detailed reviews on screening mechanisms.

3.2.1. Potential Screening

This method of screening relies on a well-chosen potential V ( ϕ ) for the scalar field. Assuming negligible time variation of ϕ and a pressureless perfect fluid, the equation of motion for the scalar field in Equation (27) is given by
2 ϕ = 1 M Pl 2 β ( ϕ ) ρ m V ( ϕ ) V eff ( ϕ ) ,
where we defined an effective potential Veff. Screening is realized when the effective potential finds a minimum, i.e., 2 ϕ = 0 . With appropriate choices of V ( ϕ ) and β ( ϕ ) (or, equivalently, A ( ϕ ) ), this can be made to happen in regions of high density such as the Solar System. With this choice, we note that moving beyond the Solar System, to, say, larger scales and lower densities, objects may begin to feel the additional force contribution, which will have an impact on structure formation. Typically such theories aim to reproduce Λ CDM at high redshift, where measurements of the cosmic microwave background (CMB) leave little room for deviation [3]. This leaves the possibility of low redshift, large-scale impacts of the scalar field (see, for example, [84,85,86,87]).
Returning to Equation (32), it would appear that we have two free functions with an infinite number of new degrees of freedom with which to match the data. While in principle this is true, these would be ruled out from a Bayesian standpoint. A simple and popular example of such a theory that promised to both exhibit screening and was able to provide acceleration without the CC 13 is the Hu–Sawicki model of f ( R ) [90]:
f ( R ) = m 2 c 1 ( R / m 2 ) n 1 + c 2 ( R / m 2 ) n .
Here, c 1 or c 2 is fixed to necessarily match a Λ CDM expansion history. m 2 is also typically fixed, leaving two free parameters, n and c 1 or c 2 , the latter usually being recast in terms of the value of the scalar field today f R 0 d f d R | z = 0 . This functional choice also fixes both V ( ϕ ) and β ( ϕ ) in the equivalent Einstein frame cast scalar tensor theory.
To see how screening works, let us consider a spherically symmetric matter distribution with radius R and total density ρ m , embedded in a background of density ρ ¯ . By the symmetry of the distribution we should be able to find some r S < R such that the effective potential V eff ( ϕ ) 0 , which defines the screened regime. Further, for r > r S , we can find a region such that ϕ = ϕ ¯ + δ ϕ ( r ) , where δ ϕ is a small perturbation. The equation of motion for this perturbation is then
2 δ ϕ = 2 ϕ 2 ϕ ¯ 1 M Pl 2 ( ρ m ρ ¯ ) β ¯ ( V ( ϕ ) V ( ϕ ¯ ) ) 1 M Pl 2 β ¯ δ ρ ¯ δ ϕ V ( ϕ ¯ ) ,
where we defined the density contrast δ ρ ¯ ρ m ρ ¯ , where a bar denotes a background quantity. In the first line, we only considered the background term in β , to stay at consistent order in the perturbations β ( ϕ ) β ( ϕ ¯ ) β ¯ , which can be restated as
d β d ϕ δ ϕ β ¯ for r > r S .
The second approximation assumes that V ( ϕ ) is slowly varying. For cosmologically relevant fields, we typically choose their background mass m 0 2 V ( ϕ ¯ ) H 0 , where H 0 is the Hubble constant, so that we only see modified gravity effects on very large scales. Near the screening regime we can then ignore the V term in Equation (34). Integrating Equation (34) for r > r S gives the acceleration associated with the fifth force:
a 5 β ¯ ϕ = β ¯ δ ϕ = 2 β ¯ 2 G M ( r ) r 2 1 M ( r S ) M ( r ) for r > r S ,
where the density contrast sourced mass within radius r is given by
M ( r ) = 4 π 0 r x 2 δ ρ ¯ ( x ) d x .
From Equation (36), we see that if r S = R , then the object is fully screened ( F 5 = 0 ), but, conversely, if r S = 0 , the object is fully unscreened. Typically, we would have 0 < r S < R , which would offer a partially screened object with the fifth force being sourced by a shell of mass between r S < r < R . The screening radius r S can be derived from the field profile once the potential and conformal transformation have been defined. These can be chosen such that r S is only slightly smaller than R so that the fifth force is only sourced by a thin shell of mass and is consequently limited in size and scales of effect. Of course, one could also find fully screened configurations, but these do not offer interesting signatures which can be looked for experimentally.
Chameleon models are a class of such models exhibiting this type of screening, and they have been extensively studied in the literature [91,92,93] (see [88,94] for reviews of fairly recent constraints). The choice for potential and conformal transformation in these models usually take the form
V ( ϕ ) = m 0 n + 4 ϕ n , A ( ϕ ) = e B ϕ ,
with B being a constant. Note that we can also obtain chameleon screening from minimally coupled matter in f ( R ) gravity (see, for example, Equation (33)) after the field redefinition necessary to find the Einstein frame. Given the freedom these models allow, data constraints are fairly restricted to specific models or parametrizations (see [94] for constraints on the model in Equation (38)).

3.2.2. Derivative Screening

To see how this type of screening works, we can proceed analogously to the self-interaction screening case. First, let us return to the modified geodesic equation, Equation (31). One way of shutting off the scalar field term on the RHS is to include additional derivative interactions in the action. As an example we can take the five-dimensional DGP braneworld model [70] in the decoupling limit, which, it has been argued, reduces to the cubic Galileon theory [72,95]
L G 3 = M Pl 2 R 2 1 2 g μ ν μ ϕ ν ϕ 1 λ 0 2 g μ ν μ ϕ ν ϕ ϕ + α V ϕ 1 M Pl 2 T ,
where λ 0 2 gives the energy scale of the theory, and we made the coupling to the trace of the matter sector EMT explicit through the dimensionless constant α V . We note that this can be achieved similarly by setting A 2 ( ϕ ) = e 2 ϕ α V in Equation (27). This action produces the following equation of motion:
ϕ + 2 λ 0 2 ( ϕ ) 2 ν μ ϕ ν μ ϕ = 1 M Pl 2 α V ρ m ,
where μ is the covariant derivative. In DGP gravity we have [96] λ 0 2 = 3 2 r c 2 β DGP and α V = 1 3 β DGP , where β DGP is associated with the cosmological background in the model and r c is the cross-over scale below which gravity becomes four-dimensional.
As in the chameleon example, if we consider small enough scales where time derivatives are small compared to spatial ones and consider a spherically symmetric mass distribution of total mass M, Equation (40) can be written as (see [97] for example)
1 r 2 d d r r 2 d ϕ d r + 4 r λ 0 2 d ϕ d r 2 = 1 M Pl 2 α V ρ m .
Integrating this equation yields
r 2 d ϕ d r + 4 r λ 0 2 d ϕ d r 2 = α V M 4 π M Pl 2 2 α V r 2 F N ,
where F N is the Newtonian gravitational force per unit mass. Taking the small r limit of Equation (42), when the second term on the LHS dominates, yields
F 5 F N = 2 α V r r V 3 / 2 ,
where F 5 = d ϕ / d r is the fifth force per unit mass, and we defined the so-called Vainshtein radius r V 8 G N M α V λ 0 2 1 / 3 , named after Arkady Vainshtein who first proposed this mechanism, to protect Fierz–Pauli theory (see Section 4.3.2) from Solar System constraints [98]. We see F 5 r 3 / 2 , and so for large Vainshtein radii, we obtain a large suppression of the fifth force when r r V . This mechanism is highly efficient, for example, in the case that M = M and α V of O ( 1 ) , we find that r V is seven orders of magnitude larger than the Solar System. On the other hand, in the large-scale limit where the first term of Equation (42) dominates, we have a modification of F 5 / F N 2 α V , which can offer interesting phenomenology to venture out and test experimentally (again, see [94] for recent constraints).
As an illustration of both types of screening mechanism, we can consider the effective modification to Newton’s gravitational constant Geff, defined using Equation (30) as
2 Φ ρ m 2 M Pl 2 G eff G N .
This quantity has been derived under various assumptions, such as a spherically symmetric density distribution, in both the chameleon screened model of Equation (33) [99] and in the Vainshtein screened DGP model [100]. In Figure 1 we plot G eff / G N against the normalized halo radius. We can clearly see that at small scales the screening mechanism kicks in and G eff G N as required. Interestingly, the DGP model (red, solid) also modifies the large-scale growth of structure, while f ( R ) (blue, dashed) experiences Yukawa suppression at large scales, giving distinct and interesting phenomenology. Further, we see the fifth force acting at some intermediate range of scales which affects the growth of structure. Detection of fifth force impact on structure and signatures of modified gravity is a prime science goal of upcoming and ongoing large-scale structure experiments [101,102,103].
We conclude this section with some general remarks on the theoretical viability of screened models. Note that screening effects kick in at scales where the low-energy theory should begin to fail. In the case of self-interaction screening, the field potential should begin to receive nonnegligible high-energy corrections in the screening regime, which poses a problem for these models. Similar considerations were made for nonlinear interaction models, such as Vainshtein screened ones [61,63], although for these models it has been argued that they may be protected against high-energy corrections as they typically satisfy nonrenormalization theorems [62,104]. We also direct the reader to Section 6 of [105] and the references therein for more on this discussion.

4. Modified Gravity Approaches

Now that we have clarified the issues surrounding the cosmological constant and reviewed the no-go theorems constraining possible modifications of gravity, we are apt to discuss some of the approaches that exploit the loopholes of these theorems, as discussed in Section 3.1.4. Including our tacit discussion about semi-classical gravity which combines GR and QFT, we will discuss nine different proposals that attempt to instantiate the self-tuning mechanism. As we already mentioned in the introduction, this is in no way exhaustive and only constitutes a representative sample of the rich literature available on modifications of gravity14.
We start in Section 4.1 by briefly recalling the challenges introduced in Section 2 and defining two additional requirements related to astrophysical and cosmological data that a successful proposal has to fulfill. Then, in Section 4.2 we discuss a general class of models that makes use of top forms to constrain the gravitational dynamics; in Section 4.3 we discuss the degravitation mechanism and its realization in the context of massive gravity, and we finish by reviewing how self-tuning can be implemented in the context of Horndeski theories in Section 4.4.

4.1. What Do We Want?

As we saw in Section 2, there are at least three problems that a model can attempt to solve: the classical and new CCP, and the DEP. However, even if a model successfully solves one or more of these problems by modifying gravity, we need to guarantee that the model does not spoil phenomenology that has been successfully accounted for with only GR.
Nowadays, we have strict constraints on the ensuing dynamics of the early Universe that led to the CMB [3] and we have tight constraints on Solar System physics [59,60,65,66,67,68,69]. Thus, to assess the statuses of the models discussed below, we need to also take into account whether their proposed solutions to any of the CCPs is able to accommodate the CMB data together with the success of the Λ CDM model, and whether these models contain a screening mechanism, as described in Section 3.2, that leaves local physics on small astrophysical scales untouched. We will call these requirements cosmic history constraints (CHCs) and astrophysical constraints (ACs).
Combining all the requirements above, we thus have a wishlist with five conditions that a successful model should satisfy. As we now shift our discussion towards how different modified gravity proposals can address the CCPs, these are the conditions that we should keep in mind while assessing how successful these different proposals are. Later, in Section 5, we summarize this wishlist and we briefly review which requirements each approach is able to tackle (see Table 1).

4.2. Constraining Gravity

There is a class of approaches that modify GR by constraining its dynamics such that the vacuum energy does not gravitate. We will discuss four of them: global and local sequestering mechanisms [106,107,108], unimodular gravity [109,110], and the nonlocal approach introduced in [111].
As we will see, local sequestering is the the most successful proposal addressing the CCPs discussed in Section 2. Thus, this section starts by discussing the global sequestering mechanism where the physical intuition is clearer as a way to further motivate the introduction of the local mechanism. Finally, since both unimodular gravity and the nonlocal approach can be easily related to the sequestering mechanism, they are discussed comparatively.
Naturally, this does not exhaust the list of proposals in the same vein. More recent proposals include [112], that promotes the Plank mass to a Lagrange multiplier, thus constraining the averaged Ricci scalar and preventing the vacuum energy from gravitating. Meanwhile, [113,114] promote the CC to an integration constant, but instead of canceling the vacuum energy, they present a mechanism that relaxes it through multiple membrane nucleations (each localized in time).

4.2.1. The Global Vacuum Energy Sequestering

One way to avoid the no-go theorems discussed in Section 3.1.2 was introduced in [106], where the vacuum energy generated by matter loops is gravitationally decoupled. This mechanism is called global vacuum energy sequestering and is based on promoting the bare cosmological constant, Λ b , and the Planck mass, M Pl , to Lagrange multipliers. It further introduces a global interaction term σ as a function of Λ b and a mass scale μ around the QFT cutoff scale. Thus, this is a minimal modification of GR as it does not introduce any new propagating degrees of freedom. For clarity, we denote the Planck mass by κ in this and the next sections to indicate that it is now a variable.
In the Jordan conformal gauge, the action reads [108]
S = d 4 x g κ 2 2 R Λ b + L m + σ Λ b μ 4 ,
where L m stands for the matter Lagrangian, and we note that σ is outside the integral. Varying it with respect to Λ b and κ 2 , respectively, yields two global constraints:
1 μ 4 d σ d Λ b = d 4 x g , d 4 x g R = 0 .
Due to the presence of the smooth function σ , the first constraint does not force the world volume to vanish. Nonetheless, in order to have a large and old universe, it is required that σ cannot be a linear function; otherwise, μ would have to be very small, lying at scales many orders of magnitude below particle physics scales and the cutoff of this effective field theory.
Meanwhile, the second equation can be recast as R = 0 , where the four-volume average is defined as ( ) d 4 x g ( ) / d 4 x g . Thus, the scalar curvature is always averaged to zero. Naturally, this is only well-defined if the four-volume is nonvanishing, imposing that d σ / d Λ b is nonzero.
Now, varying the action with respect to the metric produces the gravitational field equations, which, upon taking the trace, averaging, and using the second global constraint in Equation (46), leads to
Λ b = 1 4 T μ μ μ ,
completely fixing the bare cosmological constant in terms of the matter sources15. This is only possible because Λ b is considered as a dynamical variable whose value is determined by Equation (47). Then, we can finally rewrite the Einstein field equations as
κ 2 G μ ν μ = T μ ν μ 1 4 δ μ ν μ T α α α ,
where the bare cosmological constant is a global dynamical field fixed by the field equations. Notably, as the matter source is shifted by the last term on the RHS, both classical and quantum contributions that take the form of a cosmological constant will not gravitate16. In particular, the latter is guaranteed due to general covariance, which imposes that the loop corrections contribute in the same way at any order of the QFT loop expansion [24].
More explicitly, as is shown in [108], we can consider the matter Lagrangian at any given order in loops to be split between the renormalized quantum vacuum energy contributions, ρ vac , and local excitations, which leads to T μ ν μ = ρ vac δ μ ν μ + τ μ ν μ . Then, the field equations Equation (48) become
κ 2 G μ ν μ = τ μ ν μ 1 4 δ μ ν μ τ α α α ,
and the regularized vacuum energy drops as a source. Nonetheless, we still have a residual cosmological constant above given by the second term on the RHS which arises as the historic average of the trace of matter excitations. Thus, it is a nonlocal source both in space and time, and it is typically small in universes which grow large and old due to the spacetime averaging. Therefore, global sequestering could potentially also address the DE problem.
Something similar happens with contributions coming from phase transitions. They are not completely suppressed in the field equations Equation (49), but for large and old universes they become far smaller than the current critical density after the transition [106]. The intuition here is that the contribution is suppressed by the spacetime volume, and typically the Universe spends a relatively short time in the false vacuum [117]. Note that before the transition, despite these contributions not being sequestered, within reasonable assumptions they remain consistent with early universe cosmological data [118].
We make some additional remarks about the above construction:
  • The first constraint in Equation (46) imposes that the world volume is finite as we assume a smooth global interaction σ . This necessarily selects a universe with spatially compact sections that is temporally finite (see [119] for a model realizing this condition);
  • The action Equation (45) is not additive over spacetime due to the global interaction term σ , leading to subtleties for its quantization and its embedding in a complete UV theory.
Both these challenges can be evaded by the local version of the sequestering mechanics discussed below. For more discussions about the global proposal and its cosmological implications, see [118].

4.2.2. Local Sequestering

In the original formulation of vacuum energy sequestering, the smooth function σ ( x ) is added directly to the gravitational action rather than to its Lagrangian, while κ 2 and Λ are rigid Lagrangian multipliers. A local formulation of energy sequestering was proposed in [107], where σ ( x ) has a local form and the Lagrangian multipliers are promoted to fields. This is accomplished after modifying the action Equation (45) to
S = d 4 x g κ 2 ( x ) 2 R Λ ( x ) + L m + σ Λ ( x ) μ 4 F ( 4 ) + σ ^ κ 2 ( x ) M Pl 2 F ^ ( 4 ) ,
where F ( 4 ) and F ^ ( 4 ) are two four-forms that satisfy usual Bianchi identities17, i.e., locally F ( 4 ) = d A ( 3 ) and F ^ ( 4 ) = d A ^ ( 3 ) , and the functions σ and σ ^ are assumed to be smooth. Note that Λ ( x ) and κ 2 ( x ) are spacetime functions whose equations of motion are, respectively,
σ μ 4 F μ ν ρ σ = g ε μ ν ρ σ , σ ^ M Pl 2 F ^ μ ν ρ σ = 1 2 R g ε μ ν ρ σ ,
where ε μ ν ρ σ is the four-dimensional Levi–Civita symbol. The equation of motion for the metric is
κ 2 G ν μ = T ν μ + ( μ ν δ ν μ 2 ) κ 2 ( x ) Λ ( x ) δ ν μ ,
where T μ ν is the matter energy-momentum tensor. Due to the topological nature of the four-form actions, there are no flux contributions on the RHS in the above.
On the other hand, it follows from the equations of motion for A ( 3 ) and A ^ ( 3 ) ,
σ μ 4 μ Λ = 0 = σ ^ M Pl 2 μ κ 2 ,
that Λ and κ 2 are constants on-shell. This makes κ set the bare Planck-scale value and Λ play a similar role to the Λ b in the global case. Meanwhile, the traceful part of the metric equation gives
Λ = 1 4 T α α + 1 4 κ 2 R ,
where the last term can also be written in terms of the four-form fluxes as
1 4 κ 2 R = μ 4 κ 2 2 M Pl 2 σ ^ σ F ^ ( 4 ) F ( 4 ) Δ Λ .
Hence, compared to the global case, Equation (48) becomes modified to
κ 2 G ν μ = T ν μ 1 4 δ ν μ T α α Δ Λ δ ν μ .
The cancellation of the matter loop corrections to the cosmological constant operates as in the global case, but there is now an extra term Δ Λ in the residual cosmological constant. This is the main difference between the global and local approaches to vacuum energy sequestering. The other difference is that the first constraint in Equation (46) is absent such that infinite spacetime volumes can be considered.
Although the Λ and κ 2 dependence of Δ Λ makes it UV-sensitive, the authors of [107] argue that the smoothness of σ and σ ^ guarantees that the variation of the prefactors in Δ Λ , that depend actually on the dimensionless variables κ 2 / M Pl 2 and Λ / μ 4 , are bounded by O ( 1 ) numbers. Moreover, note that the fluxes are UV-insensitive to the choice of UV cutoff, being dominated by the IR scale in the integrals. For instance, for a constant F ( 4 ) and in a spacetime where R is bounded, both F ( 4 ) and F ^ ( 4 ) fluxes diverge with the spacetime volume, such that their ratio can be finite. In summary, the value of Δ Λ can be finite, small, it is UV-stable, and should be ultimately determined by observations. Finding support in the latter is a challenge for sequestering, due to its similarity with GR and, hence, difficulties in finding unique phenomenological signatures. Another challenge is to come up with a UV embedding (see [120,121] for attempts in the context of axion monodromy).

4.2.3. Nonlocal Approach

Another constrained approach, related to sequestering, was proposed by Carroll and Remmen [111] (CR), where a nonlocal constraint is applied to the action, and the averaged Lagrangian density is forced to vanish on-shell. The action being considered is
S CR = η d 4 x g M Pl 2 2 ( R 2 Λ b ) + L m 1 48 F ( 4 ) 2 + 1 6 · ( F ( 4 ) A ( 3 ) ) ,
where the parameter η is a constant that acts as a global Lagrange multiplier. It generalizes the conventional measure through g d 4 x g η d 4 x . In very general terms, the idea of introducing η here is to enforce a cancellation between the four-form F ( 4 ) and Λ b in the equation of motion18.
Although Equation (57) is different from Equations (45) and (50), its connection to both global and local sequestering can be seen by recasting it in an equivalent form (ignoring surface terms) as [108]
S CR = η d 4 x g M Pl 2 2 R Λ ˜ + L m + η σ Λ ˜ μ 4 F ( 4 ) ,
where Λ ˜ = θ 2 / 2 is defined as a new field variable with θ the magnetic dual of F ( 4 ) , Λ b is absorbed by L m , we introduced a function σ 2 ( z ) = 2 z , and the four-form is rescaled as F ( 4 ) F ( 4 ) / μ 2 . Comparing with the sequestering proposals, we see that Equation (58) bears resemblance with both the global and local sequestering mechanism. One of the flux terms used in local sequestering is recovered here, but η remains a global constraint and Λ ˜ is also a global variable. Although η can be made into a local parameter [122], it would not resolve the main problem of the CR model, i.e., radiative instability, that will be discussed below.
Going back to CR’s approach, we proceed from Equation (57) to see how the CC is forced to vanish. The constant scalar field η admits, upon regularization d 4 x g V , an equation of motion that forces the averaged Lagrangian to vanish, i.e.,
1 V d 4 x g L = L = 0 .
Now, introducing a constant scalar θ such that F μ ν ρ σ = θ g ε μ ν ρ σ (obviously fulfilling the flux equations), we could rewrite the last two terms in Equation (57) as
1 48 F μ ν ρ σ 2 = 1 2 θ 2 , 1 6 μ ( F μ ν ρ σ A ν ρ σ ) = θ 2 .
These can be substituted into the constraint Equation (59) to give
M Pl 2 2 R 2 Λ b + L m 1 2 θ 2 = 0 ,
which can be read as fixing θ in terms of Λ b . The Einstein equation in this model reads
0 = M Pl 2 2 R μ ν 1 2 R g μ ν + Λ b g μ ν 1 2 T μ ν + 1 4 g μ ν θ 2 .
The crucial observation is that Λ b is canceled when we use Equation (61) to substitute for θ
R μ ν 1 2 R g μ ν + 1 2 R g μ ν = M Pl 2 T μ ν L m g μ ν .
Comparing Equation (62) and Equation (63), an effective CC can be identified in terms of the averaged quantities as
Λ eff = 1 2 R + M Pl 2 L m .
We see that the CR model indeed achieves a cancellation of the tree-level CC and hence features a form of self-tuning. At the same time, it remains compatible with all known observations, since standard GR solutions including the Friedmann equations are unchanged in this model. It also circumvents Weinberg’s no-go, as the global parameter η alters the spacetime measure, similar to what is performed in sequestering. In terms of its nonlocality and the implication to the DE problem and phase transitions, it is similar to the case of global sequestering (see discussions in Section 4.2.1).
On the other hand, CR’s mechanism is not stable under radiative corrections. A simple and elegant argument was provided in [108]: as η multiplies the whole action, it can be identified with 1 / , corresponding to the 0-th order in an expansion. Now, since the loop-corrected action contains higher powers in it will also come with higher powers of 1 / η . This, however, changes the structure of the whole theory. In particular, a variation with respect to η will no longer yield Equation (63).
A remedy to this problem is to have a purely geometrical global constraint, which can be achieved, for example, by an action such as
S hybrid = d 4 x g η M Pl 2 2 R M Pl 2 Λ b + L m 1 48 F ( 4 ) 2 + 1 6 · ( F ( 4 ) A ( 3 ) ) ,
where the Lagrange multiplier is only attached to R. This however, takes us back to sequestering. To be precise, when we identify η M Pl 2 = κ 2 and perform the same steps that led to (58), we obtain a hybrid between the local and global sequestering mechanisms.

4.2.4. Unimodular Gravity

Unimodular gravity is a modification to GR already proposed by Einstein in 1919 [109]. Similar to sequestering models where nonlocal constraints are imposed, in unimodular gravity the determinant of spacetime is constrained to be constant, g = const . This condition is known as the unimodular condition and, as we will see below, leads to the same field equations as those of GR with a CC. In the following, we present the basic premises of unimodular gravity at a classical level, as well as its connection to the CCP. An extensive overview of the details of quantum unimodular gravity lies beyond the scope of the present review. For further details on the quantization of unimodular gravity, see [123,124,125,126], and for criticisms, see [124,127].
Let us start from the Einstein–Hilbert action (without a CC) in the presence of matter and including a Lagrange multiplier λ ( x ) [124],
S EH = d 4 x g M Pl 2 2 R λ ( x ) g ϵ 0 + S m .
Here, S m refers to the action containing matter fields and ϵ 0 is a constant volume element. Varying Equation (66) with respect to the metric leads to the field equations
R μ ν 1 2 g μ ν R = 1 M Pl 2 T μ ν λ ( x ) 2 g μ ν ,
while varying with respect to the Lagrange multiplier leads to g = ϵ 0 . Note that this unimodular condition breaks the full diffeomorphism invariance of the theory. In Equation (67), the EMT is given by the variation of the matter action, S m , as shown in Equation (4). Taking the divergence of Equation (67) and accounting for energy-momentum conservation, μ T μ ν = 0 , leads to
μ λ ( x ) = 0 ,
with solution
λ = constant .
In other words, the Lagrange multiplier is fixed to be a constant and the CC can be identified with Λ = λ / 2 .
An alternative derivation of the field equations is performed in the Henneaux–Teitelboim formulation of unimodular gravity [123,128,129], where full diffeomorphism invariance remains unbroken. The action of the theory is written as [129]
S HT = d 4 x g M Pl 2 2 R 2 ϕ + L m + M Pl 2 ϕ μ A ˜ μ ,
where A ˜ α = 1 6 ϵ α β γ δ A β γ δ is the vector density associated with a three-form A α β γ . In this formalism, the unimodular condition arises from the equation of motion of the scalar field ϕ ( x ) , as g = μ A ˜ μ . Additionally, the equation of motion of the vector density field A ˜ μ is μ ϕ ( x ) = 0 . The latter, again, shows that ϕ = const . can be identified as the cosmological constant Λ . Finally, varying the action with respect to the metric g μ ν leads to the Einstein equations with a CC.
Thus, it becomes clear that, within unimodular gravity, the CC arises as an integration constant in the Einstein field equations. It has been argued that this provides a conceptual exit from the CCP [110,123]. A way to understand this argument is by stating that any change in the EMT of the form T μ ν T μ ν = T μ ν + g μ ν C can be absorbed into the integration constant Λ by a shift Λ Λ + C / M Pl 2 . This is interpreted as the curvature of spacetime not coupling to quantum corrections of the form g μ ν C . The claim, then, is that such a shift can address the CCP by canceling off a given amount of vacuum energy. Here, we should note that the same logic can be applied in GR, where Λ , again, is a free constant whose value can be chosen arbitrarily.
The argument above could, indeed, constitute a solution to the CCP, if the latter was limited to its old formulation, as explained in detail in Section 2. In other words, if the issue merely consisted of the fact that accounting for the vacuum energy contribution to the EMT leads to a much larger CC than the one measured by observations today, such a problem could be resolved by using the freedom of fine-tuning an integration constant to match the observed value. However, a solution cannot be provided when one considers the new-CCP. In the latter, the renormalization prescription of the effective theory is unstable under higher-order loop corrections, requiring the bare value of the cosmological constant to be retuned every time an additional loop contribution is taken into account. Since in unimodular gravity (such as in GR) the integration constant playing the role of the CC can only be fixed once, the conceptual problem remains [124,130]. In fact, the same difficulty is encountered when attempting to solve the classical CCP in the unimodular gravity framework. This should not come as a surprise since unimodular gravity constitutes a minimal modification to GR that is also not dynamical in any way. Hence, unimodular gravity is only able to resolve the CCP problem to the same extent as standard GR itself.

4.3. Massive Gravity

Massive gravity theories are a natural extension of GR, and the formulation is of inherent theoretical interest. However, massive gravity also provides an interesting angle to the cosmological constant problems via the degravitation mechanism, where the graviton mass acts as a filtering scale in the coupling between gravity and a cosmological constant term. In that sense, the right question to ask about the CCPs might not be why the cosmological constant term is so small but, rather, why it gravitates so little. Within GR, diffeomorphism invariance guarantees that gravitation couples universally to all sources, but in more generalized theories of gravity, there is a meaningful difference. In the following, we briefly recap how this mechanism works, show how it emerges naturally in the context of massive gravity, and review the development and current status of massive gravity theories.
For extensive reviews about massive gravity, see [35,37,131,132,133].

4.3.1. Degravitation

The degravitation mechanism [134,135,136] can be seen as an IR modification of GR that prevents sources characterized by wavelengths λ S larger than a given IR scale L from gravitating. At the level of Einstein equations, it is phenomenologically implemented by promoting Newton’s constant to a scale-dependent differential operator
G μ ν = 8 π G N ( L 2 ) T μ ν .
Expanding the sources in terms of its mode functions, the gravitational coupling behaves as
G N ( L 2 λ S 2 ) = G N for λ S L 0 for λ S L ,
acting as a high-pass filter, thus preventing deep-IR modes, λ > L , from gravitating. The cosmological constant is a particular case with an infinite characteristic wavelength; thus, it degravitates. Therefore, regardless of the value of ρ vac , as long as it behaves as a cosmological constant, it does not source Einstein equations.
In the context of gravitational theories, this mechanism can be implemented by having massive gravitons or a resonance19 [135,136]. DGP, briefly discussed both in the context of derivative screening in Section 3.2.2 and of braneworld models in Section 4.5, can be considered a special case of the resonant graviton. In the rest of this section, we focus on the linear and nonlinear constructions of massive gravity.

4.3.2. Linear Massive Gravity

The search for a consistent generalization of GR that describes the behavior of massive spin-2 fields is an old problem. The Lagrangian for a linear theory of a massless spin-2 field h μ ν can be written as
L lin = h μ ν E μ ν α β h α β M Pl 1 h μ ν T μ ν ,
where the kinetic terms are given by
E μ ν α β = 1 2 η μ α η ν β η ν β μ α η μ β ν α + η μ ν α β + η α β μ ν η μ ν η α β .
Fierz and Pauli showed in 1939 [58] that the linear theory of a massive spin-2 field is uniquely given by
L FP = h μ ν E μ ν α β h α β m FP 2 h μ ν h μ ν h 2 M Pl 1 h μ ν T μ ν ,
with the Fierz–Pauli mass m FP , while other possible mass terms lead to the existence of ghosts. From the trace and the divergence of the equations of motion, we can derive
m FP 2 h μ ν = M Pl 1 T μ ν 1 3 η μ ν 1 m FP 2 μ ν T ,
which realizes the degravitation condition outlined above, where the source T μ ν is seen through a high-pass filter 1 m FP 2 1 . For sources with small associated scale λ S m FP , the filter plays a negligible role, while for sources λ S m FP , such as a cosmological constant, the coupling is suppressed. Considering nonrelativistic sources in Equation (76) leads to a Yukawa-type potential for the gravitational interaction [135]
V ( r ) exp m FP r r ,
which again suppresses the strength of the gravitational interaction on large scales. For the theory to be compatible with a standard cosmic evolution, the graviton mass should not exceed m FP H 0 10 33 eV .
Unfortunately, not all is well with the linearized theory. It is already apparent in Equation (76) that the limit m FP 0 is not smooth, as pointed out by van Dam, Veltman, and Zakharov (the vDVZ discontinuity, [137,138]). A careful accounting of the propagating degrees of freedom in the theory (e.g., by using the Stückelberg formalism [139]) demonstrates that the massless limit of Fierz–Pauli gravity is indeed not GR, but GR plus an additional attractive scalar field. This leads to a mismatch with observation: if the coupling to nonrelativistic sources is supposed to match Newtonian theory (where both the tensor and the new scalar field couple), then gravitational lensing in the new theory is weaker by a factor of 3 / 4 [35] since the scalar does not couple to the traceless energy-momentum tensor of light.
However, the additional interaction responsible for this discrepancy differs once higher orders of h μ ν are included in the action. Nonlinear interactions lead to a suppression of the additional scalar degree of freedom inside the Vainshtein radius around massive sources [98], independent of the specific nonlinear completion [140], as explained in Section 3.2.2. As a result, GR is expected to be restored close to massive sources such as the Earth (at the price of leaving the regime of applicability of the linear theory).

4.3.3. Nonlinear Theories

The question is, then, whether a nonlinear completion of the Lagrangian given by Equation (75) can be brought in agreement with observations and theoretical consistency requirements. After all, most nonlinear extensions are plagued by the Boulware–Deser ghost, since the unique structure of the mass term that guarantees the consistency of the Fierz–Pauli theory is generally spoiled by higher-order corrections [141].
A solution was found by de Rham, Gabadaze, and Tolley (dRGT) [142,143] by constructing a nonlinear generalization of Equation (75) around an arbitrary fixed background metric, and Hassan and Rosen proved the absence of ghosts in the theory to all orders [144]. Nonetheless, even for ghost-free theories, demanding unitarity and analyticity of the S-matrix severely constrains the possible graviton mass scale [131,145,146]. When combined with observational constraints from propagation of gravitational waves [147,148] and solar system measurements [149], it seems that the theory becomes nonviable.
Making the reference metric of dRGT gravity itself dynamical leads to bimetric theories of gravity with interacting massless and massive spin-2 fields [131,144] (and corresponding multigravity extensions with additional tensor fields (e.g., [150])). The physical metric responsible for gravitational effects is a mixture of the two tensor fields, and the behavior in the linearized regime is governed by a mixing angle between massive and massless modes. As the mixing becomes small, one recovers the massless spin-2 predictions of GR, while a large mixing leads to the phenomenology of linear massive gravity.
The free mixing angle allows bimetric gravity to circumvent the existing limits on the graviton mass that are so problematic for dRGT. On the other hand, a large mixing is required for degravitation in the linear regime [151], where the theory approaches the general phenomenology of linear massive gravity, with all advantages and problems connected to that: effects of a bare cosmological constant can be made to degravitate, but the additional fifth force from Fierz–Pauli theory described in Section 4.3.2 reappears. Within the Vainshtein radius, the gravitational effect of Λ does not degravitate—depending on how large the Vainshtein radius is chosen, this can slightly alleviate the problem, since upper bounds on vacuum energies from, e.g., planetary orbits are 16 orders of magnitude larger than the cosmological bound [24].
Setting constraints from cosmological structure formation remains difficult since perturbations in bigravity become nonlinear very early on. While these instabilities do not necessarily rule out the theory [152], there is no framework available yet to calculate nonlinear predictions for cosmological scales in generality. It is, however, possible to use probes of the cosmic background expansion such as BBN [153], BAO, or supernova data [154,155] to constrain the theory. There is a large parameter space available that produces a Λ CDM-like expansion history, but this parameter space cannot degravitate CC effects while also fitting the available data.
In absence of degravitation, the problem of UV sensitivity of vacuum energy contributions remains. An interesting new angle is that if the problem of UV sensitivity is solved, interaction terms between the two metrics provide an additional source that behaves as a cosmological constant, but is itself protected from quantum corrections [144].

4.4. Self-Tuning with Horndeski Theories

The essence of self-tuning is to introduce a dynamical degree of freedom to GR, typically a scalar field ϕ ( t ) , which is used to cancel the contribution of the CC’s energy density ρ vac in the evolution equations of the Universe, effectively decoupling it from gravity in the cosmological context. We note here that we have broken Poincaré invariance by enforcing a nontrivial time dependence of ϕ , which allows us to bypass Weinberg’s theorem (see Section 3.1.2). The cancellation relies on trivially satisfying the equation of motion for ϕ , or making the equation redundant. Once this is achieved, the first Friedmann equation defines the relation ϕ ˙ = f ( ρ vac ) , where f is a general function, which provides a dynamical constraint equation for ϕ in terms of ρ vac . This allows ϕ to dynamically match ρ vac across cosmic evolution and through any phase transition the Universe may undergo. Such dynamical matching is what we will refer to in this section as self-tuning. Note that this approach is not restricted to the cosmological setting, as we shall shortly see, although its applicability to multiple settings simultaneously is not guaranteed.
The models we look at here will all fall under the Horndeski scalar–tensor class of theories (see Equation (26)) which ensures the model is ghost-instability free. We note there have been significant extensions to this class [76,77,78,156,157,158,159,160,161,162,163] which we do not consider here. The Horndeski class as well as these beyond-Horndeski models have been well constrained by data, in particular by measurements of the propagation speed of gravitational waves (see [164] for a review on recent developments). Despite this, a large theory space remains open with interesting phenomenology which can be tested with cosmology.
Horndeski theory permits four free functions, G i ( ϕ , X ) , of ϕ and its canonical kinetic energy X = ( ϕ ) 2 / 2 in the Lagrangian. In this context, for a flat Friedman–Lemaître–Robertson–Walker (FLRW) background, we have the following three coupled differential equations coming from the 00-th component of the metric field equations, their trace, and the variation of the Lagrangian with respect to the scalar field, respectively,
F 1 ( H , ϕ , ϕ ˙ , ρ m , ρ vac ) = 0 ,
F 2 ( H , H ˙ , ϕ , ϕ ˙ , ϕ ¨ , ρ m , ρ vac ) = 0 ,
F 3 ( H , H ˙ , ϕ , ϕ ˙ , ϕ ¨ ) = 0 ,
where we assumed a pressureless matter component. We then assume that the scale factor has an attractor solution. Given current observations, one could choose this to be a de Sitter solution, i.e., a = e H 0 t , which describes the asymptotic future of a universe containing a cosmological constant. This being said, equally pertinent to the CCPs is the need for such self-tuning to occur in the Solar System (so, within black hole solutions [165] or the Minkowski solution [166], for example).
Considering the de Sitter attractor, a viable self-tuning model is restricted to the subset of Horndeski functions that fulfill the following (see Section 7 of [167], for example):
  • The equation of motion for ϕ Equation (78c) is trivially satisfied on the attractor or redundant with the other two equations;
  • The Friedmann equation Equation (78a) should depend on ϕ ˙ for the tuning to be dynamic;
  • The Friedmann equation must admit nontrivial expansion histories before hitting the attractor solution;
  • The theory incorporates a screening mechanism (see Section 3.2).
The first condition leaves us with only two dynamical equations on the attractor. In particular, F 1 attractor is then allowed to set the value of ϕ ˙ in terms of ρ vac , i.e., we have ϕ ˙ = f ( ρ vac ) . Moreover, any change in ρ vac , even a discontinuous one, appropriately changes the value of ϕ , allowing for self-tuning over phase transitions.
Given these basic criteria, we now look at two different models in the literature that achieve self-tuning. These essentially differ in how they achieve the first criterion.

4.4.1. The Fab-4

One of the first instances of self-tuning in this context was the Fab-4, originally proposed in [166,168] and named in homage to the Beatles. The Lagrangian functions were named after the band members
L john = g V john ( ϕ ) G μ ν μ ϕ ν ϕ ,
L paul = g V paul ( ϕ ) P μ ν α β μ ϕ ν ϕ α ϕ β ϕ ,
L george = g V george ( ϕ ) R ,
L ringo = g V ringo ( ϕ ) G ^ ,
where R is the Ricci scalar, G μ ν the Einstein tensor, P μ ν α β the double dual of the Riemann tensor, and G ^ the Gauss–Bonnet scalar. The various potentials are chosen to allow for self-tuning on the attractor, while the coupling to curvature terms is how the trivial solution to the Klein–Gordon equation is ensured.
The original setup assumed a Minkowski vacuum attractor solution, which we consider for illustration. We find F 3 = 0 (see Equation (78c)) identically on vanishing curvature backgrounds as each term is coupled to curvature (see [169] to see this explicitly). The Friedmann equation on the attractor is then found to be [169]
F 1 attractor ( H , ϕ , ϕ ˙ , ρ m , ρ vac ) = V john ( ϕ ) ( ϕ ˙ H ) 2 + V paul ( ϕ ) ( ϕ ˙ H ) 3 V george ( ϕ ) ( ϕ ˙ H ) + ρ vac = 0 ,
which places a constraint on the forms of the possible potentials V ( ϕ ) . These can also be chosen so as to reproduce the correct eras in our Universe’s expansion history as according to observations [170]. Further, the nontrivial derivative interactions of the john and paul terms leave the possibility for Vainshtein screening, although this has not been shown to be viable in tandem with self-tuning explicitly [57], and is generally problematic for light scalar fields [171].
Finally, we remark that the recent measurement of the speed of gravitational waves [147,172] has severely constrained the Fab-4 model [173], although the applicability of these constraints for such theories has been challenged in [174,175].

4.4.2. Well-Tempered Self-Tuning

The second model we present arose to address some of the initial problems of the Fab-4, in particular its struggle in consistently reproducing all of the correct eras of the expansion history together while also canceling a large value of ρ vac [176]. This being said, specific classes of scalar field potentials in Fab-4 also ameliorate this issue [170]. In [177,178,179], the authors attempt to moderate the self-tuning in order to allow for the correct expansion history. The model of [177] achieves the first self-tuning criterion by making the identification on the attractor
F 2 attractor ( H , ϕ , ϕ ˙ , ϕ ¨ , ρ m ) = F 3 attractor ( H , ϕ , ϕ ˙ , ϕ ¨ ) ,
where we removed the ρ vac dependency using Equation (78a) and assumed a de Sitter attractor with H ˙ = 0 . This identification makes Equation (78b) and Equation (78c) degenerate and imposes a constraint on the Horndeski functions allowed in this scenario. The advantage of making these two equations degenerate only on the attractor, rather than have a trivial solution to Equation (78c), as in the Fab-4, is that it implicitly requires ρ m = 0 for the degeneracy to hold, i.e., the scalar field only screens ρ vac , allowing for a nontrivial expansion history.
Further, the constraint on the Horndeski class imposed by Equation (81) does not preclude a screening mechanism. In fact, the original proposal of [177] considered the inclusion of G 2 ( ϕ , X ) and G 3 ( ϕ , X ) Lagrangian terms which give the capacity for Vainshtein and chameleon screening, although such a mechanism has not been shown to be possible explicitly. Lastly, we note that the functional forms of G 2 and G 3 should also be such that a ϕ ˙ dependency remains in Equation (78a).

4.4.3. Outlook

Horndeski self-tuning can grant a removal of the CC contribution from various spacetimes, including Minkowski and FLRW, and is robust to phase transitions. Further, it is able to recover a nontrivial expansion history in the cosmological context, which allows for late-time acceleration. The models also have the capacity to screen any deviations from GR at small physical scales in accordance with solar system tests. Lastly, the presence of a shift symmetry in some such models, for example, that of the model described in [180]
ϕ ϕ + a , Λ Λ + a c 1 M 3 λ ,
where M is a mass of ϕ and c 1 the scalar field potential’s Lagrangian coefficient, may be useful in addressing the impact of quantum corrections [41].
However, despite the Horndeski self-tuning program being an interesting approach, it leaves many things to be desired. Self-tuning requires tuning of the mass scales and Lagrangian coefficients in order to recover the correct expansion history and effectively cancel the CC. These parameters are not guaranteed to be stable against radiative corrections (see [180] for an example of such parameter tunings).
On this point, in these models, one may assume that ρ vac is the value coming from the UV-complete theory, with all corrections to the CC accounted for. The solution for ϕ can then adapt dynamically to this value. This then says that our theory is “stable” against quantum corrections to the vacuum energy. The question is then whether or not the parameters of the theory, for example, the scalar field potential or couplings to curvature, are stable against such corrections. One can force the coupling to the matter sector of the field to be weak (see, for example, [181]), which helps protect the theory against radiative corrections and, at the same time, avoids fifth-force constraints, although it may have implications for phase transitions as the response to changes in vacuum energy is slowed down.
An issue already pointed out is that the mechanism cancels ρ vac on the attractor but does little to address the large value of ρ Λ away from the attractor, at, say, earlier times, which would change the expansion history significantly (see [182] for a discussion). Similarly, self-tuning, as discussed here, assumes some specific attractor solution at which the CC is canceled by ϕ rather than at the level of the field equations, as in the sequestering approach (see Section 4.2). This leaves the issue of the CC in all solutions but the attractor, which is clearly not sufficient as we do not observe a large CC contribution both on cosmological and solar system scales. Related to this is the issue of screening. As seen in Section 3.2, we require a choice of Horndeski functions that effectively cancel ϕ ’s contribution to the Poisson equation (see Equation (30)) either through a well-chosen potential or derivative term. Horndeski functions that screen this contribution to the Poisson equation and self-tune the CC are yet to be shown explicitly.
Finally, we note that there have been many interesting proposals that take different approaches to self-tuning, such as [33,89,181,183,184,185,186], which we did not discuss here but are worth considering. They all involve some cancellation of the CC using an additional degree of freedom, with some, such as [33,112], making a connection with sequestering and others invoking higher dimensions, such as [185].

4.5. Braneworld Models

The no-go theorems in Section 3.1.2 assumed vacuum energy to be constant throughout space and time. In a nutshell, extra-dimensional models, or, more specifically, braneworld models, relax that assumption by breaking the translational invariance of vacuum energy in the direction of the extra space. This is possible because the Standard Model, and, hence, our Universe, is assumed to be confined on a spatial hypersurface, referred to as a brane, that is embedded in a higher-dimensional bulk spacetime. In fact, in the simplest models, only gravity is allowed to propagate in the bulk spacetime. As a consequence, the matter loops that give rise to vacuum energy (and its quantum corrections) are equally confined on the brane. Since vacuum energy is still constant along the intrinsic brane directions, it acts as a surface or brane tension λ . Crucially, this implies that it gravitates differently from a space-filling vacuum energy in a four-dimensional (4D) theory.
We are mostly interested in the low-energy phenomenology of braneworld models; after all, the cosmological constant is an extreme IR source. It should still be noted that branes are important building blocks of string theory, which makes this class of solutions particularly interesting from a high-energy perspective. In fact, D-branes are nonperturbative states in string theory and appear both in tentative string constructions of SM [187,188] and four-dimensional maximally symmetric spaces such as dS 4 [189,190]. Since they source the p-form massless fields in the theory, D-branes have become an important ingredient in flux compactifications and string phenomenology [191,192]. Discussing the fate of the CCP in string theory is beyond the scope of this review, and we are interested in extra-dimensional models regardless of their embedding into string theory. Moreover, the literature on braneworld models as solutions to the CC problem is rather abundant [193]. In this review, we therefore focus on six-dimensional models that rely on a simple geometric mechanism, featuring conical deficit angles, to achieve self-tuning. For other models in five dimensions (and a discussion of their shortcomings), see [194,195,196,197,198,199], and, more recently, also [183,185]20.
In Section 4.5.1, we start with a simple toy model to illustrate the general idea. In Section 4.5.2, we discuss a more complete model before providing a general outlook in Section 4.5.3.

4.5.1. Our Universe as a Cosmic String in Six Dimensions

To illustrate the six-dimensional mechanism (and to introduce the necessary notation), we consider the simple case where our 4D Universe is a pure tension brane in a six-dimensional (6D) bulk spacetime of infinite size. The corresponding action, S = S bulk + S brane , can be decomposed into a bulk and brane part, respectively,
S bulk = M 6 4 2 d 6 X g 6 R 6 ,
S brane = d 4 x g 4 λ bare + L m .
Here, M 6 is the fundamental gravity scale in the bulk and R 6 is the 6D Ricci scalar. We further distinguish the bulk and brane geometries that come with coordinates and metrics X A , ( g 6 ) A B and x μ , ( g 4 ) μ ν , respectively. The brane and bulk metric are related through the pull-back ( g 4 ) = f A x μ f B x ν ( g 6 ) A B , where f A ( x μ ) is an embedding function. To solve the vacuum system, it is sufficient to make the following ansatz for the 6D metric:
d s 2 = ( g 4 ) μ ν d x μ d x ν + d r 2 + C ( r ) 2 d φ 2 ,
where x μ and X A = ( x μ , r , φ ) are the brane and bulk coordinates, respectively. For simplicity, we also assumed a vanishing bulk cosmological constant (in the absence of massive bulk fields, this is a radiatively stable choice). Moreover, we assume the induced metric on the brane, ( g 4 ) μ ν to be maximally symmetric, i.e., we allow for a de Sitter phase on the brane. With these definitions, the full trace of the 6D Einstein equation reads
M 6 4 R 4 + 2 C ( r ) C ( r ) = 2 λ δ ( r ) 2 π C ( r ) ,
where R 4 = const . is the Ricci scalar built from the brane metric ( g 4 ) μ ν , and λ = λ bare + ρ vac contains the bare tension and the vacuum energy arising from the matter Lagrangian. In other words, the brane-induced energy momentum tensor in vacuum takes on the form ( T 4 ) μ ν = ρ vac ( g 4 ) μ ν . The localization of the brane tension at r = 0 is described in terms of a Dirac delta function. The general solution of this equation that reduces to 6D Minkowski, i.e., C ( r ) = r , in the absence of sources is
R 4 = 0 ,
C ( r ) = 1 λ 2 π M 6 4 r .
This is the straightforward generalization of the geometry of a cosmic string in four dimensions [201,202,203], with our 4D Universe playing the role of the string (for a more detailed derivation in 6D see, for example, [204,205]). In any event, the remarkable observation is that the 4D curvature is vanishing despite the presence of vacuum energy on the brane. Stated differently, an observer on the brane that, similar to us, is ignorant of the presence of two extra dimensions will conclude that vacuum energy, if present, does not gravitate (seemingly violating the equivalence principle). In contrast, a 6D observer living in the bulk will notice a geometrical response: A circular path around the brane at fixed radius r has a reduced circumference of ( 2 π δ deficit ) r , quantified in terms of a deficit angle δ deficit = λ M 6 4 . In an embedding picture, such a geometry corresponds to a cone with the brane at its tip (left panel in Figure 2). This type of extradimensional self-tuning is also fully dynamical: if there is a phase transition on the brane or we change the particle content of our EFT, the deficit angle will adapt while the 4D curvature remains zero. In fact, any value of the vacuum energy that is below the gravity cutoff M 6 will be shielded, hence preventing the need to fine-tune the bare EFT parameters21.
Of course, this model cannot be the final answer. As it stands, its gravitational sector is purely six-dimensional, meaning that the Newtonian potential of a point source on the brane scales as 1 / | x | 3 rather than 1 / | x | , in clear contradiction with observations. In addition, we want the vacuum on the brane to be a de Sitter rather than a Minkowski vacuum. Two classes of models were proposed in the past that make use of this 6D mechanism and yet give rise to a 4D Newtonian potential on the brane. First, there are generalizations of the DGP model (originally proposed in 5D [70,71], see also Section 3.2) that embed the brane in an infinite volume bulk [204,207]. Here, 4D gravity is recovered due to an induced gravity term on the brane. This model, however, suffers from instabilities [208,209,210]. Second, there are models that build on the Arkani–Hamed–Dvali–Dimopolus (ADD) proposal [211,212], where the extra space is large but has a finite volume [205,213,214,215,216]. The arguably most interesting candidate among this second class of models is the supersymmetric large extra-dimension (SLED) model with two micron-sized extra dimensions that take the form of a rugby ball [217,218,219] (right panel in Figure 2). In fact, it has been viewed for several years as the leading alternative to anthropic resolutions of the cosmological constant problem with smoking gun signatures in both tabletop tests of gravity and collider experiments. In the next section, we therefore summarize the main features and shortcomings of SLED as a potential solution to the cosmological constant problem (for more details on phenomenological aspects of the model, see also the review in [52]).

4.5.2. Supersymmetric Large Extra Dimensions

The SLED proposal generalizes the action in Equation (83) to22 [218,219]
S bulk = d 6 X g 6 M 6 4 2 R 6 + ( M ϕ ) ( M ϕ ) + 1 4 e ϕ F M N F M N + 2 M 6 8 e 2 e ϕ ,
S branes = b = N , S d 4 x g 4 λ b ( ϕ ) + 1 2 A b ( ϕ ) ϵ m n F m n .
In particular, it now contains two branes, one at the south (S) and one at the north (N) pole of the rugby ball geometry. Other new features are the presence of a dilaton ϕ and a Maxwell field strength F M N with gauge coupling parameter e. The dilaton makes this low-energy model compatible with supersymmetry. To be specific, supersymmetry manifests itself at low energies in the form of a constant scaling symmetry ( g 6 ) M N ζ ( g 6 ) M N and e ϕ ζ 1 e ϕ under which S bulk ζ 2 S bulk . This also justifies the absence of a bulk cosmological constant, which would break this symmetry. The Maxwell field, on the other hand, is needed to compactify the extra space. Its flux winds around the compact extra dimensions and provides the pressure needed to prevent the rugby ball from collapsing under its own gravitational pull. In holding with the principles of an EFT expansion, the Maxwell flux also makes an appearance as an induced term on the brane, where the indices m , n only run over the two spatial extraspace coordinates, ϵ m n is the fully antisymmetric epsilon tensor with components ± 1 / | ( g 6 ) m n | , and A b ( ϕ ) controls the strength of the brane-induced flux. Both the brane tension and the induced flux term are allowed to couple to the dilaton. In fact, for the special choice
λ b = const and A b = Φ b e ϕ ,
the brane action S brane preserves scale invariance (along with the bulk action).
The emergence of 4D gravity is ensured by construction: Due to the compactness of the extra dimension, the model admits a normalizable and massless 4D graviton in its spectrum (see, for example, [211]). There is a also a continuum of Kaluza–Klein modes, the masses of which are set by the inverse size of the extra dimension. As a consequence, their contribution to the gravitational exchange amplitude is exponentially suppressed for distances | x | > 0 where 0 is the typical size of the compact directions. Post-Cavendish experiments place an upper limit of [223] 0 45 μ m .
The vacuum equations corresponding to Equation (87) can be solved by generalizing the ansatz in Equation (84) to
d s 2 = W 2 ( r ) ( g 4 ) μ ν d x μ d x ν + d r 2 + C 2 ( r ) d φ 2 ,
A = A φ ( r ) d φ ,
ϕ = ϕ ( r ) ,
where W ( r ) is a warping function and the Maxwell field A points in the azimuthal direction φ (corresponding to the angular direction of the rugby ball).
Instead of displaying the full set of equations that follow from Equation (89), here, we focus on the crucial aspects of the vacuum solution that are relevant for the model’s potential to resolve the cosmological constant problem. The Maxwell equation can be integrated trivially, yielding
F ρ θ = e ϕ Q C W 4 + 1 2 π b = N , S A b ( ϕ ) δ ( r r b ) ,
where Q is an integration constant. Due to the compactness of the extra space, the system has to be equipped with an additional equation describing the quantization of the Maxwell flux [219,224]:
Q d r e ϕ C W 4 + 1 2 π b = N , S A b ( ϕ ) e ϕ r = r b = n e ( n N ) .
Now we come to a subtle issue. From Equation (90) we see that there is a contribution to F A B that is proportional to a Dirac delta function. Since the Einstein and dilaton equations are sourced by a term that is proportional to F M N F M N , they will contain a divergence δ ( 0 ) . Physically, this is an artifact of treating the brane as an infinitesimally thin codimension-two object. This somewhat technical point has been studied in different ways: In [225], the authors introduced a counterterm that removed the singularity and made the distributional approach applicable again; in [226,227], the brane was microscopically resolved in terms of an explicit vortex model; and, finally, in [228], the brane was described in terms of a cylindrical codimension-one object stabilized through a brane-localized angular pressure23. In all studies, it was found that extra-dimensional self-tuning, i.e., R 4 = 0 , is guaranteed in the case of scale-invariant brane–dilaton couplings, as in Equation (88). This, however, is not what we were looking for: in the scale-invariant case, ϕ drops out of the flux quantization condition in Equation (91), turning it into a tuning relation on both brane tensions [225,230]:
1 λ N 2 π M 6 4 1 / 2 1 λ S 2 π M 6 4 1 / 2 + e 2 π b Φ b = n .
At this point, one should object that observations require the vacuum to be a (quasi) de Sitter rather than a Minkowski geometry. Accordingly, the condition R 4 = 0 should be relaxed to R 4 10 33 eV to account for dark energy. This also suggests that we should move away from scale invariance (which always implies that R 4 = 0 ). This possibility was investigated in [228], where the authors allowed for a class of couplings that deviated from Equation (88) (by, for example, making the brane tension dependent on ϕ ). Unfortunately, in these cases, in order for the extra-dimensional volume and the 4D curvature to fulfill their phenomenological bounds, Equation (92) has to be satisfied with a very high accuracy (although not exactly), effectively recovering the usual fine-tuning issue within standard GR.

4.5.3. Outlook

The above finding bears close resemblance with Weinberg’s no-go theorem discussed in Section 3.1.2. There, it was found that the self-tuning condition led to a scale-invariant potential, which could only be made compatible with a vanishing curvature by imposing a tuning (to be specific, the loop-sensitive quantity V 0 had to be tuned to zero in Equation (10) for R 4 = 0 to hold). Is this the end of extra-dimensional self-tuning in higher codimensional setups? Not necessarily. The general mechanism outlined in Section 4.5.1 still has merit. The problems arose when we compactified the bulk. Intuitively, this should not come as a surprise: at very large distances 0 (or energies below the Kaluza–Klein scale), the extra dimensions cannot be resolved and the theory is equivalent to a 4D theory with a massless graviton, making the usual no-gos applicable again. However, if we were to consider a different way of recovering 4D GR, which allowed for infinite volume extra dimensions, this reasoning would not apply (as the zero mode is expected to decouple). Unfortunately, as mentioned before, a simple generalization of the DGP model to higher codimensions suffers from ghost instabilities. At this point, and in order to end this section on a more positive note, we speculate that an alternative mechanism to recover a viable 4D phenomenology could be the “volcano”-trapping proposed in [231,232] generalized to higher codimensions. It relies on the idea of deforming the near-brane geometry such that it traps the gravitational field lines close to the brane. It should also be noted that other scenarios relying on five (rather than six or more) extra dimensions are still being actively explored. One interesting direction is provided by the holographic model in [185]. It avoids problems with naked bulk singularities, which plagued earlier codimension-one self-tuning models. The model’s phenomenology is currently being investigated [184,233], and it remains to be seen if it can provide a viable phenomenology in the presence of self-tuning.

5. Conclusions—You Can’t Always Get What You Want, or Can You?

The various aspects of the cosmological constant problem have remained, for decades now, the most challenging problems in theoretical physics. Combining general relativity and quantum field theory at low energies and then fitting cosmological data from both early- and late-times yields four different puzzles (Section 2). There is an abysmal mismatch between the prediction of the vacuum energy in QFT and the current observed value in cosmology; even if this gap is fixed by a bare cosmological constant at the level of Einstein equations, the fixing is unreliable due to both phase transitions in the early Universe and the UV-sensitivity of the effective field theory employed; finally, we still lack a clear understanding of what dark energy is made of. All these puzzles have different levels of relatedness depending on the theory one is considering.
Currently, these issues bring into question the phenomenology of Einstein equations applied to large scales and they also challenge the effective field theory approach to QFT. Since the Einstein equations connect the matter content with gravity, there are two potential paths out: either modifying gravity or changing the quantum-matter sector.
As the vacuum energy couples to gravity as any other form of energy in GR, it is natural to focus on different ways that we can modify gravity to tackle one or more of the CCPs. The main idea underlining these approaches, the self-tuning mechanism, makes gravity less sensitive to the presence of vacuum energy by invoking additional gravitational degrees of freedom. This can create problems by itself, as gravitational physics is well understood by only using GR both on astrophysical and cosmological scales. Thus, in order for these approaches to be successful, they also need to typically rely on some screening mechanism operative on local scales, while at the same time being able to recover the success of the current cosmological model describing the history of our Universe. Together with solving the CCPs, these requirements formed our wishlist (introduced in Section 4.1).
However, such modifications of gravity cannot be arbitrary. Weinberg pointed out that simply introducing self-tuning fields that, once coupled to matter, will cancel out the vacuum energy still leads to fine-tuning or incorrect phenomenology (Section 3.1.1). In fact, we also showed that his theorem can be generalized by analyzing the spectral decomposition of the gravitational exchange amplitude between probes and sources. Thus, only modifications of gravity that break one or more of the underlying assumptions of these theorems will yield viable models of self-tuning.
The core assumptions of these theorems (see Section 3.1.2 and Section 3.1.3) are that the vacuum exhibits translational invariance, the theory is Lorentz-invariant, the fields’ propagators take on their canonical form, there isminimal coupling between gravity and matter, and that spacetime is four-dimensional. The models discussed here break one or more of these. For instance, sequestering leads to a noncanonical propagator by a nonstandard coupling; linear and nonlinear massive gravity depart from GR by making the graviton massive, which introduces the Stückelberg field that has a nonconstant vacuum value; the Horndeski-type proposals depart from GR and break the translational invariance of their vacuum; and SLED assumes that spacetime has more dimensions than four.
As we reviewed these models, we were able to assess their success in fulfilling our wishlist. Table 1 summarizes our analysis. We emphasize that most of these models tackle the classical CCP, while still being able to recover the gravitational physics already contained in GR at small scales. Most of these approaches can also successfully describe our cosmic history (although there may not be such an analysis present in the literature, it can be easy to see that they mostly depart from GR at late-times). Finally, the dark energy problem may or may not be solved in these approaches as that depends on whether the models allow some residual-like CC term that would induce an accelerated expansion, though in general they do not readily recover the observed value for the effective CC. Nonetheless, these models can also be combined with quintessential models (see [49,50,51] for a review) that would then drive the dark energy dynamics.
Regarding the new CCP, the analysis becomes subtle. By modifying gravity such that the vacuum energy is effectively prevented from gravitating, these proposals tacitly accept that the UV sensitivity is not necessarily an issue, in the sense that they do not attempt to fix it, leaving QFT calculations untouched. Instead, gravity is modified such that the breakdown of the UV–IR decoupling does not lead to dramatic observable effects. Insofar as the main issue with the CC is this sensitivity coming from QFT calculations alone, then modifications of gravity only fix one of its symptoms that manifests in a gravitational context. Nonetheless, without self-tuning, UV sensitivity requires fine-tuning to describe observations; while with self-tuning, UV sensitivity is “only” a conceptual problem that does not lead to any observational effects.
Table 1. Here, we summarize our assessment of the approaches introduced in Section 4 in light of the theoretical and phenomenological requirements contained in our wishlist discussed in Section 4.1 (CCP stands for cosmological constant problem, DEP for dark energy problem, CHC for cosmic history constraint, and AC for astrophysical constraint). Whenever one of the CC-problems or data constraints has been shown to be solved, we use the 🗸; for requirements that potentially can be achieved within a given approach, we use (P) (for example, Fab-4 can reproduce individual cosmological epochs, but it remains to be seen if they can be stitched together [170]); while the ones that seem to be unsuccessful given the current literature are marked as X. Note that we evaluate whether a given approach tackles any of the requirements independently. For instance, we consider that “GR + QFT” can tackle the DEP classically, thus a 🗸, despite not being able to tackle the new-CCP, which deems the solution to DEP to beunstable quantum mechanically.
Table 1. Here, we summarize our assessment of the approaches introduced in Section 4 in light of the theoretical and phenomenological requirements contained in our wishlist discussed in Section 4.1 (CCP stands for cosmological constant problem, DEP for dark energy problem, CHC for cosmic history constraint, and AC for astrophysical constraint). Whenever one of the CC-problems or data constraints has been shown to be solved, we use the 🗸; for requirements that potentially can be achieved within a given approach, we use (P) (for example, Fab-4 can reproduce individual cosmological epochs, but it remains to be seen if they can be stitched together [170]); while the ones that seem to be unsuccessful given the current literature are marked as X. Note that we evaluate whether a given approach tackles any of the requirements independently. For instance, we consider that “GR + QFT” can tackle the DEP classically, thus a 🗸, despite not being able to tackle the new-CCP, which deems the solution to DEP to beunstable quantum mechanically.
Selected Modified Gravity Approaches
CC-ProblemsData Constraints
New-CCPClass-CCPDEPCHCAC
GR + QFTXX🗸🗸🗸
Global sequestering🗸🗸(P)🗸🗸
Local sequestering🗸🗸(P)🗸🗸
Nonlocal approachX(P)(P)🗸🗸
Unimodular gravityXX🗸🗸🗸
Linear massive gravity🗸🗸(P)🗸X
Nonlinear massive gravityX🗸(P)🗸(P)
Fab-4(P)🗸🗸(P)(P)
Well-tempered self-tuning(P)🗸🗸🗸(P)
SLEDXX(P)(P)🗸
Finally, we note that the table does not provide a fully comprehensive assessment of these models. Additional requirements should be demanded as a benchmark, such as having a ghost-free theory and a UV embedding. For instance, we briefly commented that the Horndeski class of models are all ghost-free, while global sequestering faces challenges to be embedded in a UV theory, but its local counterpart could more easily achieve this. Moreover, a theoretical model which is able to address more than one of the CCPs lends itself a stronger preference, given that it does not introduce a significant level of complexity.
In spite of the different proposals presented here, there is not a consensus in the community about what the solutions to the CCPs are. Among the reasons, some of these models recover most of the phenomenology already explained by GR, while others seem to be contrived. Alternative to seeking guidance from theory, new and distinct phenomenology is an essential ingredient for any prospective model. In particular, the next decade will see a host of novel experiments probing both low-energy cosmological and high-energy astrophysical scales. The largest galaxy surveys to date have either commenced or are coming online soon, for example, the ongoing dark energy survey24 [234], the upcoming Euclid mission25 [101], and the Vera C. Rubin Observatory’s Legacy Survey of Space and Time (VRO/LSST)26 [235]. These will precisely probe a yet-untested regime of structure formation where various effects due to modifications of gravity may become apparent. There is also the unprecedentedly large radio survey, the square kilometer array27 [236], set to come online in the next few years. Further, the emerging field of gravitational wave astronomy has immense potential to probe both cosmological and astrophysical scales. The planned LISA28 [237,238] experiment will precisely measure cosmological background effects, as well as effects related to perturbations both in the scalar and tensor sectors. These experiments will, at the very least, perform null tests of GR. Any tensions that emerge in the best-fit Λ CDM model between these upcoming and existing datasets will provide essential clues for theoretical progress.

Author Contributions

All authors contributed equally and have read and agreed to the published version of the manuscript.

Funding

H.B. was supported by the Fonds de recherche du Québec (PBEEE/303549) and partially by funds from NSERC. B.B. was supported by a UK Research and Innovation Stephen Hawking Fellowship (EP/W005654/1) and a Swiss National Science Foundation (SNSF) Professorship grant, Nos. 170547 & 202671. G.F. was supported by the Göran Gustafsson Foundation for Research in Natural Sciences and Medicine. Y.H. was supported by the grant 2019-04234 from the Swedish Research Council (Vetenskapsrådet). Nordita is partially supported by Nordforsk. A.L. acknowledges support by the Swedish Research Council (Vetenskapsrådet) through contract No. 638-2013-8993 and the Oskar Klein Centre for Cosmoparticle Physics. S.H. was supported by the Excellence Cluster ORIGINS which is funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy - EXC-2094 - 390783311.

Acknowledgments

We would like to thank Fawad Hassan and Marcus Högås for valuable discussions about massive gravity, and Hiranya Peiris for the support establishing this collaboration. We also would like to thank Elisa Ferreira for discussions and early-on contributions. For the purpose of open access, the authors applied a Creative Commons Attribution (CC BY) license to any Author Accepted Manuscript version arising from this submission.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CCPCosmological Constant Problem
UVUltraviolet
EFTEffective Field Theory
GRGeneral Relativity
EMTEnergy-Momentum Tensor
CDMCold Dark Matter
CCCosmological Constant
DEDark Energy
SMStandard Model of Particle Physics
QFTQuantum Field Theory
QCDQuantum Chromodynamics
DEPDark Energy Problem
IRInfrared
DGPDvali–Gabadadze–Porrati gravity
CMBCosmic Microwave Background
CHCCosmic History Constraint
ACAstrophysical Constraints
SLEDSupersymmetric Large Extra Dimension

Notes

1
Lenz [15] was the first to make such a comparison (see, e.g., [16,17,18] and references therein for historical notes).
2
Recently, there has been evidence of an early dark energy component at slightly higher energies (at the eV rather than the milli-eV scale) [45,46]. Similar to dark energy, there is the question of what early dark energy is made of.
3
Of course, this is not the only possibility to address the CCPs. In another scenario, the value of the CC at low energies could be connected to (and ultimately explained by) unknown UV physics. In such a scenario, explaining the smallness of the CC requires knowledge of the UV sector of the theory. For example, if we assume the SM of particle physics to be valid all the way up to the Planck scale, the condition of stability of the electroweak vacuum in the deep UV (as represented by the quartic Higgs coupling) is connected to the Higgs boson and top quark masses measured at much lower energies [53,54,55]. This is a specific form of UV–IR coupling and has been used to motivate alternative scenarios for addressing the CCP (see [56] as one example). We thank an anonymous referee for bringing that to our attention.
4
An often echoed concern is that self-tuning would also spell the end for inflation. However, this assumes that sources behaving similar to a cosmological constant (such as an inflaton field in slow-roll) are equally decoupled from gravity, which does not need to be true. Sequestering in Section 4.2 provides the cleanest example of a mechanism that only affects a true cosmological constant.
5
We could also look at a general conformal coupling, but as it turns out the exponential case will be sufficiently interesting.
6
For example, the Higgs vacuum is H vac = ( 0 , v H / 2 ) T = const , with vacuum expectation value v H = 246 GeV .
7
We note that this is not the only object that transforms as a density. Instead, we can use ϵ μ ν α β A μ ν α β as the integration measure, where ϵ μ ν α β is the Levi–Civita symbol and A μ ν α β a four-form field. This is the loophole exploited by some of the approaches discussed in Section 4.2.
8
Spin-1 fields linearly coupled to a conserved source cannot contribute to A .
9
These formal relations can be defined rigorously in momentum space.
10
As an aside, screening mechanisms discussed in Section 3.2 fulfill this condition based on a violation of the weak coupling assumption (iii).
11
To derive this condition, Equation (19) needs to be slightly generalized to also account for nonlocalized sources such as vacuum energy.
12
Cosmic structure formation measurements are also constraining, but to a much lesser degree [36].
13
This theory now still requires the CC for cosmic acceleration, having been strongly constrained by data (see, for example, [88]). In the context of self-tuning, see [89].
14
In particular, the literature moves beyond proposals simply trying to tackle the CCPs, as modifications of gravity have also been extensively considered as an explanation for the late-time accelerated expansion of the Universe and an alternative to dark matter [34].
15
In contrast to unimodular gravity, for example, see Section 4.2.4.
16
There is a caveat here concerning gravitational loop corrections, which would introduce κ -dependent terms to Equation (45) that are not sequestered. This problem can be avoided if topological curvature invariants are also considered, as discussed in [115,116].
17
We denote the rank of a p-form by the subscript ( p ) .
18
On a more fundamental level, the occurrence of F ( 4 ) can be linked to the presence of membranes, but it has been argued that such an explanation might compromise the mechanism [108].
19
Superposition of small-mass massive spin-2 states.
20
These models typically generalize the Randall–Sundrum model [200], which itself cannot address the CCP due to an immediate tuning between the brane tension and the bulk cosmological constant.
21
For higher value of the brane tension, the geometry on the brane changes to de Sitter [206].
22
This model builds on the Nishino–Sezgin chiral gauged supergravity, which admits rugby ball solutions [220,221,222].
23
The last approach also addressed concerns about the general applicability of the distributional approach away from the scale-invariant case [229].
24
25
26
27
28

References

  1. Einstein, A. Kosmologische Betrachtungen zur allgemeinen Relativitätstheorie. In Sitzungsberichte der Königlich Preußischen Akademie der Wissenschaften, Berlin; Vieweg+Teubner Verlag: Wiesbaden, Germany, 1917; pp. 142–152. [Google Scholar]
  2. O’Raifeartaigh, C.; O’Keeffe, M.; Nahm, W.; Mitton, S. Einstein’s 1917 static model of the universe: A centennial review. Eur. Phys. J. H 2017, 42, 431–474. [Google Scholar] [CrossRef] [Green Version]
  3. Aghanim, N.; et al.; [Planck Collaboration] Planck 2018 results. VI. Cosmological parameters. Astron. Astrophys. 2020, 641, A6, [Erratum in Astron. Astrophys. 2021, 652, C4]. [Google Scholar] [CrossRef] [Green Version]
  4. Bautista, J.E.; Paviot, R.; Vargas Magaña, M.; de la Torre, S.; Fromenteau, S.; Gil-Marín, H.; Ross, A.J.; Burtin, E.; Dawson, K.S.; Hou, J.; et al. The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Measurement of the BAO and growth rate of structure of the luminous red galaxy sample from the anisotropic correlation function between redshifts 0.6 and 1. Mon. Not. Roy. Astron. Soc. 2020, 500, 736–762. [Google Scholar] [CrossRef]
  5. Alam, S.; Ata, M.; Bailey, S.; Beutler, F.; Bizyaev, D.; Blazek, J.A.; Bolton, A.S.; Brownstein, J.R.; Burden, A.; Chuang, C.H.; et al. The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: Cosmological analysis of the DR12 galaxy sample. Mon. Not. Roy. Astron. Soc. 2017, 470, 2617–2652. [Google Scholar] [CrossRef] [Green Version]
  6. Hildebrandt, H.; Viola, M.; Heymans, C.; Joudaki, S.; Kuijken, K.; Blake, C.; Erben, T.; Joachimi, B.; Klaes, D.; Miller, L.; et al. KiDS-450: Cosmological parameter constraints from tomographic weak gravitational lensing. Mon. Not. Roy. Astron. Soc. 2017, 465, 1454. [Google Scholar] [CrossRef] [Green Version]
  7. Doux, C.; Jain, B.; Zeurcher, D.; Lee, J.; Fang, X.; Rosenfeld, R.; Amon, A.; Camacho, H.; Choi, A.; Secco, L.F.; et al. Dark Energy Survey Year 3 results: Cosmological constraints from the analysis of cosmic shear in harmonic space. Mon. Not. Roy. Astron. Soc. 2022, 515, 1942–1972. [Google Scholar] [CrossRef]
  8. Brout, D.; Scolnic, D.; Popovic, B.; Riess, A.G.; Carr, A.; Zuntz, J.; Kessler, R.; Davis, T.M.; Hinton, S.; Jones, D.; et al. The Pantheon+ Analysis: Cosmological Constraints. Astrophys. J. 2022, 938, 110. [Google Scholar] [CrossRef]
  9. Costanzi, M.; et al.; [DES and SPT Collaborations] Cosmological constraints from DES Y1 cluster abundances and SPT multiwavelength data. Phys. Rev. D 2021, 103, 043522. [Google Scholar] [CrossRef]
  10. de Haan, T.; Benson, B.A.; Bleem, L.E.; Allen, S.W.; Applegate, D.E.; Ashby, M.L.N.; Bautz, M.; Bayliss, M.; Bocquet, S.; Brodwin, M.; et al. Cosmological Constraints from Galaxy Clusters in the 2500 square-degree SPT-SZ Survey. Astrophys. J. 2016, 832, 95. [Google Scholar] [CrossRef]
  11. Huterer, D.; Shafer, D.L. Dark energy two decades after: Observables, probes, consistency tests. Rept. Prog. Phys. 2018, 81, 016901. [Google Scholar] [CrossRef]
  12. Birrell, N.D.; Davies, P.C.W. Quantum Fields in Curved Space; Cambridge Monographs on Mathematical Physics; Cambridge Univ. Press: Cambridge, UK, 1984. [Google Scholar] [CrossRef]
  13. Nernst, W. Über einen Versuch, von quantentheoretischen Betrachtungen zur Annahme stetiger Energieänderungen zurückzukehren. Verhandlungen der Deutschen Physikalischen Gesellschaft 1916, 18, 83–116. [Google Scholar]
  14. Enz, C.P.; Thellung, A. Nullpunktsenergie und Anordnung nicht vertauschbarer Faktoren im Hamiltonoperator. Helv. Phys. Acta 1960, 33, 839. [Google Scholar]
  15. Lenz, W. Das Gleichgewicht von Materie und Strahlung in Einsteins geschlossener Welt. Phys. Zs. 1926, 27, 642–645. [Google Scholar]
  16. Peruzzi, G.; Realdi, M. The quest for the size of the universe in early relativistic cosmology (1917–1930). Arch. Hist. Exact Sci. 2011, 65, 659. [Google Scholar] [CrossRef] [Green Version]
  17. Kragh, H. Walther Nernst: Grandfather of dark energy? Astron. Geophys. 2012, 53, 24–26. [Google Scholar] [CrossRef] [Green Version]
  18. Kragh, H.S.; Overduin, J.M. The Weight of the Vacuum: A Scientific History of Dark Energy; Springer: Berlin, Germany, 2014. [Google Scholar]
  19. Zeldovich, Y.B. Cosmological Constant and Elementary Particles. JETP Lett. 1967, 6, 316. [Google Scholar] [CrossRef]
  20. Zel’dovich, Y.B.; Krasinski, A.; Zeldovich, Y.B. The Cosmological constant and the theory of elementary particles. Sov. Phys. Usp. 1968, 11, 381–393. [Google Scholar] [CrossRef]
  21. ’t Hooft, G.; Veltman, M.J.G. One loop divergencies in the theory of gravitation. Ann. Inst. H. Poincare Phys. Theor. A 1974, 20, 69–94. [Google Scholar]
  22. Goroff, M.H.; Sagnotti, A. The Ultraviolet Behavior of Einstein Gravity. Nucl. Phys. B 1986, 266, 709–736. [Google Scholar] [CrossRef]
  23. Weinberg, S. The Cosmological Constant Problem. Rev. Mod. Phys. 1989, 61, 1–23. [Google Scholar] [CrossRef]
  24. Martin, J. Everything You Always Wanted to Know About the Cosmological Constant Problem (But Were Afraid to Ask). Comptes Rendus Phys. 2012, 13, 566–665. [Google Scholar] [CrossRef] [Green Version]
  25. Copeland, E.J.; Sami, M.; Tsujikawa, S. Dynamics of dark energy. Int. J. Mod. Phys. D 2006, 15, 1753–1936. [Google Scholar] [CrossRef] [Green Version]
  26. Barrow, J.D.; Tipler, F.J. The Anthropic Cosmological Principle; Oxford U. Pr.: Oxford, UK, 1988. [Google Scholar]
  27. Gribbin, J.; Rees, M. Cosmic Coincidences: Dark Matter, Mankind, and Anthropic Cosmology; Black Swan: London, UK, 1989. [Google Scholar]
  28. Weinberg, S. Anthropic Bound on the Cosmological Constant. Phys. Rev. Lett. 1987, 59, 2607. [Google Scholar] [CrossRef] [PubMed]
  29. Efstathiou, G. An anthropic argument for a cosmological constant. Mon. Not. R. Astron. Soc. 1995, 274, L73–L76. [Google Scholar] [CrossRef] [Green Version]
  30. Martel, H.; Shapiro, P.R.; Weinberg, S. Likely values of the cosmological constant. Astrophys. J. 1998, 492, 29. [Google Scholar] [CrossRef] [Green Version]
  31. Garriga, J.; Livio, M.; Vilenkin, A. The Cosmological constant and the time of its dominance. Phys. Rev. D 2000, 61, 023503. [Google Scholar] [CrossRef] [Green Version]
  32. Peacock, J.A. Testing anthropic predictions for Lambda and the CMB temperature. Mon. Not. Roy. Astron. Soc. 2007, 379, 1067–1074. [Google Scholar] [CrossRef] [Green Version]
  33. Sobral Blanco, D.; Lombriser, L. Local self-tuning mechanism for the cosmological constant. Phys. Rev. D 2020, 102, 043506. [Google Scholar] [CrossRef]
  34. Clifton, T.; Ferreira, P.G.; Padilla, A.; Skordis, C. Modified Gravity and Cosmology. Phys. Rept. 2012, 513, 1–189. [Google Scholar] [CrossRef] [Green Version]
  35. Hinterbichler, K. Theoretical Aspects of Massive Gravity. Rev. Mod. Phys. 2012, 84, 671–710. [Google Scholar] [CrossRef] [Green Version]
  36. Koyama, K. Cosmological Tests of Modified Gravity. Rept. Prog. Phys. 2016, 79, 046902. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. de Rham, C. Massive Gravity. Living Rev. Rel. 2014, 17, 7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Nobbenhuis, S. Categorizing different approaches to the cosmological constant problem. Found. Phys. 2006, 36, 613–680. [Google Scholar] [CrossRef]
  39. Polchinski, J. The Cosmological Constant and the String Landscape. In Proceedings of the 23rd Solvay Conference in Physics: The Quantum Structure of Space and Time, Brussels, Belgium, 1–3 December 2005; pp. 216–236. [Google Scholar]
  40. Bousso, R. TASI Lectures on the Cosmological Constant. Gen. Rel. Grav. 2008, 40, 607–637. [Google Scholar] [CrossRef] [Green Version]
  41. Padilla, A. Lectures on the Cosmological Constant Problem. arXiv 2015, arXiv:1502.05296. [Google Scholar]
  42. Nojiri, S.; Odintsov, S.D.; Oikonomou, V.K. Modified Gravity Theories on a Nutshell: Inflation, Bounce and Late-time Evolution. Phys. Rept. 2017, 692, 1–104. [Google Scholar] [CrossRef] [Green Version]
  43. Riess, A.G.; Filippenko, A.V.; Challis, P.; Clocchiatti, A.; Diercks, A.; Garnavich, P.M.; Gilliland, R.L.; Hogan, C.J.; Jha, S.; Kirshner, R.P.; et al. Observational evidence from supernovae for an accelerating universe and a cosmological constant. Astron. J. 1998, 116, 1009–1038. [Google Scholar] [CrossRef] [Green Version]
  44. Perlmutter, S.; Aldering, G.; Goldhaber, G.; Knop, R.A.; Nugent, P.; Castro, P.G.; Deustua, S.; Fabbro, S.; Goobar, A.; Groom, D.E.; et al. Measurements of W and L from 42 high redshift supernovae. Astrophys. J. 1999, 517, 565–586. [Google Scholar] [CrossRef]
  45. Poulin, V.; Smith, T.L.; Karwal, T.; Kamionkowski, M. Early Dark Energy Can Resolve The Hubble Tension. Phys. Rev. Lett. 2019, 122, 221301. [Google Scholar] [CrossRef] [Green Version]
  46. Kamionkowski, M.; Riess, A.G. The Hubble Tension and Early Dark Energy. arXiv 2022, arXiv:2211.04492. [Google Scholar]
  47. Wells, J.D. Lectures on Higgs Boson Physics in the Standard Model and Beyond. In Proceedings of the 38th British Universities Summer School in Theoretical Elementary Particle Physics, Cambridge, UK, 25 August–6 September 2008. [Google Scholar]
  48. Koren, S. The Hierarchy Problem: From the Fundamentals to the Frontiers. Ph.D. Thesis, University of California, Santa Barbara, CA, USA, 2020. [Google Scholar]
  49. Wetterich, C. Cosmology and the Fate of Dilatation Symmetry. Nucl. Phys. B 1988, 302, 668–696. [Google Scholar] [CrossRef] [Green Version]
  50. Peebles, P.J.E.; Ratra, B. Cosmology with a Time Variable Cosmological Constant. Astrophys. J. Lett. 1988, 325, L17. [Google Scholar] [CrossRef]
  51. Tsujikawa, S. Quintessence: A Review. Class. Quant. Grav. 2013, 30, 214003. [Google Scholar] [CrossRef] [Green Version]
  52. Burgess, C.P. The Cosmological Constant Problem: Why it’s hard to get Dark Energy from Micro-physics. In Proceedings of the 100e Ecole d’Ete de Physique: Post-Planck Cosmology, Les Houches, France, 8 July–2 August 2013. [Google Scholar] [CrossRef] [Green Version]
  53. Bednyakov, A.V.; Kniehl, B.A.; Pikelner, A.F.; Veretin, O.L. Stability of the Electroweak Vacuum: Gauge Independence and Advanced Precision. Phys. Rev. Lett. 2015, 115, 201802. [Google Scholar] [CrossRef] [PubMed]
  54. Degrassi, G.; Di Vita, S.; Elias-Miro, J.; Espinosa, J.R.; Giudice, G.F.; Isidori, G.; Strumia, A. Higgs mass and vacuum stability in the Standard Model at NNLO. J. High Energy Phys. 2012, 8, 98. [Google Scholar] [CrossRef] [Green Version]
  55. Jegerlehner, F. The Standard model as a low-energy effective theory: What is triggering the Higgs mechanism? Acta Phys. Polon. B 2014, 45, 1167. [Google Scholar] [CrossRef] [Green Version]
  56. Jegerlehner, F. The Hierarchy Problem and the Cosmological Constant Problem Revisited—A new view on the SM of particle physics. Found. Phys. 2019, 49, 915–971. [Google Scholar] [CrossRef] [Green Version]
  57. Niedermann, F.; Padilla, A. Gravitational Mechanisms to Self-Tune the Cosmological Constant: Obstructions and Ways Forward. Phys. Rev. Lett. 2017, 119, 251306. [Google Scholar] [CrossRef] [Green Version]
  58. Fierz, M.; Pauli, W. On relativistic wave equations for particles of arbitrary spin in an electromagnetic field. Proc. Roy. Soc. Lond. A 1939, 173, 211–232. [Google Scholar] [CrossRef]
  59. Will, C.M. The Confrontation between General Relativity and Experiment. Living Rev. Rel. 2014, 17, 4. [Google Scholar] [CrossRef] [Green Version]
  60. Belgacem, E.; Finke, A.; Frassino, A.; Maggiore, M. Testing nonlocal gravity with Lunar Laser Ranging. J. Cosmol. Astropart. Phys. 2019, 2, 35. [Google Scholar] [CrossRef] [Green Version]
  61. Burrage, C.; Kaloper, N.; Padilla, A. Strong Coupling and Bounds on the Spin-2 Mass in Massive Gravity. Phys. Rev. Lett. 2013, 111, 021802. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. de Rham, C.; Ribeiro, R.H. Riding on irrelevant operators. J. Cosmol. Astropart. Phys. 2014, 11, 16. [Google Scholar] [CrossRef]
  63. Kaloper, N.; Padilla, A.; Saffin, P.; Stefanyszyn, D. Unitarity and the Vainshtein Mechanism. Phys. Rev. D 2015, 91, 045017. [Google Scholar] [CrossRef]
  64. Adams, A.; McGreevy, J.; Silverstein, E. Decapitating tadpoles. arXiv 2002, arXiv:hep-th/0209226. [Google Scholar]
  65. Bertotti, B.; Iess, L.; Tortora, P. A test of general relativity using radio links with the Cassini spacecraft. Nature 2003, 425, 374–376. [Google Scholar] [CrossRef]
  66. Williams, J.G.; Newhall, X.X.; Dickey, J.O. Relativity parameters determined from lunar laser ranging. Phys. Rev. D 1996, 53, 6730–6739. [Google Scholar] [CrossRef]
  67. Anderson, J.D.; Gross, M.; Nordtvedt, K.L.; Turyshev, S.G. The Solar test of the equivalence principle. Astrophys. J. 1996, 459, 365–370. [Google Scholar] [CrossRef] [Green Version]
  68. Dickey, J.O.; Bender, P.; Faller, J.; Newhall, X.; Ricklefs, R.; Ries, J.; Shelus, P.; Veillet, C.; Whipple, A.; Wiant, J.; et al. Lunar Laser Ranging: A Continuing Legacy of the Apollo Program. Science 1994, 265, 482–490. [Google Scholar] [CrossRef] [Green Version]
  69. Talmadge, C.; Berthias, J.P.; Hellings, R.W.; Standish, E.M. Model Independent Constraints on Possible Modifications of Newtonian Gravity. Phys. Rev. Lett. 1988, 61, 1159–1162. [Google Scholar] [CrossRef]
  70. Dvali, G.R.; Gabadadze, G.; Porrati, M. 4-D gravity on a brane in 5-D Minkowski space. Phys. Lett. B 2000, 485, 208–214. [Google Scholar] [CrossRef] [Green Version]
  71. Dvali, G.R.; Gabadadze, G. Gravity on a brane in infinite volume extra space. Phys. Rev. D 2001, 63, 065007. [Google Scholar] [CrossRef] [Green Version]
  72. Luty, M.A.; Porrati, M.; Rattazzi, R. Strong interactions and stability in the DGP model. J. High Energy Phys. 2003, 9, 029. [Google Scholar] [CrossRef] [Green Version]
  73. Horndeski, G.W. Second-order scalar-tensor field equations in a four-dimensional space. Int. J. Theor. Phys. 1974, 10, 363–384. [Google Scholar] [CrossRef]
  74. Woodard, R.P. Avoiding dark energy with 1/r modifications of gravity. Lect. Notes Phys. 2007, 720, 403–433. [Google Scholar] [CrossRef]
  75. Motohashi, H.; Suyama, T. Third order equations of motion and the Ostrogradsky instability. Phys. Rev. D 2015, 91, 085009. [Google Scholar] [CrossRef] [Green Version]
  76. Langlois, D.; Noui, K. Degenerate higher derivative theories beyond Horndeski: Evading the Ostrogradski instability. J. Cosmol. Astropart. Phys. 2016, 2, 034. [Google Scholar] [CrossRef] [Green Version]
  77. Langlois, D.; Noui, K. Hamiltonian analysis of higher derivative scalar-tensor theories. J. Cosmol. Astropart. Phys. 2016, 7, 016. [Google Scholar] [CrossRef] [Green Version]
  78. Crisostomi, M.; Koyama, K.; Tasinato, G. Extended Scalar-Tensor Theories of Gravity. J. Cosmol. Astropart. Phys. 2016, 4, 044. [Google Scholar] [CrossRef]
  79. Nicolis, A.; Rattazzi, R.; Trincherini, E. The Galileon as a local modification of gravity. Phys. Rev. D 2009, 79, 064036. [Google Scholar] [CrossRef] [Green Version]
  80. Jordan, P. The present state of Dirac’s cosmological hypothesis. Z. Phys 1959, 157, 112. [Google Scholar] [CrossRef]
  81. Brans, C.; Dicke, R.H. Mach’s principle and a relativistic theory of gravitation. Phys. Rev. 1961, 124, 925. [Google Scholar] [CrossRef]
  82. Joyce, A.; Lombriser, L.; Schmidt, F. Dark Energy Versus Modified Gravity. Ann. Rev. Nucl. Part. Sci. 2016, 66, 95–122. [Google Scholar] [CrossRef] [Green Version]
  83. Sakstein, J. Astrophysical Tests of Modified Gravity. Ph.D. Thesis, Cambridge University, DAMTP, Cambridge, UK, 2014. [Google Scholar] [CrossRef]
  84. Cataneo, M.; Rapetti, D. Tests of gravity with galaxy clusters. Int. J. Mod. Phys. D 2018, 27, 1848006. [Google Scholar] [CrossRef]
  85. Carroll, S.M.; Sawicki, I.; Silvestri, A.; Trodden, M. Modified-Source Gravity and Cosmological Structure Formation. New J. Phys. 2006, 8, 323. [Google Scholar] [CrossRef] [Green Version]
  86. Pogosian, L.; Silvestri, A. What can cosmology tell us about gravity? Constraining Horndeski gravity with S and m. Phys. Rev. D 2016, 94, 104014. [Google Scholar] [CrossRef] [Green Version]
  87. Peirone, S.; Koyama, K.; Pogosian, L.; Raveri, M.; Silvestri, A. Large-scale structure phenomenology of viable Horndeski theories. Phys. Rev. D 2018, 97, 043519. [Google Scholar] [CrossRef] [Green Version]
  88. Burrage, C.; Sakstein, J. Tests of Chameleon Gravity. Living Rev. Rel. 2018, 21, 1. [Google Scholar] [CrossRef] [Green Version]
  89. Bauer, F.; Sola, J.; Stefancic, H. The Relaxed Universe: Towards solving the cosmological constant problem dynamically from an effective action functional of gravity. Phys. Lett. B 2010, 688, 269–272. [Google Scholar] [CrossRef] [Green Version]
  90. Hu, W.; Sawicki, I. Models of f(R) cosmic acceleration that evade solar system tests. Phys. Rev. D 2007, 76, 064004. [Google Scholar] [CrossRef] [Green Version]
  91. Khoury, J.; Weltman, A. Chameleon cosmology. Phys. Rev. D 2004, 69, 044026. [Google Scholar] [CrossRef]
  92. Navarro, I.; Van Acoleyen, K. f(R) actions, cosmic acceleration and local tests of gravity. J. Cosmol. Astropart. Phys. 2007, 2, 022. [Google Scholar] [CrossRef] [Green Version]
  93. Faulkner, T.; Tegmark, M.; Bunn, E.F.; Mao, Y. Constraining f(R) Gravity as a Scalar Tensor Theory. Phys. Rev. D 2007, 76, 063505. [Google Scholar] [CrossRef] [Green Version]
  94. Brax, P.; Casas, S.; Desmond, H.; Elder, B. Testing Screened Modified Gravity. Universe 2021, 8, 11. [Google Scholar] [CrossRef]
  95. Nicolis, A.; Rattazzi, R. Classical and quantum consistency of the DGP model. J. High Energy Phys. 2004, 6, 059. [Google Scholar] [CrossRef]
  96. Koyama, K.; Silva, F.P. Non-linear interactions in a cosmological background in the DGP braneworld. Phys. Rev. D 2007, 75, 084040. [Google Scholar] [CrossRef] [Green Version]
  97. Barreira, A.; Li, B.; Hellwing, W.A.; Baugh, C.M.; Pascoli, S. Nonlinear structure formation in the Cubic Galileon gravity model. J. Cosmol. Astropart. Phys. 2013, 10, 027. [Google Scholar] [CrossRef] [Green Version]
  98. Vainshtein, A.I. To the problem of nonvanishing gravitation mass. Phys. Lett. B 1972, 39, 393–394. [Google Scholar] [CrossRef]
  99. Lombriser, L.; Li, B.; Koyama, K.; Zhao, G.B. Modeling halo mass functions in chameleon f(R) gravity. Phys. Rev. D 2013, 87, 123511. [Google Scholar] [CrossRef] [Green Version]
  100. Schmidt, F.; Hu, W.; Lima, M. Spherical Collapse and the Halo Model in Braneworld Gravity. Phys. Rev. D 2010, 81, 063005. [Google Scholar] [CrossRef] [Green Version]
  101. Amendola, L.; Appleby, S.; Avgoustidis, A.; Bacon, D.; Baker, T.; Baldi, M.; Bartolo, N.; Blanchard, A.; Bonvin, C.; Borgani, S.; et al. Cosmology and fundamental physics with the Euclid satellite. Living Rev. Rel. 2018, 21, 2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  102. Lee, S.; Huff, E.M.; Choi, A.; Elvin-Poole, J.; Hirata, C.; Honscheid, K.; MacCrann, N.; Ross, A.J.; Troxel, M.A.; Eifler, T.F.; et al. Probing gravity with the DES-CMASS sample and BOSS spectroscopy. Mon. Not. Roy. Astron. Soc. 2021, 509, 4982–4996. [Google Scholar] [CrossRef]
  103. Alam, S.; Arnold, C.; Aviles, A.; Bean, R.; Cai, Y.C.; Cautun, M.; Cervantes-Cota, J.L.; Cuesta-Lazaro, C.; Devi, N.C.; Eggemeier, A.; et al. Towards testing the theory of gravity with DESI: Summary statistics, model predictions and future simulation requirements. J. Cosmol. Astropart. Phys. 2021, 11, 050. [Google Scholar] [CrossRef]
  104. Brax, P.; Valageas, P. Quantum field theory of K-mouflage. Phys. Rev. D 2016, 94, 043529. [Google Scholar] [CrossRef]
  105. Joyce, A.; Jain, B.; Khoury, J.; Trodden, M. Beyond the Cosmological Standard Model. Phys. Rept. 2015, 568, 1–98. [Google Scholar] [CrossRef] [Green Version]
  106. Kaloper, N.; Padilla, A. Sequestering the Standard Model Vacuum Energy. Phys. Rev. Lett. 2014, 112, 091304. [Google Scholar] [CrossRef] [Green Version]
  107. Kaloper, N.; Padilla, A.; Stefanyszyn, D.; Zahariade, G. Manifestly Local Theory of Vacuum Energy Sequestering. Phys. Rev. Lett. 2016, 116, 051302. [Google Scholar] [CrossRef] [Green Version]
  108. D’Amico, G.; Kaloper, N.; Padilla, A.; Stefanyszyn, D.; Westphal, A.; Zahariade, G. An étude on global vacuum energy sequester. J. High Energy Phys. 2017, 9, 074. [Google Scholar] [CrossRef] [Green Version]
  109. Einstein, A.; Mayer, W. Albert Einstein: Akademie-Vorträge: Sitzungsberichte der Preußischen Akademie der Wissenschaften 1914–1932; Wiley: Weinheim, Germany, 2006; pp. 375–392. [Google Scholar] [CrossRef]
  110. Ng, Y.J.; van Dam, H. Unimodular Theory of Gravity and the Cosmological Constant. J. Math. Phys. 1991, 32, 1337–1340. [Google Scholar] [CrossRef]
  111. Carroll, S.M.; Remmen, G.N. A Nonlocal Approach to the Cosmological Constant Problem. Phys. Rev. D 2017, 95, 123504. [Google Scholar] [CrossRef] [Green Version]
  112. Lombriser, L. On the cosmological constant problem. Phys. Lett. B 2019, 797, 134804. [Google Scholar] [CrossRef]
  113. Kaloper, N.; Westphal, A. A Quantum-Mechanical Mechanism for Reducing the Cosmological Constant. arXiv 2022, arXiv:2204.13124. [Google Scholar] [CrossRef]
  114. Kaloper, N. Pancosmic Relativity and Nature’s Hierarchies. arXiv 2022, arXiv:2202.08860. [Google Scholar]
  115. Kaloper, N.; Padilla, A. Vacuum Energy Sequestering and Graviton Loops. Phys. Rev. Lett. 2017, 118, 061303. [Google Scholar] [CrossRef] [PubMed]
  116. El-Menoufi, B.K.; Nagy, S.; Niedermann, F.; Padilla, A. Quantum corrections to vacuum energy sequestering (with monodromy). Class. Quant. Grav. 2019, 36, 215014. [Google Scholar] [CrossRef] [Green Version]
  117. Kaloper, N.; Padilla, A.; Stefanyszyn, D. Sequestering effects on and of vacuum decay. Phys. Rev. D 2016, 94, 025022. [Google Scholar] [CrossRef] [Green Version]
  118. Kaloper, N.; Padilla, A. Vacuum Energy Sequestering: The Framework and Its Cosmological Consequences. Phys. Rev. D 2014, 90, 084023, Addendum in Phys. Rev. D 2014, 90, 109901. [Google Scholar] [CrossRef] [Green Version]
  119. Kaloper, N.; Padilla, A. Sequestration of Vacuum Energy and the End of the Universe. Phys. Rev. Lett. 2015, 114, 101302. [Google Scholar] [CrossRef] [Green Version]
  120. Kaloper, N. Irrational Monodromies of Vacuum Energy. J. High Energy Phys. 2019, 11, 106. [Google Scholar] [CrossRef] [Green Version]
  121. Padilla, A. Monodromy inflation and an emergent mechanism for stabilising the cosmological constant. J. High Energy Phys. 2019, 1, 175. [Google Scholar] [CrossRef] [Green Version]
  122. Oda, I. Manifestly Local Formulation of Nonlocal Approach to the Cosmological Constant Problem. Phys. Rev. D 2017, 95, 104020. [Google Scholar] [CrossRef] [Green Version]
  123. Smolin, L. The Quantization of unimodular gravity and the cosmological constant problems. Phys. Rev. D 2009, 80, 084003. [Google Scholar] [CrossRef] [Green Version]
  124. Padilla, A.; Saltas, I.D. A note on classical and quantum unimodular gravity. Eur. Phys. J. C 2015, 75, 561. [Google Scholar] [CrossRef] [Green Version]
  125. Percacci, R. Unimodular quantum gravity and the cosmological constant. Found. Phys. 2018, 48, 1364–1379. [Google Scholar] [CrossRef]
  126. Eichhorn, A. On unimodular quantum gravity. Class. Quant. Grav. 2013, 30, 115016. [Google Scholar] [CrossRef]
  127. Kuchar, K.V. Does an unspecified cosmological constant solve the problem of time in quantum gravity? Phys. Rev. D 1991, 43, 3332–3344. [Google Scholar] [CrossRef]
  128. Henneaux, M.; Teitelboim, C. The Cosmological Constant and General Covariance. Phys. Lett. B 1989, 222, 195–199. [Google Scholar] [CrossRef]
  129. Fiol, B.; Garriga, J. Semiclassical Unimodular Gravity. J. Cosmol. Astropart. Phys. 2010, 8, 015. [Google Scholar] [CrossRef] [Green Version]
  130. Nojiri, S.; Odintsov, S.D.; Oikonomou, V.K. Unimodular F(R) Gravity. J. Cosmol. Astropart. Phys. 2016, 5, 046. [Google Scholar] [CrossRef] [Green Version]
  131. Schmidt-May, A.; von Strauss, M. Recent developments in bimetric theory. J. Phys. A 2016, 49, 183001. [Google Scholar] [CrossRef] [Green Version]
  132. de Rham, C.; Deskins, J.T.; Tolley, A.J.; Zhou, S.Y. Graviton Mass Bounds. Rev. Mod. Phys. 2017, 89, 025004. [Google Scholar] [CrossRef] [Green Version]
  133. Heisenberg, L. A systematic approach to generalisations of General Relativity and their cosmological implications. Phys. Rept. 2019, 796, 1–113. [Google Scholar] [CrossRef] [Green Version]
  134. Arkani-Hamed, N.; Dimopoulos, S.; Dvali, G.; Gabadadze, G. Nonlocal modification of gravity and the cosmological constant problem. arXiv 2002, arXiv:hep-th/0209227. [Google Scholar]
  135. Dvali, G.; Hofmann, S.; Khoury, J. Degravitation of the cosmological constant and graviton width. Phys. Rev. D 2007, 76, 084006. [Google Scholar] [CrossRef] [Green Version]
  136. de Rham, C.; Hofmann, S.; Khoury, J.; Tolley, A.J. Cascading Gravity and Degravitation. J. Cosmol. Astropart. Phys. 2008, 2, 011. [Google Scholar] [CrossRef]
  137. van Dam, H.; Veltman, M.J.G. Massive and massless Yang-Mills and gravitational fields. Nucl. Phys. B 1970, 22, 397–411. [Google Scholar] [CrossRef] [Green Version]
  138. Zakharov, V.I. Linearized gravitation theory and the graviton mass. JETP Lett. 1970, 12, 312. [Google Scholar]
  139. Stueckelberg, E.C.G. Theory of the radiation of photons of small arbitrary mass. Helv. Phys. Acta 1957, 30, 209–215. [Google Scholar]
  140. Arkani-Hamed, N.; Georgi, H.; Schwartz, M.D. Effective field theory for massive gravitons and gravity in theory space. Annals Phys. 2003, 305, 96–118. [Google Scholar] [CrossRef] [Green Version]
  141. Boulware, D.G.; Deser, S. Can gravitation have a finite range? Phys. Rev. D 1972, 6, 3368–3382. [Google Scholar] [CrossRef] [Green Version]
  142. de Rham, C.; Gabadadze, G.; Tolley, A.J. Resummation of Massive Gravity. Phys. Rev. Lett. 2011, 106, 231101. [Google Scholar] [CrossRef] [Green Version]
  143. de Rham, C.; Gabadadze, G. Generalization of the Fierz-Pauli Action. Phys. Rev. D 2010, 82, 044020. [Google Scholar] [CrossRef] [Green Version]
  144. Hassan, S.F.; Rosen, R.A. On Non-Linear Actions for Massive Gravity. J. High Energy Phys. 2011, 7, 009. [Google Scholar] [CrossRef] [Green Version]
  145. Bellazzini, B.; Riva, F.; Serra, J.; Sgarlata, F. Beyond Positivity Bounds and the Fate of Massive Gravity. Phys. Rev. Lett. 2018, 120, 161101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  146. de Rham, C.; Melville, S.; Tolley, A.J.; Zhou, S.Y. Positivity Bounds for Massive Spin-1 and Spin-2 Fields. J. High Energy Phys. 2019, 3, 182. [Google Scholar] [CrossRef] [Green Version]
  147. LIGO Scientific Collaboration and Virgo Collaboration; Fermi GBM; INTEGRAL; IceCube Collaboration; AstroSat Cadmium Zinc Telluride Imager Team; IPN Collaboration; The Insight-HXMT Collaboration; ANTARES Collaboration; The Swift Collaboration; AGILE Team; et al. Multi-messenger Observations of a Binary Neutron Star Merger. Astrophys. J. Lett. 2017, 848, L12. [Google Scholar] [CrossRef]
  148. Baker, T.; Bellini, E.; Ferreira, P.G.; Lagos, M.; Noller, J.; Sawicki, I. Strong constraints on cosmological gravity from GW170817 and GRB 170817A. Phys. Rev. Lett. 2017, 119, 251301. [Google Scholar] [CrossRef]
  149. Bernus, L.; Minazzoli, O.; Fienga, A.; Gastineau, M.; Laskar, J.; Deram, P. Constraining the mass of the graviton with the planetary ephemeris INPOP. Phys. Rev. Lett. 2019, 123, 161103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  150. Khosravi, N.; Rahmanpour, N.; Sepangi, H.R.; Shahidi, S. Multi-Metric Gravity via Massive Gravity. Phys. Rev. D 2012, 85, 024049. [Google Scholar] [CrossRef] [Green Version]
  151. Platscher, M.; Smirnov, J. Degravitation of the Cosmological Constant in Bigravity. J. Cosmol. Astropart. Phys. 2017, 3, 051. [Google Scholar] [CrossRef] [Green Version]
  152. Högås, M.; Torsello, F.; Mörtsell, E. On the stability of bimetric structure formation. J. Cosmol. Astropart. Phys. 2020, 4, 046. [Google Scholar] [CrossRef]
  153. Högås, M.; Mörtsell, E. Constraints on bimetric gravity from Big Bang nucleosynthesis. J. Cosmol. Astropart. Phys. 2021, 11, 1. [Google Scholar] [CrossRef]
  154. Caravano, A.; Lüben, M.; Weller, J. Combining cosmological and local bounds on bimetric theory. J. Cosmol. Astropart. Phys. 2021, 9, 035. [Google Scholar] [CrossRef]
  155. Högås, M.; Mörtsell, E. Constraints on bimetric gravity. Part II. Observational constraints. J. Cosmol. Astropart. Phys. 2021, 5, 002. [Google Scholar] [CrossRef]
  156. Gleyzes, J.; Langlois, D.; Piazza, F.; Vernizzi, F. Healthy theories beyond Horndeski. Phys. Rev. Lett. 2015, 114, 211101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  157. Lin, C.; Mukohyama, S.; Namba, R.; Saitou, R. Hamiltonian structure of scalar-tensor theories beyond Horndeski. J. Cosmol. Astropart. Phys. 2014, 10, 071. [Google Scholar] [CrossRef]
  158. Gleyzes, J.; Langlois, D.; Piazza, F.; Vernizzi, F. Exploring gravitational theories beyond Horndeski. J. Cosmol. Astropart. Phys. 2015, 2, 018. [Google Scholar] [CrossRef]
  159. Gao, X. Hamiltonian analysis of spatially covariant gravity. Phys. Rev. D 2014, 90, 104033. [Google Scholar] [CrossRef] [Green Version]
  160. Kase, R.; Tsujikawa, S. Effective field theory approach to modified gravity including Horndeski theory and Hořava–Lifshitz gravity. Int. J. Mod. Phys. D 2015, 23, 1443008. [Google Scholar] [CrossRef] [Green Version]
  161. Frusciante, N.; Raveri, M.; Vernieri, D.; Hu, B.; Silvestri, A. Hořava Gravity in the Effective Field Theory formalism: From cosmology to observational constraints. Phys. Dark Univ. 2016, 13, 7–24. [Google Scholar] [CrossRef] [Green Version]
  162. Horava, P. Quantum Gravity at a Lifshitz Point. Phys. Rev. D 2009, 79, 084008. [Google Scholar] [CrossRef] [Green Version]
  163. Motohashi, H.; Noui, K.; Suyama, T.; Yamaguchi, M.; Langlois, D. Healthy degenerate theories with higher derivatives. J. Cosmol. Astropart. Phys. 2016, 7, 033. [Google Scholar] [CrossRef] [Green Version]
  164. Kobayashi, T. Horndeski theory and beyond: A review. Rept. Prog. Phys. 2019, 82, 086901. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  165. Babichev, E.; Charmousis, C. Dressing a black hole with a time-dependent Galileon. J. High Energy Phys. 2014, 8, 106. [Google Scholar] [CrossRef] [Green Version]
  166. Charmousis, C.; Copeland, E.J.; Padilla, A.; Saffin, P.M. General second order scalar-tensor theory, self tuning, and the Fab Four. Phys. Rev. Lett. 2012, 108, 051101. [Google Scholar] [CrossRef] [Green Version]
  167. Brax, P. What makes the Universe accelerate? A review on what dark energy could be and how to test it. Rept. Prog. Phys. 2018, 81, 016902. [Google Scholar] [CrossRef]
  168. Charmousis, C.; Copeland, E.J.; Padilla, A.; Saffin, P.M. Self-tuning and the derivation of a class of scalar-tensor theories. Phys. Rev. D 2012, 85, 104040. [Google Scholar] [CrossRef]
  169. Charmousis, C. From Lovelock to Horndeski’s Generalized Scalar Tensor Theory. Lect. Notes Phys. 2015, 892, 25–56. [Google Scholar] [CrossRef]
  170. Copeland, E.J.; Padilla, A.; Saffin, P.M. The cosmology of the Fab-Four. J. Cosmol. Astropart. Phys. 2012, 12, 026. [Google Scholar] [CrossRef] [Green Version]
  171. Chiba, T.; Smith, T.L.; Erickcek, A.L. Solar System constraints to general f(R) gravity. Phys. Rev. D 2007, 75, 124014. [Google Scholar] [CrossRef] [Green Version]
  172. Abbott, B.P.; et al.; [LIGO Scientific Collaboration and Virgo Collaboration] GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral. Phys. Rev. Lett. 2017, 119, 161101. [Google Scholar] [CrossRef] [Green Version]
  173. Ezquiaga, J.M.; Zumalacárregui, M. Dark Energy After GW170817: Dead Ends and the Road Ahead. Phys. Rev. Lett. 2017, 119, 251304. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  174. de Rham, C.; Melville, S. Gravitational Rainbows: LIGO and Dark Energy at its Cutoff. Phys. Rev. Lett. 2018, 121, 221101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  175. Baker, T.; Calcagni, G.; Chen, A.; Fasiello, M.; Lombriser, L.; Martinovic, K.; Pieroni, M.; Sakellariadou, M.; Tasinato, G.; Bertacca, D.; et al. Measuring the propagation speed of gravitational waves with LISA. J. Cosmol. Astropart. Phys. 2022, 8, 031. [Google Scholar] [CrossRef]
  176. Linder, E.V. How Fabulous Is Fab 5 Cosmology? J. Cosmol. Astropart. Phys. 2013, 12, 032. [Google Scholar] [CrossRef] [Green Version]
  177. Appleby, S.; Linder, E.V. The Well-Tempered Cosmological Constant. J. Cosmol. Astropart. Phys. 2018, 7, 034. [Google Scholar] [CrossRef] [Green Version]
  178. Appleby, S.; Linder, E.V. The Well-Tempered Cosmological Constant: The Horndeski Variations. J. Cosmol. Astropart. Phys. 2020, 12, 036. [Google Scholar] [CrossRef]
  179. Bernardo, R.C.; Said, J.L.; Caruana, M.; Appleby, S. Well-tempered teleparallel Horndeski cosmology: A teleparallel variation to the cosmological constant problem. J. Cosmol. Astropart. Phys. 2021, 10, 078. [Google Scholar] [CrossRef]
  180. Khan, A.; Taylor, A. A minimal self tuning model to solve the cosmological constant problem. arXiv 2022, arXiv:2201.09016. [Google Scholar] [CrossRef]
  181. Copeland, E.J.; Ghataore, S.; Niedermann, F.; Padilla, A. Generalised scalar-tensor theories and self-tuning. J. Cosmol. Astropart. Phys. 2022, 3, 004. [Google Scholar] [CrossRef]
  182. Appleby, S.; Bernardo, R.C. Tadpole cosmology: Self tuning without degeneracy. J. Cosmol. Astropart. Phys. 2022, 7, 035. [Google Scholar] [CrossRef]
  183. Lacombe, O.; Mukohyama, S. Self-tuning of the cosmological constant in brane-worlds with P(X,f). arXiv 2022, arXiv:2203.16322. [Google Scholar]
  184. Amariti, A.; Charmousis, C.; Forcella, D.; Kiritsis, E.; Nitti, F. Brane cosmology and the self-tuning of the cosmological constant. J. Cosmol. Astropart. Phys. 2019, 10, 007. [Google Scholar] [CrossRef] [Green Version]
  185. Charmousis, C.; Kiritsis, E.; Nitti, F. Holographic self-tuning of the cosmological constant. J. High Energy Phys. 2017, 9, 031. [Google Scholar] [CrossRef] [Green Version]
  186. Bernardo, R.C.; Said, J.L.; Caruana, M.; Appleby, S. Well-tempered Minkowski solutions in teleparallel Horndeski theory. Class. Quant. Grav. 2022, 39, 015013. [Google Scholar] [CrossRef]
  187. Blumenhagen, R.; Cvetic, M.; Langacker, P.; Shiu, G. Toward realistic intersecting D-brane models. Ann. Rev. Nucl. Part. Sci. 2005, 55, 71–139. [Google Scholar] [CrossRef] [Green Version]
  188. Maharana, A.; Palti, E. Models of Particle Physics from Type IIB String Theory and F-theory: A Review. Int. J. Mod. Phys. A 2013, 28, 1330005. [Google Scholar] [CrossRef] [Green Version]
  189. Kachru, S.; Kallosh, R.; Linde, A.D.; Trivedi, S.P. De Sitter vacua in string theory. Phys. Rev. D 2003, 68, 046005. [Google Scholar] [CrossRef]
  190. Balasubramanian, V.; Berglund, P.; Conlon, J.P.; Quevedo, F. Systematics of moduli stabilisation in Calabi-Yau flux compactifications. J. High Energy Phys. 2005, 03, 007. [Google Scholar] [CrossRef] [Green Version]
  191. Grana, M. Flux compactifications in string theory: A Comprehensive review. Phys. Rept. 2006, 423, 91–158. [Google Scholar] [CrossRef] [Green Version]
  192. Blumenhagen, R.; Kors, B.; Lust, D.; Stieberger, S. Four-dimensional String Compactifications with D-Branes, Orientifolds and Fluxes. Phys. Rept. 2007, 445, 1–193. [Google Scholar] [CrossRef] [Green Version]
  193. Koyama, K. The cosmological constant and dark energy in braneworlds. Gen. Rel. Grav. 2008, 40, 421–450. [Google Scholar] [CrossRef] [Green Version]
  194. Arkani-Hamed, N.; Dimopoulos, S.; Kaloper, N.; Sundrum, R. A Small cosmological constant from a large extra dimension. Phys. Lett. B 2000, 480, 193–199. [Google Scholar] [CrossRef] [Green Version]
  195. Kachru, S.; Schulz, M.B.; Silverstein, E. Selftuning flat domain walls in 5-D gravity and string theory. Phys. Rev. D 2000, 62, 045021. [Google Scholar] [CrossRef] [Green Version]
  196. Forste, S.; Lalak, Z.; Lavignac, S.; Nilles, H.P. A Comment on selftuning and vanishing cosmological constant in the brane world. Phys. Lett. B 2000, 481, 360–364. [Google Scholar] [CrossRef] [Green Version]
  197. Forste, S.; Lalak, Z.; Lavignac, S.; Nilles, H.P. The Cosmological constant problem from a brane world perspective. J. High Energy Phys. 2000, 09, 034. [Google Scholar] [CrossRef] [Green Version]
  198. Csaki, C.; Erlich, J.; Grojean, C.; Hollowood, T.J. General properties of the selftuning domain wall approach to the cosmological constant problem. Nucl. Phys. B 2000, 584, 359–386. [Google Scholar] [CrossRef] [Green Version]
  199. Binetruy, P.; Cline, J.M.; Grojean, C. Dynamical instability of brane solutions with a self-tuning cosmological constant. Phys. Lett. B 2000, 489, 403–410. [Google Scholar] [CrossRef]
  200. Randall, L.; Sundrum, R. An Alternative to compactification. Phys. Rev. Lett. 1999, 83, 4690–4693. [Google Scholar] [CrossRef] [Green Version]
  201. Vilenkin, A. Gravitational Field of Vacuum Domain Walls and Strings. Phys. Rev. D 1981, 23, 852–857. [Google Scholar] [CrossRef]
  202. Gott, J.R., III. Gravitational lensing effects of vacuum strings: Exact solutions. Astrophys. J. 1985, 288, 422–427. [Google Scholar] [CrossRef]
  203. Hiscock, W.A. Exact Gravitational Field of a String. Phys. Rev. D 1985, 31, 3288–3290. [Google Scholar] [CrossRef] [PubMed]
  204. Kaloper, N.; Kiley, D. Charting the landscape of modified gravity. J. High Energy Phys. 2007, 05, 045. [Google Scholar] [CrossRef] [Green Version]
  205. Sundrum, R. Compactification for a three-brane universe. Phys. Rev. D 1999, 59, 085010. [Google Scholar] [CrossRef] [Green Version]
  206. Niedermann, F.; Schneider, R. Radially stabilized inflating cosmic strings. Phys. Rev. D 2015, 91, 064010. [Google Scholar] [CrossRef] [Green Version]
  207. Niedermann, F.; Schneider, R.; Hofmann, S.; Khoury, J. Universe as a cosmic string. Phys. Rev. D 2015, 91, 024002. [Google Scholar] [CrossRef] [Green Version]
  208. Dubovsky, S.L.; Rubakov, V.A. Brane induced gravity in more than one extra dimensions: Violation of equivalence principle and ghost. Phys. Rev. D 2003, 67, 104014. [Google Scholar] [CrossRef] [Green Version]
  209. Hassan, S.F.; Hofmann, S.; von Strauss, M. Brane Induced Gravity, its Ghost and the Cosmological Constant Problem. J. Cosmol. Astropart. Phys. 2011, 01, 020. [Google Scholar] [CrossRef]
  210. Eglseer, L.; Niedermann, F.; Schneider, R. Brane induced gravity: Ghosts and naturalness. Phys. Rev. D 2015, 92, 084029. [Google Scholar] [CrossRef] [Green Version]
  211. Arkani-Hamed, N.; Dimopoulos, S.; Dvali, G.R. The Hierarchy problem and new dimensions at a millimeter. Phys. Lett. B 1998, 429, 263–272. [Google Scholar] [CrossRef] [Green Version]
  212. Arkani-Hamed, N.; Dimopoulos, S.; Dvali, G.R. Phenomenology, astrophysics and cosmology of theories with submillimeter dimensions and TeV scale quantum gravity. Phys. Rev. D 1999, 59, 086004. [Google Scholar] [CrossRef] [Green Version]
  213. Chen, J.W.; Luty, M.A.; Ponton, E. A Critical cosmological constant from millimeter extra dimensions. J. High Energy Phys. 2000, 9, 012. [Google Scholar] [CrossRef]
  214. Carroll, S.M.; Guica, M.M. Sidestepping the cosmological constant with football shaped extra dimensions. arXiv 2003, arXiv:hep-th/0302067. [Google Scholar]
  215. Navarro, I. Codimension two compactifications and the cosmological constant problem. J. Cosmol. Astropart. Phys. 2003, 9, 004. [Google Scholar] [CrossRef]
  216. Cline, J.M.; Descheneau, J.; Giovannini, M.; Vinet, J. Cosmology of codimension two brane worlds. J. High Energy Phys. 2003, 6, 048. [Google Scholar] [CrossRef] [Green Version]
  217. Aghababaie, Y.; Burgess, C.P.; Parameswaran, S.L.; Quevedo, F. Towards a naturally small cosmological constant from branes in 6-D supergravity. Nucl. Phys. B 2004, 680, 389–414. [Google Scholar] [CrossRef] [Green Version]
  218. Burgess, C.P.; van Nierop, L. Large Dimensions and Small Curvatures from Supersymmetric Brane Back-reaction. J. High Energy Phys. 2011, 4, 078. [Google Scholar] [CrossRef] [Green Version]
  219. Burgess, C.P.; van Nierop, L. Technically Natural Cosmological Constant From Supersymmetric 6D Brane Backreaction. Phys. Dark Univ. 2013, 2, 1–16. [Google Scholar] [CrossRef] [Green Version]
  220. Salam, A.; Sezgin, E. Chiral Compactification on Minkowski x S**2 of N = 2 Einstein-Maxwell Supergravity in Six-Dimensions. Phys. Lett. B 1984, 147, 47. [Google Scholar] [CrossRef]
  221. Randjbar-Daemi, S.; Salam, A.; Sezgin, E.; Strathdee, J.A. An Anomaly Free Model in Six-Dimensions. Phys. Lett. B 1985, 151, 351–356. [Google Scholar] [CrossRef]
  222. Nishino, H.; Sezgin, E. The Complete N=2, d=6 Supergravity With Matter and Yang-Mills Couplings. Nucl. Phys. B 1986, 278, 353–379. [Google Scholar] [CrossRef]
  223. Kapner, D.J.; Cook, T.S.; Adelberger, E.G.; Gundlach, J.H.; Heckel, B.R.; Hoyle, C.D.; Swanson, H.E. Tests of the gravitational inverse-square law below the dark-energy length scale. Phys. Rev. Lett. 2007, 98, 021101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  224. Randjbar-Daemi, S.; Salam, A.; Strathdee, J.A. Spontaneous Compactification in Six-Dimensional Einstein-Maxwell Theory. Nucl. Phys. B 1983, 214, 491–512. [Google Scholar] [CrossRef]
  225. Niedermann, F.; Schneider, R. Fine-tuning with Brane-Localized Flux in 6D Supergravity. J. High Energy Phys. 2016, 2, 025. [Google Scholar] [CrossRef] [Green Version]
  226. Burgess, C.P.; Diener, R.; Williams, M. EFT for Vortices with Dilaton-dependent Localized Flux. J. High Energy Phys. 2015, 11, 054. [Google Scholar] [CrossRef] [Green Version]
  227. Burgess, C.P.; Diener, R.; Williams, M. Self-Tuning at Large (Distances): 4D Description of Runaway Dilaton Capture. J. High Energy Phys. 2015, 10, 177. [Google Scholar] [CrossRef] [Green Version]
  228. Niedermann, F.; Schneider, R. SLED Phenomenology: Curvature vs. Volume. J. High Energy Phys. 2016, 3, 130. [Google Scholar] [CrossRef] [Green Version]
  229. Burgess, C.P.; Diener, R.; Williams, M. A problem with d-functions: Stress-energy constraints on bulk-brane matching (with comments on arXiv:1508.01124). J. High Energy Phys. 2016, 1, 017. [Google Scholar] [CrossRef] [Green Version]
  230. Gibbons, G.W.; Gueven, R.; Pope, C.N. 3-branes and uniqueness of the Salam-Sezgin vacuum. Phys. Lett. B 2004, 595, 498–504. [Google Scholar] [CrossRef]
  231. Karch, A.; Randall, L. Locally localized gravity. J. High Energy Phys. 2001, 5, 008. [Google Scholar] [CrossRef] [Green Version]
  232. Kaloper, N.; Sorbo, L. Locally localized gravity: The Inside story. J. High Energy Phys. 2005, 8, 070. [Google Scholar] [CrossRef] [Green Version]
  233. Ghosh, J.K.; Kiritsis, E.; Nitti, F.; Witkowski, L.T. De Sitter and Anti-de Sitter branes in self-tuning models. J. High Energy Phys. 2018, 11, 128. [Google Scholar] [CrossRef] [Green Version]
  234. Abbott, T.M.C.; et al.; [DES Collaboration] Dark Energy Survey Year 3 Results: Constraints on extensions to LCDM with weak lensing and galaxy clustering. arXiv 2022, arXiv:2207.05766. [Google Scholar]
  235. Ishak, M.; Baker, T.; Bull, P.; Pedersen, E.M.; Blazek, J.; Ferreira, P.G.; Leonard, C.D.; Lin, W.; Linder, E.; Pardo, K.; et al. Modified Gravity and Dark Energy models Beyond w(z)CDM Testable by LSST. arXiv 2019, arXiv:1905.09687. [Google Scholar]
  236. Weltman, A.; Bull, P.; Camera, S.; Kelley, K.; Padmanabhan, H.; Pritchard, J.; Raccanelli, A.; Riemer-Sørensen, S.; Shao, L.; Andrianomena, S.; et al. Fundamental physics with the Square Kilometre Array. Publ. Astron. Soc. Aust. 2020, 37, e002. [Google Scholar] [CrossRef] [Green Version]
  237. Belgacem, E.; Calcagni, G.; Crisostomi, M.; Dalang, C.; Dirian, Y.; Ezquiaga, J.M.; Fasiello, M.; Foffa, S.; Ganz, A.; García-Bellido, J.; et al. Testing modified gravity at cosmological distances with LISA standard sirens. J. Cosmol. Astropart. Phys. 2019, 7, 024. [Google Scholar] [CrossRef] [Green Version]
  238. Arun, K.G.; Belgacem, E.; Benkel, R.; Bernard, L.; Berti, E.; Bertone, G.; Besancon, M.; Blas, D.; Böhmer, C.G.; Brito, R.; et al. New horizons for fundamental physics with LISA. Living Rev. Rel. 2022, 25, 4. [Google Scholar] [CrossRef]
Figure 1. The normalized effective Newton’s constant as a function of dark matter halo radius R h in the chameleon screened Hu–Sawicki f ( R ) gravity (blue, dashed) and the Vainshtein screened DGP gravity (red, solid). R h , i is the initial size of the halo and Ω rc 1 / ( 4 H 0 2 r c 2 ) . The amplitudes in the unscreened regime reflect the specific values of f R 0 and r c chosen.
Figure 1. The normalized effective Newton’s constant as a function of dark matter halo radius R h in the chameleon screened Hu–Sawicki f ( R ) gravity (blue, dashed) and the Vainshtein screened DGP gravity (red, solid). R h , i is the initial size of the halo and Ω rc 1 / ( 4 H 0 2 r c 2 ) . The amplitudes in the unscreened regime reflect the specific values of f R 0 and r c chosen.
Universe 09 00063 g001
Figure 2. Geometry (a) and (b) are employed in higher-dimensional generalizations of the DGP model (infinite extra space volume) and the SLED model (finite extra space volume), respectively. The red dot corresponds to our Universe.
Figure 2. Geometry (a) and (b) are employed in higher-dimensional generalizations of the DGP model (infinite extra space volume) and the SLED model (finite extra space volume), respectively. The red dot corresponds to our Universe.
Universe 09 00063 g002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Foundational Aspects of Dark Energy (FADE) Collaboration; Bernardo, H.; Bose, B.; Franzmann, G.; Hagstotz, S.; He, Y.; Litsa, A.; Niedermann, F. Modified Gravity Approaches to the Cosmological Constant Problem. Universe 2023, 9, 63. https://doi.org/10.3390/universe9020063

AMA Style

Foundational Aspects of Dark Energy (FADE) Collaboration, Bernardo H, Bose B, Franzmann G, Hagstotz S, He Y, Litsa A, Niedermann F. Modified Gravity Approaches to the Cosmological Constant Problem. Universe. 2023; 9(2):63. https://doi.org/10.3390/universe9020063

Chicago/Turabian Style

Foundational Aspects of Dark Energy (FADE) Collaboration, Heliudson Bernardo, Benjamin Bose, Guilherme Franzmann, Steffen Hagstotz, Yutong He, Aliki Litsa, and Florian Niedermann. 2023. "Modified Gravity Approaches to the Cosmological Constant Problem" Universe 9, no. 2: 63. https://doi.org/10.3390/universe9020063

APA Style

Foundational Aspects of Dark Energy (FADE) Collaboration, Bernardo, H., Bose, B., Franzmann, G., Hagstotz, S., He, Y., Litsa, A., & Niedermann, F. (2023). Modified Gravity Approaches to the Cosmological Constant Problem. Universe, 9(2), 63. https://doi.org/10.3390/universe9020063

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop