Quantum corrections to vacuum energy sequestering (with monodromy)

Field theory models of axion monodromy have been shown to exhibit vacuum energy sequestering as an emergent phenomenon for cancelling radiative corrections to the cosmological constant. We study one loop corrections to this class of models coming from virtual axions using a heat kernel expansion. We find that the structure of the original sequestering proposals is no longer preserved at low energies. Nevertheless, the cancellation of radiative corrections to the cosmological constant remains robust, even with the new structures required by quantum corrections.

theory: they are radiatively stable thanks to (approximate) chiral symmetry but their values are given by a measurement [20]. The cancellation mechanism in sequestering models can also be understood in terms of so-called decapitation [21,22].
It is natural to ask whether or not the sequestering mechanism can emerge in a low energy effective theory arising from string compactifications. To this end, it was very recently shown [16,17] how sequestering models can emerge from a pair of (deformed) field theory monodromies [23][24][25][26][27][28][29][30]. In field theory monodromy, there is a bilinear mixing between a scalar and a four-form field strength. The dynamics is that of a very massive scalar field and the dual of the four form. The latter is constant in spacetime but can play the role of a rigid degree of freedom familiar to sequestering. Field theory monodromy set-ups have been proposed as the low energy description of monodromy inflation in string theory [31,32], so we can be optimistic about finding a realisation of sequestering in a fundamental theory.
To gain further insight into what could emerge from string theory, we adopt a bottom up approach and consider quantum corrections to monodromy models with an emergent sequestering mechanism. In particular, we integrate out fluctuations in the two axion fields at one-loop using heat kernel techniques to expand the effective action in powers of derivatives. We find that the original constructions are not closed under renormalisation, even to a fixed order in the expansion. To carry out the renormalization program, we henceforth modify the original theory by including all necessary operators. With the renormalized theory in hand, we ask whether or not the sequestering mechanism is still operative. Remarkably, we find that under very reasonable assumptions on the parameters of the theory, the sequestering mechanism survives. The rest of this paper is organised as follows. In section 2, we review the emergence of the sequestering mechanism in the original deformations of a pair of field theory models of axion monodromy. We generalise the results of [16,17], slightly, allowing both axions to fluctuate and including different versions of low energy sequestering. In section 3, we compute the quantum corrections to these theories due to virtual axions, at one loop, using a heat kernel expansion. We also elucidate the role of the global degrees of freedom at the level of the path integral. In the end we find corrections to the original models of [16,17]. In section 4, we explicitly demonstrate how these corrections still allow for an emergent cancellation of radiative corrections to vacuum energy, via a generalisation of the sequestering mechanism. We conclude in section 5.

