Axions are blind to anomalies

TL;DR

This paper demonstrates that axion couplings to gauge bosons originate from non-anomalous fermion loops, challenging the traditional anomaly-based understanding, especially impacting chiral gauge theories and axion-like particle searches.

Contribution

It proves that axion couplings are independent of anomaly structures and provides explicit calculations showing deviations from expected anomaly patterns in chiral gauge theories.

Findings

Axion couplings arise from non-anomalous fermion loops.

Couplings to electroweak gauge bosons do not follow anomaly patterns.

Implications for axion-like particle searches are significant.

Abstract

The axion couplings to SM gauge bosons are derived in various models, and shown to always arise entirely from non-anomalous fermion loops. They are thus independent of the anomaly structure of the model. This fact is without consequence for vector gauge interactions like QCD and QED, but has a major impact for chiral gauge theories. For example, in the DFSZ axion model, the couplings of axions to electroweak gauge bosons do not follow the pattern expected from chiral anomalies, as we prove by an explicit calculation. The reason for this mismatch is traced back to triangle Feynman diagrams sensitive to the anomalous breaking of the vector Ward identity, and is ultimately related to the conservation of baryon and lepton numbers. Though our analyses are entirely done for true axion models, this observation could have important consequences for axion-like particle searches.

Tables

Table 1. Table 1: Coefficients of the effective axion to gauge boson couplings of Eq. ( 95 ) in the m u , d , e → ∞ → subscript 𝑚 𝑢 𝑑 𝑒 m_{u,d,e}\rightarrow\infty limit. For the linear representation, those are found directly from the THDM amplitudes of Ref.  [ 19 ] . For the polar representation, the contributions of the local anomalous terms and that of the triangle amplitudes built on the axial ( A 𝐴 A ) and vector ( V 𝑉 V ) derivative interactions have to be added together. The fact that only the three independent coefficients 𝒩 g ​ g superscript 𝒩 𝑔 𝑔 \mathcal{N}^{gg} , 𝒩 γ ​ γ superscript 𝒩 𝛾 𝛾 \mathcal{N}^{\gamma\gamma} , and 𝒩 L subscript 𝒩 𝐿 \mathcal{N}_{L} occur for the local anomalous interactions comes from their S ​ U ​ ( 2 ) L ⊗ U ​ ( 1 ) Y tensor-product 𝑆 𝑈 subscript 2 𝐿 𝑈 subscript 1 𝑌 SU(2)_{L}\otimes U(1)_{Y} invariance. The coefficient 𝒩 L = − 1 / 2 ​ ( 3 ​ α + β ) subscript 𝒩 𝐿 1 2 3 𝛼 𝛽 \mathcal{N}_{L}=-1/2(3\alpha+\beta) , with α 𝛼 \alpha and β 𝛽 \beta being the free parameters tuning the U ​ ( 1 ) ℬ 𝑈 subscript 1 ℬ U(1)_{\mathcal{B}} and U ​ ( 1 ) ℒ 𝑈 subscript 1 ℒ U(1)_{\mathcal{L}} components of U ​ ( 1 ) P ​ Q 𝑈 subscript 1 𝑃 𝑄 U(1)_{PQ} . The explicit presence of β 𝛽 \beta in the a 0 → Z ​ Z , W ​ W → superscript 𝑎 0 𝑍 𝑍 𝑊 𝑊 a^{0}\rightarrow ZZ,WW triangles is due to the peculiar nature of the neutrinos, which are kept massless. They do not contribute in the linear representation but have to explicitly appear in the anomalous interactions and derivative terms since those are S ​ U ​ ( 2 ) L ⊗ U ​ ( 1 ) Y tensor-product 𝑆 𝑈 subscript 2 𝐿 𝑈 subscript 1 𝑌 SU(2)_{L}\otimes U(1)_{Y} invariant.

Linear	Polar
$a^0\overline{\psi }{\gamma }_5\psi$	Anomalous	${\partial }_{\mu }{a}^0\overline{\psi }{\gamma }^{\mu }{\gamma }_5\psi$		${\partial }_{\mu }{a}^0\overline{\psi }{\gamma }^{\mu }\psi$
	interactions	$AVV$	$AAA$	$VAV$
${𝒩}^{gg}=\frac{1}{2}\left(x+\frac{1}{x}\right)$	${𝒩}^{gg}$	0	$-$
${𝒩}^{\gamma \gamma }=\frac{4}{3}\left(x+\frac{1}{x}\right)$	${𝒩}^{\gamma \gamma }$	0	$-$
${𝒩}_1^{\gamma Z}=\frac{1}{2}\left(x+\frac{1}{x}\right)$	${𝒩}_L$	0	$-$
${𝒩}_2^{\gamma Z}={𝒩}^{\gamma \gamma }$	${𝒩}^{\gamma \gamma }$	0	$-$	0
${𝒩}_1^{ZZ}=\frac{1}{4}x+\frac{1}{3x}$	${𝒩}_L$	$\frac{\beta }{16}$	$-\frac{1}{2}{𝒩}_1^{ZZ}+\frac{\beta }{16}$
${𝒩}_2^{ZZ}={𝒩}_1^{\gamma Z}$	${𝒩}_L$	0	0
${𝒩}_3^{ZZ}={𝒩}^{\gamma \gamma }$	${𝒩}^{\gamma \gamma }$	0	0	0
${𝒩}^{WW}=\frac{x}{4}+\frac{3}{8x}$	${𝒩}_L$	$\frac{3}{2}{𝒩}^{WW}-\frac{3}{2}{𝒩}_1^{\gamma Z}+\frac{\beta }{16}$	$-\frac{1}{2}{𝒩}^{WW}+\frac{\beta }{16}$

Equations

$${\mathcal{L}}_{eff}=\frac{a^0}{16{\pi }^{2}v}\left(g_s^2{\mathcal{N}}_C{G}_{\mu \nu }^a{\tilde{G}}^{a,\mu \nu }+e^2{\mathcal{N}}_{em}{F}_{\mu \nu }{\tilde{F}}^{\mu \nu }\right),$$

$${\mathcal{L}}_{eff}=\frac{a^0}{16{\pi }^{2}v}\left(g_s^2{\mathcal{N}}_C{G}_{\mu \nu }^a{\tilde{G}}^{a,\mu \nu }+g^2{\mathcal{N}}_L{W}_{\mu \nu }^i{\tilde{W}}^{i,\mu \nu }+g^{\prime 2}{\mathcal{N}}_Y{B}_{\mu \nu }{\tilde{B}}^{\mu \nu }\right),$$

$$\mathcal{L}=-\frac{1}{4}{F}_{\mu \nu }{F}^{\mu \nu }+{\bar{\psi }}_L(i/D){\psi }_L+{\bar{\psi }}_R(i/D){\psi }_R+(y\varphi {\bar{\psi }}_L{\psi }_R+h.c.)+{\partial }_{\mu }{\varphi }^{†}{\partial }^{\mu }\varphi -V(\varphi ),$$

$$\varphi \to \exp (-i\theta )\varphi ,{\psi }_L\to \exp (i\alpha \theta ){\psi }_L,{\psi }_R\to \exp (i(\alpha +1)\theta ){\psi }_R.$$

$$\varphi (x)={\sigma }^0(x)+i{a}^0(x)+v,$$

$${\mathcal{L}}_{\mathrm{Linear}}$$

$$+\frac{1}{2}({\partial }_{\mu }{\sigma }^0{\partial }^{\mu }{\sigma }^0-2\lambda v^2({\sigma }^0{)}^2)-\lambda v{\sigma }^0(({\sigma }^0{)}^2+(a^0{)}^2)-\frac{\lambda }{4}(({\sigma }^0{)}^2+(a^0{)}^2{)}^2-\frac{\lambda v^4}{4}.$$

$${\mathcal{T}}_{PVV}^{\alpha \beta }$$

$$=-i\frac{1}{2{\pi }^2}{\epsilon }^{\alpha \beta \rho \sigma }{q}_{1,\rho }{q}_{2,\sigma }(m{C}_0(m^2)-M{C}_0(M^2)).$$

$$\underset{m\to \infty }{\lim }{C}_0(m^2)=\frac{-1}{2m^2},\underset{m\to 0}{\lim }{m}^2{C}_0(m^2)=0.$$

$$\mathcal{M}(a^0\to \gamma \gamma {)}_{\mathrm{Linear}}=-i\frac{m}{v}{e}^2{\mathcal{T}}_{PVV}^{\alpha \beta }\epsilon (q_1{)}_{\alpha }^{\ast }\epsilon (q_2{)}_{\beta }^{\ast }=-\frac{e^2}{2{\pi }^{2}v}{m}^2{C}_0(m^2){\epsilon }^{\alpha \beta \rho \sigma }\epsilon (q_1{)}_{\alpha }^{\ast }\epsilon (q_2{)}_{\beta }^{\ast }{q}_{1,\rho }{q}_{2,\sigma }.$$

$${\mathcal{L}}_{\mathrm{Linear}}^{\mathrm{eff}}=-\frac{e^2}{16{\pi }^{2}v}{a}^0{F}_{\mu \nu }{\tilde{F}}^{\mu \nu }.$$

$$\varphi (x)=\frac{1}{\sqrt{2}}(v+{\sigma }^0(x))e^{-i{a}^0(x)/v}.$$

$${\psi }_L(x)\to \exp (i\alpha a^0(x)/v){\psi }_L(x),{\psi }_R(x)\to \exp (i(\alpha +1)a^0(x)/v){\psi }_R(x).$$

$$\delta {\mathcal{L}}_{\mathrm{Der}}=-\frac{{\partial }_{\mu }{a}^0}{v}(\alpha {\bar{\psi }}_L{\gamma }^{\mu }{\psi }_L+(\alpha +1){\bar{\psi }}_R{\gamma }^{\mu }{\psi }_R)=-\frac{{\partial }_{\mu }{a}^0}{2v}((2\alpha +1)\bar{\psi }{\gamma }^{\mu }\psi +\bar{\psi }{\gamma }^{\mu }{\gamma }_5\psi ),$$

$$\delta {\mathcal{L}}_{\mathrm{Jac}}=\frac{e^2}{16{\pi }^{2}v}{a}^0(\alpha -(\alpha +1))F_{\mu \nu }{\tilde{F}}^{\mu \nu }=-\frac{e^2}{16{\pi }^{2}v}{a}^0{F}_{\mu \nu }{\tilde{F}}^{\mu \nu }.$$

$${\mathcal{L}}_{{}_{\mathrm{Polar}}}$$

