This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

UV and IR Effects in Axion Quality Control

C.P. Burgess,1,2 Gongjun Choi3 and F. Quevedo4
1 Department of Physics & Astronomy, McMaster University
  1280 Main Street West, Hamilton ON, Canada.
2 Perimeter Institute for Theoretical Physics
  31 Caroline Street North, Waterloo ON, Canada.
3 CERN, Theoretical Physics Department, Genève 23, Switzerland.
4 DAMTP, University of Cambridge, Wilberforce Road, Cambridge, CB3 0WA, UK.
Abstract

Motivated by recent discussions and the absence of exact global symmetries in UV completions of gravity we re-examine the axion quality problem (and naturalness issues more generally) using antisymmetric Kalb-Ramond (KR) fields rather than their pseudoscalar duals, as suggested by string and higher dimensional theories. Two types of axions can be identified: a model independent SS-type axion dual to a two form BμνB_{\mu\nu} in 4D and a TT-type axion coming directly as 4D scalar Kaluza-Klein (KK) components of higher-dimensional tensor fields. For TT-type axions our conclusions largely agree with earlier workers for the axion quality problem, but we also reconcile why TT-type axions can couple to matter localized on 3-branes with Planck suppressed strength even when the axion decay constants are of order the KK scale. For SS-type axions, we review the duality between form fields and massive scalars and show how duality impacts naturalness arguments about the UV sensitivity of the scalar potential. In particular UV contributions on the KR side suppress contributions on the scalar side by powers of m/Mm/M with mm the axion mass and MM the UV scale. We re-examine how the axion quality problem is formulated on the dual side and compare to recent treatments. We study how axion quality is affected by the ubiquity of pp-form gauge potentials (for both p=2p=2 and p=3p=3) in string vacua and identify two criteria that can potentially lead to a problem. We also show why most fields do not satisfy these criteria, but when they do the existence of multiple fields also provides mechanisms for resolving it. We conclude that the quality problem is easily evaded.

preprint: CERN-TH-2022-176

1 Introduction

String theory giveth and string theory taketh away, at least where axions111We follow the string literature and broadly refer to any low-energy Goldstone boson enjoying a rigid compact shift symmetry as an ‘axion’ (as opposed to the ‘dilatons’ associated with rigid scaling symmetries), something that would be called an ALP (axion-like particle) by particle phenomenologists. Our later focus is on those Goldstone bosons whose symmetries have a QCD anomaly and so can take part in the strong-CP problem StrongCP ; Weinberg:1977ma ; Wilczek:1977pj (which is what a particle physicist would usually mean by an ‘axion’). are concerned. On one hand axions are said to be ubiquitous in the spectrum of particles predicted around most string vacua StringUbiquity ; stringaxions . This observation motivates the study of their phenomenological consequences Axiverse , with a particular focus of late on their possible role as a light form of dark matter AxionReviews .

On the other hand, string theory equally generally forbids222Although there are known ways out Burgess:2008ri the conclusion is nonetheless broadly true and global symmetries tend to be both rare and approximate. the existence of exact rigid (or global) symmetries NoGlobal , in principle including the rigid shift symmetries on which low-energy axion properties are founded. For Goldstone bosons this breaking can keep them from being light, and can interfere with any mechanisms that rely on the survival of axions down to the low-energy theory. As applied to the QCD axion this has come to be known as the axion ‘quality’ problem QualityProblem .

So which is it? Are axions as abundant as dirt or as diamonds in low-energy string vacua? The resolution (which has long been known) is that there is a sense they are both. The absence of global symmetries really does mean that one never really directly finds scalar axions aa with shift symmetries in string vacua. Instead these scalars arise indirectly as Kaluza-Klein (KK) modes from fields that not themselves scalars; commonly arising3334D axions can also arise as KK modes from other types of extra-dimensional fields, but we focus on the Kalb-Ramond field because it allows a unified treatment of two different types of 4D axion. as components of 2-form Kalb-Ramond gauge fields KRGauge , B=12BMNdzMdzNB=\frac{1}{2}\,B_{{\scriptscriptstyle M\hbox{\kern-0.50003pt}N}}\,{\rm d}z^{\scriptscriptstyle M}\wedge{\rm d}z^{\scriptscriptstyle N}, subject to the gauge symmetries BdλB\to{\rm d}\lambda for some arbitrary field λM(x)\lambda_{\scriptscriptstyle M}(x). Fields like BMNB_{{\scriptscriptstyle M\hbox{\kern-0.50003pt}N}} arise so frequently in string vacua because they are related to other fields (notably the metric) by supersymmetry in higher dimensions.

1.1 Types of UV axion pedigree

Low-energy scalars typically emerge in the 4D effective theory from such fields in one of two ways:

  • TT-type axions: b(x)b(x) are specific cases of Kaluza-Klein (KK) modes arising when dimensionally reducing the extra-dimensional components Bmn(x,y)=b(x)ωmn(y)B_{mn}(x,y)=b(x)\,\omega_{mn}(y), where xμx^{\mu} denote the observed 4 dimensions, ymy^{m} are extra-dimensional coordinates and ωmn(y)\omega_{mn}(y) is a harmonic 2-form field within the extra dimensions.

  • SS-type axions: a(x){a}(x) arise directly as the 4-dimensional components Bμν(x,y)=bμν(x)ω(y)B_{\mu\nu}(x,y)=b_{\mu\nu}(x)\,\omega(y), which in four dimensions are known to be dual to scalar fields with shift symmetries Savit:1979ny through relations of the form μaϵμνλρνBλρ\partial^{\mu}{a}\propto\epsilon^{\mu\nu\lambda\rho}\partial_{\nu}B_{\lambda\rho} (much more about which below). Here ω(y)\omega(y) is a harmonic 0-form field – typically a yy-independent constant that can depend on extra-dimensional moduli.

This UV provenance is of course relevant to the axion quality problem, which is in essence an issue of UV sensitivity. One of our goals with this paper is to explore the ways that it helps, for both TT- and SS-type axions. Some of our conclusions are similar to earlier discussions of this issue Kallosh:1995hi ; DualStrongCP , in particular that the problem gets rephrased in dual form (for SS-type axions) in terms of the existence of multiple 3-form gauge potentials.

Since these issues have recently been revisited anew Sakhelashvili:2021eid ; Dvali:2022fdv we clarify what properties these fields must have to actually cause a quality problem and use this to argue why gravitational examples specifically (and the great abundance of such potentials in string vacua more generally) need not pose a problem in themselves. The dual formulation also suggests how the presence of multiple axions (as is common in string theory) can help alleviate the quality problem. The upshot is that the UV can, but need not, cause a quality problem. Whether or not it does cannot be decided purely at low energies because it depends on what happens in the UV.444The same is is also true of other naturalness problems; they arise because of strong dependence on physical masses for states that actually appear in the UV theory and not a dependence on cutoffs, as is sometimes mistakenly asserted (for a summary of these issues see e.g. Burgess:2013ara ).

But our discussion has implications that apply more broadly than just to the quality problem for the QCD axion. Along the way we identify more generally how dimensional ‘naturalness’ arguments for the scalar potential give very different estimates depending on whether they are made directly for the scalar or are first done for its dual and then mapped to the scalar using duality. In particular terms involving nn powers of the canonically normalized scalar arise additionally suppressed by powers of (m/M)n(m/M)^{n} where MM is the UV scale and mm is the axion mass (an observation also made in the past for inflationary models NaturalnessForms ).

We find a number of other ways that axion properties suggested by string-motivated extra-dimensional physics can be informative. For instance we describe a simple model for which TT-type axions have physical axion-matter couplings gaffg_{{a\hbox{\kern-0.70004pt}f\hbox{\kern-0.70004pt}f}} that are dramatically smaller than the naive value 1/f1/f read off from the axion kinetic term. In the example given (motivated by the models of YogaDE ) gaffg_{{a\hbox{\kern-0.70004pt}f\hbox{\kern-0.70004pt}f}} is order 1/Mp1/M_{p} despite ff being an ordinary particle-physics scale. Decoupling these scales from one another could have practical implications for axion phenomenology.

We show why the same hierarchy does not arise in these models for SS-type axions and we clarify why not. Physical couplings of SS-type axions really are of order 1/f1/f and we identify which interactions in the UV completion are responsible for the breakdown of the E/fE/f expansion at energies E>fE\mathrel{\raise 1.29167pt\hbox{$>$\kern-8.50006pt\lower 4.30554pt\hbox{$\sim$}}}f. SS-type axions illustrate how the scalar and dual representations can provide instances of weak/strong coupling duality, for which both the scalar and the dual cannot be within the weakly coupled regime. In the extra-dimensional example studied it is the Kalb-Ramond formulation that is weakly coupled. This could also have phenomenological implications to the extent that an axion that is dual to a weakly coupled system is unlikely to be well-described by the semiclassical methods that are universally used when exploring its physical implications.

Some of these observations imply that the use of the scalar (rather than Kalb-Ramond) variable can be misleading in some circumstances. This can seem surprising at first sight because the duality between axions and Kalb-Ramond fields is in essence a field redefinition and so scalar and dual formulations should be completely equivalent; it shouldn’t matter that string theory hands you Kalb-Ramond fields if scalar axions are equivalent and are much simpler to work with. Why should one care that a more complicated framework exists if it only obscures implications drawn using more transparent methods? We argue here that phenomena like weak/strong coupling duality are special cases of Weinberg’s Third Law of Progress in Theoretical Physics Weinberg:1981qq : You can use any degrees of freedom you like to describe a physical system, but if you use the wrong ones you’ll be sorry.

1.2 Non-propagating low-energy forms

These duality arguments touch on a related rich vein of physics with broader significance: the importance of keeping non-propagating entities like auxiliary and/or topological fields when formulating Wilsonian effective theories. These are fields that can be integrated out without changing the types of particles that propagate, and so it is tempting to think one should do so once and for all and simply ignore them thereafter. However such fields bring to the low-energy effective theory information about how its UV completion responds, e.g. to environments with nontrivial topology. They arise in concrete situations (such as in EFTs for 3-dimensional Quantum Hall systems, where the presence of emergent non-propagating gauge fields is essential for capturing the fractional quantization of Hall plateaux and the unusual charge and statistics of some excitations QHE ; EFTBook ).

Evidence is building that a similar role is played more widely by 3-form gauge potentials in four spacetime dimensions, C:=16CμνλdxμdxνdxλC:=\frac{1}{6}\,C_{\mu\nu\lambda}\,{\rm d}x^{\mu}{\rm d}x^{\nu}{\rm d}x^{\lambda} subject to the gauge freedom CC+dΛC\to C+{\rm d}\Lambda where Λμν(x)=Λνμ(x)\Lambda_{\mu\nu}(x)=-\Lambda_{\nu\mu}(x) is an arbitrary 2-form field. These are known to bring to the low-energy 4D effective theory topological information coming from integrated out extra dimensions Bousso:2000xa ; Burgess:2015lda , and more generally provide the origin for the auxiliary fields that appear in the 4D supergravities that are the low-energy limits of string vacua Bielleman:2015ina ; Herraez:2018vae . They appear in the QCD quality story because they can give masses to Kalb-Ramond fields Quevedo:1996uu through a Higgs mechanism that is dual to more mundane methods of axion mass generation. Because the field strength H=dCH={\rm d}C often appears in the action with a definite sign (often as a square), its presence can alter the implications of naturalness arguments for the scalar potential Burgess:2021juk . Indeed such terms provide the 4D understanding of why 6D SLED models SLED can in some circumstances suppress the 4D vacuum energy, but also why they struggle to do so enough to solve the cosmological constant problem Burgess:2015gba ; Burgess:2015lda ; Burgess:2015kda ; Niedermann:2015via . Their interplay with accidental scaling symmetries lies behind a recent attempt to find a dynamical relaxation mechanism for vacuum energies in four dimensions YogaDE .

In what follows we build our case for the above story using concrete examples. We first, in §2, briefly review the duality construction – in particular its extension to massive axions Quevedo:1996uu , which provides a way to think about scalar masses arising through a Higgs mechanism. §3 then briefly reviews and clarifies its use to dualize the axion solution to the strong-CP problem DualStrongCP , highlighting in particular how the quality problem gets rephrased in the dual language and how the apparent UV sensitivity of terms in the EFT differs between the axion formulation and its dual. Finally §4 provides a concrete extra-dimensional example – inspired by a UV completion of YogaDE – that illustrates both how axion/Kalb-Ramond duality can map weak to strong couplings, and how enormous hierarchies can arise with gaff,gaγγ1/Mpg_{{a\hbox{\kern-0.70004pt}f\hbox{\kern-0.70004pt}f}},g_{a\gamma\gamma}\sim 1/M_{p} even with ff as low as eV scales.

2 Axions and duality

We start with a review of why 2-form gauge potentials like BμνB_{\mu\nu} are dual Savit:1979ny ; Buscher:1987qj to scalar fields, both in the standard shift-symmetric massless case and for massive scalars, following the discussion of Quevedo:1996uu (which in turn generalizes earlier arguments JuliaToulouse aimed at describing particle/vortex duality in Kosterlitz-Thouless transitions KTTransitions ).

2.1 Axion/2-form duality

Consider the following path integral

Ξ[J]=𝒟BeiS1[B]\Xi[J]=\int{\cal D}B\;e^{iS_{1}[B]} (1)

where S1=d4x1S_{1}=\int{\rm d}^{4}x\;{\cal L}_{1} with 1{\cal L}_{1} chosen (at least to start) to be

1(B)=𝒵23!GμνλGμνλ13!ϵμνλρGμνλJρ,{\cal L}_{1}(B)=-\frac{{\cal Z}}{2\cdot 3!}\,G_{\mu\nu\lambda}G^{\mu\nu\lambda}-\frac{1}{3!}\,\epsilon^{\mu\nu\lambda\rho}G_{\mu\nu\lambda}J_{\rho}\,, (2)

with G=dBG={\rm d}B the exterior derivative of a 2-form field BμνB_{\mu\nu} and 𝒵{\cal Z} and JρJ_{\rho} possibly depending on other fields (collectively denoted ψ\psi). B=12BμνdxμdxνB=\frac{1}{2}\,B_{\mu\nu}\,{\rm d}x^{\mu}\wedge{\rm d}x^{\nu} is only defined up to the gauge redundancy BB+dλB\to B+{\rm d}\lambda for an arbitrary 1-form λ=λμdxμ\lambda=\lambda_{\mu}\,{\rm d}x^{\mu}.

The duality starts by trading the integration over BμνB_{\mu\nu} for an integral over GμνλG_{\mu\nu\lambda} subject to a constraint that imposes the Bianchi identity dG=0{\rm d}G=0. These are equivalent because the Bianchi identity is sufficient to guarantee the local existence of a field BμνB_{\mu\nu} with G=dBG={\rm d}B. The constraint is imposed by integrating over a scalar Lagrange-multiplier field a{a}, and so writing

Ξ[J]=𝒟G𝒟aeiS0\Xi[J]=\int{\cal D}G\,{\cal D}{a}\;e^{iS_{0}} (3)

where S0=d4x0S_{0}=\int{\rm d}^{4}x\;{\cal L}_{0} with

0(G,a)=𝒵23!GμνλGμνλ13!aϵμνλρμGνλρ13!ϵμνλρGμνλJρ,{\cal L}_{0}(G,{a})=-\frac{{\cal Z}}{2\cdot 3!}\,G_{\mu\nu\lambda}G^{\mu\nu\lambda}-\frac{1}{3!}{a}\,\epsilon^{\mu\nu\lambda\rho}\partial_{\mu}G_{\nu\lambda\rho}-\frac{1}{3!}\,\epsilon^{\mu\nu\lambda\rho}G_{\mu\nu\lambda}J_{\rho}\,, (4)

Integrating out a{a} imposes the Bianchi identity dG=0{\rm d}G=0 and allows the integral over GG to be replaced with the integral over BB, leading back to (2).