Vacuum energy sequestering and monodromy
In [16,17], it was shown how the sequestering mechanism could emerge from a pair of (deformed) monodromy constructions in field theory. Generalising slightly, these models could be described by the following action, (2.1) The first line is just the Einstein-Hilbert action with matter fields, Ψ, minimally coupled to the metric, g µν . The second line contains a pair of field theory monodromy constructions [23][24][25][26][27][28][29][30], with axions σ i and four-form field strengths, F i µναβ = 4∂ [µ A i ναβ] whose bilinear mixing is controlled by the mass scale µ i . The final line contains the deformations that break the (discrete) shift symmetry for the axions, with independent mass scalesm and λ. Note that in the interests of brevity we have neglected to include boundary terms dependent on g µν and σ i , required in order to have a well defined variational principle. Their precise form is non-trivial and specified in [13,18], but not required for our analysis.
If we identify σ 2 with the inflaton, as in [16], one can generate the potential λ 4 F λ mσ 2 λ 2 from loops of matter coupled to the inflaton 1 . The other deformation is assumed to be of gravitational origin, consistent with a breaking of the shift symmetry induced by gravitational instantons [35,36]. Its precise form depends on the model of sequestering we wish to emerge at low energies. For earlier models of sequestering (see e.g. [12]) we take O[g] = R. For the more sophisticated later model [14], we take O[g] = R GB def = R µναβ R µναβ −4R µν R µν +R 2 to be the Gauss-Bonnet combination. This model was recently dubbed Omnia Sequestra [18] on account of its ability to sequester radiative corrections to vacuum energy from both matter and graviton loops. Note that our convention for the metric is (− + ++), the Riemann tensor is R µ ναβ = ∂ β Γ µ να − · · · and the Ricci tensor is R µν = R λ µνλ . We also denote the mass dimension of the operator O by [O]. To see how sequestering emerges at low energies, it is convenient to switch to a dual description by integrating out the four-form field strengths and writing the action in terms of their magnetic duals. To do this, we add a pair of continuous Lagrange multipliers, , and include functional variation over both Q i and F i µναβ . The four-forms now enter algebraically, and at quadratic order, so they can be integrated out exactly, yielding The Q i are constrained to be constant on-shell by the three-form variation. Furthermore, their expectation values are quantised in units of the membrane charge, q i , as in Q i = 2πN q i [38,39]. The second line of the action is manifestly invariant under a simultaneous 1 Such couplings can be motivated by reheating [16]. For example, suppose there is a coupling ησ 2 2 h 2 to some massive scalar h of mass m h (possibly the Higgs), then loops of h will generate a term of the form m 2 σ 2 2 wherem ∼ ηm h . For an effective theory cut-off at some scale M ∼ 4πλ, we can motivate the higher order operators using naive dimensional analysis [33,34], treating the mass,m, as a spurion [16].
shift in σ i and Q i . Non-perturbative effects break this down to a discrete symmetry [40], consistent with the quantisation condition on the Q i , where the periodicity of the axions is given by f i = q i /µ i . We now see how the axion, σ i has gained a mass, µ i , from its bilinear mixing with the four-form field strength. To connect directly to sequestering models, we re-define A i → A i /µ i and take the limit where µ i → ∞, Q i → ∞ holding s i = −Q i /µ i fixed. In this limit the axion settles into its minimum at s i so that the action now reads This is precisely the form of the actions presented in [12,14], suggesting that the low energy dynamics of our deformed monodromies (2.1) is equivalent to a sequestering theory. Furthermore, given that the cancellation of vacuum energy occurs at infinite wavelength in sequestering, we might expect it to still occur even in the full monodromy theory (2.1).
We shall now show that this is indeed the case. We start with the action (2.2) and go to the global limit by integrating out the threeforms (see appendix A). This yields where the Q i are assumed to be global Lagrange multipliers constant in spacetime and C i def = M dA i = ∂M A i is the integral of the three-form flux across the boundary, ∂M, of the spacetime, M. Our notation reflects the fact that this action is a functional of the fields σ i , g, Ψ, but a function of the global variables Q i . Variation with respect to the axions and the metric yields the following local equations δg µν is the energy-momentum tensor in the matter sector and is the variation of the gravity deformation (see e.g. [41]). Here P µανβ is the double dual of the Riemann tensor. The energy-momentum tensor for the axions appears in the last terms of (2.6c). Variation with respect to the global Lagrange multipliers yields the following global constraints, √ g σ i )/µ i , and we can rewrite these global constraints as Integrating the scalar equations of motion (2.6a) and (2.6b) over spacetime and taking the ratio yields a global geometric constant which can be written as where we have introduced the weighted spacetime average Q w = −gw , and we are taking w = F ′ O . This constraint is crucial. As long as the operator O depends on the scale dependent part of the curvature, the constraint (2.10) fixes the large wavelength mode of the scalar curvature in a radiatively stable way [14,15]. At the level of the gravity equation (2.6c), we identify the long wavelength mode of λ 4 F λ + i 1 2 (µ i σ i + Q i ) 2 with the cosmological constant counterterm, whose value is dynamically fixed by the constraint (2.10). Indeed, we can eliminate λ 4 F λ + i 1 2 (µ i σ i + Q i ) 2 w from (2.6c) by taking the trace and (weighted) spacetime average, allowing us to rewrite the gravity equation as where δχ = χ − χ w is the localised fluctuation in a generic field χ and is the effective Planck mass. Notice that we have split the effective potential into two parts: the localised part, corresponding to δ λ 4 F λ + . . . , and the global, infinite wavelength part, Λ global . The latter is given by (2.14) Derivation of the last relation uses the fact we can write the Gauss-Bonnet term as where W µναβ is the Weyl tensor [14]. Let us first consider radiative corrections from the matter sector and their backreaction on to the geometry. To see this cleanly, it is convenient to go to the global limit and consider maximally symmetric configurations of constant curvature, R µναβ =R 12 (g µα g νβ − g µβ g να ), constant scalars, σ i =σ i , and a vacuum energy source, T µν = −V vac g µν . Plugging this ansatz into the effective gravity equation (2.11) yields It follows that the curvature is stable against radiative corrections from the matter sector provided O w is radiatively stable, in both cases. However, O w is constrained by the fluxes according to equation (2.10), so the question is now whether or not the right hand side of (2.10) is radiatively stable. The fluxes C i are set as as infra-red boundary condition and have no ultra-violet sensitivity. In contrast, radiative corrections will affect the values for the scalars,σ i , and this could contaminate O w through the averaged prefactor, going in our global limit. The radiative stability of this prefactor can be inferred For a theory cut-off at some ultra-violet scale, M, these go as [42] where m is a typical mass scale in the effective field theory governing the matter sector.
which is expected to be the case provided M M P . Similarly, radiative stability of F λ requires it to go as O(1) M 4 16π 2 λ 4 . This is also expected to be the case precisely because λ 4 F λ plays the role of the counter-term cancelling off radiative corrections to vacuum energy going as O M 4 16π 2 . For generic smooth functions F O and F λ , this stability is transferred to the prefactor We therefore see that O w is radiatively stable against matter loops, and by association the same must be true for the long wavelength mode of the spacetime curvature,R.
Although we have focussed on the global limit, this is just for convenience of presentation. The local fluctuations are irrelevant when considering the gravitational effects of vacuum energy since they cannot be sourced by it. That is not to say the local fluctuations do not contribute to the geometry on large scales. They can, through their (weighted) spacetime averages entering in (2.13). However, these contributions are immune from the radiative instabilities contained in the vacuum energy. For the monodromy set-up studied here, the local fluctuations can, of course, backreact on to the geometry on shorter scales. For the case of σ 2 , such fluctuations can even give rise to inflation, as suggested in [16]. Fluctuations in σ 1 are potentially more exotic since they represent a departure from local General Relativity [43]. As long as we make the mass of the fluctuations, µ 1 , sufficiently large, beyond the scale of inflation, there will be no conflict with observation.
Let us now turn our attention to graviton loops. This is where our two models go their separate ways. For O[g] = R, as in earlier sequestering proposals [12], graviton loops can and will induce a UV sensitive potential for σ 1 . This is problematic because it will source equation (2.6a) and contaminate the constraint (2.10) [14]. In contrast, when O[g] = R GB , there is an approximate shift symmetry F O → F O + c, becoming exact in the limit µ 1 , Q 1 → ∞, µ 1 C 1 → 0, holding Q 1 /µ 1 fixed, as we can easily see by staring at the limiting action (2.4). The shift symmetry is also exact even without taking these limits, provided we take F O to be linear. In any event, this shift symmetry, exact or otherwise, helps protect us from a perturbatively generated UV sensitive potential for σ 1 and contamination of the constraint (2.10). Non-perturbative effects will break the symmetry down to discrete shifts, of course, although these corrections will be suppressed for large instanton action [44].
The UV completion of sequestering via monodromy brings in another source of radiative corrections, from axion loops. The main concern, here, is that the action (2.1) is not closed under quantum corrections, and that to fix this we spoil the cancellation mechanism that renders the observed cosmological constant radiatively stable. We will investigate this in detail in the next two sections.