$$+\frac{1}{2}{\partial }_{\mu }{a}^0{\partial }^{\mu }{a}^0{\left(1+\frac{{\sigma }^0}{v}\right)}^2-\frac{{\partial }_{\mu }{a}^0}{2v}((2\alpha +1)\bar{\psi }{\gamma }^{\mu }\psi +\bar{\psi }{\gamma }^{\mu }{\gamma }_5\psi )-\frac{e^2}{16{\pi }^{2}v}{a}^0{F}_{\mu \nu }{\tilde{F}}^{\mu \nu }$$

$$+\frac{1}{2}({\partial }_{\mu }{\sigma }^0{\partial }^{\mu }{\sigma }^0-2\lambda v^2({\sigma }^0{)}^2)-\lambda v({\sigma }^0{)}^3-\frac{\lambda }{4}({\sigma }^0{)}^4-\frac{\lambda v^4}{4}.$$

$${\mathcal{T}}_{AVV}^{\gamma \alpha \beta }=\int \frac{d^{4}k}{(2\pi {)}^4}(-1)\mathrm{Tr}\left[\frac{i}{/k-/q_1-m}{\gamma }^{\alpha }\frac{i}{/k-m}{\gamma }^{\beta }\frac{i}{/k+/q_2-m}{\gamma }^{\gamma }{\gamma }_5\right]+(1,\alpha \leftrightarrow 2,\beta ).$$

$$\mathcal{M}(a^0\to \gamma \gamma {)}_{\mathrm{Der}}$$

$$=\frac{-e^2}{4{\pi }^{2}v}(2m^2{C}_0(m^2)-2M^2{C}_0(M^2)){\epsilon }^{\alpha \beta \rho \sigma }\epsilon (q_1{)}_{\alpha }^{\ast }\epsilon (q_2{)}_{\beta }^{\ast }{q}_{1,\rho }{q}_{2,\sigma }.$$

$$\mathcal{M}(a^0\to \gamma \gamma {)}_{{}_{\mathrm{Polar}}}$$

$$=-\frac{e^2}{2{\pi }^{2}v}{m}^2{C}_0(m^2){\epsilon }^{\alpha \beta \rho \sigma }\epsilon (q_1{)}_{\alpha }^{\ast }\epsilon (q_2{)}_{\beta }^{\ast }{q}_{1,\rho }{q}_{2,\sigma }=\mathcal{M}(a^0\to \gamma \gamma {)}_{\mathrm{Linear}}.$$

$${\partial }_{\mu }{A}^{\mu }-\frac{e^2}{8{\pi }^2}{F}_{\mu \nu }{\tilde{F}}^{\mu \nu }=2imP,$$

$$\mathcal{M}(a^0\to \gamma \gamma {)}_{\mathrm{Linear}}=\mathcal{M}(a^0\to \gamma \gamma {)}_{{}_{\mathrm{Polar}}}\stackrel{m\to 0}{=}0.$$

$$\mathcal{M}(a^0\to \gamma \gamma {)}_{\mathrm{Der}}\stackrel{m\to \infty }{=}0,$$

$${\mathcal{L}}_{\mathrm{KSVZ}}={\mathcal{L}}_{\mathrm{SM}}+{\bar{\psi }}_L(i/D){\psi }_L+{\bar{\psi }}_R(i/D){\psi }_R+(y\varphi {\bar{\psi }}_L{\psi }_R+h.c.)+{\partial }_{\mu }{\varphi }^{†}{\partial }^{\mu }\varphi -V(\varphi ),$$

$$D^{\mu }={\partial }^{\mu }-i{g}_s{G}_a^{\mu }{T}^a-ig{W}_i^{\mu }{T}^i-i{g}^{\prime }\frac{Y}{2}{B}^{\mu },$$

$${\mathcal{L}}_{\mathrm{Linear}}^{\mathrm{eff},\mathrm{KSVZ}}=-\frac{1}{16{\pi }^{2}v}{a}^0(g_s^2{d}_L^{\psi }{C}_C^{\psi }{G}_{\mu \nu }^a{\tilde{G}}^{a,\mu \nu }+g^2{d}_C^{\psi }{C}_L^{\psi }{W}_{\mu \nu }^i{\tilde{W}}^{i,\mu \nu }+g^{\prime 2}{d}_L^{\psi }{d}_C^{\psi }{C}_Y^{\psi }{B}_{\mu \nu }{\tilde{B}}^{\mu \nu }),$$

$${\mathcal{T}}_{LLL}^{\alpha \beta \gamma }=\int \frac{d^{4}k}{(2\pi {)}^4}(-1)\mathrm{Tr}\left[\frac{i}{/k-/q_1}{\gamma }^{\beta }{P}_L\frac{i}{/k}{\gamma }^{\gamma }{P}_L\frac{i}{/k+/q_2}{\gamma }^{\alpha }{P}_L\right]+(1,\beta \leftrightarrow 2,\gamma ),$$

$$i(q_1+q_2{)}_{\alpha }{\mathcal{T}}_{LLL}^{\alpha \beta \gamma }$$

$$-i(q_1{)}_{\beta }{\mathcal{T}}_{LLL}^{\alpha \beta \gamma }$$

$$-i(q_2{)}_{\gamma }{\mathcal{T}}_{LLL}^{\alpha \beta \gamma }$$

Full text

March 2019

Axions are blind to anomalies

Jérémie Quevillon1 and Christopher Smith2

Laboratoire de Physique Subatomique et de Cosmologie, Université Grenoble-Alpes, CNRS/IN2P3, Grenoble INP, 38000 Grenoble, France.

Abstract

The axion couplings to SM gauge bosons are derived in various models, and shown to always arise entirely from non-anomalous fermion loops. They are thus independent of the anomaly structure of the model. This fact is without consequence for vector gauge interactions like QCD and QED, but has a major impact for chiral gauge theories. For example, in the DFSZ axion model, the couplings of axions to electroweak gauge bosons do not follow the pattern expected from chiral anomalies, as we prove by an explicit calculation. The reason for this mismatch is traced back to triangle Feynman diagrams sensitive to the anomalous breaking of the vector Ward identity, and is ultimately related to the conservation of baryon and lepton numbers. Though our analyses are entirely done for true axion models, this observation could have important consequences for axion-like particle searches.

††footnotetext: ${}^{1}\;$jeremie.quevillon@lpsc.in2p3.fr††footnotetext: ${}^{2}\;$chsmith@lpsc.in2p3.fr

1 Introduction

Axions have been around for more than 40 years, and even if the simplest models have been ruled out, they still remain probably the best solution for the strong $CP$ problem of the Standard Model. They come under many guises, but the basic principle is always the same: design a spontaneously broken global $U(1)_{PQ}$ symmetry and assign chiral charges to some colored fermions [1]. This ensures the associated Goldstone boson $a^{0}$, the axion [2, 3], has a direct coupling to gluons, and possibly also to photons, of the form

[TABLE]

where $v$ is the vacuum expectation value (VEV) of the scalar breaking the $U(1)_{PQ}$ symmetry, and $\mathcal{N}_{C,em}$ some constants related to the fermion $U(1)_{PQ}$ and gauge charges. Naively, the $CP$-violating QCD coupling $\theta_{QCD}G_{\mu\nu}\tilde{G}^{\mu\nu}$ can then be rotated away no matter the value of $\theta_{QCD}$ since the Goldstone boson is shifted as $a^{0}\rightarrow a^{0}+v\theta$ under a $U(1)_{PQ}$ transformation of parameter $\theta$. In reality, what happens is that non-perturbative QCD effects dressing $G_{\mu\nu}\tilde{G}^{\mu\nu}$ create an effective potential for the axion field [1] (see also e.g. Ref. [5]), whose minimum is attained precisely when $\left\langle a^{0}\right\rangle+v\theta_{QCD}=0$. In the process, the axion acquires a small QCD-induced mass, typically well below the eV scale [6, 7]. For a recent review, see for example Ref. [8].

Though this picture is correct, the nature of the couplings in Eq. (1) is often wrongly ascribed to the anomaly of the $U(1)_{PQ}$ fermionic current. This idea comes from Noether’s theorem: the Goldstone boson $a^{0}$ is coupled to its symmetry current, $vp^{\mu}=\langle 0|J_{PQ}^{\mu}|a^{0}(p)\rangle$, and the current is anomalous, $\partial_{\mu}J_{PQ}^{\mu}\sim g_{s}^{2}\mathcal{N}_{C}G_{\mu\nu}^{a}\tilde{G}^{a,\mu\nu}+e^{2}\mathcal{N}_{em}F_{\mu\nu}\tilde{F}^{\mu\nu}$, hence it seems Eq. (1) is immediately recovered. It has been known for a long time that this derivation of the axion couplings is not correct but sufficient for practical purpose [9, 10]. It permits to identify the couplings of axions to gluons and photons as induced by heavy fermions. So, by a common abuse of language, the $a\rightarrow\gamma\gamma$ and $a\rightarrow gg$ processes are said to be induced by the anomaly in the $U(1)_{PQ}$ current.

However, as we will explicitly demonstrate in this paper, using this same procedure to derive the couplings of the axion to electroweak gauge bosons fails whenever chiral fermions are charged under $U(1)_{PQ}$. In that case, the true gauge couplings do not follow the pattern expected from the anomalies in the $U(1)_{PQ}$ current. More precisely, for the $SU(3)_{C}\otimes SU(2)_{L}\otimes U(1)_{Y}$ gauge group, one expects from chiral anomalies [9] (see also Ref. [11])

[TABLE]

where $\mathcal{N}_{C,L,Y}$ are derived from the $U(1)_{PQ}$ charges and gauge quantum numbers of the SM fermions. While for QCD and QED, which are vector theories, the coefficients $\mathcal{N}_{C}$ and $\mathcal{N}_{em}=\mathcal{N}_{L}+\mathcal{N}_{Y}$ indeed tune the $a^{0}\rightarrow gg$ and $\gamma\gamma$ decays, we will prove that the $a^{0}\rightarrow\gamma Z$, $ZZ$, and $W^{+}W^{-}$ decay amplitudes are independent of $\mathcal{N}_{L}$ and $\mathcal{N}_{Y}$! The true processes driving these decays are not anomalous and their relative strengths do not follow the pattern expected from Eq. (2). The reason for this mismatch between the anomalous couplings and the true axion couplings will be traced to the existence of additional anomalies, occurring only for chiral theories. These additional anomalies cancel exactly the contributions in Eq. (2), leaving as remainder only non-local, non-anomalous processes.