The dual version is obtained from (4) by instead integrating out GμνλG_{\mu\nu\lambda} so that a{a} is the remaining field. The result inherits a shift symmetry aa+{a}\to{a}\;+ constant because 0{\cal L}_{0} transforms into a total derivative. The GG integration is gaussian, whose saddle point is Gμνλ=𝒢μνλG_{\mu\nu\lambda}={\cal G}_{\mu\nu\lambda} where

𝒢μνλ=𝒵1ϵμνλρ(ρa+Jρ),{\cal G}_{\mu\nu\lambda}=-{\cal Z}^{-1}\epsilon_{\mu\nu\lambda\rho}\Bigl{(}\partial^{\rho}{a}+J^{\rho}\Bigr{)}\,, (5)

and so the integration gives the new lagrangian density

2(a)=12𝒵(μa+Jμ)(μa+Jμ).{\cal L}_{2}({a})=-\frac{1}{2{\cal Z}}(\partial_{\mu}{a}+J_{\mu})(\partial^{\mu}{a}+J^{\mu})\,. (6)

If 𝒵=1{\cal Z}=1 then a{a} is a canonically normalized massless scalar derivatively coupled to the same local current JμJ_{\mu} as in the original formulation. Because (2) and (6) are both obtained from (4) they must describe equivalent physics. Although the implied field redefinition from BμνB_{\mu\nu} to a{a} is in principle nonlocal the physics on both sides is nonetheless local because this is true of the relation between the field strengths given in (5).

Significance of 𝒵𝒵1{\cal Z}\leftrightarrow{\cal Z}^{-1}

In reality the above gaussian action is always supplemented by other non-gaussian interactions int{\cal L}_{\rm int} within a low-energy Wilsonian effective field theory (EFT). To the extent that both BμνB_{\mu\nu} and a{a} are derivatively coupled perturbative semiclassical methods in the presence of nongaussian terms like (GμνλGμνλ)2int(G_{\mu\nu\lambda}G^{\mu\nu\lambda})^{2}\in{\cal L}_{\rm int} are ultimately justified by a low-energy derivative expansion that applies equally well on both sides of a duality relationship because relationships like (5) involve equal numbers of derivatives on both sides.

The inversion of 𝒵𝒵1{\cal Z}\to{\cal Z}^{-1} as one passes from (2) to (6) is a noteworthy feature of duality. When 𝒵1{\cal Z}\gg 1 this implies 2-point correlators of GμνλG_{\mu\nu\lambda} are order 𝒵1{\cal Z}^{-1} in size while those of μa\partial_{\mu}{a} are instead order 𝒵{\cal Z}. The significance of the change 𝒵𝒵1{\cal Z}\to{\cal Z}^{-1} depends on whether or not BμνB_{\mu\nu} and a{a} can be freely rescaled to remove 𝒵{\cal Z} by going to canonically normalized variables. If this is so then 𝒵{\cal Z} in any case drops out of observables. For instance, when Jμ0J_{\mu}\neq 0 this rescaling shows that Ξ\Xi is really only a function of J~μ:=𝒵1/2Jμ\widetilde{J}_{\mu}:={\cal Z}^{-1/2}J_{\mu} rather than depending on 𝒵{\cal Z} and JμJ_{\mu} separately. Although 𝒵𝒵1{\cal Z}\leftrightarrow{\cal Z}^{-1} is sometimes called weak/strong coupling duality, Ξ[J~]\Xi[\widetilde{J}] is the same on both sides of the duality and so both sides agree on its functional dependence if expanded order-by-order in powers of 𝒵1{\cal Z}^{-1} (say).

One situation where this kind of rescaling is not possible is when 𝒵{\cal Z} depends on other fields and the target-space metric in field space is not flat. Another case where physics can depend explicitly on 𝒵{\cal Z} is when the field BμνB_{\mu\nu} or a{a} is quantized555This is generic the case in string theory for which the symmetries associated to antisymmetric tensors and axions are compact (meaning there always exist magnetic-like branes). For a general discussion see Banks:2010zn ., perhaps satisfying a boundary condition like Wdxμμa=2πnf\oint_{\scriptscriptstyle W}{\rm d}x^{\mu}\partial_{\mu}{a}=2\pi nf for some curve 𝒲{\cal W}, integer nn and mass scale ff, or perhaps CB=2πn~f1\oint_{\scriptscriptstyle C}B=2\pi\tilde{n}f^{-1} for some 2-cycle CC and possibly different integer n~\tilde{n} and mass scale f~\tilde{f}. In these situations physical results can depend on 𝒵{\cal Z} (i.e. on ff and/or f~\tilde{f}) and JμJ_{\mu} separately, and the relation 𝒵𝒵1{\cal Z}\to{\cal Z}^{-1} can carry physical significance.

2.2 A Higgs mechanism for scalar masses

Although the above makes the shift symmetry (and so also masslessness) of a{a} seem automatic, we next summarize how duality extends to massive scalars, following Quevedo:1996uu . A scalar potential is achieved in the dual framing through a Higgs mechanism in which the field BμνB_{\mu\nu} ‘eats’ (or is eaten by) a non-propagating gauge potential666Known string vacua can also contain a large number of these 3-form gauge potentials. CμνλC_{\mu\nu\lambda}. Because CμνλC_{\mu\nu\lambda} does not propagate this meal does not change the number of propagating degrees of freedom.

To this end consider the following gaussian path integral

Ξ[J]=𝒟C𝒟BeiS1\Xi[J]=\int{\cal D}C\,{\cal D}B\;e^{iS_{1}} (7)

where S1=d4x1S_{1}=\int{\rm d}^{4}x\;{\cal L}_{1} and

1(C,B)\displaystyle{\cal L}_{1}(C,B) =\displaystyle= 124!HμνλρHμνλρ123!(Gμνλ+mCμνλ)(Gμνλ+mCμνλ)\displaystyle-\frac{1}{2\cdot 4!}H_{\mu\nu\lambda\rho}H^{\mu\nu\lambda\rho}-\frac{1}{2\cdot 3!}(G_{\mu\nu\lambda}+mC_{\mu\nu\lambda})(G^{\mu\nu\lambda}+mC^{\mu\nu\lambda}) (8)
13!ϵμνλρ(Gμνλ+mCμνλ)Jρ.\displaystyle\qquad\qquad-\frac{1}{3!}\,\epsilon^{\mu\nu\lambda\rho}(G_{\mu\nu\lambda}+mC_{\mu\nu\lambda})J_{\rho}\,.

Here CμνλC_{\mu\nu\lambda} is a 3-form gauge potential with field strength H=dCH={\rm d}C and BμνB_{\mu\nu} is a 2-form gauge potential with G=dBG={\rm d}B while mm is a parameter with dimension mass.

This lagrangian has the gauge symmetry CC+dΛC\to C+{\rm d}\Lambda and BBmΛB\to B-m\,\Lambda for an arbitrary 2-form Λ\Lambda. So when m0m\neq 0 we can choose a gauge B=0B=0. The field equation for CC that follows from this action then is

DμHμνλρ+m2Cνλρ+mϵνλρμJμ=0.D_{\mu}H^{\mu\nu\lambda\rho}+m^{2}C^{\nu\lambda\rho}+m\,\epsilon^{\nu\lambda\rho\mu}J_{\mu}=0\,. (9)

This describes a single spin state propagating with mass mm once all the gauge symmetries are used, as can be seen by counting the massless states from which it is built. (In 4D BμνB_{\mu\nu} is shown above to be equivalent to a massless scalar and CμνλC_{\mu\nu\lambda} contains no propagating degrees of freedom at all because one can always write Hμνλρ=hϵμνλρH_{\mu\nu\lambda\rho}=h\,\epsilon_{\mu\nu\lambda\rho} with field equation μHμνλρ=0\partial_{\mu}H^{\mu\nu\lambda\rho}=0 in the massless limit, which implies hh is a constant and so does not propagate.)

The dual should therefore be a massive scalar and this can be verified by trading the integral over BB for an integral over GG and introducing (as before) a lagrange multiplier a{a} to impose the Bianchi identity777One can equivalently omit the mCμνλmC_{\mu\nu\lambda} terms everywhere and instead impose the modified Bianchi identity dG=mH{\rm d}G=mH. dG=0{\rm d}G=0, leading to the lagrangian density

0(C,G,a)\displaystyle{\cal L}_{0}(C,G,{a}) =\displaystyle= 124!HμνλρHμνλρ123!(Gμνλ+mCμνλ)(Gμνλ+mCμνλ)\displaystyle-\frac{1}{2\cdot 4!}H_{\mu\nu\lambda\rho}H^{\mu\nu\lambda\rho}-\frac{1}{2\cdot 3!}(G_{\mu\nu\lambda}+mC_{\mu\nu\lambda})(G^{\mu\nu\lambda}+mC^{\mu\nu\lambda}) (10)
13!aϵμνλρμGνλρ13!ϵμνλρ(Gμνλ+mCμνλ)Jρ.\displaystyle\qquad\qquad-\frac{1}{3!}\,{a}\,\epsilon^{\mu\nu\lambda\rho}\partial_{\mu}G_{\nu\lambda\rho}-\frac{1}{3!}\,\epsilon^{\mu\nu\lambda\rho}(G_{\mu\nu\lambda}+mC_{\mu\nu\lambda})J_{\rho}\,.

Integrating out a{a} returns us to the above formulation, but instead performing the integration over GG leads to the saddle point

𝒢μνλ=mCμνλϵμνλρ(ρa+Jρ),{\cal G}_{\mu\nu\lambda}=-mC_{\mu\nu\lambda}-\epsilon_{\mu\nu\lambda\rho}\Bigl{(}\partial^{\rho}{a}+J^{\rho}\Bigr{)}\,, (11)

and so to the lagrangian

2(C,a)=124!HμνλρHμνλρm4!aϵμνλρHμνλρ12μaμaJμμa12JμJμ.{\cal L}_{2}(C,{a})=-\frac{1}{2\cdot 4!}H_{\mu\nu\lambda\rho}H^{\mu\nu\lambda\rho}-\frac{m}{4!}\,{a}\,\epsilon^{\mu\nu\lambda\rho}H_{\mu\nu\lambda\rho}-\frac{1}{2}\partial_{\mu}{a}\,\partial^{\mu}{a}-J^{\mu}\partial_{\mu}{a}-\frac{1}{2}J_{\mu}J^{\mu}\,. (12)

Next we perform the integral over CμνλC_{\mu\nu\lambda}, and this is equivalent to simply performing the gaussian integral over HμνλρH_{\mu\nu\lambda\rho} because the integrability condition for writing H=dCH={\rm d}C is dH=0{\rm d}H=0 which is always true (in 4D). The saddle point for the HH integral occurs for Hμνλρ=μνλρH_{\mu\nu\lambda\rho}={\cal H}_{\mu\nu\lambda\rho} where

μνλρ=maϵμνλρ{\cal H}_{\mu\nu\lambda\rho}=-m\,{a}\,\epsilon_{\mu\nu\lambda\rho} (13)

and so leads to the scalar lagrangian

2(a)=12(a)2m22a2Jμμa12JμJμ.{\cal L}_{2}({a})=-\frac{1}{2}(\partial{a})^{2}-\frac{m^{2}}{2}{a}^{2}-J^{\mu}\partial_{\mu}{a}-\frac{1}{2}J_{\mu}J^{\mu}\,. (14)

This is the expected massive scalar.

2.2.1 Scalar potential

For future reference notice that it is only this last step that would differ if we’d had higher-dimension terms like δ=W(X)\delta{\cal L}=W(X) in the lagrangian with X=14!ϵμνλρHμνλρX=\frac{1}{4!}\epsilon^{\mu\nu\lambda\rho}H_{\mu\nu\lambda\rho} and so X2=14!HμνλρHμνλρX^{2}=-\frac{1}{4!}H_{\mu\nu\lambda\rho}H^{\mu\nu\lambda\rho} and so on. The above discussion is the special case W=12X2W=\frac{1}{2}X^{2} but one could entertain, for example,

W=c1M2X+12X2+2c33M2X3+c44M4X4+W=c_{1}M^{2}X+\frac{1}{2}\,X^{2}+\frac{2c_{3}}{3M^{2}}\,X^{3}+\frac{c_{4}}{4M^{4}}\,X^{4}+\cdots (15)

where the coefficients cic_{i} are dimensionless and MM is a UV scale inserted everywhere on dimensional grounds (with HμνλρH_{\mu\nu\lambda\rho} canonically normalized).

For non-quadratic WW the integral over HH is no longer gaussian, but we can proceed assuming a semiclassical saddle-point approximation is valid, in which case the saddle point (13) is modified to

(WX)H==ma,\left(\frac{\partial W}{\partial X}\right)_{H={\cal H}}=m\,{a}\,, (16)

which agrees with (13) when W=12X2W=\frac{1}{2}X^{2}. For example, for the choice (15) this becomes

c1M2+X(1+2c3M2X+c4M4X2+)mac_{1}M^{2}+X\left(1+\frac{2c_{3}}{M^{2}}X+\frac{c_{4}}{M^{4}}X^{2}+\cdots\right)\simeq m\,{a} (17)

and so

Xmac1M22c3M2(mac1M2)2+𝒪[(mac1M2)3/M4].X\simeq m\,{a}-c_{1}M^{2}-\frac{2c_{3}}{M^{2}}\left(m\,{a}-c_{1}M^{2}\right)^{2}+{\cal O}\left[\left(m\,{a}-c_{1}M^{2}\right)^{3}/M^{4}\right]\,. (18)

Once used in the lagrangian this shows how non-quadratic pieces of WW map over to non-quadratic contributions to the scalar potential for a{a} in the dual lagrangian 2{\cal L}_{2}. In particular the axion potential becomes

V(a)=W(X)+maX=12(mac1M2)22c33M2(mac1M2)3+.V({a})=-W(X)+m{a}X=\frac{1}{2}\left(m{a}-c_{1}M^{2}\right)^{2}-\frac{2c_{3}}{3M^{2}}\left(m{a}-c_{1}M^{2}\right)^{3}+\cdots\,. (19)

Two features are noteworthy about this potential:

  • First, notice it shares the usual Legendre property

    Va=mX+(WX+ma)Xa=mX,\frac{\partial V}{\partial{a}}=mX+\left(-\frac{\partial W}{\partial X}+m{a}\right)\frac{\partial X}{\partial{a}}=mX\,, (20)

    where the last equality uses (16). Even if new non-quadratic terms introduce new stationary points for V(a)V({a}) (or shifts the positions of old ones) eq. (20) ensures X=0X=0 for all of them.

  • Second, once a{a} is shifted so that the minimum is at a=0{a}=0 the potential depends on mm and a{a} only through the combination mam{a}. Consequently, a term proportional to an{a}^{n} comes suppressed by a power of (m/M)n(m/M)^{n} relative to what would naively be expected on dimensional grounds for V(a)V({a}). This is how the dual theory reproduces the same MM-dependence as found for higher powers of HμνλρH_{\mu\nu\lambda\rho} given that a{a} has canonical dimension mass while HH has dimension (mass)2. This shows how a dimensional assessment of how UV scales appear in the low-energy theory can care about the existence of a dual formulation.

3 Naturalness issues for dual systems

This section examines how naturalness arguments look for TT- and SS-type axions, and for SS-type axions how they depend on which side of the duality relation they are made. We do so using the axion quality problem as a representative example.

3.1 QCD and the dual PQ mechanism

To this end we extend the above reasoning to the main event: QCD and the θ\theta-term. The idea is to dualize the coupling of the axion to QCD to see how the strong-CP problem gets formulated, along the general lines of DualStrongCP . We then ask how UV physics might complicate the story in the dual theory. Consider then adding a gauge potential AμA_{\mu} (with field strength FμνF_{\mu\nu}) to represent the QCD gauge sector888We do not write quarks explicitly but flag the few places where their implicit presence affects what is written. and this time consider the path integral

Ξ[J]=𝒟G𝒟A𝒟aeiS0\Xi[J]=\int{\cal D}G\,{\cal D}A\,{\cal D}{a}\;e^{iS_{0}} (21)

where S0=d4x0S_{0}=\int{\rm d}^{4}x\;{\cal L}_{0} and