The quantum effective action
To what extent does the sequestering mechanism survive radiative corrections from the gravitational sector? A complete analysis requires us to consider the effects of quantum fluctuations of all fields, including gravitons, axions, 3-forms and, of course, the matter fields. As a first step, here we will integrate out the fluctuations in the axions at one loop and compute the corresponding one-loop effective action for sequestering. The metric will be treated as a classical background field while quantum fluctuations in the matter sector will only be considered through their contributions to the vacuum energy. We will neglect sources for the 3-forms allowing us to integrate them out exactly.
Our starting point will be the action (2.1). The path integral is given by the functional integration over the dynamical fields where N is a constant normalisation factor. We have focussed on the integration over the axions and the 3-forms, although as stated above, here we only consider sources for the axions. The action S[A i , σ i , . . .] is given by (2.1). It also includes boundary terms in order to have a well defined variational principle. As explained in [13,18], for the sequestering cancellation to occur we need to retain the global variations in the (pseudo) scalars. This means we must impose Neumann, rather than Dirichlet, boundary conditions on them, which in turn requires us to adjust the boundary condition on the metric 2 . In what follows, we assume the action to be supplemented with boundary terms necessary for the heat kernel method to be well defined given our choice of boundary conditions. To proceed, we pass to the dual description described in the previous section. To do so, recall that we promote F i µναβ to a field, introduce a Lagrange multiplier, Q i , to impose the constraint, 1 4! ǫ µναβ (F i µναβ − 4∂ [µ A i ναβ] ) = 0, and include functional variation over both Q i and F i µναβ . Since the F i µναβ now enters the action up to quadratic order we are able to integrate it out exactly, giving the action (2.2) and an overall constant normalisation in the path integral, which we specify in the appendix A. The result is where N ′ is, again, a constant normalisation factor and S[A i , σ i , Q i . . .] is given by (2.2). Now, as explained in detail in the appendix A, performing the 3-form integration yields a delta functional, forcing the gradient, ∂ µ Q i , to vanish. This suppresses local, but not global, off-shell fluctuations in the Q i , reducing the path integration over fields, Q i (x), to an ordinary integration over the global variables, Q i . As a result, we now work with the action of the form (2.5) and a path integral This structure is reminiscent of a gas of "Universes" whose coupling constants are random variables to be summed over (see e. g. [45,46]). Focussing on one such "Universe" (with Q i fixed) we now integrate out the fluctuations in the axions at one loop. To this end, we derive the resulting quantum effective action using a heat kernel expansion. Rather than discussing both sequestering models in turn as done in the previous section, we will study a generalised version of action (2.5) which couples If we want to make contact with the original models, we can project on the corresponding parameter subspace whenever deemed necessary. Anticipating the final result, in order to renormalize the theory we are also forced to generalize the action to include a coupling between σ 2 and curvature and thus we have instead Here we treat M P , λ,m as fixed scales and capture the RG flow in terms of dimensionless couplings, such as Z, or coefficients in the functional definition of the F R , F GB and F λ . Note that both F R and F GB retain full dependence on both axions, except that we shall still assume the absence of any direct mixing (see the discussion in section 5), so that For a generic scalar field, σ, in d dimensional curved space, the effective action can be expressed as the trace of the heat kernel where the heat kernel is a solution to a heat-diffusion equation and as such is defined as a bi-tensor density of weight-1/2 and D is a second-order differential operator whose eigenmodes are scalar densities of weight-1/2.
where Q(x) is a space-time dependent function of the metric and all background fields.
In each case we have a candidate mass parameter given by the µ i . Formally, ultraviolet divergences in the effective action (3.6) are controlled by the small s behaviour ofH(x, y; s), and we can use the asymptotic expansion, as s → 0, of the heat kernel; see [37] and references therein. Explicitly, the trace takes the form where the ∼ sign reflects the fact that this is an asymptotic expansion. The central problem then reduces to computing the coefficient functions E k for a given operator. Let us first note that where {R 3 , ∇∇R 2 , R} denotes the set of all purely geometric invariants containing six derivatives.
Given the form of the action (3.4), for each axion, σ i , we obtain the differential operators where where the superscript 0 denotes background values, nevertheless, we will drop it from now on to simplify notation. It is essential to notice that, regardless of the boundary conditions imposed on quantum fluctuations, the differential operator in Eq. (3.15) is self-adjoint with respect to the natural inner product on the space of real functions [37]. Moreover, any boundary terms resulting due to integration by parts can be made to vanish under the relevant boundary conditions. With the above ingredients, we re-write the effective action as follows In what follows, we shall use the following replacement to pass to 4 dimensions where µ is the scale of dimensional regularization. We shall also use the MS renormalization scheme and define where γ E is the Euler-Mascheroni constant. As stated above, the k = 0, 1, 2 contributions in (3.18) are divergent and in general we obtain, Now, we have to specify the exact forms of the functions F R,i , F GB,i and F λ . We elucidated in the previous section that these functions must not be linear, but otherwise arbitrary; accordingly, In the following we set c (i) k≥3 = 0 and β k≥5 = 0. 4 We will see that within our analysis this is a self-consistent choice as it is not spoiled by axion loops. Let us stress though that these operators (as well as those describing a mixing between both axion sectors or an axion coupling to higher curvature invariants) are expected to be produced by quantum gravity corrections and are therefore generically present in an effective theory of sequestering. Here, on the other hand, we are interested in a first nontrivial check of the model's stability under axion loops rather than a exhaustive analysis of all EFT operators. Correspondingly, we only discuss the simplified version of (3.22) and speculate later about the relevance of the neglected operators (see section 5 for a more detailed discussion of this point).
We further amend the effective action to include terms quadratic in curvature, where u, v are the bare Wilson coefficients. We can now plug back in (3.21) and simply determine the renormalization constants which are given in appendix B by Eqs. (B.2). The corresponding beta functions of the various couplings are listed in the appendix as well.
We take as our initial conditions for the RGEs the values of the couplings at some scale µ ⋆ . It is indeed possible to solve the coupled set of RGEs, nevertheless, we choose not to perform the explicit computation since it is not essential for our purposes. Hence, we finally arrive at the renormalised sequestering action correct up to one loop (suppressing boundary terms)