The paper is organized as follow. To set the stage, we start in the next section by presenting in details an axion toy model [1]. Though simple, this model illustrates many important physical features of more realistic axion models. In particular, it will be clear that the axion couplings to gauge bosons are not anomalous, but that this has no quantitative consequence for photons or gluons. Then, before turning our attention to full-fledged axion models, we derive a number of important results for anomalies in Section 3. Specifically, to be able to treat chiral gauge theories, the Ward identities applied to $AVV$ or $AAA$ triangle graphs have to be properly calculated, where $V_{\mu}=\bar{\psi}\gamma_{\mu}\psi$ and $A_{\mu}=\bar{\psi}\gamma_{\mu}\gamma_{5}\psi$ denote the vector and axial currents. Indeed, it is necessary to go beyond simple regularization procedures to be able to locate the anomaly of the $AVV$ triangle in one of the vector Ward identities or to break explicitly the Bose symmetry of the $AAA$ triangle. Equipped with these results, we turn to the Peccei-Quinn axion model in Section 4, derive the correct axion couplings to gauge bosons, and identify precisely where the naive procedure leading to Eq. (2) fails. This analysis is then trivially extended to invisible axion models. Finally, our results are summarized in Section 5, along with their implications for axion-like particle searches.

2 An axionic toy model

Let us consider the following simple extension of QED. Starting with a massless Dirac fermion $\psi$, charged under a local symmetry $U(1)_{em}$, we add a scalar field $\phi$ charged under a $U(1)_{PQ}$ global symmetry. The complete Lagrangian is

[TABLE]

where $D^{\mu}=\partial^{\mu}+ieA^{\mu}$ and $e$ is the electric charge of $\psi$, and $\psi_{L,R}=1/2(1\mp\gamma_{5})\psi$. When the Yukawa coupling is non-zero, $y\neq 0$, the fermion fields have to be charged under $U(1)_{PQ}$ also. If we normalize the scalar charge to $-1$, then under a $U(1)_{PQ}$ transformation of parameter $\theta$,

[TABLE]

At the classical level, this clearly leaves the Lagrangian invariant whatever the free parameter $\alpha$. Actually, this parameter corresponds to the conservation of the global fermion number, which is aligned with the electric charge in this simple one-fermion model.

When the global $U(1)_{PQ}$ symmetry is broken spontaneously, a Goldstone boson is left behind and will be identified with the axion [2, 3]. So, we choose the potential as $V(\phi^{\dagger}\phi)=\mu^{2}\phi^{\dagger}\phi+\lambda(\phi^{\dagger}\phi)^{2}$ with $\mu^{2}<0$. Let us analyze the resulting theory using two different representations for the scalar field.

2.1 Linear representation

Shifted to the true vacuum, the usual linear representation for the scalar field is

[TABLE]

with $v=-\mu^{2}/\lambda$. The Goldstone boson $a^{0}$ remains massless and is the axion, while the $\sigma^{0}$ field acquires a mass $m_{\sigma}^{2}=2\lambda v^{2}$. Assuming $y$ is real and denoting the fermion mass as $m\equiv vy$, the whole Lagrangian when $\phi$ is expended around the true vacuum becomes

[TABLE]

Since $\phi$ is not charged under $U(1)_{em}$, $a^{0}$ does not directly couple to photons. However, this coupling arises through one-loop triangle graphs, see Fig. 1. Adopting the Pauli-Villars procedure to regulate the loop amplitude in intermediate steps, we compute

[TABLE]

The three-point scalar function $C_{0}(m^{2})\equiv C_{0}(q_{1}^{2},q_{2}^{2},(q_{1}+q_{2})^{2};m^{2},m^{2},m^{2})$ obeys111We use the notations and conventions of Ref. [4]

[TABLE]

The regulator $M$ can thus safely be sent to infinity, and the decay amplitude in the linear representation for the scalar field is

[TABLE]

In the $m\rightarrow\infty$ limit, this amplitude corresponds to the local interaction

[TABLE]

It is interesting to remark that this term is not decoupling 222This non-decoupling is to be understood in the formal sense, since strictly speaking the $m\rightarrow\infty$ limit is not compatible with our perturbative treatment of the Yukawa couplings., even though $\mathcal{T}_{PVV}^{\alpha\beta}\rightarrow 0$ as $m\rightarrow\infty$, thanks to the $m$ factor coming from the $a^{0}$ coupling to the fermion. Still, at the level of the linearly-realized theory, anomalies do not show up and the loop amplitude can be safely computed using a naive regularization procedure, even if we know from the fermion charges in Eq. (4) that $U(1)_{PQ}$ is anomalous.

2.2 Polar representation

Instead of the linear representation for the scalar field, there is another representation more suited to the circular geometry of the vacuum333Note that the polar representation still allows the pseudoscalar of having a small mass which could be treated as a perturbation as this is usually done in axion phenomenology.. The polar or exponential representation is

[TABLE]

Then, under the action of a $U(1)_{PQ}$ rotation of parameter $\theta$, the component fields have simple (though non-linear) transformation properties: the $\sigma$ field is constant, $\sigma^{0}\rightarrow\sigma^{0}$, while the axion is shifted by a constant, $a^{0}\rightarrow a^{0}+v\theta$. Importantly, and even if we are using the same notations $a^{0}$ and $\sigma^{0}$, these fields are not the same ones as in Eq. (5). In particular, the $a^{0}$ couples differently to fermions. Clearly, expanding the exponential, couplings of arbitrarily high dimensions are present, and the theory is no longer manifestly renormalizable.

As a next step, after plugging the polar representation of $\phi$ in the Lagrangian, we perform a reparametrization of the fermion fields to make them invariant under $U(1)_{PQ}$:

[TABLE]

This has three effects. First, the Goldstone field disappears from the Yukawa coupling. Second, the fermion kinetic terms induce derivative interactions,

[TABLE]

where one can recognize the fermionic current associated to the $U(1)_{PQ}$ symmetry. Third, the fermionic path integral measure is not invariant, and the Jacobian of the transformation has to be included as a new local interaction:

[TABLE]

The whole non-linear Lagrangian is thus

[TABLE]

Under this form, the Lagrangian is manifestly $U(1)_{PQ}$-symmetric but for the local Jacobian term $a^{0}F_{\mu\nu}\tilde{F}^{\mu\nu}$ since $\sigma^{0}$, $\psi_{L,R}$, and $A^{\mu}$ are invariant but $a^{0}$ is shifted by a constant. This permits to get rid of any pre-existing $\theta_{em}F_{\mu\nu}\tilde{F}^{\mu\nu}$ term in the Lagrangian, by shifting $a^{0}\rightarrow a^{0}+v\theta_{em}$. Of course, for an abelian theory, $F_{\mu\nu}\tilde{F}^{\mu\nu}$ is a total derivative that can be safely discarded, but this toy model can easily be adapted to the non-abelian case, and closely resembles the KSVZ invisible axion model [12] when applied to QCD.

Let us now compute $a^{0}\rightarrow\gamma\gamma$ within the non-linearly realized theory. Together with the local amplitude from $\delta\mathcal{L}_{\mathrm{Jac}}$, we must include the triangle graphs arising from $\delta\mathcal{L}_{\mathrm{Der}}$, which have both a vector and axial component. The vector current $V^{\mu}=\bar{\psi}\gamma^{\mu}\psi$ does not contribute since the photon couplings are also vectorial, and the corresponding triangle graph vanishes thanks to Furry’s theorem. The axial current $A^{\mu}=i\bar{\psi}\gamma^{\mu}\gamma_{5}\psi$, on the other hand, gives a non-zero contribution. We recognize the well-known $AVV$ triangle graph:

[TABLE]

Adopting a Pauli-Villars regularization, we add to this a second term with $m\rightarrow M$, and proceed with the calculation in $d=4$ dimensions. With the derivative from $\partial_{\mu}a^{0}$ amounting to taking the divergence, we find

[TABLE]

The regulator term does not decouple in the $M\rightarrow\infty$ limit but gives a finite local contribution. Interestingly, this contribution precisely cancel with the local term from $\delta\mathcal{L}_{\mathrm{Jac}}$, so altogether, we recover precisely the result obtained in the linear case, see Eq. (9):

[TABLE]

This exercise illustrates a number of important features:

• •

In the polar representation, the appearance of the anomalous local term $a^{0}F_{\mu\nu}\tilde{F}^{\mu\nu}$ is spurious. When computing specific amplitudes, it only serves to cancel out the anomalous term arising from the derivative interaction $\partial_{\mu}a^{0}\bar{\psi}\gamma^{\mu}\gamma_{5}\psi$. At the end, this is nothing but an application of the well-known axial current Ward identity:

[TABLE]

where $P=\bar{\psi}\gamma_{5}\psi$, $A^{\mu}=\bar{\psi}\gamma^{\mu}\gamma_{5}\psi$. The right-hand side corresponds to the amplitude in the linear representation, and the left-hand side to the two contributions arising in the polar representation. It is thus clear that it would be wrong to understand the $a^{0}\rightarrow\gamma\gamma$ coupling as induced by the anomaly. The anomalous interaction of Eq. (14) only arises because of the reparametrization in Eq. (12), and necessarily comes together with the appropriate derivative interactions.

• •

In the $m\rightarrow 0$ limit, the $a^{0}\rightarrow\gamma\gamma$ amplitude vanishes exactly,

[TABLE]

In the linear representation, this trivially follows from the vanishing of the coupling of $a^{0}$ to fermions, see Eq. (6). In the polar representation, it requires an exact cancellation of the local anomalous contribution with the non-local triangle amplitudes. Again, this can be understood in terms of the axial current Ward identity. In this limit, the spurious nature of the contact interaction is manifest.

• •

The reason why $a^{0}\rightarrow\gamma\gamma$ is often misinterpreted as induced by the anomaly can be understood looking at the $m\rightarrow\infty$ limit. Indeed, the contribution of the derivative term vanishes,

[TABLE]

as can be trivially seen in Eq. (17). Since then all that remains in the non-linear theory is the local term from $\delta\mathcal{L}_{\mathrm{Jac}}$, it necessarily corresponds to the contribution surviving in the $m\rightarrow\infty$ limit in the linear theory, Eq. (10). Yet, this is only a parametric correspondence, and certainly not a physical identification. The anomaly still cancels in $\mathcal{M}(a^{0}\rightarrow\gamma\gamma)_{\mathrm{Der}}+\mathcal{M}(a^{0}\rightarrow\gamma\gamma)_{\mathrm{Jac}}$, and the only surviving contribution is actually the first term of Eq. (17). Ultimately, the local anomalous term $\delta\mathcal{L}_{\mathrm{Jac}}$ in Eq. (14) is no more than a convenient book-keeping device tracking all the fields that have been integrated out [9, 10].

• •

Furry’s theorem together with the vector current Ward identity ensures the absence of any dependence on the free parameter $\alpha$ in the non-linear theory, in agreement with the absence of this parameter in the corresponding linear representation.

• •