0(G,A,a)\displaystyle{\cal L}_{0}(G,A,{a}) =\displaystyle= 123!GμνλGμνλa3!ϵμνλρ(μGνλρ14Ωμνλρ)13!ϵμνλρGμνλJρ\displaystyle-\frac{1}{2\cdot 3!}G_{\mu\nu\lambda}G^{\mu\nu\lambda}-\frac{{a}}{3!}\,\epsilon^{\mu\nu\lambda\rho}\left(\partial_{\mu}G_{\nu\lambda\rho}-\frac{1}{4}\,\Omega_{\mu\nu\lambda\rho}\right)-\frac{1}{3!}\,\epsilon^{\mu\nu\lambda\rho}G_{\mu\nu\lambda}J_{\rho} (22)
14FμνFμνθ2ϵμνλρFμνFλρ.\displaystyle\qquad\qquad-\frac{1}{4}F_{\mu\nu}F^{\mu\nu}-\frac{\theta}{2}\,\epsilon^{\mu\nu\lambda\rho}F_{\mu\nu}F_{\lambda\rho}\,.

We suppress both gauge-group indices and traces over them to avoid notational clutter. FμνF_{\mu\nu} is the field strength for the gauge potential AμA_{\mu} but GμνλG_{\mu\nu\lambda} is an arbitrary 3-form until the integral over a{a} is performed.

Integrating out a{a} imposes the Bianchi identity dG=Ω{\rm d}G=\Omega where Ω\Omega is a gauge-invariant quantity built from the gauge field that on grounds of consistency must satisfy dΩ=0{\rm d}\Omega=0, for which we take

112ϵμνλρΩμνλρ=1fϵμνλρFμνFλρ\frac{1}{12}\,\epsilon^{\mu\nu\lambda\rho}\Omega_{\mu\nu\lambda\rho}=\frac{1}{f}\,\epsilon^{\mu\nu\lambda\rho}F_{\mu\nu}F_{\lambda\rho} (23)

The mass scale ff is here required on dimensional grounds. Doing this allows the GG integral to be traded for one over BB as before and gives the lagrangian

1(B,A)=123!GμνλGμνλ13!ϵμνλρGμνλJρ14FμνFμνθ2ϵμνλρFμνFλρ.{\cal L}_{1}(B,A)=-\frac{1}{2\cdot 3!}G_{\mu\nu\lambda}G^{\mu\nu\lambda}-\frac{1}{3!}\,\epsilon^{\mu\nu\lambda\rho}G_{\mu\nu\lambda}J_{\rho}-\frac{1}{4}F_{\mu\nu}F^{\mu\nu}-\frac{\theta}{2}\,\epsilon^{\mu\nu\lambda\rho}F_{\mu\nu}F_{\lambda\rho}\,. (24)

where G=dB+SG={\rm d}B+S where dΩ=0{\rm d}\Omega=0 implies there locally exists an SμνλS_{\mu\nu\lambda} – the Chern-Simons 3-form – that satisfies Ω=dS\Omega={\rm d}S.

The dual formulation instead integrates out GG and leaves a{a} as the dual variable. Integrating out GG leads to the lagrangian density

2(A,a)\displaystyle{\cal L}_{2}(A,{a}) =\displaystyle= 12(a)2Jμμa12JμJμ+a4!ϵμνλρΩμνλρ14FμνFμνθ2ϵμνλρFμνFλρ\displaystyle-\frac{1}{2}(\partial{a})^{2}-J^{\mu}\partial_{\mu}{a}-\frac{1}{2}J_{\mu}J^{\mu}+\frac{{a}}{4!}\,\epsilon^{\mu\nu\lambda\rho}\Omega_{\mu\nu\lambda\rho}-\frac{1}{4}F_{\mu\nu}F^{\mu\nu}-\frac{\theta}{2}\,\epsilon^{\mu\nu\lambda\rho}F_{\mu\nu}F_{\lambda\rho} (25)
=\displaystyle= 12(a)2Jμμa12JμJμ14FμνFμν+12(afθ)ϵμνλρFμνFλρ.\displaystyle-\frac{1}{2}(\partial{a})^{2}-J^{\mu}\partial_{\mu}{a}-\frac{1}{2}J_{\mu}J^{\mu}-\frac{1}{4}F_{\mu\nu}F^{\mu\nu}+\frac{1}{2}\left(\frac{{a}}{f}-\theta\right)\,\epsilon^{\mu\nu\lambda\rho}F_{\mu\nu}F_{\lambda\rho}\,.

This shows that the standard axion-gauge coupling is the dual of the 2-form/QCD coupling given in 1{\cal L}_{1} and that ff can be interpreted as its decay constant.

Below the QCD scale

In the standard axion-QCD story integrating out QCD leaves a residual axion potential due its anomalous coupling to FFF\wedge F. This minimum is argued to be minimized where a=θ¯f{a}={\overline{\theta}}\,f (where θ¯{\overline{\theta}} is the usual combination of θ\theta and phases in the quark mass matrices) which ensures that the CP-odd contribution turns off. We seek to express how physics below the QCD scale works in the dual language involving BμνB_{\mu\nu}.

Below ΛQCD\Lambda_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D} the gauge degrees of freedom are integrated out, naively leaving only hadrons coupled to BμνB_{\mu\nu}. The key thought is that this is not quite right: the QCD EFT below ΛQCD\Lambda_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D} contains a path integral over low-energy hadrons and an integration over a low-energy field CμνλC_{\mu\nu\lambda}, whose emergent presence the strongly coupled vacuum of QCD makes mandatory. The field CμνλSμνλC_{\mu\nu\lambda}\propto\langle S_{\mu\nu\lambda}\rangle is the low-energy counterpart of the Chern-Simons field appearing in the topological susceptibility Luscher:1978rn above the QCD scale, where dS=FF{\rm d}S=F\wedge F.

Having this field in the low-energy theory below the QCD scale does not affect the existence of a gap or the spectrum of the known hadrons because CμνλC_{\mu\nu\lambda} does not propagate. It is an auxiliary field that is required in order for the low-energy theory to capture properly the response of QCD to any topology in its environment. Similar fields are known to arise in this way in other concrete systems like the EFTs describing Quantum Hall systems QHE ; EFTBook . This 3-form potential differs from many of the others that often arise in string vacua because it arises from the IR properties of QCD rather than from the physics of UV compactification.

On dimensional grounds we write H=dCH={\rm d}C with

112Λ~QCD2ϵμνλρHμνλρ=ϵμνλρFμνFλρ,\frac{1}{12}\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2}\epsilon^{\mu\nu\lambda\rho}H_{\mu\nu\lambda\rho}=\epsilon^{\mu\nu\lambda\rho}\langle F_{\mu\nu}F_{\lambda\rho}\rangle\,, (26)

where Λ~QCD\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D} denotes a parameter of order the QCD scale that ensures that HH has canonical dimension (mass)2. The lagrangian (24) above the QCD scale is then replaced with its low-energy counterpart

1(C,B)=123!GμνλGμνλ13!ϵμνλρGμνλJρθ¯4!Λ~QCD2ϵμνλρHμνλρ124!HμνλρHμνλρ+,{\cal L}_{1}(C,B)=-\frac{1}{2\cdot 3!}G_{\mu\nu\lambda}G^{\mu\nu\lambda}-\frac{1}{3!}\,\epsilon^{\mu\nu\lambda\rho}G_{\mu\nu\lambda}J_{\rho}-\frac{{\overline{\theta}}}{4!}\,\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2}\epsilon^{\mu\nu\lambda\rho}H_{\mu\nu\lambda\rho}-\frac{1}{2\cdot 4!}H_{\mu\nu\lambda\rho}H^{\mu\nu\lambda\rho}+\cdots\,, (27)

where the explicit term proportional to θX\theta X combines with quark mass phases – that also enter as terms linear in XX, as in the c1c_{1} term of (15) – to produce θ¯X{\overline{\theta}}X. The ellipses in (27) are at least cubic in XX (or involve derivatives of XX).

Combining eq. (23) (and the discussion just above it) with (26) implies

dG=Ω=Λ~QCD2fH,{\rm d}G=\langle\Omega\rangle=\frac{\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2}}{f}\,H\,, (28)

and so comparing this to dG=mH{\rm d}G=mH (as would follow from G=dB+mCG={\rm d}B+mC) allows us to read off the mass relation m=Λ~QCD2/fm=\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2}/f. We see that the mCmC term captures the expectation value S/f\langle S\rangle/f of the Chern-Simons term in the UV theory above the QCD scale if mm scales with ff in the same way that the usual axion mass depends on its decay constant.

We expect the low-energy presence of such a 4-form field HH to give BB a nonzero mass, as we check by introducing the lagrange multiplier a{a} in the usual way and integrating out GG. This leads to the result

2(C,a)=12(a)2Jμμa12JμJμ+14!(maθ¯Λ~QCD2)ϵμνλρHμνλρ124!HμνλρHμνλρ+.{\cal L}_{2}(C,{a})=-\frac{1}{2}(\partial{a})^{2}-J^{\mu}\partial_{\mu}{a}-\frac{1}{2}J_{\mu}J^{\mu}+\frac{1}{4!}(m{a}-{\overline{\theta}}\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2})\epsilon^{\mu\nu\lambda\rho}H_{\mu\nu\lambda\rho}-\frac{1}{2\cdot 4!}H_{\mu\nu\lambda\rho}H^{\mu\nu\lambda\rho}+\cdots\,. (29)

Integrating out HH leads to the saddle point Hμνλρ=μνλρH_{\mu\nu\lambda\rho}={\cal H}_{\mu\nu\lambda\rho} with

μνλρ=(maθ¯Λ~QCD2)ϵμνλρ,{\cal H}_{\mu\nu\lambda\rho}=\left(m{a}-{\overline{\theta}}\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2}\right)\,\epsilon_{\mu\nu\lambda\rho}\,, (30)

and so gives the axion lagrangian

2(a)=12(a)2Jμμa12JμJμ12(maθ¯Λ~QCD2)2,{\cal L}_{2}({a})=-\frac{1}{2}(\partial{a})^{2}-J^{\mu}\partial_{\mu}{a}-\frac{1}{2}J_{\mu}J^{\mu}-\frac{1}{2}\left(m{a}-{\overline{\theta}}\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2}\right)^{2}\,, (31)

showing that the minimum indeed occurs where a=θ¯Λ~QCD2/m=θ¯f{a}={\overline{\theta}}\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2}/m={\overline{\theta}}f, which turns off the CP-violating term of (29).

In general integrating out the UV QCD sector also generates more complicated low-energy interactions involving CC, such as the function W(X)W(X) of X=14!ϵμνλρHμνλρX=\frac{1}{4!}\epsilon^{\mu\nu\lambda\rho}H_{\mu\nu\lambda\rho}. As above, such terms semiclassically change the saddle point to

(WX)H==maθ¯Λ~QCD2,\left(\frac{\partial W}{\partial X}\right)_{H={\cal H}}=m\,{a}-{\overline{\theta}}\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2}\,, (32)

and so leads to the axion potential

V(a)=W(X)+(maθ¯Λ~QCD2)X.V({a})=-W(X)+(m{a}-{\overline{\theta}}\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2})X\,. (33)

This satisfies

Va=mX+(WX+maθ¯Λ~QCD2)Xa=mX,\frac{\partial V}{\partial{a}}=mX+\left(-\frac{\partial W}{\partial X}+m{a}-{\overline{\theta}}\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2}\right)\frac{\partial X}{\partial{a}}=mX\,, (34)

and so again ensures that X=0X=0 at any of the stationary points of VV. We see that the presence of interactions like W(X)W(X) show that VV is minimized at ma=θ¯Λ~QCD2m{a}={\overline{\theta}}\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2} if W/X\partial W/\partial X vanishes when X=0X=0.

3.2 The Quality Problem

We now have the tools required to explore UV sensitivity and the axion quality problem. We start by restating the original formulaton of the quality problem and then how it is rephrased in 2-form language for both TT-type (this section) and SS-type (next section) axions.

The axion quality problem asks two related questions QualityProblem :

  1. 1.

    Do corrections to the QCD axion potential change its minimum in a way that preserves a sufficiently small effective vacuum angle: θ¯eff1010\bar{\theta}_{\rm eff}\lesssim 10^{-10}?

  2. 2.

    Do corrections to the QCD axion potential change the usual expression for the axion mass (that assumes it is dominantly generated by the ‘IR-dominated’ QCD instanton with size ρΛQCD1\rho\sim\Lambda_{{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}}^{-1})?

The first of these essentially asks if the QCD axion remains a good solution to the strong CP problem when perturbed by new physics, whereas the second asks the same of our understanding of axion mass. The axion mass question can apply more generally to ALPs as well, whereas the first one is specific to the QCD axion.

Any UV completion must decide what happens at energies above the axion decay constant ff above which the low-energy expansion in powers of E/fE/f breaks down. We consider in turn the original formulation and the TT- and SS-type axions that arise within an extra-dimensional context.

3.2.1 Original formulation

In the initial formulation the UV completion for scales above ff was assumed to involve a second scalar that combines with the axion to linearly realize the PQ symmetry as a complex scalar Φ\Phi. In this picture the modulus of Φ\Phi acquires a mass proportional to fΦf\sim\langle\Phi\rangle and the axion starts life as the phase of Φeia/f\Phi\propto e^{i{a}/f}.

Motivated by string theory and black-hole thought experiments it is then assumed that UV physics cannot support an unbroken global symmetry, and so at some large scale MM the form of the scalar potential for Φ\Phi cannot be assumed to be invariant under re-phasings of Φ\Phi. As an expansion in powers of Φ\Phi, the generic potential form would be

VUV(Φ)=M42n=1(cnΦnMn+h.c.),V_{\scriptscriptstyle U\hbox{\kern-0.50003pt}V}(\Phi)=\frac{M^{4}}{2}\sum_{n=1}^{\infty}\left(c_{n}\frac{\Phi^{n}}{M^{n}}+\hbox{h.c.}\right)\,, (35)

where the cnc_{n}’s are in general complex. This is true even if the UV physics is assumed to be CP-invariant because cnc_{n} will inherit the phase of the fermion mass matrix after chiral PQ rotations. In the initial formulation MM is assumed to be the Planck mass MpM_{p}, and although we can see that such a choice would dominate smaller MM for the terms with n<4n<4 it is likely that M<MpM<M_{p} would be more dangerous for n>4n>4. Early workers typically assumed that the renormalizable part of the potential would be tuned to make the axion potential sufficiently shallow and so effectively started the sum in (35) at n=5n=5.

Freezing the modulus field at Φ=f\langle\Phi\rangle=f and integrating it out at the classical level leads to the following effective axion potential

VUV(a)=M42n=1|cn|fnMn(eiδneina/f+h.c.)=M4n=5|cn|fnMncos(naf+δn),V_{\scriptscriptstyle U\hbox{\kern-0.50003pt}V}({a})=\frac{M^{4}}{2}\sum_{n=1}|c_{n}|\frac{f^{n}}{M^{n}}\left(e^{i\delta_{n}}\,e^{in{{a}}/{f}}+\hbox{h.c.}\right)=M^{4}\sum_{n=5}|c_{n}|\frac{f^{n}}{M^{n}}\cos\left(\frac{n{a}}{f}+\delta_{n}\right)\,, (36)

where we shift fields so that the standard QCD solution is a=0{a}=0. The QCD minimum therefore remains unchanged if VUV(0)=0V_{\scriptscriptstyle U\hbox{\kern-0.50003pt}V}^{\prime}(0)=0 and this would be true if all of the δn\delta_{n}’s were to vanish. Although the axion potential height (and therefore possibly axion mass) might still change due to the presence of VUV(a)V_{\scriptscriptstyle U\hbox{\kern-0.50003pt}V}({a}), evasion of the strong CP problem requires only that the minimum for a{a} remains unmoved.


Stability of the minimum:

For δn,|cn|𝒪(1)\delta_{n}\,,|c_{n}|\sim{\cal O}(1) we can estimate the size of the effective value of θ¯eff\bar{\theta}_{\rm eff} by perturbing around the QCD minimum at a=aQCD{a}={a}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}:

θ¯effVUV(aQCD)fVQCD′′(aQCD)VUV(aQCD)VQCD(aQCD)M4ΛQCD4(fM)n0,\bar{\theta}_{\rm eff}\simeq-\frac{V^{\prime}_{\scriptscriptstyle U\hbox{\kern-0.50003pt}V}({a}_{{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}})}{fV^{\prime\prime}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}({a}_{{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}})}\sim\frac{V_{\scriptscriptstyle U\hbox{\kern-0.50003pt}V}({a}_{{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}})}{V_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}({a}_{{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}})}\sim\frac{M^{4}}{\Lambda_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{4}}\left(\frac{f}{M}\right)^{n_{0}}\,, (37)

where n0n_{0} represents the first power appearing in the sum. For example, requiring θ¯eff<1010\bar{\theta}_{\rm eff}<10^{-10} for the example f=1012f=10^{12} GeV, M=Mp=1018M=M_{p}=10^{18} GeV and ΛQCD0.2\Lambda_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}\simeq 0.2 GeV in (37) requires n0>15n_{0}\mathrel{\raise 1.29167pt\hbox{$>$\kern-8.50006pt\lower 4.30554pt\hbox{$\sim$}}}15.


Stability of the axion mass:

The change to the axion mass induced by the UV axion potential is given by

δma2=2VUV(a)a2|a=0=M2n=1n2|cn|(fM)n2cosδn,\delta m_{a}^{2}=\left.\frac{\partial^{2}V_{\scriptscriptstyle U\hbox{\kern-0.50003pt}V}({a})}{\partial{a}^{2}}\right|_{{a}=0}=M^{2}\sum_{n=1}n^{2}|c_{n}|\left(\frac{f}{M}\right)^{n-2}\cos\delta_{n}\,, (38)

which can be significant unless the coefficients |cn||c_{n}|’s are extremely small even if all the δn\delta_{n}’s could be contrived to vanish. When significant such contributions spoil the relation mafmπFπm_{a}f\sim m_{\pi}F_{\pi} that holds for the low-energy QCD contribution and on which most axion phenomenology is based. Because the mass is not inversely proportional to ff this expression shows that the relation between mam_{a} and ff need not be inversely proportional to one other, for example allowing a very heavy axion to be still very weakly coupled to matter – a drastic change relative to standard axion phenomenology.

3.2.2 TT-type axions

The story is similar for TT-type axions, at least below the Kaluza-Klein scale where they are 4D scalars. No quality issue arises above the KK scale because here the relevant fields are higher-dimensional form fields HMNPH_{{\scriptscriptstyle M\hbox{\kern-0.50003pt}N\hbox{\kern-0.50003pt}P}} and the only symmetries involved are gauge symmetries like BB+dλB\to B+{\rm d}\lambda Kallosh:1995hi .

Recalling that TT-type axions, bb, arise as extra-dimensional reductions of the form Bmn(x,y)=b(x)ωmn(y)B_{mn}(x,y)=b(x)\,\omega_{mn}(y), with ωmn\omega_{mn} a harmonic form in the extra dimensions, the origin of the low-energy shift symmetry bb+cb\to b+c (for constant cc) has its origins as the extra-dimensional transformation BmnBmn+cωmnB_{mn}\to B_{mn}+c\,\omega_{mn}. This is a symmetry of H=dBH={\rm d}B because harmonic forms are closed: dω=0{\rm d}\omega=0. It is strictly speaking a ‘large’ gauge transformation because harmonic forms are not exact: there does not globally999The situation resembles a gauge field Am(x,y)A_{m}(x,y) dimensionally reduced on a circle, so Am(x,y+L)=Am(x,y)A_{m}(x,y+L)=A_{m}(x,y). In this case the massless scalar would be Am(x,y)=a(x)ω(y)A_{m}(x,y)=a(x)\omega(y) where ω(y)\omega(y) is independent of yy, for which the shift symmetry aa+ca\to a+c locally corresponds to a gauge transformation AmAm+mζA_{m}\to A_{m}+\partial_{m}\zeta if ζ/y=c\partial\zeta/\partial y=c, but this cannot be done globally because the solution cannot satisfy ζ(y+L)=ζ(y)\zeta(y+L)=\zeta(y). exist a λm\lambda_{m} such that ω=dλ\omega={\rm d}\lambda.

The quality problem arises because the shift symmetry in the low-energy 4D theory is not a local gauge symmetry and so it in principle need not be respected by UV effects. One consequently cannot completely preclude the generation of a scalar potential,

VUV(b)M4ncn(bM)n,V_{\scriptscriptstyle U\hbox{\kern-0.50003pt}V}(b)\sim M^{4}\sum_{n}c_{n}\left(\frac{b}{M}\right)^{n}\,, (39)

where cnc_{n} are dimensionless order-unity coefficients. But its failure to be a local gauge symmetry is a global obstruction rather than a local one and this means that UV effects cannot generate VUV(b)V_{\scriptscriptstyle U\hbox{\kern-0.50003pt}V}(b) until scales are integrated out that ‘see’ the topology that can distinguish ω\omega from dλ{\rm d}\lambda. This implies two sorts of changes to the standard quality-problem argument. First, the scale MM where problems first arise cannot be higher than the KK scale M1/RKKM\sim 1/R_{{\scriptscriptstyle K\hbox{\kern-0.50003pt}K}} corresponding to the size of the 2D cycle in the extra dimensions whose presence is associated with the existence of the harmonic form ωmn(y)\omega_{mn}(y). Second, the physics at scale MM that generates the potential must itself be sensitive to the nontrivial topology, often leading to additional suppressions.

For instance, an example of physics that can generate PQ-violating operators in (35) identified in Kallosh:1995hi ; Holman:1992ah is wormhole Giddings:1987cg . For these the coefficients cnc_{n} in (36) are exponentially suppressed, given by Kallosh:1995hi

cneSe(MpL)2c_{n}\sim e^{-S}\sim e^{-(M_{p}L)^{2}} (40)

where SS is a wormhole action and LL the size of its throat. Maintaining the success of the PQ mechanism requires S>190S\mathrel{\raise 1.29167pt\hbox{$>$\kern-8.50006pt\lower 4.30554pt\hbox{$\sim$}}}190. More complicated configurations are possible for extra-dimensional theories, for which MpM_{p} can be replaced by another UV gravity scale MgM_{g}, that might be the string scale or the extra-dimensional Planck scale MgM_{g} in specific examples. Similarly LL can be one of the geometric scales of the background, that could (but need not) be approximately a compactification scale RKKR_{{\scriptscriptstyle K\hbox{\kern-0.50003pt}K}}. All known semiclassical arguments of this type must assume MgL1M_{g}L\gg 1 for the calculation to be under control, because semiclassical methods are justified within an expansion in powers of (MgL)1(M_{g}L)^{-1} within any gravitational EFT. MgL14M_{g}L\sim 14 suffices to ensure S>190S\mathrel{\raise 1.29167pt\hbox{$>$\kern-8.50006pt\lower 4.30554pt\hbox{$\sim$}}}190 and so satisfying this constraint seems not that difficult within the semiclassical regime. These kinds of arguments were used in NaturalnessForms to argue for the absence of large gravitational correction to the inflaton potential.

3.3 The dual Quality Problem

For SS-type axions the representation directly obtained from UV physics is the field bμνb_{\mu\nu} dual to the scalar axion. And as alluded to earlier – c.f. §2.2.1 – issues of UV sensitivity can look very different in dual formulations to scalar theories, with for example the existence of a dual implying that the effective couplings for terms like anVUV(a){a}^{n}\in V_{\scriptscriptstyle U\hbox{\kern-0.50003pt}V}({a}) come suppressed by powers of the axion mass (m/M)n(m/M)^{n} relative to generic scalar estimates. Such suppressions can be enormous given the small size of mm relative to UV scales.

We therefore revisit earlier discussions of how the axion quality problem arises in the dual formulation, partly motivated by recent discussions Sakhelashvili:2021eid ; Dvali:2022fdv that argue that gravity causes new problems. Although we confirm the important role played by multiple 3-form potentials DualStrongCP in the framing of the dual quality problem, we also show that the many 3-forms found in string vacua do not generically pose a problem. Problems are only caused where strongly interacting systems make instanton-like effects important and this is not the case for the many ‘elementary’ 3-forms that descend from extra dimensional vacua. We argue that for similar reasons 4D gravitational Chern-Simons forms also need not cause problems (such as for string vacua where the UV completion of gravity is described by weakly coupled physics).

To the extent that the shape of the axion potential V(a)V({a}) is dual to interactions like W(X)W(X) involving the 4-form field strength X=14!HμνλρϵμνλρX=\frac{1}{4!}H_{\mu\nu\lambda\rho}\epsilon^{\mu\nu\lambda\rho}, one might think that the dual version of the axion quality issue should hinge on the detailed form of UV contributions to W(X)W(X). This proves not to be right, as we now argue. The central point turns on the Legendre transformation relating V(a)V({a}) to W(X)W(X); in particular on (32) and (34), that state

(WX)H==maθ¯Λ~QCD2andVa=mX.\left(\frac{\partial W}{\partial X}\right)_{H={\cal H}}=m\,{a}-{\overline{\theta}}\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2}\quad\hbox{and}\quad\frac{\partial V}{\partial{a}}=mX\,. (41)

On the scalar side the strong-CP problem is not solved unless ma=θ¯Λ~2m{a}={\overline{\theta}}\tilde{\Lambda}^{2} at the minimum of VV, and the quality problem is the statement that corrections to VV can perturb the minimum so that this relation fails. Although XX always vanishes at a minimum for VV, eq. (41) suggests that on the dual side the criterion for satisfying the strong-CP problem is that W/X=0\partial W/\partial X=0 is satisfied when X=0X=0. So the quality problem seems to hinge on whether or not UV physics can introduce a linear term δW=ηX\delta W=\eta X whose inclusion would modify (41) in a way that obstructs having ma=θ¯Λ~QCD2m\,{a}={\overline{\theta}}\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2} be a solution to V/a=0\partial V/\partial{a}=0.

Suppose, then, that one finds after integrating out the UV physics an EFT below the QCD scale of the form (27), but with a linear term in XX whose coefficient is not proportional to the CP violating parameter θ¯{\overline{\theta}}:

1(C,B)\displaystyle{\cal L}_{1}(C,B) =\displaystyle= 123!(Gμνλ+mCμνλ)(Gμνλ+mCμνλ)13!ϵμνλρ(Gμνλ+mCμνλ)Jρ\displaystyle-\frac{1}{2\cdot 3!}(G_{\mu\nu\lambda}+mC_{\mu\nu\lambda})(G^{\mu\nu\lambda}+mC_{\mu\nu\lambda})-\frac{1}{3!}\,\epsilon^{\mu\nu\lambda\rho}(G_{\mu\nu\lambda}+mC_{\mu\nu\lambda})J_{\rho} (42)
14!(θ¯+η)Λ~QCD2ϵμνλρHμνλρ124!HμνλρHμνλρ+,\displaystyle\qquad-\frac{1}{4!}({\overline{\theta}}+\eta)\,\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2}\epsilon^{\mu\nu\lambda\rho}H_{\mu\nu\lambda\rho}-\frac{1}{2\cdot 4!}H_{\mu\nu\lambda\rho}H^{\mu\nu\lambda\rho}+\cdots\,,

with two low-energy CP-violating parameters θ¯{\overline{\theta}} and η\eta. Dualizing this system as above then shows that scalar potential on the scalar side is given by a function of ma(θ¯+η)Λ~2m{a}-({\overline{\theta}}+\eta)\tilde{\Lambda}^{2}, in which θ¯{\overline{\theta}} and η\eta only appear as a sum. The arguments of §3.1 now show that this potential is minimized when ma(θ¯+η)Λ~2=0m{a}-({\overline{\theta}}+\eta)\tilde{\Lambda}^{2}=0. Repeating the calculation of the neutron electric dipole moment (edm) in this case – for a recent review, see for example Hook:2018dlk – then shows that the neutron edm also depends only on the sum θ¯+η{\overline{\theta}}+\eta and so would continue to vanish when a{a} is evaluated at the potential’s minimum. Interestingly, just introducing new terms linear in HμνλρH_{\mu\nu\lambda\rho} in (27) appears not to cause a quality problem.

3.3.1 A second strong sector

Just introducing a linear term in HμνλρH_{\mu\nu\lambda\rho} in (27) does not cause a quality problem because doing so below the QCD scale is like introducing the new CP-violating parameter η\eta only in the FFF\wedge F term of (24) above the QCD scale (i.e. shifting θθ+η\theta\to\theta+\eta). This also would not cause a quality problem on the scalar side. For there to be a problem requires there to be a CP-violating contribution to V(a)V({a}) that is independent of the CP-violation in the θ\theta-term.

What might this look like on the dual side? One way to proceed is to imagine a specific type of CP-violating UV completion and ask what happens in this case. One such an example would add another strongly interacting nonabelian gauge sector that also contributes to the axion anomaly. In this case VUV(a)V_{\scriptscriptstyle U\hbox{\kern-0.50003pt}V}({a}) is obtained by integrating out the new gauge sector and this is by construction independent of the QCD-generated part. A dual formulation of this type of system would involve a new Chern-Simons form EμνλE_{\mu\nu\lambda} for the new sector in addition to the QCD field CμνλC_{\mu\nu\lambda}, since both gauge sectors have their own Chern-Simons fields and either of these can be the field that is eaten by BμνB_{\mu\nu}. Instead of (42) below the QCD scale one would find the following low-energy action

1(C,E,B)\displaystyle{\cal L}_{1}(C,E,B) =\displaystyle= 123!GμνλGμνλ13!ϵμνλρGμνλJρ14!ϵμνλρ(θ¯Λ~QCD2Hμνλρ+ηΛ~X2Kμνλρ)\displaystyle-\frac{1}{2\cdot 3!}G_{\mu\nu\lambda}G^{\mu\nu\lambda}-\frac{1}{3!}\,\epsilon^{\mu\nu\lambda\rho}G_{\mu\nu\lambda}J_{\rho}-\frac{1}{4!}\epsilon^{\mu\nu\lambda\rho}\Bigl{(}{\overline{\theta}}\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2}H_{\mu\nu\lambda\rho}+\eta\tilde{\Lambda}_{\scriptscriptstyle X}^{2}K_{\mu\nu\lambda\rho}\Bigr{)} (43)
124!(HμνλρHμνλρ+KμνλρKμνλρ)+,\displaystyle\qquad\qquad-\frac{1}{2\cdot 4!}\Bigl{(}H_{\mu\nu\lambda\rho}H^{\mu\nu\lambda\rho}+K_{\mu\nu\lambda\rho}K^{\mu\nu\lambda\rho}\Bigr{)}+\cdots\,,

where K=dEK={\rm d}E and H=dCH={\rm d}C and G=dB+mC+m~EG={\rm d}B+mC+\tilde{m}E.

Proceeding as before we introduce a Lagrange multiplier a{a} to enforce the GG Bianchi identity and then semiclassically integrate out GG, HH and KK to find

2(a)=12(a)2Jμμa12JμJμV(a),{\cal L}_{2}({a})=-\frac{1}{2}(\partial{a})^{2}-J^{\mu}\partial_{\mu}{a}-\frac{1}{2}J_{\mu}J^{\mu}-V({a})\,, (44)

where defining X=14!ϵμνλρHμνλρX=\frac{1}{4!}\epsilon^{\mu\nu\lambda\rho}H_{\mu\nu\lambda\rho} and Y=14!ϵμνλρKμνλρY=\frac{1}{4!}\epsilon^{\mu\nu\lambda\rho}K_{\mu\nu\lambda\rho} we find

V(a)=W(X,Y)+(maθ¯Λ~QCD2)X+(m~aηΛ~X2)Y,V({a})=-W(X,Y)+(m{a}-{\overline{\theta}}\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2})X+(\tilde{m}{a}-\eta\tilde{\Lambda}_{\scriptscriptstyle X}^{2})Y\,, (45)

where W=12(X2+Y2)+W=\frac{1}{2}(X^{2}+Y^{2})+(higher powers). At the saddle point (H,K)=(,𝒦)(H,K)=({\cal H},{\cal K}) we have