25c)
and we defined effective coupling constants as functions of the physically measurable quantities. Explicitly, Let us comment finally on several important issues: • Radiative stability: We can now investigate the issue of radiative stability of the various couplings in the theory. Looking at the beta functions in the appendix, we easily notice that provided we take all beta functions remain small in the perturbative regime. In particular, all potentially dangerous mass ratios are naturally O(1). This remains true even if we consider the contribution of ordinary matter fields to the running of α, provided λ is as high as the cutoff of the effective theory. Moreover, from (3.26b) and (3.26c) it can be inferred that the choice β k≥5 = 0 and c (i) k≥3 = 0 is stable as advertised. Also note that condition (3.27b) disappears if we truncate F R GB ,2 at quadratic order (corresponding to d (2) j≥3 = 0) as follows from (3.26d).
• Effective action truncation: The expansion of the effective action we utilized originates from the asymptotic expansion of the heat kernel in Eq. (3.10), which necessarily does not converge. One may nevertheless view the result as a derivative expansion, with the cutoff being the axion masses µ i . In order that the expansion is meaningful, we need ask how good it is to ignore E 3 (x), Eq. (3.14), from our analysis of sequestering in the quantum-corrected theory. In curved space, however, one should be precise about the meaning of the derivative expansion because general covariance makes it difficult to define an expansion parameter. Nevertheless, we can define a local relational measure as follows which dictates that in units of axion masses the curvature invariants forming E 3 (x) are subdominant to those in E 2 (x). Still this is not the full story, since the background values of the axions must be constrained as well to ensure the validity of the truncation. Explicitly which then translates into the following constraints where we have assumed that R ≪ µ 2 i consistent with (3.28).
• In order to make contact with both sequestering models discussed in section 2, we specify two subclasses of (3.4) that are both stable under axion loops. To be precise, we recover the old sequestering proposal by setting d