As an aside, it should be clear that though the anomalous coupling in Eq. (14) is not affected by radiative corrections [13], neither is the anomalous part of the triangle graph in Eq. (17). So, the two exactly cancel at all orders, but the $a^{0}\rightarrow\gamma\gamma$ amplitude does get corrected by higher order effects since it is actually not induced by the anomaly. Further, note that the theory in Eq. (3) is obviously renormalizable, so radiative corrections can be calculated perturbatively using standard techniques. In this respect, if the photons are also coupled to some other fermions $\chi$, the two-loop process $a^{0}\rightarrow\gamma\gamma\rightarrow\bar{\chi}\chi$ is even UV finite in the linear representation [14]. This fact would clearly be difficult to guess using the polar representation, where the local $a^{0}\rightarrow\gamma\gamma$ vertex from $\delta\mathcal{L}_{\mathrm{Jac}}$ leads to a UV divergent diagram.

• •

As said earlier, the present toy model can be generalized to more complicated, and more realistic KSVZ-like axion models [12]. Consider for instance the SM, to which a gauge-singlet scalar $\phi$ and a set of vector-like fermions $\psi$ are added:

[TABLE]

with the same scalar potential as in Eq. (3). A priori, the covariant derivative can include all three SM interactions,

[TABLE]

where $T^{a}$ and $T^{i}$ are the $SU(3)_{C}$ and $SU(2)_{L}$ generators in the representation carried by $\psi$, and $Y$ its hypercharge. Then, if $yv\gg v_{EW}$ and $v\gg v_{EW}$, with $v_{EW}$ the electroweak vacuum expectation value, the effect of the heavy fermion is to induce (see Eq. (10))

[TABLE]

with the quadratic invariants $C_{C}^{\psi}\delta^{ab}=\operatorname*{Tr}[T^{a}T^{b}]$, $C_{L}^{\psi}\delta^{ij}=\operatorname*{Tr}[T^{i}T^{j}]$, and $C_{Y}^{\psi}=Y^{2}/4$, and $d_{C,L}^{\psi}$ the corresponding dimensions of the $SU(3)_{C}$ and $SU(2)_{L}$ representations. Yet, several points must be clear: (1) These couplings are not anomalous, but result from the appropriately generalized Eq. (9) in the $y\rightarrow\infty$ limit. (2) The integrated fermion must have vector couplings to gauge fields, otherwise the free parameter $\alpha$ corresponding to fermion number does not cancel in Eq. (14). We will see later on how to deal with chiral theories. (3) Even if vector-like, no mass term is allowed for $\psi$ because it would explicitly break the $U(1)_{PQ}$ symmetry. (4) The couplings in $\mathcal{L}_{\mathrm{Linear}}^{\mathrm{eff,KSVZ}}$ are not protected from radiative corrections since they are not anomalous. (5) Finally, the nature of the divergence arising when using this effective Lagrangian to compute e.g. the couplings of $a^{0}$ to SM fermions is clear. In the UV complete theory, Eq. (23), the $a^{0}\rightarrow\gamma\gamma$, $WW$, and $gg$ vertices are never local. Instead, the pseudoscalar triangle of Eq. (9) acts as a form factor and is sufficient to regulate the UV behavior of the vector boson loops, making them finite.

3 On the consistent use of anomalies

In the course of the previous section, the anomaly appeared several times, and there is a number of peculiarities that are worth detailing. First, consider the calculation of the local term $\delta\mathcal{L}_{\mathrm{Jac}}$, see Eq. (14). The simplest way to compute the Jacobian of the transformation is to compute a triangle graph, this time with three left-(or right-)handed currents (the momentum flow is defined as in Fig. 1):

[TABLE]

with $P_{L}=(1-\gamma_{5})/2$. Following Ref. [15], the divergences of this amplitude are calculated by first keeping track of the ambiguities in the loop momentum routing, and then from the surviving surface terms, and one finds:

[TABLE]

Yet, it is worth stressing that this choice $a=-b=1$ is not the one usually adopted, as it breaks the Bose symmetry of the amplitude under the exchange of the three currents. The consistent anomaly is that with $a=-b=1/3$, which restores the Bose symmetry as

[TABLE]

Importantly, this choice is not always explicit. For example, using the Pauli-Villars regularization, one immediately arrives at the Bose symmetric result (see Ref. [16]).

In terms of vector and axial currents, the $VVV$ and $VAA$ triangle amplitudes cancel out because of Furry’s theorem. This leaves the $AVV$ and $AAA$ diagrams, for which we obtain using the same method as in Ref. [15]:

[TABLE]

These results are valid in the massless limit for the fermions. As for the consistent anomaly, the usual conventions are $a=-b=1$ for the $AVV$ triangle to keep the vector currents conserved, and $a=-b=1/3$ for the $AAA$ triangle to enforce the Bose symmetry. The former case is the basis for Eq. (17). These conventions are automatically enforced when using simple regularization procedures like Pauli-Villars or dimensional regularization. However, for the following, the freedom to assign the anomaly in one of the vector currents or to explicitly break the Bose symmetry will be essential, so these regularization procedures are inadequate.

At the level of the $AVV$ and $AAA$ triangles, there is no obstruction to consider massive fermions. These contributions can be separated from the anomalous pieces, and do not depend on the chosen momentum routing. Actually, a simple calculation leads to

[TABLE]

and

[TABLE]

with the pseudoscalar triangle amplitudes (see Eq. (7))

[TABLE]

where $C_{i}(m^{2})=C_{i}(q_{1}^{2},q_{2}^{2},(q_{1}+q_{2})^{2},m^{2},m^{2},m^{2})$.

In some sense, the usual conventions are automatic for the mass terms: the vector currents are always conserved and Bose symmetry must be manifest for the $AAA$ triangle. Indeed, once the anomalous terms are separated, the rest obeys the naive Ward identities $\partial_{\mu}V^{\mu}=0$ and $\partial_{\mu}A^{\mu}=2imP$, with $P$ the pseudoscalar current. This means that the mass-dependent term must be absent from $-i(q_{1})_{\beta}\mathcal{T}_{AVV}^{\alpha\beta\gamma}$ and $-i(q_{2})_{\gamma}\mathcal{T}_{AVV}^{\alpha\beta\gamma}$, and it must be equal for all three divergences of $\mathcal{T}_{AAA}^{\alpha\beta\gamma}$. This picture is trivially confirmed using dispersion relations [17], see Fig. 2: cutting out a vector current necessarily gives zero since the intermediate fermions are on their mass shell and $\partial_{\alpha}\bar{\psi}\gamma^{\alpha}\psi=0$, while cutting out any one of the axial currents, the identity $\partial_{\alpha}\bar{\psi}\gamma^{\alpha}\gamma_{5}\psi=2im\bar{\psi}\gamma_{5}\psi$ must hold.

Finally, in the $m\rightarrow\infty$ limit, using $C_{0}\rightarrow-1/(2m^{2})$ and $C_{1}(m^{2})\rightarrow 1/(6m^{2})$, the various divergences are

[TABLE]

and

[TABLE]

Crucially for the following, the vanishing of all the divergences in Eq. (33) in the $m\rightarrow\infty$ limit is only true when the anomaly is carried by the axial current in the $AVV$ triangle, which implies the parametric choice $a=-b=1$. Similarly, the vanishing of all the divergences in Eq. (34) in the $m\rightarrow\infty$ limit is only true when the anomaly is carried equally by the three axial currents in the $AAA$ triangle, which implies the choice $a=-b=1/3$. But evidently, there is no reason for this to hold in general and we will see in the next section how this plays a role in the derivation of the axion couplings to electroweak gauge bosons.

4 The Peccei-Quinn axion and its couplings

Historically, the first axion model was proposed by Peccei and Quinn [1]. It is based on a Two Higgs Doublet Model (THDM), but with an enlarged symmetry. Denoting the Higgs doublets as

[TABLE]

the scalar potential is assumed to allow for the independent rephasing

[TABLE]

In the standard notation (see Ref. [18]), this forces the parameters $m_{12}$, $\lambda_{5,6,7}$ to vanish. One combination of these $U(1)$s can be identified with the hypercharge, since both doublets have the same gauge quantum numbers, while another combination is denoted $U(1)_{PQ}$. The vacuum expectation values

[TABLE]

break $SU(2)_{L}\otimes U(1)_{Y}\otimes U(1)_{PQ}\rightarrow U(1)_{em}$ spontaneously, where $v_{1}^{2}+v_{2}^{2}\equiv v^{2}\approx\left(246\,\text{GeV}\right)^{2}$ and $\tan\beta\equiv v_{1}/v_{2}$. The $U(1)_{Y}$ symmetry is local and its Goldstone boson is eaten up by the $Z$ boson, but the $U(1)_{PQ}$ is only global and the massless left over state is the axion. Indeed, assuming Yukawa couplings of Type II,

[TABLE]

the $U(1)_{PQ}$ symmetry is necessarily chiral and anomalous. Note that at this level, there is a lot of freedom in how to assign $U(1)_{1}\otimes U(1)_{2}$ charges for the fermions, and we will see later on how to (partially) fix them. Remark also that we did not include a phase difference between $v_{1}$ and $v_{2}$ in Eq. (37) because it can be removed by a $U(1)_{PQ}$ transformation. Actually, the $U(1)_{PQ}$ symmetry being anomalous, this rather moves the phase difference between $v_{1}$ and $v_{2}$ into the strong and electroweak $\theta$ terms. This is precisely the interplay that permits to solve the strong $CP$ puzzle once the $U(1)_{PQ}$ symmetry is spontaneously broken, since it entangles the phase difference with the axion field.

To further explore this simple model, we follow the same strategy as for the toy model of Section 2, and perform the analysis using either a linear or a polar representation for the scalar fields $\Phi_{1,2}$.

4.1 Linear representation

The most convenient basis for the scalar states is obtained by performing a rotation of angle $\beta$,

[TABLE]

so that the VEV is carried by one doublet only

[TABLE]

where only the SM Goldstone bosons have been exponentiated. The pseudoscalar $A^{0}$ is the massless Goldstone boson corresponding to the $U(1)_{PQ}$ breaking, and it is parametrized linearly. The scalar states $(\phi_{h},\phi_{H})$ are related to the mass eigenstates $(h^{0},H^{0})$ through a rotation of angle $\alpha$ which depends on the scalar potential parameters, and $H^{+}$ is the charged scalar boson. We note in passing that the charge assignments corresponding to $U(1)_{1}\otimes U(1)_{2}\leftrightarrow U(1)_{Y}\otimes U(1)_{PQ}$ thus require orthogonality of the Goldstone states $G^{0}$ and $A^{0}$, and cannot be guessed a priori from the $U(1)$ charges of the doublets in Eq. (36). Actually, since $A^{0}$ is a $\beta$-dependent combination of $\operatorname{Im}H_{1}^{0}$ and $\operatorname{Im}H_{2}^{0}$, the $U(1)_{PQ}$ charges of the doublets are also function of $\beta$. We will come back to this point later, and for now, we simply proceed by plugging Eq. (40) in the THDM Lagrangian.