(WX)Y=maθ¯Λ~QCD2and(WY)X=m~aηΛ~X2,\left(\frac{\partial W}{\partial X}\right)_{\scriptscriptstyle Y}=m\,{a}-{\overline{\theta}}\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2}\quad\hbox{and}\quad\left(\frac{\partial W}{\partial Y}\right)_{\scriptscriptstyle X}=\tilde{m}\,{a}-\eta\tilde{\Lambda}_{\scriptscriptstyle X}^{2}\,, (46)

where the subscripts indicate what is held fixed in the derivative. Differentiating (45) implies

Va=mX+m~Y.\frac{\partial V}{\partial{a}}=mX+\tilde{m}Y\,. (47)

This does have a quality problem because the competition between the two gauge sectors drives the axion away from the minimum for which the neutron electric dipole moment vanishes. For the simplest example – where W=12(X2+Y2)W=\frac{1}{2}(X^{2}+Y^{2}) – we can see explicitly how the shift of the global minimum of the axion potential is induced. From (47) we learn that V(a)/a=0\partial V({a})/\partial{a}=0 takes place at Y=(m/m~)XY=-(m/\tilde{m})X. From (46), we obtain

X=maθ¯Λ~QCD2andY=(mm~)X=m~aηΛ~X2X=m{a}-\bar{\theta}\tilde{\Lambda}_{QCD}^{2}\quad\hbox{and}\quad Y=-\left(\frac{m}{\tilde{m}}\right)X=\tilde{m}{a}-\eta\tilde{\Lambda}_{X}^{2} (48)

Equating these two expressions for XX and solving for a{a}, we obtain

amin=mθ¯Λ~QCD2+m~ηΛ~X2m2+m~2=aQCD+(m~ηΛ~X2/m2)1+(m~/m)2aQCD+m~ηΛ~X2m2,{a}_{\rm min}=\frac{m\bar{\theta}\tilde{\Lambda}_{{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}}^{2}+\tilde{m}\eta\tilde{\Lambda}_{{\scriptscriptstyle X}}^{2}}{m^{2}+\tilde{m}^{2}}=\frac{{a}_{{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}}+({\tilde{m}\eta\tilde{\Lambda}_{X}^{2}}/{m^{2}})}{1+\left({\tilde{m}}/{m}\right)^{2}}\simeq{a}_{{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}}+\frac{\tilde{m}\eta\tilde{\Lambda}_{X}^{2}}{m^{2}}\,, (49)

which denotes the global minimum before introducing an extra three form gauge field by aQCD=θ¯Λ~QCD2/m{a}_{{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}}=\bar{\theta}\tilde{\Lambda}_{QCD}^{2}/m. The approximate equality assumes mm~m\gg\tilde{m} so as not to spoil the QCD axion solution the strong CP problem.

Finally, defining the UV contribution to the effective vacuum angle by θ¯eff:=(aminaQCD)/f{\overline{\theta}}_{\rm eff}:=({a}_{\rm min}-{a}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D})/f where mfΛ~QCD2mf\simeq\tilde{\Lambda}_{QCD}^{2}, we obtain the constraint

θ¯effη(m~m)(Λ~XΛ~QCD)2<1010.\bar{\theta}_{\rm eff}\sim\eta\left(\frac{\tilde{m}}{m}\right)\left(\frac{\tilde{\Lambda}_{X}}{\tilde{\Lambda}_{QCD}}\right)^{2}\mathrel{\raise 1.29167pt\hbox{$<$\kern-8.50006pt\lower 4.30554pt\hbox{$\sim$}}}10^{-10}\,. (50)

Although this derivation assumed the simplest form W=12(X2+Y2)W=\frac{1}{2}(X^{2}+Y^{2}), the reasoning presented here can be applied to a more complicated W(X,Y)W(X,Y). In such a case (47) remains unchanged while (46) and (48) are modified. But amin{a}_{\rm min} remains connected to the value for (X,Y)(X,Y) that makes V/a\partial V/\partial{a} vanish via (46) and (47). Once amin{a}_{\rm min} is expressed in terms of aQCD{a}_{{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}}, one can always infer θ¯eff\bar{\theta}_{\rm eff} as above and impose the constraint θ¯eff<1010\bar{\theta}_{\rm eff}<10^{-10}.

The upshot is this: the requirement of multiple strongly coupled sectors on the dual side to generate a quality problem is much more explicit because the contribution of each sector is described by a separate 3-form potential, rather than having everything all be rolled into the same scalar potential.

3.3.2 Multiple fundamental 3-forms

At first sight the previous section makes it sound like string theory should typically have a huge quality problem, because of the generic appearance there of multiple 3-form potentials. We identify the circumstances under which these potentials could cause a quality problem and argue why such a problem generically does not happen. We also discuss how these criteria bear on a recent realization of these issues Sakhelashvili:2021eid ; Dvali:2022fdv .

To start consider how the EFT (43) above the QCD scale would be modified by the presence of many 3-form potentials 𝒞μνλA{\cal C}^{\scriptscriptstyle A}_{\mu\nu\lambda} (where A=1,,NA=1,\dots,N distinguishes the different UV potentials):

1(B,A,𝒞)\displaystyle{\cal L}_{1}(B,A,{\cal C}) =\displaystyle= 123!GμνλGμνλ13!ϵμνλρGμνλJρ14FμνFμνθ2ϵμνλρFμνFλρ\displaystyle-\frac{1}{2\cdot 3!}G_{\mu\nu\lambda}G^{\mu\nu\lambda}-\frac{1}{3!}\,\epsilon^{\mu\nu\lambda\rho}G_{\mu\nu\lambda}J_{\rho}-\frac{1}{4}F_{\mu\nu}F^{\mu\nu}-\frac{\theta}{2}\,\epsilon^{\mu\nu\lambda\rho}F_{\mu\nu}F_{\lambda\rho} (51)
14!ηAμνλρAϵμνλρ124!μνλρAAμνλρ+.\displaystyle\qquad\qquad-\frac{1}{4!}\,\eta_{\scriptscriptstyle A}{\cal H}^{\scriptscriptstyle A}_{\mu\nu\lambda\rho}\epsilon^{\mu\nu\lambda\rho}-\frac{1}{2\cdot 4!}{\cal H}^{\scriptscriptstyle A}_{\mu\nu\lambda\rho}{\cal H}_{\scriptscriptstyle A}^{\mu\nu\lambda\rho}+\cdots\,.

where A=d𝒞A{\cal H}^{\scriptscriptstyle A}={\rm d}{\cal C}^{\scriptscriptstyle A} and G=dB+SG={\rm d}B+S for the QCD Chern-Simons 3-form that satisfies Ω=dS\Omega={\rm d}S with Ω\Omega as given in (23). To the extent that none of the new fields μνλρA{\cal H}^{\scriptscriptstyle A}_{\mu\nu\lambda\rho} appear in the Bianchi identity dG=Ω{\rm d}G=\Omega they do not couple to QCD or to BμνB_{\mu\nu} and so play no role in the duality transformation from BμνB_{\mu\nu} to a{a}. One then arrives below the QCD scale with the lagrangian

1(𝒞,B)\displaystyle{\cal L}_{1}({\cal C},B) =\displaystyle= 123!GμνλGμνλ13!ϵμνλρGμνλJρθ¯4!Λ~QCD2ϵμνλρHμνλρ124!HμνλρHμνλρ\displaystyle-\frac{1}{2\cdot 3!}G_{\mu\nu\lambda}G^{\mu\nu\lambda}-\frac{1}{3!}\,\epsilon^{\mu\nu\lambda\rho}G_{\mu\nu\lambda}J_{\rho}-\frac{{\overline{\theta}}}{4!}\,\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2}\epsilon^{\mu\nu\lambda\rho}H_{\mu\nu\lambda\rho}-\frac{1}{2\cdot 4!}H_{\mu\nu\lambda\rho}H^{\mu\nu\lambda\rho} (52)
14!ηAμνλρAϵμνλρ124!μνλρAAμνλρ+.\displaystyle\qquad\qquad-\frac{1}{4!}\,\eta_{\scriptscriptstyle A}{\cal H}^{\scriptscriptstyle A}_{\mu\nu\lambda\rho}\epsilon^{\mu\nu\lambda\rho}-\frac{1}{2\cdot 4!}{\cal H}^{\scriptscriptstyle A}_{\mu\nu\lambda\rho}{\cal H}_{\scriptscriptstyle A}^{\mu\nu\lambda\rho}+\cdots\,.

Dualization proceeds as before, with the introduction of the scalar a{a} to enforce dG=Ω{\rm d}G=\Omega, and the saddle point in the integral over the 3-form potentials becomes

(WX)Y=maθ¯Λ~QCD2and(WYA)X=ηA,\left(\frac{\partial W}{\partial X}\right)_{{\scriptscriptstyle Y}}=m\,{a}-{\overline{\theta}}\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2}\quad\hbox{and}\quad\left(\frac{\partial W}{\partial Y^{\scriptscriptstyle A}}\right)_{\scriptscriptstyle X}=-\eta_{\scriptscriptstyle A}\,, (53)

where

W=12X2+12YAYA+(higher powers),W=\frac{1}{2}X^{2}+\frac{1}{2}Y^{\scriptscriptstyle A}Y_{\scriptscriptstyle A}+\hbox{(higher powers)}\,, (54)

and we define as before X=14!ϵμνλρHμνλρX=\frac{1}{4!}\epsilon^{\mu\nu\lambda\rho}H_{\mu\nu\lambda\rho} and YA=14!ϵμνλρμνλρAY^{\scriptscriptstyle A}=\frac{1}{4!}\epsilon^{\mu\nu\lambda\rho}{\cal H}^{\scriptscriptstyle A}_{\mu\nu\lambda\rho}. The dual lagrangian is

2(a)=12(a)2Jμμa12JμJμV(a),{\cal L}_{2}({a})=-\frac{1}{2}(\partial{a})^{2}-J^{\mu}\partial_{\mu}{a}-\frac{1}{2}J_{\mu}J^{\mu}-V({a})\,, (55)

where

V(a)\displaystyle V({a}) =\displaystyle= W(X,YA)+(maθ¯Λ~QCD2)XηAYA\displaystyle-W(X,Y^{\scriptscriptstyle A})+(m{a}-{\overline{\theta}}\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2})X-\eta_{\scriptscriptstyle A}Y^{\scriptscriptstyle A} (56)
=\displaystyle= 12X2+(maθ¯Λ~QCD2)X+12ηAηA,\displaystyle-\frac{1}{2}\,X^{2}+(m{a}-{\overline{\theta}}\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2})X+\frac{1}{2}\eta_{\scriptscriptstyle A}\eta^{\scriptscriptstyle A}\,,

and so

Va=mX.\frac{\partial V}{\partial{a}}=mX\,. (57)

We see that X=0X=0 in the vacuum and this implies from (53) and (54) that the strong-CP problem remains solved.

These arguments also show that two ingredients are required for additional 3-form potentials to cause a problem:

  1. 1.

    The additional 3-form potential 𝒞A0{\cal C}^{{\scriptscriptstyle A}_{0}} must contribute to the Bianchi identity for GG, and so κA00\kappa_{{\scriptscriptstyle A}_{0}}\neq 0 in the expression dG=Ω+κAA{\rm d}G=\Omega+\kappa_{\scriptscriptstyle A}{\cal H}^{\scriptscriptstyle A}, where A=d𝒞A{\cal H}^{\scriptscriptstyle A}={\rm d}{\cal C}^{\scriptscriptstyle A}; and

  2. 2.

    The additional 3-form potential must appear linearly in WW, so ηA00\eta_{{\scriptscriptstyle A}_{0}}\neq 0 in (52).

When both of these are satisfied then a{a} couples to A{\cal H}^{\scriptscriptstyle A} and leads to the competition of minima as in (46) along the lines described in §3.3.1. The need for both of these conditions to be true is why the bound (50) is proportional to both η\eta and Λ~X2\tilde{\Lambda}_{\scriptscriptstyle X}^{2}. The good news is that the vanishing of κA\kappa_{\scriptscriptstyle A} can be enforced by a gauge symmetry, since κA\kappa_{\scriptscriptstyle A} can only be nonzero if BB transforms as BBκAΛAB\to B-\kappa_{\scriptscriptstyle A}\Lambda^{\scriptscriptstyle A} under the 3-form gauge transformations CACA+dΛAC^{\scriptscriptstyle A}\to C^{\scriptscriptstyle A}+{\rm d}\Lambda^{\scriptscriptstyle A}.

There is at least one example of a 3-form potential which we know must exist and which also contributes to the Bianchi identity dG{\rm d}G: the gravitational Chern Simons 3-form, SgS_{g}. The existence of a PQ-Lorentz-Lorentz anomaly requires this form to appear in GG and so have a nonzero coefficient κg\kappa_{g} in the same way that the PQ-QCD-QCD anomaly requires the QCD Chern Simons form to appear there. Ref. DualStrongCP argues that this is real trouble whose evasion requires model-building, such as that done in Dvali:2022fdv .

Whether the existence of this form is a problem or not depends on whether it also satisfies item 2 above: i.e. whether or not it appears linearly in the lagrangian with coefficient ηg0\eta_{g}\neq 0. How big should ηg\eta_{g} be expected to be? Because any 4-form field strength =d𝒞{\cal H}={\rm d}{\cal C} is locally a total derivative it wants to drop out of perturbative physics when it appears linearly in the action (much as does FFF\wedge F). Consequently its appearance in a low-energy action requires some sort of nonperturbative process (like an instanton) to contribute to physical processes. This is indeed what happens for QCD for which the linear term in Ω\Omega appears with coefficient

Λ~QCD2M2e2πb/α\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2}\propto M^{2}\,e^{-2\pi b/\alpha} (58)

with MM a UV scale, bb a pure number and α=g2/4π\alpha=g^{2}/4\pi the QCD coupling. The tell-tale nonperturbative dependence on α\alpha is a semiclassical consequence of the topological character of FF\int F\wedge F and H\int H.

This suggests that for gravity a linear term in g{\cal H}^{g} should similarly be of size

ηM2e(ML)2\eta\propto M^{2}\,e^{-(ML)^{2}} (59)

for a characteristic instanton length scale LL and gravitational UV scale MM given that (ML)2(ML)^{-2} plays the role of the semiclassical expansion parameter (compare to (40)). This can be extremely small within the domain of validity of semiclassical reasoning, for which ML1ML\gg 1 (as would presumably apply when the UV completion is weakly coupled, such as for perturbative string vacua).

Examples of three forms characterized by η\eta in (59) include Eguchi-Hanson instantons Gibbons:1978tef ; Eguchi:1978gw and the gravitational Chern-Simons 3-form made up of gravitational connection. For the Chern-Simons 3-form ref. DualStrongCP argues that gravity indeed becomes strong in the UV, as would be required for η\eta to be significant. This could well be true, but the evidence for there being a problem hinges on how convinced one is about gravitational interactions becoming strong in the UV.

3.3.3 Multiple-axion solution

We close this section by remarking that having multiple axion candidates (as is often true for string vacua) can alleviate the above problem associated with multiple 3-form fields, even if the above two conditions are satisfied.101010The use of multiple axions to solve the quality problem is mentioned also in Heidenreich:2020 , who have different but related motivations for there being a plethora of form fields present in the UV. This observation points to an equally general quality control mechanism on the scalar side of the duality as well.

To see why, we introduce a second Kalb-Ramond field 𝔅μν{\mathfrak{B}}_{\mu\nu} to the model of §3.3.1, and supplementing the lagrangian of (43) with the appropriate additional kinetic term gives