The robustness of sequestering under quantum corrections
We are now in a position to test the robustness of the sequestering mechanism after correcting the action by axion loops. For both classes of model, the quantum corrected action can be written in the compact form (4.2c) and the couplings are understood to be the effective couplings, although we drop the superscript "eff" for brevity. We will now proceed in a way similar to section 2. The field equations now yield and where T µν is the energy momentum tensors for matter and axions, defined in section 2.
Recall that we also defined a linear operator of the form The quadratic curvature corrections contribute a term Variation with respect to the global Lagrange multipliers yields the global constraints (2.8), given in terms of the flux of the 3-forms through the boundary, C i = ∂M A i . It is convenient to adjust these constraints to the form (2.9) by a redefinition of the 3-forms, Integrating the scalar equations of motion (4.3) and (4.4) over spacetime and using the adjusted form of the constraint (2.9), we obtain the following algebraic condition on the (weighted) spacetime average of the Ricci scalar, and we have introduced the shorthand The choice of weighting factor, w, is not especially important, as long as it is built from radiatively stable quantities. To make contact with section 2, a natural choice would be thereby eliminating some of the terms linear in δR = R − R w . The constraint equation (4.8) should still be enough to fix the large wavelength mode of the scalar curvature in a radiatively stable way, provided the coefficients A, B, C are themselves radiatively stable. Further, it is clear that it can be fixed to an arbitrarily small value with a judicious, and radiatively stable, choice of the flux ratio, C 1 /C 2 . The long wavelength mode of λ 4 F λ + i 1 2 (µ i σ i + Q i ) 2 will again play the role of the cosmological constant counter term, forced by (4.8) to take on precisely the right value to cancel off radiative corrections to vacuum energy. Let us see this explicitly by calculating the effective gravity equation. By taking traces and (weighted) spacetime averages of (4.5), we find that where we define δχ = χ − χ w as the localised fluctuation in a generic field χ, in terms of the weighted averages with weight w. The effective Planck mass is As before, the potential contains a localised part, corresponding to δ λ 4 F λ + . . . , and a global, infinite wavelength part, given by where Λ * = 1 4 M 2 eff R w with the corresponding formula for R w now given by the solution to (4.8). To examine the impact of radiative corrections to the vacuum energy most succinctly, it is convenient to go to the global limit in which the spacetime geometry, the scalars and the matter source are all maximally symmetric, as we did in section 2. In other words, we have constant curvature, R µναβ =R 12 (g µα g νβ −g µβ g να ), constant scalars, σ i =σ i , and a vacuum energy source, T µν = −V vac g µν . Since H µν vanishes on these configurations the effective gravity equation will once again yield an expression for the curvature given asR = R w , with R w constrained by the fluxes according to equation (4.8). The only question is whether or not the solution to this equation is guaranteed to be radiatively stable. This will indeed be the case provided the coefficients A, B, C are themselves radiatively stable. Nowm, λ and M P are fixed scales so they don't receive any corrections 5 . The fluxes, C i , are set as infra-red boundary conditions, so the only possible source of ultra-violet sensitivity in (4.8) is through the spacetime averages of the combinations of (derivatives of) F O,i and F λ . In the global limit, the spacetime averages just reduce to their global values. Now, as argued in the previous section, the couplings c j and β j , corresponding to the coefficients of the defining functions (4.2), are at most order one and stable under axion loops if conditions (3.27) are fulfilled. In particular λ µ 1 µ 2 , which we can also take as an indication of the stability under matter loops, provided that the scale λ lies at or above the field theory cut-off.
What about the size and sensitivity to changes in the vacuum energy of the arguments To study this it is convenient to consider the original equations (4.3) to (4.5) in the global limit, Given the smallness of the observed value of M 2 PR in terms of V vac we can infer good approximations for theσ i by studying this equation in the limit of vanishingR. Assuminḡ mσ 2 λ 2 1, one can straightforwardly show that Our approximation is self-consistent as long as λ lies at or above the field theory cutoff, ensuring that V vac λ 4 . We have also used the fact thatm ≪ µ 2 [required for the convergence of the heat kernel expansion, cf. (3.30a)] and that the functions contain order one couplings in their definitions (4.2). We now see how the axion field values have no UV sensitivity to leading order. Of course, the subleading corrections toσ i do depend on the vacuum energy V vac but the point is that the sensitivity is weak thanks to the assumptions we have made on the parameters of the theory. An order one rescaling of the vacuum energy, ∆V vac = O(1)V vac , will only induces changes in the axion that go as ∆σ 2 = O(1)σ 2 ∆V vac /λ 4 O(1)σ 2 and leaveσ 1 unchanged.
We conclude that the quantum corrected theory (2.1) retains the ability to sequester radiative corrections to the vacuum energy arising from loops of matter. This is a nontrivial result since the quantum corrections have introduced a new structure. As in section 2, we have simplified our analysis by focussing on the global limit, in which we consider maximally symmetric configurations. As explained then, the local fluctuations that we have ignored are irrelevant for considering the effects of vacuum energy since they cannot be sourced by it. Their only contribution to the geometry on large scales is indirect, through their contribution to the spacetime averages entering in (4.14).
There is one last thing to check. In deriving the quantum corrected action (4.1), we made an assumption that the second derivatives satisfied F ′′ required for a convergent heat kernel expansion. The condition on F ′′ R,i holds trivially given its quadratic form, with order one couplings, since F ′′ . For F ′′ λ we take the form of (4.2a), with order one couplings, and plug in the solution (4.19) , we generically require λ 2 mM P in accordance with (3.27b). However, since M P , λ µ 1 µ 2 ≫m this condition is tricky to satisfy. For example, if we were to takem to be at the GUT scale (∼ 10 16 GeV), consistent with high scale inflation, the condition becomes λ 10 17 GeV and so we can only marginally satisfy λ µ i ≫m. 6 The alternative, of course, is to assume that F R GB ,i is at most quadratic or to sacrifice the inflation scenario by choosinḡ m below the GUT scale. In any event, we conclude that the assumptions made in deriving the effective action in the previous section can be made self consistent.