Because $\Phi_{H}$ has no vacuum expectation value, the $A^{0}$ boson does not couple to pairs of gauge bosons at tree level. On the other hand, it does couple to fermions. From Eq. (38), we derive

[TABLE]

To reach this form, the mass terms are identified as $\sin\beta v\mathbf{Y}_{u}\equiv\sqrt{2}\mathbf{m}_{u}$ and $\cos\beta v\mathbf{Y}_{d,e}\equiv\sqrt{2}\mathbf{m}_{d,e}$ and the fermions are rotated to their mass basis. At the one-loop level, the $A^{0}$ then couples to pairs of gauge bosons through the $PVV$ and $PAA$ triangle graphs. Denoting the fermion currents coupled to gauge bosons as $\bar{\psi}_{f^{\prime}}(g_{V}^{f}\gamma^{\mu}-g_{A}^{f}\gamma^{\mu}\gamma_{5})\psi_{f}$ with, using the short-hand $c_{W}=\cos\theta_{W},s_{W}=\sin\theta_{W}$:

[TABLE]

where $SU(3)$ generators $T_{a}^{f}$ are normalized such that $\operatorname*{Tr}(T_{a}^{f}T_{b}^{f})=1/2\delta^{ab}$, it is a simple exercise to find (see Eq. (32)):

[TABLE]

where $\varepsilon(q_{1})_{\alpha}^{\ast}$ and $\varepsilon(q_{2})_{\beta}^{\ast}$ are the polarizations of $V_{1}$ and $V_{2}$, carrying respectively the momentum $q_{1}$ and $q_{2}$. To account for the absence of crossed diagrams for $W^{+}W^{-}$, we also define $\sigma=1/2$ for $V_{1,2}=W$, $\sigma=1$ for $V_{1,2}=g,\gamma,Z$. Also, the trace over color is understood for quarks, bringing a factor $\delta^{ab}/2$ for $V_{1,2}=g$, and $N_{C}$ for $V_{1,2}=W,Z,\gamma$. Finally, the loop functions in $\mathcal{T}_{PVV}^{\alpha\beta}(m_{f})$ and $\mathcal{T}_{AVV}^{\alpha\beta}(m_{f})$ are the same as in Eq. (32) for $V_{1,2}=g,\gamma,Z$, and generalized as $C_{i}(m_{f}^{2})=C_{i}(q_{1}^{2},q_{2}^{2},(q_{1}+q_{2})^{2},m_{f}^{2},m_{f^{\prime}}^{2},m_{f}^{2})$ with $f^{\prime}$ the $SU(2)$ partner of $f$ for $W^{+}W^{-}$. In that case, the CKM entries sum up to one when the three generations are included. The expression Eq. (43) for the decay amplitudes agrees with the standard result of Ref. [19].

A few limiting cases are interesting. First, remark that

[TABLE]

since the $A^{0}$ coupling to fermions is proportional to their masses. Second, in the opposite limit of infinite fermion masses, the loop functions behave as $C_{0}\rightarrow-1/(2m^{2})$ and $C_{1}(m^{2})\rightarrow 1/(6m^{2})$, so

[TABLE]

The only special case is that of the lepton contribution to $A^{0}\rightarrow W^{+}W^{-}$, since neutrinos remain massless. In that case, $C_{0}(m_{e}^{2})\rightarrow-1/m_{e}^{2}$ and $C_{1}(m_{e}^{2})\rightarrow 1/(2m_{e}^{2})$ in the $m_{e}\rightarrow\infty$ limit, and $\mathcal{T}_{PAA}^{\alpha\beta}$ cancels out. Plugging all these limits in the amplitudes, they take the simple form

[TABLE]

for some mass-independent coefficients $\mathcal{N}_{V_{1},V_{2}}$. Working out their values, to these amplitudes correspond the effective interactions

[TABLE]

with the coefficients (to simplify the comparison, we adopt the standard convention $x\equiv 1/\tan\beta$):

[TABLE]

Often, one interprets the non-decoupling of these amplitudes when $m_{u,d,e}\rightarrow\infty$ as an indirect manifestation of the underlying anomalies, but it is important to understand this statement correctly. The pseudoscalar triangles do decouple in that limit, $\mathcal{T}_{PVV}^{\alpha\beta}(m_{f})\rightarrow 0$ and $\mathcal{T}_{PAA}^{\alpha\beta}(m_{f})\rightarrow 0$ when $m_{f}\rightarrow\infty$, but this is compensated by the mass-dependent couplings of the $A^{0}$ to fermions. The fact that this compensation works, leaving a constant remainder in the $m_{f}\rightarrow\infty$, is guaranteed by the anomalous axial Ward identity of Eq. (20). Yet, there is no anomalous contribution of any kind to the $A^{0}\rightarrow V_{1}V_{2}$ amplitudes, and those are driven entirely by non-local and well-behaved fermion loops. The situation is quite similar for the loop-induced Higgs decay $h\rightarrow\gamma\gamma$ in the SM, whose non-decoupling behavior in the $m_{u,d,e}\rightarrow\infty$ limit is interpreted as a manifestation of the trace anomaly [20]. Yet, there is of course no direct contribution of that anomaly to $h\rightarrow\gamma\gamma$, which can be safely computed perturbatively.

4.2 Polar representation

The Higgs boson fields are written in the polar representation555Note that one cannot start from the rotated basis Eq. (40) to exponentiate the Goldstone field $A^{0}$, because such a representation would not be canonical for $\Phi_{H}$ since it has no vacuum expectation value [21]. as (see Eq. (11)):

[TABLE]

where $v_{1}=v\sin\beta$, $v_{2}=v\cos\beta$. The $\eta_{1}$ and $\eta_{2}$ states are the Goldstone bosons associated to the $U(1)_{1}$ and $U(1)_{2}$ symmetries. They are related to the physical Goldstone bosons $a^{0}$ and $G^{0}$ by the same rotation as in Eq. (39):

[TABLE]

Note that the $G^{0}$ components of $\Phi_{1}$ and $\Phi_{2}$ are the same since they have the same hypercharge, and the $a^{0}$ and $G^{0}$ states are orthogonal to each other. Now, the $G^{0}$ can be gauged away via a $U(1)_{Y}$ rotation, leaving

[TABLE]

where we introduce $x=v_{2}/v_{1}=1/\tan\beta$ to stick to common conventions. The remaining neutral Goldstone field, here denoted $a^{0}$ instead of $A^{0}$, is the axion. Under this form, the true $U(1)_{PQ}$ charges of the Higgs doublets are finally apparent, since under such a transformation, $a^{0}\rightarrow a^{0}+v\theta$. Then, the Yukawa Lagrangian can also be made invariant under $U(1)_{PQ}$ provided the fermion charges are set as

[TABLE]

with

[TABLE]

where $\alpha$ and $\beta$ are free parameters, corresponding to the conservation of baryon and lepton numbers $\mathcal{B}$ and $\mathcal{L}$.

Following the same steps as in Section 2, we now perform a reparametrization of the fermion fields to force them to be invariant under $U(1)_{PQ}$. This is achieved by the field-dependent rotations $\psi\rightarrow\exp(i\chi_{\psi}a^{0}(x)/v)\psi$. The axion field disappears from the Yukawa Lagrangian, but reappears in derivative interactions from the fermion kinetic terms

[TABLE]

with

[TABLE]

The vector current couplings are $\alpha$ and $\beta$-dependent, since those parameters reflect the conservation of $\mathcal{B}$ and $\mathcal{L}$, which are vectorial [9]. The $\beta$ dependence of $\chi_{A}^{\nu}$ is related to the peculiar nature of neutrinos, which are kept massless. As we will see in Section 4.3.4, it disappears if right handed neutrinos $\nu_{R}$ together with the $\mathcal{L}$-invariant Yukawa interaction $\bar{\nu}_{R}Y_{\nu}\ell_{L}\Phi_{1}$ are added.

In addition, given the $U(1)_{PQ}$ charges in Eq. (52), the fermion reparametrization is anomalous and its Jacobian has to be included

[TABLE]

Given the SM fermion quantum numbers, the coefficients are

[TABLE]

where $d_{C,L}^{\psi}$, $C_{C,L}^{\psi}$ are the $SU(3)_{C}$ and $SU(2)_{L}$ dimensions and quadratic Casimir invariant of the representation carried by the field $\psi$, respectively, and by extension, $C_{Y}^{\psi}=Y(\psi)^{2}/4$ with $Y(q_{L},u_{R},d_{R},\ell_{L},e_{R})=(1/3,4/3,-2/3,-1,-2)$. These expressions correspond to Eq. (27), choosing $a=-b=1$ to reflect the absence of anomalies for the SM gauge symmetries (so no Bose symmetry here).

The parameters $\alpha$ and $\beta$ enter in a combination reflecting the anomaly in the $\mathcal{B}+\mathcal{L}$ current. Interestingly, neither the QED nor the QCD coefficients depend on them, but those to electroweak gauge bosons do. To see this, it sufficies to plug $W_{\mu}^{3}=c_{W}Z_{\mu}+s_{W}A_{\mu}$, $B_{\mu}=-s_{W}Z_{\mu}+c_{W}A_{\mu}$ and $g^{\prime}c_{W}=e=gs_{W}$ in Eq. (56):

[TABLE]

and the free parameters cancel out in

[TABLE]

Concerning the strong $CP$ puzzle, the polar representation makes the shift symmetry explicit. Under a $U(1)_{PQ}$ transformation of parameter $\theta$, all the fields are invariant but for $a^{0}\rightarrow a^{0}+v\theta$. So, if there is a $\theta_{QCD}G_{\mu\nu}^{a}\tilde{G}^{a,\mu\nu}$ term in the gauge sector, it can be eliminated by an adequate shift of the $a^{0}$ field. As said in the introduction, what happens in reality is that one is forced by non-perturbative QCD effects to shift the $a^{0}$ field by precisely that amount [1]. Once the $a^{0}$ shift is set by QCD, the electroweak $\theta_{L}W_{\mu\nu}^{i}\tilde{W}^{i,\mu\nu}$ can still be removed by an appropriate choice of $N_{C}\alpha+\beta$. All that remains then is the abelian $\theta_{Y}B_{\mu\nu}\tilde{B}^{\mu\nu}$, which is an inoccuous total derivative. We thus recover the well-known fact that the electroweak $\theta$ term can be rotated away thanks to the $\mathcal{B}$ and $\mathcal{L}$ invariance of the model.