1(C,E,B,)\displaystyle{\cal L}_{1}(C,E,B,{\cal B}) =\displaystyle= 123!𝔊μνλ𝔊μνλ123!GμνλGμνλ13!ϵμνλρGμνλJρ\displaystyle-\frac{1}{2\cdot 3!}{\mathfrak{G}}_{\mu\nu\lambda}{\mathfrak{G}}^{\mu\nu\lambda}-\frac{1}{2\cdot 3!}G_{\mu\nu\lambda}G^{\mu\nu\lambda}-\frac{1}{3!}\,\epsilon^{\mu\nu\lambda\rho}G_{\mu\nu\lambda}J_{\rho}
14!ϵμνλρ(θ¯Λ~QCD2Hμνλρ+ηΛ~X2Kμνλρ)124!(HμνλρHμνλρ+KμνλρKμνλρ)+,\displaystyle\qquad-\frac{1}{4!}\epsilon^{\mu\nu\lambda\rho}\Bigl{(}{\overline{\theta}}\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2}H_{\mu\nu\lambda\rho}+\eta\tilde{\Lambda}_{\scriptscriptstyle X}^{2}K_{\mu\nu\lambda\rho}\Bigr{)}-\frac{1}{2\cdot 4!}\Bigl{(}H_{\mu\nu\lambda\rho}H^{\mu\nu\lambda\rho}+K_{\mu\nu\lambda\rho}K^{\mu\nu\lambda\rho}\Bigr{)}+\cdots\,,

where as before K=dEK={\rm d}E and H=dCH={\rm d}C and G=dB+mC+m~EG={\rm d}B+mC+\tilde{m}E, but now also

𝔊:=d𝔅+mE.{\mathfrak{G}}:={\rm d}{\mathfrak{B}}+m_{\star}E\,. (61)

This system dualizes much as before: we introduce Lagrange multipliers a{a} and b{b} to enforce the GG and 𝔊{\mathfrak{G}} Bianchi identities dG=mH+m~K{\rm d}G=mH+\tilde{m}K and d𝔊=mK{\rm d}{\mathfrak{G}}=m_{\star}K and then integrate out GG, 𝔊{\mathfrak{G}}, HH and KK to find

2(a)=12(b)212(a+J)2V(a,b),{\cal L}_{2}({a})=-\frac{1}{2}(\partial{b})^{2}-\frac{1}{2}(\partial{a}+J)^{2}-V({a},{b})\,, (62)

with

V(a,b)=W(X,Y)+(maθ¯Λ~QCD2)X+(m~a+mbηΛ~X2)Y,V({a},{b})=-W(X,Y)+(m{a}-{\overline{\theta}}\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2})X+(\tilde{m}{a}+m_{\star}{b}-\eta\tilde{\Lambda}_{\scriptscriptstyle X}^{2})Y\,, (63)

and we define as before X=14!ϵμνλρHμνλρX=\frac{1}{4!}\epsilon^{\mu\nu\lambda\rho}H_{\mu\nu\lambda\rho} and Y=14!ϵμνλρKμνλρY=\frac{1}{4!}\epsilon^{\mu\nu\lambda\rho}K_{\mu\nu\lambda\rho}. For the simplest example of W=12(X2+Y2)W=\frac{1}{2}(X^{2}+Y^{2}), at the saddle point (H,K)=(,𝒦)(H,K)=({\cal H},{\cal K}) gives the following relation between (X,Y)(X,Y) and (a,b)({a},{b}):

(WX)Y=maθ¯Λ~QCD2and(WY)X=m~a+mbηΛ~X2.\left(\frac{\partial W}{\partial X}\right)_{\scriptscriptstyle Y}=m\,{a}-{\overline{\theta}}\tilde{\Lambda}_{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}^{2}\quad\hbox{and}\quad\left(\frac{\partial W}{\partial Y}\right)_{\scriptscriptstyle X}=\tilde{m}\,{a}+m_{\star}\,{b}-\eta\tilde{\Lambda}_{\scriptscriptstyle X}^{2}\,. (64)

Differentiating (63) with respect to a{a} and b{b} implies

Va=mX+m~Y,Vb=mY.\frac{\partial V}{\partial{a}}=mX+\tilde{m}Y,\qquad\frac{\partial V}{\partial{b}}=m_{\star}Y\,. (65)

and so shows that all extrema of the potential satisfy X=Y=0X=Y=0 (provided m,m~m,\tilde{m} and mm_{\star} are nonzero). Because W/X\partial W/\partial X vanishes at X=0X=0 it follows that the dynamics chooses amin{a}_{\rm min} to satisfy θ¯Λ~QCD2/m=θ¯f{\overline{\theta}}\tilde{\Lambda}_{{\scriptscriptstyle Q\hbox{\kern-0.50003pt}C\hbox{\kern-0.50003pt}D}}^{2}/m={\overline{\theta}}f through (64); the axion quality problem essentially disappears.

What happened? Why does introducing another axion resolve the quality problem? The crux of the mechanism lies in the difference between eq. (65) and (47). The derivative of the potential always sets a linear combination of 4-form field strengths to zero and if there are as many equations as there are fields the only solution is generically to have all 4-form field strengths vanish. Once this is true then the first of eqs. (64) ensures that this solution solves the strong-CP problem. Trouble only arises – as it did in §3.3.1 – when there are fewer equations than unknowns (i.e. fewer axions than 3-form potentials), since then XX need not vanish and eqs. (64) become competing conditions on the same axion variable.

A similar mechanism also exists on the scalar side of the duality. If two sectors generate contributions to the QCD axion potential then the problem arises because these compete in the value they imply for the axion expectation value. Introducing a second anomalous U(1)U(1) symmetry that also has anomalies with the same two sectors provides enough latitude to minimize each sector’s potential separately, thereby removing the troublesome competition.

For instance, suppose there was a new non-Abelian gauge sector 𝒢\mathcal{G} and suppose the usual PQ symmetry has both a QCD anomaly and an anomaly in the 𝒢{\cal G} sector. This is the kind of thing that can cause a quality problem because of the contradictory conditions the two sectors impose on the QCD axion. But also introducing another global U(1)U(1) with only a 𝒢{\cal G}-sector anomaly can help because there is a linear combination of the PQ symmetry and the new U(1)U(1) that is anomaly free in the 𝒢{\cal G} sector and the PQ mechanism then goes through using this new symmetry.111111Ref. Dvali:2022fdv uses a special case of this general mechanism by introducing an extra U(1)U(1) symmetry in the leptonic sector to resolve the problem raised by the assumption that gravity is strongly coupled.

4 UV completion and matter couplings

Since the motivations both for considering Kalb-Ramond fields and for the absence of global symmetries come from the UV, it is useful to ask whether there are other potential surprises for axion physics having their roots in the UV. This section examines two such examples; one each for TT-type and for SS-type axions. For TT-type axions we provide simple examples for which physical axion-matter couplings like gaffg_{{a\hbox{\kern-0.70004pt}f\hbox{\kern-0.70004pt}f}} can be much smaller than the naive value 1/fb1/f_{b} read off from the axion kinetic term. In the example shown here gauge invariant matter couplings like gaffg_{{a\hbox{\kern-0.70004pt}f\hbox{\kern-0.70004pt}f}} are order 1/Mp1/M_{p} despite fbf_{b} being an ordinary particle physics scale, while anomalous gauge couplings remain order 1/fb1/f_{b} in size (if they exist at all).

For SS-type axions we show that the corresponding physical couplings indeed are of order 1/fa1/f_{a} and we identify the UV physics to which couplings of size E/faE/f_{a} match at energies E>faE\mathrel{\raise 1.29167pt\hbox{$>$\kern-8.50006pt\lower 4.30554pt\hbox{$\sim$}}}f_{a}. We also show how SS-type axions can be examples of weak/strong duality, and that it is the Kalb-Ramond side of the duality that is usually weakly coupled. Weak/strong coupling interchange due to duality could be relevant to applications for which the effects of scalar axions are explored using semiclassical reasoning, and if so would provide a further motivation for taking the Kalb-Ramond formulation as primary.

4.1 Extra-dimensional UV completion

To this end suppose that both Kalb-Ramond field and the standard model arise within an extra-dimensional model. For concreteness’ sake we take the higher-dimensional kinetic term for the 2-form field and the Einstein-Hilbert part of the action to be121212For simplicity we ignore extra-dimensional warping in this discussion. We also do not canonically normalize BMNB_{{\scriptscriptstyle M\hbox{\kern-0.50003pt}N}}, which here is taken to be dimensionless.

Skin=12M2+dd4xddyg~(D)(~+13!eλϕGMNPGMNP),S_{\rm kin}=-\frac{1}{2}\,M^{2+d}\int{\rm d}^{4}x\,{\rm d}^{d}y\;\sqrt{-\tilde{g}_{({\scriptscriptstyle D})}}\left(\widetilde{\cal R}+\frac{1}{3!}\,e^{-\lambda\phi}\,G_{{\scriptscriptstyle M\hbox{\kern-0.50003pt}N\hbox{\kern-0.50003pt}P}}G^{{\scriptscriptstyle M\hbox{\kern-0.50003pt}N\hbox{\kern-0.50003pt}P}}\right)\,, (66)

where there are D=4+dD=4+d spacetime dimensions and MM is a UV scale – the higher-dimensional Planck scale. ~\widetilde{\cal R} here denotes the Ricci scalar and g~(D)\tilde{g}_{({\scriptscriptstyle D})} is the determinant of the full DD-dimensional metric g~MN\tilde{g}_{{\scriptscriptstyle M\hbox{\kern-0.50003pt}N}}. As above H=dB+H={\rm d}B+\cdots is the Kalb-Ramond field strength and ϕ\phi is the extra-dimensional dilaton that often arises within the higher-dimensional gravity supermultiplet. The parameter λ\lambda depends on higher-dimensional details, with (d,λ)=(2,2)(d,\lambda)=(2,2) for chiral 6D supergravity 6Dsugra , (d,λ)=(6,1)(d,\lambda)=(6,1) for Neveu-Schwarz 2-forms in 10D supergravity and (d,λ=1)(d,\lambda=-1) for Ramond 2-forms in 10D supergravity 10Dsugra (for example).

The derivation of this type of lagrangian as the low-energy limit of a string vacuum usually relies on two approximations: the low-energy approximation (or α\alpha^{\prime} expansion) where energies are well below the string scale EMsE\ll M_{s}; and the weak string coupling approximation, which involves expanding in powers of eϕ1e^{\phi}\ll 1. For simplicity we restrict ourselves to this limit as well, and specialize to the simplest case (d,λ)=(2,2)(d,\lambda)=(2,2) corresponding to 6D chiral supergravity.

Dimensional reduction to 4D proceeds by integrating out the two extra dimensions and putting the 4D Einstein-Hilbert term into standard form (4D Einstein frame) by appropriately rescaling the 4D part of the metric

g~μν=(𝒱2𝒱2)gμν=1𝒱2(Mp2M2)gμνwhere𝒱d:=Mdddyg~(d)\tilde{g}_{\mu\nu}=\left(\frac{{\cal V}_{2\star}}{{\cal V}_{2}}\right)g_{\mu\nu}=\frac{1}{{\cal V}_{2}}\left(\frac{M_{p}^{2}}{M^{2}}\right)g_{\mu\nu}\quad\hbox{where}\quad{\cal V}_{d}:=M^{d}\int{\rm d}^{d}y\sqrt{\tilde{g}_{(d)}} (67)

is the dimensionless extra-dimensional volume and the subscript ‘\star’ on a field denotes its present-day value131313The factor of 𝒱2{\cal V}_{2\star} ensures the rescaling is trivial at present, as required to not change present-day units of length. and the 4D Planck massis is defined by

Mp2=𝒱2M2.M_{p}^{2}={\cal V}_{2\star}\,M^{2}\,. (68)

SS-type axion

The kinetic term for bμνb_{\mu\nu} in 4D Einstein frame that is obtained by dimensional reduction is

kin=112M2𝒱2g~(4)e2ϕg~μνg~βρg~ξζHμβξHνρζ=M412Mp2ge2ϕ𝒱22hμνβhμνβ,{\cal L}_{\rm kin}=-\frac{1}{12}\,M^{2}{\cal V}_{2}\sqrt{-\tilde{g}_{(4)}}\;e^{-2\phi}\tilde{g}^{\mu\nu}\tilde{g}^{\beta\rho}\tilde{g}^{\xi\zeta}\,H_{\mu\beta\xi}H_{\nu\rho\zeta}=-\frac{M^{4}}{12M_{p}^{2}}\sqrt{-g}\;e^{-2\phi}\,{\cal V}_{2}^{2}\,h_{\mu\nu\beta}h^{\mu\nu\beta}\,, (69)

where hμνλ=μbνλh_{\mu\nu\lambda}=\partial_{\mu}b_{\nu\lambda} + (cyclic). This last form can be written in terms of a scalar by dualizing as in earlier sections, imposing the Bianchi identity141414The factors of MM here are chosen so that Ω\Omega has dimension (mass)4. dh=Ω/M2{\rm d}h=\Omega/M^{2}, leading to the dual result

dual=g[Mp2e2ϕ𝒱22μ𝔞μ𝔞+13!𝔞ϵμνβρΩμνβρ]{\cal L}_{\rm dual}=-\sqrt{-g}\left[\frac{M_{p}^{2}e^{2\phi}}{{\cal V}_{2}^{2}}\,\partial_{\mu}{\mathfrak{a}}\,\partial^{\mu}{\mathfrak{a}}+\frac{1}{3!}\,{\mathfrak{a}}\,\epsilon^{\mu\nu\beta\rho}\Omega_{\mu\nu\beta\rho}\right] (70)

which suggests its decay constant can be written fa=(Mp/𝒱d)eϕf_{a}=(M_{p}/{\cal V}_{d\star})\,e^{\phi_{\star}}.

Two things are noteworthy here. First, notice that the volume dependence means that faf_{a} can be very much smaller than Planckian size. In the extreme case of two large extra dimensions (and working in the weak-coupling regime for which eϕe^{\phi} is moderately small) the size of the extra dimensions can be as large as MRKK<1014MR_{\scriptscriptstyle K\hbox{\kern-0.50003pt}K}\mathrel{\raise 1.29167pt\hbox{$<$\kern-8.50006pt\lower 4.30554pt\hbox{$\sim$}}}10^{14} and so 𝒱2(MRKK)2<1028{\cal V}_{2}\sim(MR_{\scriptscriptstyle K\hbox{\kern-0.50003pt}K})^{2}\mathrel{\raise 1.29167pt\hbox{$<$\kern-8.50006pt\lower 4.30554pt\hbox{$\sim$}}}10^{28} can be enormous (potentially allowing faMpf_{a}\lll M_{p} to be as small as eV energies).

Second, notice that although (69) has large coefficients when e2ϕ1e^{2\phi}\ll 1, the same is not true of the kinetic term in (70). This reflects how Kalb-Ramond/axion duality is a weak-strong coupling duality from the point of view of the string coupling gseϕg_{s}\sim e^{\phi}. To the extent that semiclassical expansions rely on the leading action being proportional to the inverse of a small coupling151515When this is true then powers of gs2g_{s}^{2} and powers of \hbar are equivalent when evaluating a path integral over eiS0e^{iS_{0}}.0=L0/gs2{\cal L}_{0}=L_{0}/g_{s}^{2} – semiclassical methods should fail for the scalar representation but hold for its dual.

TT-type axion

For TT-type axions we use Bmn(x,y)=𝔟(x)ωmn(y)B_{mn}(x,y)={\mathfrak{b}}(x)\,\omega_{mn}(y) where in six dimensions the harmonic form can be taken to be proportional to the extra-dimensional volume form εmn(y)\varepsilon_{mn}(y). Typically ωmn\omega_{mn} satisfies a quantization condition that states the integral of ωmn\omega_{mn} over the two extra dimensions Cω\oint_{\scriptscriptstyle C}\omega is a pure number, proportional to an integer. Because this result is volume independent it follows that ωmn=𝒱21εmn\omega_{mn}={\cal V}_{2}^{-1}\,\varepsilon_{mn}.

The kinetic term for the TT-type scalar 𝔟{\mathfrak{b}} obtained in this way is therefore proportional to