Conclusions
The extension of vacuum energy sequestering models [9][10][11][12][13][14][15][16][17][18]to field theory models of monodromy [16,17] was an important step in the long term goal of realising the sequestering mechanism as an emergent phenomenon in a fundamental UV complete theory of Nature. However, this development also brought in additional moving parts, and it is important to ask whether or not they really are compatible with the original low energy models of sequestering, and more importantly, the vacuum energy cancellation mechanism itself. In this paper we have answered both of these questions, having calculated the effect of loop corrections due to those additional moving parts, the axions. The result is that the original low energy models of sequestering are incomplete because the axion loops necessarily generate new structures not seen in those original theories. However, as we have seen, for a new class of Lagrangians, closed under renormalisation up to a particular order in power counting, the cancellation of vacuum energy remains intact. As a result, the sequestering mechanism does seem to be robust against quantum corrections from axion loops.
The details of our analysis relied on a number of reasonable assumptions regarding the parameters in the theory. For the quantum corrected actions of the form (4.1), we ought to assume the following hierarchy of scales where M is the field theory cut-off and k is the momentum scale. We also assume that the dimensionless couplings in the definitions (4.2) of F O,i and F λ are at most order one. Some of these conditions are only required to ensure the convergence of the heat kernel expansion (in particular, µ i ≫m, k), while the rest play some role in ensuring the robustness of the sequestering mechanism. For example, the scale λ is always chosen to lie at or above the cut-off M in order to avoid large contributions in the effective gravity equation that go as V vac /λ 4 with V vac ∼ M 4 . Although the choice of scales is important, our analysis in section 4 does seem to suggest that the precise structure of the monodromy deformation is less important from the point of view of a successful sequestering mechanism. Indeed, as long as the axion dependence contains suitable scalings, we anticipate the sequestering mechanism to survive a variety of deformations to the specific examples we have studied here. These may include (i) couplings between the axion sectors; (ii) higher powers of the axion fields (β j≥5 = 0 and c j≥2 = 0); and even (iii) a completely new type of coupling between σ i and higher curvature invariants (like σ i R 2 ). Although our analysis was not quite as general as this, it still provided a non-trivial challenge to the classical mechanism having included a consistent subset of loop-induced operators. These terms, which have not been discussed before in a classical context, necessarily arise from integrating out axion loops and we have shown that they leave sequestering intact.
When, then, do we expect the sequestering mechanism to fail? The answer to this is suggested by the arguments in favour of Omnia Sequestra [14,18] over the original proposals [9,10,12]. Problems occur when the form of the couplings themselves become too UV sensitive. For example, in [14], it was argued that loop corrections involving virtual gravitons alongside virtual matter fields would introduce UV sensitivity in the form of F λ . This introduced UV sensitivity into the global geometric constraint so that the residual cosmological constant present in the effective gravity equations itself became UV sensitive. In this paper, this problem never arose because we only considered pure axion loops, and there was no direct coupling between the axions and the matter fields. Of course, virtual gravitons can mediate such an interaction but they were not considered. A complete analysis including the impact of graviton loops will certainly reduce the space of monodromy models for which sequestering remains robust. This will be left for future work, but for now we speculate that possibly only the Omnia Sequestra model with a linear axion-Gauss-Bonnet coupling will survive, thanks to the exact shift symmetry such a theory would possess. If such a symmetry is indeed crucial, it would then be important to ask how it can emerge from a consistent quantum theory of gravity.