4.3 Matching the polar and linear representations

The physics must not depend on the representation chosen for the scalar fields or on the parametrization of the fermion fields. So, the one-loop $A^{0}\rightarrow V_{1}V_{2}$ amplitudes of Eq. (43) must match the $a^{0}\rightarrow V_{1}V_{2}$ amplitudes computed using the local terms $\delta\mathcal{L}_{\text{{Jac}}}$ together with the triangle graphs arising from the derivative interactions in $\delta\mathcal{L}_{\text{{Der}}}$ (these two contributions are of the same order in the coupling constants). At this stage, it is immediately clear that the triangle graphs will play a crucial role. Indeed, taken alone, the contributions from the anomalous local terms of $\delta\mathcal{L}_{\text{{Jac}}}$ match those found in the linear case in the $m_{u,d,e}\rightarrow\infty$ limit, Eq. (47), only for the $\gamma\gamma$ and $gg$ final states. For the others, there is no dependence on $\alpha$ and $\beta$ in Eq. (47), but there is one through $\mathcal{N}_{L}$ in $\delta\mathcal{L}_{\text{{Jac}}}$. Even worse, no choice of these parameters could make the two results compatible because the relative strengths of the $\gamma\gamma$, $\gamma Z$, $ZZ$, and $W^{+}W^{-}$ decays are irremediably different. As we will now detail, it is only once the anomalous contributions from $\delta\mathcal{L}_{\text{{Jac}}}$ precisely cancel out with those hidden in the triangle graphs that the $a^{0}\rightarrow V_{1}V_{2}$ amplitudes match the $A^{0}\rightarrow V_{1}V_{2}$ amplitudes.

Throughout this section, the momentum flow is defined as $a^{0}(q_{1}+q_{2})\rightarrow V_{1}(q_{1},\alpha)V_{2}(q_{2},\beta)$, with $V_{1,2}=g,\gamma,Z,W$, and neither the final gauge bosons nor the initial axion are necessarily on-shell. The final amplitudes thus depend on $q_{1}^{2}$, $q_{2}^{2}$, and $(q_{1}+q_{2})^{2}$, and are to be compared to those obtained in the linear case in Eq. (43).

4.3.1 The $a^{0}\rightarrow\gamma\gamma$ and $a^{0}\rightarrow gg$ decays

For the decays into photons or gluons, only the $AVV$ triangle contributes since $SU(3)$ and $U(1)_{em}$ are vector symmetries. Then, the situation is exactly the same as in the toy model of Section 2. Setting $a=-b=1$ in Eq. (30), we find

[TABLE]

and

[TABLE]

where $d_{C}^{u,d}=N_{C}$, $d_{C}^{e,\nu}=1$, $Q^{f}$ is the electric charge of $f$, and $C_{C}^{u,d}=1/2$, $C_{C}^{e,\nu}=0$. The amplitudes in the polar representation thus trivially match that in the linear case, see Fig. 3,

[TABLE]

In other words, the anomalous contact interactions do cancel out systematically with the anomalous part of the triangle graphs.

As for the toy model, the polar representation is thus interesting only to make the shift symmetry manifest, and because the contact $a^{0}\gamma\gamma$ and $a^{0}gg$ interactions read off $\delta\mathcal{L}_{\text{{Jac}}}$ are reliable book-keeping of the effects of heavy fermions. Specifically, $\mathcal{M}(a^{0}\rightarrow\gamma\gamma,gg)_{\text{{Der}}}^{AVV}\overset{m\rightarrow\infty}{=}0$ implies that $\mathcal{M}(A^{0}\overset{}{\rightarrow}\gamma\gamma,gg)_{\mathrm{Linear}}\overset{m\rightarrow\infty}{=}\mathcal{M}(a^{0}\overset{}{\rightarrow}\gamma\gamma,gg)_{\text{{Jac}}}$. Finally, remark that the cancellation of the local anomalous terms ensures $\mathcal{M}(a^{0}\overset{}{\rightarrow}\gamma\gamma,gg)_{{}_{\mathrm{Polar}}}=0$ in the $m_{u,d,e}\rightarrow 0$ limit. So, though interpreting the axion coupling to photons or gluons as induced by the anomaly is incorrect, this misidentification does not lead to serious consequences for those final states. For heavy fermions, the coupling to gluons is tuned by $\mathcal{N}_{C}$, and that to photons by $\mathcal{N}_{em}$, and their ratio, when restricted to quarks, give back the usual $\mathcal{N}_{em}^{q}/\mathcal{N}_{C}=8/3$. However, as we will see in the next subsection, interpreting the axion coupling involving at least one electroweak gauge boson as induced by the anomaly is not only wrong in principle but also leads to incorrect couplings.

4.3.2 The $a^{0}\rightarrow\gamma Z$ decay

For electroweak gauge bosons in the final state, the situation is less simple. Consider first the $\gamma Z$ final state. The derivative interactions induce again only the $AVV$ triangle graphs since the photon coupling is vectorial. However, this time we have two possible contributions, depending on which current is carrying the anomaly, see Fig. 4. First, there are the usual $A(\partial_{\mu}a^{0})-V(\gamma)-V(Z)$ triangles, for which the anomaly is in the axial current. Using Eq. (30) with $a=-b=1$, we find

[TABLE]

where $\varepsilon^{\varepsilon_{1}^{\ast}\varepsilon_{2}^{\ast}q_{1}q_{2}}\equiv\varepsilon^{\alpha\beta\mu\nu}\varepsilon(q_{1})_{\alpha}^{\ast}\varepsilon(q_{2})_{\beta}^{\ast}q_{1\mu}q_{2\nu}$. But contrary to the $\gamma\gamma$ and $gg$ final state, there are now new contributions from the $V(\partial_{\mu}a^{0})-V(\gamma)-A(Z)$ triangles, with the axion vector couplings of Eq. (54). Using again Eq. (30) but this time with $a=b=1$ to preserve $SU(2)_{L}\otimes U(1)_{Y}$, and noting that $-i(q_{1})_{\gamma}\mathcal{T}_{AVV}^{\alpha\gamma\beta}=-i(q_{1}+q_{2})_{\gamma}\mathcal{T}_{VAV}^{\gamma\alpha\beta}$,

[TABLE]

As explained in Section 3, this contribution is free of any mass-dependent term because the naive vector current conservation must hold, up to the anomaly.

These two contributions from the derivative interactions combine with the Jacobian term from Eq. (60),

[TABLE]

to give the total decay amplitude in the polar representation:

[TABLE]

as it should. We can now understand why $\mathcal{M}(a^{0}\overset{}{\rightarrow}\gamma Z)_{\text{{Jac}}}$ does not match $\mathcal{M}(A^{0}\overset{}{\rightarrow}\gamma Z)_{\mathrm{Linear}}$ in the $m_{u,d,e}\rightarrow\infty$ limit. Indeed, while $\mathcal{M}(a^{0}\rightarrow\gamma Z)_{\text{{Der}}}^{AVV}$ vanishes in the $m_{u,d,e}\rightarrow\infty$ limit, $\mathcal{M}(a^{0}\overset{}{\rightarrow}\gamma Z)_{\text{{Der}}}^{VAV}$ obviously does not since it is independent of $m_{u,d,e}$. In that limit, we should write

[TABLE]

in which the parameters $\alpha$ and $\beta$ cancel out. So, for a chiral gauge theory, the local terms coming from $\delta\mathcal{L}_{\text{{Jac}}}$ are no longer reliable book-keeping of the effect of heavy fermions, because part of the anomaly is hidden in the $VAV$ triangle.

Actually, the reason why the $VAV$ triangle plays a role can be understood directly from the fermion reparametrization, given the charges in Eq. (52). When $\alpha$ and/or $\beta$ are different from zero, the fermionic current associated to $U(1)_{PQ}$ has a component aligned with baryon and/or lepton number, respectively. The fermionic reparametrization thus generates the anomalous $\mathcal{B}+\mathcal{L}$ interactions: those correspond to the $N_{C}\alpha+\beta$ terms of $\mathcal{N}_{L}$ and $\mathcal{N}_{Y}$ in Eq. (56). But as for the toy model, these interactions are spurious and must cancel with the anomalies present in the triangle graphs induced by the derivative interactions $\delta\mathcal{L}_{\text{{Der}}}$. This must necessarily come from a breakdown of the vector Ward identity for the $VAV$ triangle since $\mathcal{B}$ and $\mathcal{L}$ are purely vectorial symmetries [22].

4.3.3 The $a^{0}\rightarrow ZZ$ and $a^{0}\rightarrow W^{+}W^{-}$ decays

Turning first to the $a^{0}\rightarrow ZZ$ amplitude, the derivative interactions induce new types of diagrams: the $AAA$ triangles and graphs with neutrinos circulating in the loop, see Fig. 5.

Proceeding with the calculation, the three contributions from the derivative terms are

[TABLE]

where $\mathcal{M}(A^{0}\overset{}{\rightarrow}ZZ)_{\mathrm{Linear}}^{PVV}$ is the part of Eq. (43) proportional to $\mathcal{T}_{PVV}^{\alpha\beta}$. The neutrino contribution is absent from that term since the loop function in Eq. (32) vanishes for massless fermions. The next piece comes from the divergence of the axionic vector current,

[TABLE]

There is no mass-dependent terms here because of the vector Ward identity. Finally, the $AAA$ triangle gives, using Eq. (31) with $a=-b=1$

[TABLE]

where $\mathcal{M}(A^{0}\overset{}{\rightarrow}ZZ)_{\mathrm{Linear}}^{PAA}$ is the part of Eq. (43) proportional to $\mathcal{T}_{PAA}^{\alpha\beta}$. Again, the neutrino contribution to that part vanishes. Combining the three contributions of the derivative interactions with that of the Jacobian

[TABLE]

all the local anomalous contributions cancel exactly, leaving:

[TABLE]

Though the anomalous terms again sum up to zero, they are now completely entangled in the three types of triangle graphs, which all depend on the free parameters $\alpha$ and $\beta$. As a result, none of the derivative amplitudes vanish in the $m_{u,d,e}\rightarrow\infty$ limit. Specifically, because neutrinos remain massless,

[TABLE]

where $\beta$ comes from the neutrino coupling $\chi_{A}^{\nu}$ in Eq. (55). Also,

[TABLE]