kin=12M2𝒱2g~(4)e2ϕg~μνμ𝔟ν𝔟𝒱22=gMp2e2ϕ𝒱22gμνμ𝔟ν𝔟.{\cal L}_{\rm kin}=-\frac{1}{2}\,M^{2}{\cal V}_{2}\sqrt{-\tilde{g}_{(4)}}\;e^{-2\phi}\tilde{g}^{\mu\nu}\partial_{\mu}{\mathfrak{b}}\,\partial_{\nu}{\mathfrak{b}}\,{\cal V}_{2}^{-2}=-\sqrt{-g}\;M_{p}^{2}e^{-2\phi}{\cal V}_{2}^{-2}g^{\mu\nu}\partial_{\mu}{\mathfrak{b}}\,\partial_{\nu}{\mathfrak{b}}\,. (71)

Notice that the kinetic term, both here and in (70), takes the form

kin=12gMp2[(𝔟)2τ2+(𝔞)2σ2]{\cal L}_{\rm kin}=-\frac{1}{2}\sqrt{-g}\,M_{p}^{2}\;\left[\frac{(\partial{\mathfrak{b}})^{2}}{\tau^{2}}+\frac{(\partial{\mathfrak{a}})^{2}}{\sigma^{2}}\right] (72)

with τ=𝒱2eϕ\tau={\cal V}_{2}\,e^{\phi} in (71), and σ=𝒱2eϕ\sigma={\cal V}_{2}\,e^{-\phi} in (70).

4.2 Coupling strengths

What matters for phenomenology is the couplings of the fields 𝔟{\mathfrak{b}} and 𝔞{\mathfrak{a}} to matter. This is controlled by the size of {\cal F} for axion couplings of the form

ax=12μaμa1μaJμ,{\cal L}_{\rm ax}=-\frac{1}{2}\,\partial_{\mu}{a}\,\partial^{\mu}{a}-\frac{1}{{\cal F}}\,\partial_{\mu}{a}J^{\mu}\,, (73)

where a{a} is the canonically normalized axion field and JμJ^{\mu} is a matter current. 1=gaff{\cal F}^{-1}=g_{{a\hbox{\kern-0.70004pt}f\hbox{\kern-0.70004pt}f}} is the axion-fermion current if JμJ^{\mu} is built from fermion bilinears and 1gagg{\cal F}^{-1}\simeq g_{agg} or gaγγg_{a\gamma\gamma} if JμJ^{\mu} is the Hodge dual of the QCD or QED Chern-Simons 3-form.

For concreteness’ sake we evaluate the size of this coupling in the perturbative semiclassical regime where 𝒱2{\cal V}_{2} is large and the UV physics is weakly coupled (and so eϕe^{\phi} small). In this limit we have fb>Mp/𝒱2>faf_{b}\mathrel{\raise 1.29167pt\hbox{$>$\kern-8.50006pt\lower 4.30554pt\hbox{$\sim$}}}M_{p}/{\cal V}_{2}\mathrel{\raise 1.29167pt\hbox{$>$\kern-8.50006pt\lower 4.30554pt\hbox{$\sim$}}}f_{a} and both are much smaller than MpM_{p}. In both cases we will see that {\cal F} can (but need not) be simply given by the corresponding decay constant faf_{a} or fbf_{b}.

In higher dimensional constructions very often ordinary matter is localized on a space-filling brane, Σ\Sigma, within the extra dimensions. Σ\Sigma could be a four-dimensional 3-brane or a higher-dimensional pp-brane with 3p3+d3\leq p\leq 3+d. If p>3p>3 then the extra-dimensional part of the brane typically wraps some topological cycle within the extra dimensions, and if this were a two-cycle (e.g. if p=5p=5) it would also have an associated harmonic 2-form ωmn(y)\omega_{mn}(y) required to ensure that TT-type axions appear in the low-energy 4D theory. We here explore the simplest case p=3p=3.

SS-type axion

A generally covariant low-dimension interaction between HMNPH_{{\scriptscriptstyle M\hbox{\kern-0.50003pt}N\hbox{\kern-0.50003pt}P}} and matter fields living on the brane, that is linear in BMNB_{{\scriptscriptstyle M\hbox{\kern-0.50003pt}N}} is161616The first equality shows that this interaction is independent of the metric, and this can also be seen after the second equality from the observation that ϵμνλρ=±g\epsilon_{\mu\nu\lambda\rho}=\pm\sqrt{-g} and so ϵμνλρ=±(g)1/2\epsilon^{\mu\nu\lambda\rho}=\pm(-g)^{-1/2}.

Sint=c^ΣeβϕHJ=c3!d4xgeβϕϵμνλρhμνλJρ,S_{\rm int}=-\hat{c}\int_{\Sigma}e^{-\beta\phi}H\wedge J=-\frac{c}{3!}\int{\rm d}^{4}x\sqrt{-g}\;e^{-\beta\phi}\,\epsilon^{\mu\nu\lambda\rho}h_{\mu\nu\lambda}J_{\rho}\,, (74)

where JρJ_{\rho} is a current built from brane-localized matter fields and β\beta is a parameter – like λ\lambda in (66) – that is predicted by any specific extra-dimensional UV completion. The matter current JρJ_{\rho} has dimension (mass)3 – making the coupling parameters c^\hat{c} and cc dimensionless – and so could be a fermion bilinear or the Hodge dual of a gauge boson Chern-Simons term (though for the Chern Simons term gauge invariance would require β=0\beta=0). Because this term is covariant without use of the metric it does not acquire factors of 𝒱2{\cal V}_{2} or Mp/MM_{p}/M when going to 4D Einstein frame.

The dual effective theory for 𝔞{\mathfrak{a}} is then found by adding (74) to the kinetic term (69), imposing the Bianchi identity dG=Ω/M2{\rm d}G=\Omega/M^{2} and integrating out hμνλh_{\mu\nu\lambda}, modifying (70) to

dual=g[Mp2e2ϕ2𝒱22Dμ𝔞Dμ𝔞+13!𝔞ϵμνβρΩμνβρ]{\cal L}_{\rm dual}=-\sqrt{-g}\left[\frac{M_{p}^{2}e^{2\phi}}{2{\cal V}_{2}^{2}}\,D_{\mu}{\mathfrak{a}}\,D^{\mu}{\mathfrak{a}}+\frac{1}{3!}\,{\mathfrak{a}}\,\epsilon^{\mu\nu\beta\rho}\Omega_{\mu\nu\beta\rho}\right] (75)

where

Dμ𝔞:=μ𝔞+cMp2e(β+2)ϕ𝒱22Jμ=μ𝔞+cfa2eβϕJμ.D_{\mu}{\mathfrak{a}}:=\partial_{\mu}{\mathfrak{a}}+\frac{c}{M_{p}^{2}}\,e^{-(\beta+2)\phi}\,{\cal V}_{2}^{2}\,J_{\mu}=\partial_{\mu}{\mathfrak{a}}+\frac{c}{f_{a}^{2}}\,e^{-\beta\phi}\,J_{\mu}\,. (76)

As before we use the kinetic term to identify fa=(Mp/𝒱2)eϕf_{a}=(M_{p}/{\cal V}_{2\star})\,e^{\phi_{\star}}. Using Mp2=M2𝒱2M_{p}^{2}=M^{2}{\cal V}_{2\star} with 𝒱2=(MRKK)2{\cal V}_{2\star}=(MR_{\scriptscriptstyle K\hbox{\kern-0.50003pt}K})^{2} for a Kaluza-Klein length scale RKKR_{\scriptscriptstyle K\hbox{\kern-0.50003pt}K}, this implies faM𝒱21/2eϕ(1/RKK)eϕf_{a}\sim M{\cal V}_{2\star}^{-1/2}e^{\phi_{\star}}\sim(1/R_{\scriptscriptstyle K\hbox{\kern-0.50003pt}K})\,e^{\phi_{\star}}, and so famKK1/RKKf_{a}\sim m_{\scriptscriptstyle K\hbox{\kern-0.50003pt}K}\sim 1/R_{\scriptscriptstyle K\hbox{\kern-0.50003pt}K} when eϕe^{\phi_{\star}} is not that much smaller than order unity.

The physical coupling that comes from comparing the kinetic and μ𝔞Jμ\partial_{\mu}{\mathfrak{a}}\,J^{\mu} term to (73) is

gagggaff=1affc𝒱2e(β+1)ϕMpceβϕfafor couplings to J.g_{agg}\sim g_{{a\hbox{\kern-0.70004pt}f\hbox{\kern-0.70004pt}f}}=\frac{1}{{\cal F}_{{a\hbox{\kern-0.70004pt}f\hbox{\kern-0.70004pt}f}}}\sim\frac{c{\cal V}_{2\star}e^{-(\beta+1)\phi_{\star}}}{M_{p}}\sim\frac{c\,e^{-\beta\phi_{\star}}}{f_{a}}\quad\hbox{for couplings to $J$}\,. (77)

In the special case where JμJ_{\mu} is the Hodge dual of a gauge-field Chern Simons term, gauge invariance also requires we take β=0\beta=0, and once this is done the coupling in (77) agrees (up to numerical factors) with the physical coupling to Ω\Omega implied by (75). This coupling becomes strong when EfaE\sim f_{a}, which we’ve seen is of order the Kaluza-Klein scale in the special case of two extra dimensions.

TT-type axion

The lowest-dimension generally covariant and gauge invariant interaction that couples HMNPH_{{\scriptscriptstyle M\hbox{\kern-0.50003pt}N\hbox{\kern-0.50003pt}P}} to matter localized on a space-filling 3-brane and that is linear in the components HμmnH_{\mu mn} has the form

Sint=Σe2ϕHJMp2M2d4xge2ϕ𝒱22gμνμ𝔟(x)Jν(x),S_{\rm int}=\int_{\Sigma}\;e^{-2\phi}\;{}^{\star}H\wedge J\ni\frac{M_{p}^{2}}{M^{2}}\int{\rm d}^{4}x\sqrt{-g}\;\frac{e^{-2\phi}}{{\cal V}_{2}^{2}}g^{\mu\nu}\partial_{\mu}{\mathfrak{b}}(x)J_{\nu}(x)\,, (78)

where H{}^{\star}H denotes the 6D Hodge dual and we choose the ϕ\phi coupling to be the same as the kinetic term.

The kinetic and interactions terms combine to give the effective action (in 4D Einstein frame)

Seff=d4xgMp2τ2[(𝔟)2+μ𝔟JμM2]S_{\rm eff}=\int{\rm d}^{4}x\,\sqrt{-g}\;\frac{M_{p}^{2}}{\tau^{2}}\left[(\partial{\mathfrak{b}})^{2}+\frac{\partial_{\mu}{\mathfrak{b}}J^{\mu}}{M^{2}}\right] (79)

where τ:=𝒱2eϕ\tau:={\cal V}_{2}\,e^{\phi} as before. Inspection of the kinetic term identifies the decay constant as fbMp/τf_{b}\simeq M_{p}/\tau_{\star} with τ𝒱2\tau_{\star}\propto{\cal V}_{2\star} denoting the present value of τ\tau. Because Mp2M2𝒱2M_{p}^{2}\simeq M^{2}{\cal V}_{2\star} we see that fbM𝒱21/2RKK1f_{b}\sim M{\cal V}^{-1/2}_{2\star}\sim R_{\scriptscriptstyle K\hbox{\kern-0.50003pt}K}^{-1} is of order the KK scale in size.

Canonically normalizing by rescaling b=Mp𝔟/τb=M_{p}\,{\mathfrak{b}}/\tau_{\star} —- then produces a lagrangian of the form (73) with

M2τMpMpfbMpτ.{\cal F}\sim\frac{M^{2}\tau_{\star}}{M_{p}}\sim M_{p}\gg f_{b}\sim\frac{M_{p}}{\tau_{\star}}\,. (80)

As is typical for KK modes the field 𝔟Bmn{\mathfrak{b}}\in B_{mn} couples with gravitational strength. Notice that the ratio /fb𝒱2{\cal F}/f_{b}\propto{\cal V}_{2\star} can be enormous, since 𝒱2{\cal V}_{2\star} can be as large as 102810^{28} in the extreme case of two large eV-scale dimensions. In this case gauge invariance precludes choosing JJ to be the dual of the Chern Simons form of a brane localized gauge sector, even in the absence of any ϕ\phi-dependence in (78).

It is not that surprising to have a breakdown of 4D EFT methods at the KK scale, but the above discussion shows there is a difference between what happens at this scale for TT- and SS-type axions. For TT-type axions the coupling to matter is order 1/Mp1/M_{p} and this remains true above the KK scale. The breakdown of the 4D EFT is about the appearance of a multitude of new KK modes, all of which couple with gravitational strength. But the SS-type axion’s coupling to matter is proportional to E/faE/f_{a} and so actually grows to become order unity at the KK scale; what does this order unity coupling match to in the UV theory? It matches to a dimensionless extra-dimensional coupling in the UV theory: either to the coupling cc appearing in (74) or to the coupling gcsg_{cs} of BMNB_{{\scriptscriptstyle M\hbox{\kern-0.50003pt}N}} to the Chern-Simons term SMNPS_{{\scriptscriptstyle M\hbox{\kern-0.50003pt}N\hbox{\kern-0.50003pt}P}} that is implied by the field strength G=dB+gcsSG={\rm d}B+g_{cs}S. Although gcsg_{cs} is order 1/M21/M^{2} when BMNB_{{\scriptscriptstyle M\hbox{\kern-0.50003pt}N}} is dimensionless (as above), it is dimensionless once BMNB_{{\scriptscriptstyle M\hbox{\kern-0.50003pt}N}} is canonically normalized in six dimensions.

Since gauge invariance prevents the coupling (78) from containing a coupling between TT-type axions and a gauge sector localized on a 3-brane, one can ask whether such couplings are more generally forbidden. The answer to this is ‘no’ if we allow ourselves to consider gauge sectors localized on higher-dimensional branes. For instance for a six-dimensional 5-brane Σ6\Sigma_{6} they can arise from an interaction of the form

Sint,g=cM2Σ6BFFcd4xg𝔟ϵμνλρFμνFλρ,S_{\rm int,g}=cM^{2}\int_{\Sigma_{6}}B\wedge F\wedge F\propto c\int{\rm d}^{4}x\sqrt{-g}\;{\mathfrak{b}}\,\epsilon^{\mu\nu\lambda\rho}F_{\mu\nu}F_{\lambda\rho}\,, (81)

where the explicit factor of 𝒱2{\cal V}_{2} coming from the integration over the additional two dimensions cancels the normalization of the harmonic form ωmn𝒱21εmn\omega_{mn}\propto{\cal V}_{2}^{-1}\varepsilon_{mn}. For a canonically normalized scalar this would imply gaγγ1/fbg_{a\gamma\gamma}\sim 1/f_{b}.

The upshot is this: the model-dependent TT-type axions can couple surprisingly weakly to non-gauge matter compared to the scale set by their decay constant: 1/1/Mp1/fbRKK1/{\cal F}\sim 1/M_{p}\ll 1/f_{b}\sim R_{\scriptscriptstyle K\hbox{\kern-0.50003pt}K}. By contrast, the model-independent SS-type axion couples to matter with strength 1/1/fa1/fb1/{\cal F}\sim 1/f_{a}\sim 1/f_{b} and the same is true of TT-type couplings to gauge fields on higher-dimensional branes. From the point of view of the underlying string coupling eϕe^{\phi} the duality that maps bμνb_{\mu\nu} to 𝔞{\mathfrak{a}} is also a strong/weak coupling duality.

5 Conclusions

Axions (or ALPs) are often motivated by appealing to string theory, which seems to provide them with abundance. But string theory also provides strong concrete evidence for the assertion that exact global symmetries cannot survive contact with quantum gravity; the observation that underlies the UV quality problem for attempts to solve the strong-CP problem using a global PQ symmetry.

We here reconsider some of the implications that follow from the observation that axions arise as antisymmetric tensor fields in higher dimensions and that Standard Model fields usually live in localised objects like D-branes within the extra dimensions. Axions arise in two general types in this way: the model independent SS-type axion originating from a two-form field in 4-dimensions after compactification; and the model dependent TT-type axion such as arises as a Kaluza-Klein mode for an extra-dimensional tensor field (of which we focus for simplicity on two-form potentials in two extra dimensions).171717Axions may come from other forms such as three or four forms in ten dimensions depending of which type of string theory.