A Integrating out 3-forms and 4-forms
In deriving the action (2.2), we introduced the Lagrange multiplier fields Q i (x) and promoted the four-form field strengths to fields in their own right. The four-forms, only entering algabraically, and at quadratic order, were then integrated out, resulting in a mass term for σ i . At the level of the path integral this gave rise to an irrelevant overall factor To arrive at the global form of the action given in (2.5), we must also integrate out the three-form fields, A i . Focussing on the part of the relevant part of the action (2.2), we have wherec i βγδ = A i βγδ | ∂M corresponding to Dirichlet boundary conditions for the 3-form fields 7 . The 3-form field can be integrated out which yields a functional delta function, explicitly where the last expression is somewhat formal. We can give it a precise meaning by discretising the Q i (x) path integration in the next step. For simplicity we focus on the one dimensional case, which can be generalised to four dimensions. Dropping the index i for brevity, note that the path integral measure is defined via where each integration runs from −∞ to +∞ and Q a def = Q(x a ) is evaluated on a discrete grid with spacing ∆x = (x N − x 0 ) /N . For the functional delta we get Putting both expressions together we find that The functional integration therefore collapses to an ordinary integration over Q ≡ Q 0 = Q 1 = . . . = Q N corresponding to a constant Q(x) = Q. The factor in squared brackets does not depend on the fields and therefore is an irrelevant constant which drops out if we normalise our amplitudes. We have thus fully integrated out the 3-form fields from the original action (2.5). Going back to (A.3), we see that only the boundary term keeps a memory of their original presence due to their non-vanishing boundary valuesc i βγδ . Using the constancy of Q i (x), we can write where we defined C i = 1 3! ∂M d 3x n α ǫ αβγδci βγδ . Putting everything together, the final form of the sequestering action reads as (2.5), with Q i constant. These fields play the role of global variables reminiscent of the very first sequestering model [9,10]. Note that at the level of path integral we retain the ordinary integration, over these global variables, dQ i , as discussed at the beginning of section 3.