because in addition to massless neutrinos, the mass term respects Bose symmetry but not the anomaly which resides only in the axion axial current (there is a relative factor of 1/3, see Eq. (34) with $a=-b=1$). Finally, $\mathcal{M}(a^{0}\rightarrow ZZ)_{\text{{Der}}}^{VAV}$ does not vanish since it is mass-independent. Because of all this, $\mathcal{M}(a^{0}\overset{}{\rightarrow}ZZ)_{\text{{Jac}}}$ is clearly not identical to $\mathcal{M}(A^{0}\overset{}{\rightarrow}ZZ)_{\mathrm{Linear}}$ in the $m_{u,d,e}\rightarrow\infty$ limit, and cannot be used as an estimate of the electroweak couplings of the axion.

The $a^{0}\rightarrow W^{+}W^{-}$ amplitude shares many similarities with that for $a^{0}\rightarrow ZZ$, so we will be very brief. Following the same steps, the various contributions to the decay in the polar representation are

[TABLE]

Note well that it is $\mathcal{N}_{1}$ that occurs in these amplitudes, and not the $\mathcal{N}_{2}$ coefficients to which $\mathcal{M}(A^{0}\overset{}{\rightarrow}W^{+}W^{-})_{\mathrm{Linear}}$ tends to in the $m_{u,d,e}\rightarrow\infty$ limit. Still, the sum of these four contributions reproduces again the linear result, $\mathcal{M}(a^{0}\rightarrow W^{+}W^{-})_{\text{{Polar}}}=\mathcal{M}(A^{0}\rightarrow W^{+}W^{-})_{\mathrm{Linear}}$. As none of the triangle amplitudes vanish in the $m_{u,d,e}\rightarrow\infty$ limit, $\mathcal{M}(a^{0}\rightarrow W^{+}W^{-})_{\text{{Jac}}}$ is clearly not identical to $\mathcal{M}(A^{0}\rightarrow W^{+}W^{-})_{\mathrm{Linear}}$ in the $m_{u,d,e}\rightarrow\infty$ limit.

4.3.4 Impact of heavy neutrinos

The effective interactions in Eq. (47) involve five $\mathcal{N}_{i}$ coefficients, but this is partly due to the non-decoupling of neutrinos. To see this, let us add right handed neutrinos together with the Yukawa interaction $\bar{\nu}_{R}Y_{\nu}\ell_{L}\Phi_{1}$. Then, the effective amplitude in the linear representation is trivially obtained by extending the sum in Eq. (43) to $f=u,d,e,\nu$. If, for the sake of the argument, we now take the infinite mass limit for all fermions, $m_{u,d,e,\nu}\rightarrow\infty$, the effective interactions in Eq. (47) are recovered but with

[TABLE]

By contrast, it is amusing to remark that in the polar representation, the anomalous interactions of Eq. (56) are not modified by the presence of right-handed neutrinos since they are neutral under the whole SM gauge group. Yet, right-handed neutrinos do modify the derivative interactions in Eq. (55) to $\chi_{V}^{\nu}=2\beta+x$ and $\chi_{A}^{\nu}=x$. The various triangle amplitudes are then obtained from those quoted in the previous section by replacing $\mathcal{N}_{1}\rightarrow\mathcal{N}_{\nu}$ and discarding the $\beta$ terms. While $\mathcal{M}(a^{0}\rightarrow ZZ,W^{+}W^{-})_{\text{{Der}}}^{AVV}$ now vanish in the $m_{u,d,e,\nu}\rightarrow\infty$ limit, this does not change our conclusions and $\mathcal{M}(a^{0}\rightarrow ZZ,W^{+}W^{-})_{\text{{Jac}}}$ remain different from $\mathcal{M}(A^{0}\rightarrow ZZ,W^{+}W^{-})_{\mathrm{Linear}}$ in that limit.

As a final note, it is interesting to remark that once certain that $\alpha$ and $\beta$ disappear from the physical amplitudes, one can chose to fix them as one wishes. In particular, as said before, they can be set to eliminate the electroweak $\theta$ term. Alternatively, one can set $\alpha=\beta=0$, thereby getting rid of the $a^{0}W_{\mu\nu}^{i}\tilde{W}^{i,\mu\nu}$ coupling altogether. As we have seen, this does not forbids $a^{0}\rightarrow W^{+}W^{-}$ since the amplitude is independent of $\alpha$ and $\beta$. Finally, this freedom can also be used to allow for the presence of a lepton number violating Majorana mass term, as well as an effective dimension-five $(\bar{\ell}_{L}^{C}H_{1})(\ell_{L}H_{1})/\Lambda$ operator. We will not explore such settings here, see Ref. [23] for studies involving axions together with neutrino Majorana mass terms.

4.4 Application to the DFSZ axion model

The Peccei-Quinn model is ruled out experimentally because the coupling of the axion to fermions, tuned by the electroweak vacuum expectation value, is too large. One way to extend the model and render the axion invisible was proposed not long after the original idea and is called the DFSZ axion model [24]. Our goal here is to derive the electroweak couplings for that model.

The starting point is to extend the THDM by adding a complex scalar field $\phi$, singlet under the SM gauge symmetries, and constrain the scalar potential to

[TABLE]

In this way, the whole Lagrangian remains invariant under independent rephasing of the doublet provided the singlet has the charge:

[TABLE]

The Type II Yukawa couplings of Eq. (38) then imply chiral charges for the fermions when $\alpha_{1}\neq\alpha_{2}$, and part of the $U(1)_{1}\otimes U(1)_{2}$ symmetry is anomalous.

Let us take the linear representation (compare with Eq. (40))

[TABLE]

where $\Phi_{h,H}$ are given in terms of $\Phi_{1,2}$ in Eq. (39). Plugging these representations in the potential, and extracting the quadratic terms, the $G^{0}$ field is immediately massless and corresponds to the Would-be Goldstone of the $Z$ boson. On the other hand, the $\phi_{A}$ and $\pi_{S}$ fields mix,

[TABLE]

with $M_{A^{0}}^{2}=0$ and $M_{\pi^{0}}^{2}=2\lambda_{12}(v^{2}\sin 2\beta+v_{s}^{2}/\sin 2\beta)$. The $A^{0}$ field is the axion, and the $\pi^{0}$ field is a heavy pseudoscalar Higgs boson. While $A^{0}\rightarrow\phi_{A}$ when $v_{s}\rightarrow 0$, reproducing the THDM discussed previously, the opposite holds when $v_{s}\gg v$ [24]:

[TABLE]

Then, since $\pi_{S}$ is not coupled to fermions, it is easy to see that the $A^{0}\rightarrow V_{1}V_{2}$ amplitudes are simply rescaled by $v\sin 2\beta/v_{s}$, i.e.,

[TABLE]

with the same coefficients as in Eq. (48). In particular, the axion couplings to $\gamma Z$, $ZZ$, and $W^{+}W^{-}$ are not aligned with the purely anomalous interactions one derives looking only at the $PQ$ charges of the fermion fields.

4.5 Anomaly cancellation in generic axion models

In the previous sections, we showed how the anomalies present in the non-linearly realized axion models necessarily cancel out, since they are spuriously generated by reparametrizing an anomaly-free model. Here, we want to follow the opposite route. Starting from a Non Renormalizable (NR) effective Lagrangian valid above the electroweak scale:

[TABLE]

where $v_{a}$ is some high scale, $\mathcal{N}_{C,L,Y}^{eff}$ as well as $\chi_{V,A}^{f}$ some a priori free Wilson coefficients, we want to show that whenever the anomalies cancel, this model collapses to a simple renormalizable theory.

To do this, starting from $\mathcal{L}_{NR}^{\mathrm{eff}}$, we require the anomalous interactions $\mathcal{N}_{C,L,Y}^{eff}$ to cancel with the local anomalous terms present in the $AVV$, $VAV$, and $AAA$ triangle graphs. Imposing this for the five channels $a^{0}\rightarrow V_{1}V_{2}$, $V_{1,2}=g,\gamma,Z,W$, four independent constraints are derived:

[TABLE]

together with

[TABLE]

Once these conditions are imposed, all the $a^{0}\rightarrow V_{1}V_{2}$ processes are entirely driven by the non-anomalous parts of the $AVV$ and $AAA$ triangle graphs (see Figs. 3, 4, and 5). The $VAV$ triangles disappear since they have no mass-dependent terms. But then, these non-anomalous amplitudes satisfy the naive Ward identities, and are equivalently obtained from $PVV$ and $PAA$ triangles, provided

[TABLE]

Phenomenologically, $\mathcal{L}_{R}^{\mathrm{eff}}$ and $\mathcal{L}_{NR}^{\mathrm{eff}}$ are indistinguishable once anomaly cancellation is enforced, even if $\mathcal{L}_{R}^{\mathrm{eff}}$ has far fewer physical parameters. This can be understood as follows. Out of the initial ten parameters in $\mathcal{L}_{NR}^{\mathrm{eff}}$ for massless neutrino ($\chi_{V}^{\nu}=-\chi_{A}^{\nu}$), four are eliminated by the anomaly cancellation constraints. Then, the naive Ward identities are valid again. The three $\chi_{V}^{u,d,e}$ parameters disappear thanks to $\partial_{\mu}V^{\mu}=0$, and $\partial_{\mu}A^{\mu}=2imP$ leaves only the three $\chi_{P}^{f}=\chi_{A}^{f}$ parameters as truly physical.

Now, starting from $\mathcal{L}_{R}^{\mathrm{eff}}$, all the $a^{0}\rightarrow V_{1}V_{2}$ decay amplitudes are easily computed, and given by Eq. (43). In the $m_{u,d,e}\rightarrow\infty$ limit, keeping neutrino massless, these amplitudes match onto the effective couplings in Eq. (47), with

[TABLE]

In this limit, we thus do not recover the manifestly $SU(2)_{L}\otimes U(1)_{Y}$ invariant pattern of anomalous interactions of $\mathcal{L}_{NR}^{\mathrm{eff}}$. Finally, note that if one takes $\chi_{A}^{u}=\chi_{P}^{u}=x$, $\chi_{A}^{d}=\chi_{P}^{d}=\chi_{A}^{e}=\chi_{P}^{e}=1/x$, the results in Eq. (48) are recovered. This shows that $\mathcal{L}_{R}^{\mathrm{eff}}$ is the most generic implementation of the PQ (or DFSZ) axion model, where the $U(1)_{PQ}$ charges of the SM fermions can take any value.

5 Conclusion and perspective

In this paper, we derived the couplings of axions to gauge bosons. We started from a simplified toy model, and then extended the discussions to the KSVZ axion [12], the original PQ axion [1], and the DFSZ axion [24]. In each model, our strategy has been to match the axion decay modes computed using either a linear or a polar representation for the scalar field breaking the $U(1)_{PQ}$ symmetry. In this way, we were able to unravel the physics at play, and uncovered a number of interesting features:

• •