These allow for a rich structure of axion phenomenology and each type of axion can be adapted to realize the PQ solution to the strong CP problem. It has been known for a while that UV effects can affect the original PQ proposal by generating effective interactions that violate the global PQ symmetry and modify the prediction for the axion mass: the axion quality problem. We revisit how this problem arises for the two types of axion using the UV tools at hand.

We find that for TT-type axions the quality problem resembles the form originally studied, since the UV theory directly provides a pseudoscalar field once compactified to four dimensions. Our framework differs from early versions of the quality problem that imagine the PQ symmetry to be linearly realized by a complex scalar at energies E>fE>f, but generally agree with estimates based on the contributions due to wormholes or gravitational instantons below a compactification scale. To the extent that these contributions are exponentially suppressed their constraints are mild.

The SS-type case is more interesting since both the strong CP and axion quality problems must first be reformulated in terms of the two-form field and its field strength. The PQ mechanism involves giving a mass to the axion and so on the dual side involves the ‘eating’ of a 3-form potential along the lines proposed in Quevedo:1996uu . The required 3-form potential is generated by the QCD sector itself as a non-propagating topological field in the EFT below the QCD scale. As applied to QCD our re-analysis broadly agrees with that of DualStrongCP in concluding that the quality problem gets recast as an issue that arises when there are multiple 3-form fields present in the low-energy theory. This might be imagined to be a problem for string theory, for which 3-form potentials are as ubiquitous as axions.

Prompted by recent discussions of this problem Sakhelashvili:2021eid ; Dvali:2022fdv we formulate the two properties which new 3-form fields must have if they are to threaten the PQ solution to the strong CP problem, arguing why string-generated 3-form fields are not generically a problem, largely because these fields need not couple to QCD (in string theory it depends on how bulk fields couple to brane fields and usually only one couples to QCD). The gravitational Chern-Simons term does couple to QCD but whether or not it sparks a new strong CP problem depends on whether or not gravity is strongly coupled in the UV. The discussions of DualStrongCP ; Dvali:2022fdv assume that it is, but we argue that if it is not (such as if the UV completion is a weakly coupled string vacuum) then the estimates for the size of the problem are again exponentially suppressed and so would not pose a quality problem.

Finally we also explore other UV implications for axion physics. We found that depending on the brane configuration hosting the Standard Model, extra dimensions can dramatically suppress physical couplings between the axion and Standard Model sector relative to the axion decay constant appearing in the axion kinetic term, especially if the volume of the extra dimensions is very large. This is possible for TT-type axions but in the the examples examined does so only for non-gauge couplings (making this observation more pertinent for ALPs, whose properties would tend to be ‘fermiophobic’).

For SS-type fields both kinds of couplings have similar size.181818For supersymmetric realizations this can be seen because both the Kähler potential and gauge kinetic function depend directly on the SS field. For this case though, we argue that the duality relating the 2-form to the axion field swaps weak and strong couplings, and suggests a semiclassical description of 2-form response need not correspond to the usual semiclassical description of a scalar axion. This again motivates better exploring the 2-form side of the theory.

It is an old argument that UV information can have important implications for low-energy naturalness questions such as the strong CP problem. The observation that this could be informative in situations where the questions are solved using features like global symmetries that apparently should not be present at very high energies has sparked a revival of studies of generalised and non-invertible symmetries. Many of these ideas resonate well with string-motivated constructions, such as those we explore here.

Acknowledgements

We thank Philippe Brax and Junwu Huang for helpful conversations. CB’s research was partially supported by funds from the Natural Sciences and Engineering Research Council (NSERC) of Canada. Research at the Perimeter Institute is supported in part by the Government of Canada through NSERC and by the Province of Ontario through MRI. The work of FQ has been partially supported by STFC consolidated grants ST/P000681/1, ST/T000694/1.

References

  • (1) R. D. Peccei and H. R. Quinn, “CP Conservation in the Presence of Instantons,” Phys. Rev. Lett. 38 (1977), 1440-1443 “Constraints Imposed by CP Conservation in the Presence of Instantons,” Phys. Rev. D 16 (1977), 1791-1797.
  • (2) S. Weinberg, “A New Light Boson?,” Phys. Rev. Lett. 40 (1978), 223-226
  • (3) F. Wilczek, “Problem of Strong PP and TT Invariance in the Presence of Instantons,” Phys. Rev. Lett. 40 (1978), 279-282
  • (4) E. Witten, “Some Properties of O(32) Superstrings,” Phys. Lett. 149B (1984) 351–356. S. M. Barr, “Harmless Axions in Superstring Theories,” Phys. Lett. 158B (1985) 397–400. K. Choi and J. E. Kim, “Harmful Axions in Superstring Models,” Phys. Lett. 154B (1985) 393. [Erratum: Phys. Lett.156B,452(1985)].
  • (5) P. Svrcek and E. Witten, “Axions In String Theory,” JHEP 06 (2006), 051 [arXiv:hep-th/0605206 [hep-th]]; J. P. Conlon, “The QCD axion and moduli stabilisation,” JHEP 05 (2006), 078 [arXiv:hep-th/0602233 [hep-th]].
  • (6) A. Arvanitaki, S. Dimopoulos, S. Dubovsky, N. Kaloper and J. March-Russell, “String Axiverse,” Phys. Rev. D 81 (2010), 123530 [arXiv:0905.4720 [hep-th]]. M. Cicoli, M. Goodsell and A. Ringwald, “The type IIB string axiverse and its low-energy phenomenology,” JHEP 10 (2012), 146 [arXiv:1206.0819 [hep-th]].
  • (7) J. Preskill, M. Wise and F. Wilczek, “Cosmology of the Invisible Axion”, Phys.Lett.B 120 (1983) 127-132 M. Dine and W. Fischler, “The Not So Harmless Axion”, Phys. Lett. 120B (1983) 137–141. L. F. Abbott and P. Sikivie, “A Cosmological Bound on the Invisible Axion,” Phys. Lett. 120B (1983) 133–136.
  • (8) T. Banks and L. J. Dixon, “Constraints on String Vacua with Space-Time Supersymmetry,” Nucl. Phys. B 307 (1988), 93-108 D. Harlow and H. Ooguri, “Symmetries in quantum field theory and quantum gravity,” Commun. Math. Phys. 383 (2021) no.3, 1669-1804 [arXiv:1810.05338 [hep-th]].
  • (9) C. P. Burgess, J. P. Conlon, L. Y. Hung, C. H. Kom, A. Maharana and F. Quevedo, “Continuous Global Symmetries and Hyperweak Interactions in String Compactifications,” JHEP 07 (2008), 073 [arXiv:0805.4037 [hep-th]].
  • (10) M. Kamionkowski and J. March-Russell, “Planck scale physics and the Peccei-Quinn mechanism,” Phys. Lett. B 282 (1992), 137-141 [arXiv:hep-th/9202003 [hep-th]]. S. M. Barr and D. Seckel, “Planck scale corrections to axion models,” Phys. Rev. D 46 (1992), 539-549 S. Ghigna, M. Lusignoli and M. Roncadelli, “Instability of the invisible axion,” Phys. Lett. B 283 (1992), 278-281 R. Holman, S. D. H. Hsu, T. W. Kephart, E. W. Kolb, R. Watkins and L. M. Widrow, “Solutions to the strong CP problem in a world with gravity,” Phys. Lett. B 282 (1992), 132-136 [arXiv:hep-ph/9203206 [hep-ph]].
  • (11) M. Kalb and P. Ramond, “Classical direct interstring action,” Phys. Rev. D 9 (1974), 2273-2284
  • (12) R. Savit, “Duality in Field Theory and Statistical Systems,” Rev. Mod. Phys. 52 (1980), 453
  • (13) R. Kallosh, A. D. Linde, D. A. Linde and L. Susskind, “Gravity and global symmetries,” Phys. Rev. D 52 (1995), 912-935 [arXiv:hep-th/9502069 [hep-th]].
  • (14) G. Dvali, “Three-form gauging of axion symmetries and gravity,” [arXiv:hep-th/0507215 [hep-th]].
  • (15) O. Sakhelashvili, “Consistency of the dual formulation of axion solutions to the strong CP problem,” Phys. Rev. D 105 (2022) no.8, 085020 [arXiv:2110.03386 [hep-th]].
  • (16) G. Dvali, “Strong-CPCP with and without gravity,” [arXiv:2209.14219 [hep-ph]].
  • (17) C. P. Burgess, “The Cosmological Constant Problem: Why it’s hard to get Dark Energy from Micro-physics,” [arXiv:1309.4133 [hep-th]].
  • (18) N. Kaloper and L. Sorbo, “A Natural Framework for Chaotic Inflation,” Phys. Rev. Lett.  102 (2009) 121301 [arXiv:0811.1989 [hep-th]]; N. Kaloper, A. Lawrence and L. Sorbo, “An Ignoble Approach to Large Field Inflation,” JCAP 1103 (2011) 023 [arXiv:1101.0026 [hep-th]].
  • (19) C. P. Burgess, D. Dineen and F. Quevedo, “Yoga Dark Energy: Natural Relaxation and Other Dark Implications of a Supersymmetric Gravity Sector,” [arXiv:2111.07286 [hep-th]].
  • (20) S. Weinberg, “Why the Renormalization Group is a Good Thing,” contribution to Asymptotic Realms of Physics festschrift for Francis Low (Cambridge 1981).
  • (21) D. Tong, “Lectures on the Quantum Hall Effect,” [arXiv:1606.06687 [hep-th]]. E. Witten, “Two Lectures on Gauge Theory and Khovanov Homology,” [arXiv:1603.03854 [math.GT]].
  • (22) C. P. Burgess, Introduction to Effective Field Theory, Cambridge University Press (2020) ISBN 978-1-139-04804-0, 978-0-521-19547-8
  • (23) R. Bousso and J. Polchinski, “Quantization of four form fluxes and dynamical neutralization of the cosmological constant,” JHEP 0006, 006 (2000) [hep-th/0004134].
  • (24) C. P. Burgess, R. Diener and M. Williams, “Self-Tuning at Large (Distances): 4D Description of Runaway Dilaton Capture,” JHEP 1510 (2015) 177 [arXiv:1509.04209 [hep-th]].
  • (25) S. Bielleman, L. E. Ibanez and I. Valenzuela, “Minkowski 3-forms, Flux String Vacua, Axion Stability and Naturalness,” JHEP 1512 (2015) 119 [arXiv:1507.06793 [hep-th]].
  • (26) A. Herraez, L. E. Ibanez, F. Marchesano and G. Zoccarato, “The Type IIA Flux Potential, 4-forms and Freed-Witten anomalies,” JHEP 1809 (2018) 018 [arXiv:1802.05771 [hep-th]].
  • (27) F. Quevedo and C. A. Trugenberger, “Phases of antisymmetric tensor field theories,” Nucl. Phys. B 501 (1997), 143-172 [arXiv:hep-th/9604196 [hep-th]]; F. Quevedo and C. A. Trugenberger, “Condensation of p-branes and generalized Higgs / confinement duality,” Int. J. Mod. Phys. A 12 (1997), 1227-1236 [arXiv:hep-th/9608171 [hep-th]].
  • (28) C. P. Burgess and F. Quevedo, “Who’s Afraid of the Supersymmetric Dark? The Standard Model vs Low-Energy Supergravity,” Fortsch. Phys. 70 (2022) no.7-8, 2200077 [arXiv:2110.13275 [hep-th]].
  • (29) Y. Aghababaie, C. P. Burgess, S. L. Parameswaran and F. Quevedo, “Towards a naturally small cosmological constant from branes in 6-D supergravity,” Nucl. Phys. B680 (2004) 389 [arXiv:hep-th/0304256]; C.P. Burgess, “Supersymmetric Large Extra Dimensions and the Cosmological Constant: An Update,” Ann. Phys. 313 (2004) 283-401 [arXiv:hep-th/0402200]; “Towards a natural theory of dark energy: Supersymmetric large extra dimensions,” in the proceedings of the Texas A&M Workshop on String Cosmology, [arXiv:hep-th/0411140]; C. P. Burgess and L. van Nierop, “Technically Natural Cosmological Constant From Supersymmetric 6D Brane Backreaction,” Phys. Dark Univ.  2 (2013) 1 [arXiv:1108.0345 [hep-th]].
  • (30) C. P. Burgess, R. Diener and M. Williams, “EFT for Vortices with Dilaton-dependent Localized Flux,” JHEP 1511 (2015) 054 [arXiv:1508.00856 [hep-th]].
  • (31) C. P. Burgess, R. Diener and M. Williams, “A problem with δ\delta-functions: stress-energy constraints on bulk-brane matching (with comments on arXiv:1508.01124),” JHEP 01 (2016), 017 [arXiv:1509.04201 [hep-th]].
  • (32) F. Niedermann and R. Schneider, “Fine-tuning with Brane-Localized Flux in 6D Supergravity,” JHEP 02 (2016), 025 [arXiv:1508.01124 [hep-th]].
  • (33) T. H. Buscher, “Path Integral Derivation of Quantum Duality in Nonlinear Sigma Models,” Phys. Lett. B 201 (1988), 466-472
  • (34) B. Julia and G. Toulouse, J. Physique, Lett. 40 (1979) 396.
  • (35) M. Kosterlitz and D. Thouless, J. Phys. C6 (1973) 1181; V.L. Berezinski, JETP 34 (1972) 610; J. Villain, J. Phys. C36 (19975) 581.
  • (36) T. Banks and N. Seiberg, “Symmetries and Strings in Field Theory and Gravity,” Phys. Rev. D 83 (2011), 084019 [arXiv:1011.5120 [hep-th]].
  • (37) M. Luscher, “The Secret Long Range Force in Quantum Field Theories With Instantons,” Phys. Lett. B 78 (1978), 465-467
  • (38) Holman, R. and Kephart, T. W. and Rey, Soo-Jong, PhysRevLett.71.320, [arXiv:hep-ph/9207208].
  • (39) G. W. Gibbons and S. W. Hawking, “Gravitational Multi - Instantons,” Phys. Lett. B 78 (1978), 430
  • (40) T. Eguchi and A. J. Hanson, “Selfdual Solutions to Euclidean Gravity,” Annals Phys. 120 (1979), 82
  • (41) S. B. Giddings and A. Strominger, “Axion Induced Topology Change in Quantum Gravity and String Theory,” Nucl. Phys. B 306 (1988), 890-907
  • (42) A. Hook, “TASI Lectures on the Strong CP Problem and Axions,” PoS TASI2018 (2019), 004 [arXiv:1812.02669 [hep-ph]].
  • (43) B. Heidenreich, J. McNamara, M. Montero, M. Reece, T. Rudelius, and I. Valenzuela, Non-invertible global symmetries and completeness of the spectrum, JHEP 09 (2021) 203, [arXiv:2104.07036].
  • (44) H. Nishino and E. Sezgin, “Matter and Gauge Couplings of N=2 Supergravity in Six-Dimensions,” Phys. Lett. B 144 (1984), 187-192; H. Nishino and E. Sezgin, “The Complete N=2N=2, d=6d=6 Supergravity With Matter and Yang-Mills Couplings,” Nucl. Phys. B 278 (1986), 353-379.
  • (45) A.H. Chamseddine, “Interacting Supergravity in 10 Dimensions: the role of the six-index gauge field,” Phys. Rev. D24 (1981) 3065; E. Bergshoeff, M. de Roo, B. de Wit and P. van Nieuwenhuizen, “Ten dimensional Maxwell-Einstein Supergravity, its currents and the issue of its auxiliary fields,” Nucl. Phys. B195 (1982) 97; G. Chapline and N.S. Manton, “Unification of Yang-Mills Theory and Supergravity in 10 Dimensions,” Phys. Lett. 120B (1983) 105; J.H. Schwarz, “Covariant Field Equations of Chiral N=2N=2 and D=10D=10 Supergravity,” Nucl. Phys. B226 (1983) 269. Supergravities in Diverse Dimensions, Vols. 1 & 2 ed. by A. Salam and E. Sezgin, World Scientific (Singapore) 1989.