In the linear representation, the axion couplings to gauge bosons are not induced by the anomaly, but by non-anomalous pseudoscalar triangle graphs with fermions circulating in the loop. Though this has to our knowledge not been exploited before, these amplitudes can be identified with those of the pseudoscalar Higgs to gauge bosons calculated in the THDM, which have been known for a long time [19]. In the $m_{u,d,e}\rightarrow\infty$ limit, they match onto the effective interactions

[TABLE]

with the coefficients shown in the first column of Table 1. In the opposite limit $m_{u,d,e}\rightarrow 0$, all these amplitudes vanish since the axion couplings to fermions are proportional to the fermion masses.

• •

In the context of axion models, it is customary to adopt a polar representation for the scalar fields. The pseudoscalar axion couplings to fermions are replaced by contact anomalous interactions to the gauge bosons and axion derivative interactions to the fermions. All these interactions are entirely fixed in terms of the assigned $PQ$ charges of the SM fermions. As shown in Table 1, this parametrization is particularly convenient for $a^{0}\rightarrow\gamma\gamma$ and $a^{0}\rightarrow gg$ because the derivative interactions do not contribute in the $m_{u,d,e}\rightarrow\infty$ limit [9, 10], and thus the strengths of the $a^{0}\rightarrow\gamma\gamma$ and $a^{0}\rightarrow gg$ processes can be read off the anomalous couplings $a^{0}F_{\mu\nu}\tilde{F}^{\mu\nu}$ and $a^{0}G_{\mu\nu}^{a}\tilde{G}^{a,\mu\nu}$. By contrast, and contrary to what is conjectured in Ref. [9], the central result of this paper is that for chiral gauge theories in which chiral fermions have non-trivial $PQ$ charges, the strengths of the $a^{0}\rightarrow\gamma Z$, $ZZ$, and $W^{+}W^{-}$ processes do not match the anomalous couplings, even in the $m_{u,d,e}\rightarrow\infty$ limit, as evident comparing the first and second columns of Table 1. They thus cannot be encoded into the $SU(2)_{L}\otimes U(1)_{Y}$ invariant effective interactions a^{0}B_{\mu\nu}\tilde{B}^{\mu\nu}\and $a^{0}W_{\mu\nu}^{i}\tilde{W}^{i,\mu\nu}$.

• •

The main reason for this mismatch is due to the presence of triangle graphs arising from the derivative interactions $\partial_{\mu}a^{0}\bar{\psi}\gamma^{\mu}\gamma_{5}\psi$ and $\partial_{\mu}a^{0}\bar{\psi}\gamma^{\mu}\psi$, which do not vanish in the $m_{u,d,e}\rightarrow\infty$ limit for chiral gauge theories. Importantly, even the vector couplings play a role since the anomalous breaking of the axionic vector current conservation enters through the $VAV$ triangle graphs (see Figs. 4 and 5). One important result of this paper is the proof that once all these triangle contributions (last three columns of Table 1) are summed with the local anomalous amplitudes (second column of Table 1), the THDM results are recovered. Without surprise, physical observables do not depend on the chosen parametrization, and this further confirms that the THDM calculation is correct.

• •

Our analysis also proves that the axion to gauge boson processes are not induced by the anomaly. This is obvious in the linear representation since the fermion loops driving these processes are not anomalous. In the polar representation, what happens is that the local anomalous interactions precisely cancel with the anomalies present in triangle graphs induced by the vector and axial derivative interactions, leaving the non-anomalous pseudoscalar triangles of the THDM as the only surviving mechanism driving all the axion to gauge bosons processes. Ultimately, this precise cancellation is ensured by the anomalous vector and axial Ward identities. It is required because all the amplitudes must vanish in the $m_{u,d,e,\nu}\rightarrow 0$ limit, but the anomalous contributions are independent of the fermion mass666This interpretation is compatible with the matching observed in Table 1 in the $m_{u,d,e}\rightarrow\infty$ limit for $a^{0}\rightarrow\gamma\gamma,gg$. Indeed, in the polar representation, these amplitudes receive two contributions, $\mathcal{M}_{ano}$ from the local anomalous term, and $\mathcal{M}_{der}$ from the triangle amplitudes built from $\partial_{\mu}a^{0}\bar{\psi}\gamma^{\mu}\gamma_{5}\psi$. But in this case, the axial Ward identity translates as $\mathcal{M}_{der}=\mathcal{M}_{lin}-\mathcal{M}_{ano}$, with $\mathcal{M}_{lin}$ the THDM amplitude. So, even if parametrically, $\mathcal{M}_{lin}=\mathcal{M}_{ano}$ and $\mathcal{M}_{der}=0$ when $m_{u,d,e}\rightarrow\infty$, $a^{0}\rightarrow\gamma\gamma,gg$ are not induced by the anomaly since it cancels out in $\mathcal{M}_{ano}+\mathcal{M}_{der}$..

• •

In all axion models where chiral fermions have $PQ$ charges but conserve baryon and lepton numbers, both the vector couplings $\partial_{\mu}a^{0}\bar{\psi}\gamma^{\mu}\psi$ and the $a^{0}W_{\mu\nu}^{i}\tilde{W}^{i,\mu\nu}$, $a^{0}B_{\mu\nu}\tilde{B}^{\mu\nu}$ anomalous contact interactions depend on some free parameters [9]. In Table 1, these free parameters enter into $\mathcal{N}_{L}$, and are related to the freedom to choose the $PQ$ charge of left-handed fermions (there are also explicit occurrence of the free parameter $\beta$, but this is an artifact of keeping the neutrino massless). It is only once the failure of the naive vector Ward identity in the $VAV$ triangle graphs is properly accounted for that these free parameters cancel out, as they should [22]. To understand how this cancellation is connected with the conservation of baryon and lepton numbers $\mathcal{B}$ and $\mathcal{L}$, notice first that the freedom to choose the $PQ$ charges of the fermions is actually due to the invariance under $U(1)_{\mathcal{B}}$ and $U(1)_{\mathcal{L}}$. In general, the $U(1)_{PQ}$ symmetry has two arbitrary components aligned with $U(1)_{\mathcal{B}}$ and $U(1)_{\mathcal{L}}$. But then, these components translate into derivative couplings of the axion to the $\mathcal{B}$ and $\mathcal{L}$ fermionic vectorial currents together with a local $a^{0}W_{\mu\nu}^{i}\tilde{W}^{i,\mu\nu}$ coupling from the $\mathcal{B}+\mathcal{L}$ anomaly. Clearly, all these couplings are spurious when the axion is a Goldstone boson living in a $\mathcal{B}$ and $\mathcal{L}$-invariant vacuum. Since $\mathcal{B}$ and $\mathcal{L}$ are not spontaneously broken, $\mathcal{N}_{L}$ must systematically drop out of physical observables, and Table 1 shows that this indeed occurs.

• •

In the linear representation, the calculation is done in the electroweak broken phase throughout. In the $m_{u,d,e}\rightarrow\infty$ limit, the amplitudes do not match onto $SU(2)_{L}\otimes U(1)_{Y}$ invariant operators, as shown in Table 1. Though this is somewhat expected, what is not obvious is the very peculiar way in which the $SU(2)_{L}\otimes U(1)_{Y}$ symmetry breaking seeps in. Indeed, in the non-linear representation, the anomalous interactions and the derivative interactions are manifestly $SU(2)_{L}\otimes U(1)_{Y}$ invariant once the SM fermions are assigned $PQ$ charges in an $SU(2)_{L}\otimes U(1)_{Y}$ invariant way. Further, evidently, the anomalies present in the triangle graphs built on the derivative interactions do not break $SU(2)_{L}\otimes U(1)_{Y}$ since that would put the SM gauge invariance itself in jeopardy. So, in that representation, the fact that the chiral $SU(2)_{L}\otimes U(1)_{Y}$ symmetry is spontaneously broken enters only and entirely through the fermion masses. The triangle graph explicitly break the $SU(2)_{L}\otimes U(1)_{Y}$ symmetry because the $m_{u,d,e}\rightarrow\infty$ limit is incompatible with that symmetry. In the opposite limit $m_{u,d,e}\rightarrow 0$, the $SU(2)_{L}\otimes U(1)_{Y}$ symmetry is recovered but in a very trivial way: the triangle graphs precisely cancel the local anomalous term and all the amplitudes simply vanish.

• •

We derived generalized forms for all the triangle anomalies, including the terms proportional to the fermion masses, see Eqs. (30) and (31). This requires going beyond usual regularization procedures. To our knowledge, and except those based on dispersion relations [17], only that described by Weinberg in Ref. [15] offers sufficient freedom to choose which current in the $AVV$ and $AAA$ triangles is to carry the anomaly. This is crucial to deal with chiral gauge theories, in which some axial currents have to remain anomaly-free to preserve gauge invariance, as well as for treating the $\mathcal{B}$ and $\mathcal{L}$ vectorial currents [22].

In the present paper, we consider only true axion models, those designed to solve the strong $CP$ puzzle and for which the axion ends up extremely light and extremely weakly coupled to SM fields. However, our results could have consequences for so-called Axion-Like Particle (ALP) searches (see e.g. Refs. [25, 26, 27, 28, 29, 30, 31]). Specifically, our analysis shows that an adequate effective parametrization is

[TABLE]

with from Eq. (93),

[TABLE]

for some free parameters $\chi_{P}^{f}$, and from Eq. (25)

[TABLE]

for some free parameters $\mathcal{N}_{C,L,Y}^{eff}$. Typically, $\mathcal{L}_{\mathrm{KSVZ-like}}$ encodes the effects of heavy vector-like fermions, while $\mathcal{L}_{\mathrm{DFSZ-like}}$ encode the anomaly-free impact of charging the SM fermions under some global $U(1)_{PQ}$ symmetry whose current is coupled to the $a^{0}$ field. The impact of some new heavy chiral fermions (like a sequential fourth generation or a heavy neutrino) is not included, but could be by extending the sum in $\mathcal{L}_{\mathrm{DFSZ-like}}$ to those states also. This parametrization of $\mathcal{L}_{\mathrm{DFSZ-like}}$ ensures $a^{0}$ is truly axion-like: no $\mathcal{B}$ or $\mathcal{L}$ violating effects are implicit, only physical free parameters are introduced, and UV divergences should be under control. Still, it is important to stress that the $a^{0}\rightarrow\gamma Z$, $ZZ$, and $W^{+}W^{-}$ amplitudes induced by $\mathcal{L}_{\mathrm{DFSZ-like}}$ are entirely tuned by triangle graphs, and in the $m_{f}\rightarrow\infty$ limit, they do not match the pattern of $\mathcal{L}_{\mathrm{KSVZ-like}}$ but rather that in Eq. (47).