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

The Non-triviality of Dynamical Chern-Simons Gravity and the Standard Model

Stephon Alexander [email protected]    Heliudson Bernardo [email protected] Brown Theoretical Physics Center and Department of Physics, Brown University,
Barus Building, 340 Brook Street, Providence, RI 02912, USA
   Cyril Creque-Sarbinowski [email protected] Center for Computational Astrophysics, Flatiron Institute
162 5th Avenue, New York, New York 10010, USA
Abstract

Given the growing interest in gravitational-wave and cosmological parity-violating effects in dynamical Chern-Simons (dCS) gravity, it is crucial to investigate whether the scalar-gravitational Pontryagin term in dCS persists when formulated in the context of the U(1)BL\text{U(1)}_{\text{B}-\text{L}} anomaly in the Standard Model (SM). In particular, it has been argued that dCS gravity can be reduced to Einstein gravity after “rotating away” the gravitational-Pontryagin coupling into the phase of the Weinberg operator – analogous to the rotation of the axion zero-mode into the quark mass matrix. We find that dCS is nontrivial if the scalar field ϕ\phi has significant space-time dependence from dynamics. We provide a comprehensive consideration of the dCS classical and quantum symmetries relevant for embedding a dCS sector in the SM. We find that, because of the B-L chiral gravitational anomaly, the scalar-Pontryagin term cannot be absorbed by a field redefinition. Assuming a minimal extension of the SM, we also find that a coupling of the dCS scalar with right-handed neutrinos induces both the scalar-Pontryagin coupling and an axion-like phase in the dimension-five Weinberg operator. We comment on the issue of gauging the U(1)BL\text{U(1)}_{\text{B}-\text{L}}, the observational effects with these two operators present for upcoming experiments, and the origin of dCS gravity in string theory.

preprint: APS/123-QED

I Introduction

Given that the electroweak sector of the Standard Model (SM) maximally violates parity, it is natural to investigate whether parity-violation occurs in sectors beyond the SM. In particular, the gravitational sector has been under such observational and theoretical investigation over the past several years [1, 2, 3, 4, 5, 6, 7, 8]. The primary phenomenological motivation is that parity-violating modifications of general relativity (GR) lead to a myriad of distinct astrophysical and cosmological observables.

As originally demonstrated in Ref. [9], the cosmic microwave background (CMB) is sensitive to parity violation through its EB and TB angular power spectra. In fact, utilizing Galactic foreground emission as a polarization-angle calibrator [10, 11, 12], there are recent hints that Planck’s EB angular power spectrum shows such a signature [13, 14]. In addition, there are further indications that the four-point galaxy correlation function is parity violating [15, 16], although the CMB’s four-point correlators yield no analogous signal [17, 18]. Finally, it is possible that parity-violation may be encoded in the large-scale correlation of galaxy spins, as explored in Refs. [19, 20, 21, 22]. Clearly, as such parity-sensitive data accumulates, the need for well-motivated parity-breaking GR theories increases.

The simplest and most theoretically motivated effective theory of gravity that encodes parity violation is dynamical Chern-Simons (dCS) gravity111An earlier non-dynamical version that violates the strong equivalence principle of Chern-Simons gravity was first put forth in Ref. [23]. See also [24, 25, 26] and references therein for works on the Lorentz CS form in supersymmetric theories.There are different motivations to study such a theory. The most basic one is the effective theory perspective, which states one could add all possible terms compatible with symmetries. Thus, dCS gravity is a valid candidate for a classical scalar-tensor modification of GR. This approach to dCS has been explored in the literature before (see Ref. [1] and references therein).

DCS gravity was the framework within which the cosmic baryon/lepton number (B-L) asymmetry was instantiated in inflationary models through the production of parity asymmetry of chiral gravitational waves which sourced a B-L asymmetry by the end of inflation [27]. It was later demonstrated that these gravitational waves can induce a scalar four-point (trispectrum) that contributes to the primordial large-scale parity violation in galactic distributions [28]. In the strong gravity regime, such as compact binary systems, dCS gravity leads to distinct waveforms in the propagation and sourcing of gravitational waves as well as the phenomena of scalarization [29, 30].

It is therefore important to consider how dCS connects with the standard model and its logical extensions. From the perspective of modified theories of gravity, such as f(R)f(R) theories, the principle of general covariance is enough to formulate it. This is the case of versions of dCS that emerge as an effective theory of 4D4D Heterotic string theory, and there is no need to relate it directly to the standard model. However, dCS gravity has a term where a pseudoscalar couples to the gravitational Chern-Simons form, which also encodes the B-L global anomaly222We shall call a classical symmetry anomalous if the path integral measure is not invariant under the symmetry transformation. in the standard model. This structure has a close semblance to the axion coupling to the Yang-Mills Pontryagin density. For the latter case, in the presence of massless quarks and when the axion is stabilized at a value θ\theta, a field redefinition of the quark fields can eliminate the θFF~\theta F\tilde{F} coupling, rendering the CP-violating effects ignorable. In this paper, we investigate to what extent and how these axion features apply to the dCS pseudoscalar. This is relevant for any SM extension in which the dCS scalar-Pontryagin coupling appears after integrating out heavy fermions.

The paper is organized as follows. In the next section, we discuss the global shift symmetry of dCS gravity and use analogous axion results to comment on what to expect in the quantum gravitational regime. In Section III, we show that dCS gravity is non-trivial even after including fermions, highlighting some features relevant to the coupling with SM fermions. In Section IV, we propose an extension of the SM where a dCS sector appears after integrating out right-handed neutrinos. We conclude in Section V with a discussion on the combined shift and B-L symmetry of the model and its gauging. Some aspects of the baryonic and leptonic number anomalies in the SM are summarized in Appendix A.

Conventions: We set =c=1\hbar=c=1. Moreover, we work with the mostly plus metric-signature convention and the Weyl (chiral) basis for the Dirac matrices.

II dCS Gravity and Axions

II.1 Continuous Shift Symmetry

The dCS action is

SdCS[gμν,ϕ]=d4xg\displaystyle S_{\text{dCS}}[g_{\mu\nu},\phi]=\int d^{4}x\sqrt{-g} (12κ2R12μϕμϕ\displaystyle\left(\frac{1}{2\kappa^{2}}R-\frac{1}{2}\partial_{\mu}\phi\partial^{\mu}\phi-\right.
λϕR~μνρσRμνρσ),\displaystyle\left.-\lambda\phi\tilde{R}_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma}\right), (1)

where gμνg_{\mu\nu} is the spacetime metric with gg its determinant, ϕ\phi is the dCS pseudoscalar, RR is the Ricci scalar, RμνρσR^{\mu\nu\rho\sigma} is the Riemann tensor with R~μνρσ\tilde{R}_{\mu\nu\rho\sigma} its Hodge dual, and λ\lambda the dCS coupling constant.

Under the global symmetry ϕ(xμ)ϕ(xμ)+χ\phi(x^{\mu})\to\phi(x^{\mu})+\chi, Eq. (II.1) changes by a constant shift since the last term can be written as the derivative of the gravitational Chern-Simons (CS) current KμK^{\mu},

R~μνρσRμνρσ\displaystyle\tilde{R}_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma} =2ϵμνρσμ[tr(ω[νρωσ]+23ω[νωρωσ])]\displaystyle=2\epsilon^{\mu\nu\rho\sigma}\partial_{\mu}\left[\text{tr}\left(\omega_{[\nu}\partial_{\rho}\omega_{\sigma]}+\frac{2}{3}\omega_{[\nu}\omega_{\rho}\omega_{\sigma]}\right)\right]
μKμ,\displaystyle\equiv\nabla_{\mu}K^{\mu}, (2)

where ϵμνρσ\epsilon^{\mu\nu\rho\sigma} are the contravariant components of the Levi-Civita tensor, and the trace is taken over the indices of the Lorentz-algebra valued spin connection ωμba{\omega_{\mu}}^{a}_{\;\;b}. The equations of motion will be invariant under constant scalar shifts. We will call this shift symmetry U(1)dCS\text{U}(1)_{\text{dCS}}. Its associated Noether current is

jdCSμ=μϕ+λKμ,j^{\mu}_{\text{dCS}}=-\partial^{\mu}\phi+\lambda K^{\mu}, (3)

which satisfy

μjdCSμ=0\nabla_{\mu}j^{\mu}_{\text{dCS}}=0 (4)

upon using ϕ\phi’s equation of motion ϕλR~R=0\Box\phi-\lambda\tilde{R}R=0.333We note that the SdCSS_{\text{dCS}} is not strictly invariant under global ϕ\phi shifts because the dCS Lagrangian changes by a total derivative. Symmetries satisfying this property are often referred to as quasi-symmetries, so a UdCS(1)\text{U}_{\text{dCS}}(1) transformation is more properly referred to as a quasi-symmetry of the action, although we will make no such distinction in what follows.

For field theories in a flat spacetime with vanishing boundary conditions, Eq. (4) guarantees the existence of a conserved charge. In curved spacetime, a conserved Noether charge is obtained if we assume a globally hyperbolic and asymptotically flat spacetime. In this case, we then integrate Eq. (4) over a spacetime volume VV foliated by spacelike surfaces Σt\Sigma_{t} whose boundaries Σt\partial\Sigma_{t} are 2-spheres with a very large radius (compared to the maximum curvature scale inside VV), where tt is the parameter of the timelike curves perpendicular to Σt\Sigma_{t}. We have

0\displaystyle 0 =Vd4xgμjμ=Vd3x|h|nμjμ\displaystyle=\int_{V}d^{4}x\sqrt{-g}\nabla_{\mu}j^{\mu}=\int_{\partial V}d^{3}x\sqrt{|h|}n_{\mu}j^{\mu}
=Σt1𝑑Sj0Σt2𝑑Sj0+×S2𝑑Sniji.\displaystyle=\int_{\Sigma_{t_{1}}}dSj^{0}-\int_{\Sigma_{t_{2}}}dSj^{0}+\int_{\mathbb{R}\times S^{2}}dSn_{i}j^{i}. (5)

Note that jμj^{\mu} is defined up to a term Δjμ\Delta j^{\mu} that vanishes at the boundary of VV. If the last integral in the second line of Eq. (II.1) vanishes, then we have that

Q=Σt𝑑Sj0Q=\int_{\Sigma_{t}}dSj^{0} (6)

is independent of Σt\Sigma_{t} and thus is conserved. Typically it is assumed that niji0n_{i}j^{i}\to 0 asymptotically such that the last integral in Eq. (II.1) vanishes and we have charge conservation.

Let us now apply these facts to dCS gravity. Assuming that iϕ\partial^{i}\phi vanishes at the boundary, only the gravitational CS current will contribute to the integral of nijin_{i}j^{i}. However, the spin-connection will be pure gauge on the cylinder at infinity, and we can choose a gauge where its time component ω0\omega_{0} vanishes there [31]. In such a case, the CS contribution to nijdCSin_{i}j^{i}_{\text{dCS}} is null. Then, the conserved charge is

QdCS=Σt𝑑S(ϕ˙+λK0),Q_{\text{dCS}}=\int_{\Sigma_{t}}dS\left(\dot{\phi}+\lambda K^{0}\right), (7)

This charge is precisely the momentum of ϕ\phi, with the contribution K0K^{0} coming from the fact that the Lagrangian contains a term proportional to λϕ˙K0\lambda\dot{\phi}K^{0}. If QdCSQ_{\text{dCS}} were gauge invariant, the expression above would be independent of the gauge choice ω0|Σt=0\omega_{0}|_{\partial\Sigma_{t}}=0. As shown below, QdCSQ_{\text{dCS}} changes under large Lorentz gauge transformations. Under a Lorentz transformation Λ\Lambda, we have that the gravitational CS current transforms as KμKμK^{\mu}\to K^{\prime\mu}, where

Kμ=Kμ2ϵμνρσtr\displaystyle K^{\prime\mu}=K^{\mu}-2\epsilon^{\mu\nu\rho\sigma}\text{tr} [ν(ρΛΛ1ωσ)+\displaystyle\left[\partial_{\nu}\left(\partial_{\rho}\Lambda\Lambda^{-1}\omega_{\sigma}\right)+\right.
+13νΛΛ1ρΛΛ1σΛΛ1].\displaystyle\left.+\frac{1}{3}\partial_{\nu}\Lambda\Lambda^{-1}\partial_{\rho}\Lambda\Lambda^{-1}\partial_{\sigma}\Lambda\Lambda^{-1}\right]. (8)

It follows that [32]

Q=Q23λΣt𝑑Sϵ0ijktr[iΛΛ1jΛΛ1kΛΛ1]\displaystyle Q^{\prime}=Q-\frac{2}{3}\lambda\int_{\Sigma_{t}}dS\epsilon^{0ijk}\text{tr}\left[\partial_{i}\Lambda\Lambda^{-1}\partial_{j}\Lambda\Lambda^{-1}\partial_{k}\Lambda\Lambda^{-1}\right] (9)

and thus

QdCS=QdCS16π2λν(Λ),Q^{\prime}_{\text{dCS}}=Q_{\text{dCS}}-16\pi^{2}\lambda\nu(\Lambda), (10)

where ν(𝚲)\nu(\bf\Lambda) is the winding number of the corresponding Lorentz transformation (that might be non-vanishing even with Λ(xi)𝟏\Lambda(x^{i})\to\mathbf{1} asymptotically at infinity). In particular, we have that the Noether charge is only invariant under transformations such that ν(Λ)=0\nu(\Lambda)=0, i.e. for transformations continuously connected to the identity. But the charge transforms non-trivially for large Lorentz transformations. Eq. (10) then shows that the charge QdCSQ_{\text{dCS}} associated with global shifts in the dCS scalar transform non-trivially under large gauge transformation while still being conserved.

Some of the above results are analogous to axions coupled to non-Abelian gauge fields [33, 34, 35, 36]. To better understand what can be applied to the gravitational case, we shall review some aspects of axions in the next section.

II.2 Non-Abelian Axions

We wish to understand what happens with the shift symmetry and the conserved charge when we quantize the scalar ϕ\phi in theory. Because all the results above have analogs for axions and their coupling to gauge fields, we will first consider the system

S=d4x(12μaμaλatrF~μνFμν),S=\int d^{4}x\left(-\frac{1}{2}\partial_{\mu}a\partial^{\mu}a-\lambda a\text{tr}\tilde{F}_{\mu\nu}F^{\mu\nu}\right), (11)

which includes an axion-like field aa coupled to a SU(2)(2) gauge field AμaA_{\mu}^{a} via the term proportional to λ\lambda. We can consider this system in Minkowski space and Cartesian coordinates for simplicity since curvature is unimportant for the axion discussion in this section. Classically, the action above has the global shift symmetry aa+ca\to a+c with conserved current

jμ=μa+λCμ,j^{\mu}=-\partial^{\mu}a+\lambda C^{\mu}, (12)

where CμC^{\mu} is the Chern-Simons current of the gauge field

Cμ=2ϵμνρσtr(AνρAσ+23AνAρAσ),C^{\mu}=2\epsilon^{\mu\nu\rho\sigma}\text{tr}\left(A_{\nu}\partial_{\rho}A_{\sigma}+\frac{2}{3}A_{\nu}A_{\rho}A_{\sigma}\right), (13)

and the trace is now taken over the adjoint of the gauge group. The Noether charge is

Qa=Σ𝑑S(a˙+λC0),Q_{a}=\int_{\Sigma}dS\left(\dot{a}+\lambda C^{0}\right), (14)

where we chose the gauge A0=0A^{0}=0 without loss of generality. Under a gauge transformation Aμg1Aμg+g1μgA_{\mu}\to g^{-1}A_{\mu}g+g^{-1}\partial_{\mu}g where

g(xi)=eif(r)r^τg(x^{i})=e^{if(r)\hat{r}\cdot\tau} (15)

with τa\tau^{a} the generators of SU(2)(2) and f(r)2πnf(r)\to 2\pi n as rr\to\infty, the charge QaQ_{a} transforms as [32, 37]

Qa=Qa16π2λn(g).Q_{a}^{\prime}=Q_{a}-16\pi^{2}\lambda n(g). (16)

This gauge dependence does not introduce any issues because once we fix the background gauge field configuration, the system will not be invariant under large gauge transformations with a non-trivial winding number. To see that explicitly, note that a given background configuration Aμ(t,xi)A^{\mu}(t,x^{i}) should approach a pure gauge one as rr\to\infty, and we can choose it to be of the form (15) with f(r)2πmf(r)\to 2\pi m at infinity. In other words, the homotopy class of the gauge configuration is part of the system’s definition so that large gauge transformations would change the system altogether. Thus, by definition, large gauge transformations cannot be a symmetry of the system, and it is no surprise that such a gauge transformation changes QaQ_{a}.

The conservation of the homotopy-class dependent QaQ_{a} implies that the dynamical evolution cannot change the winding number of the background configuration A(ν)μA_{(\nu)}^{\mu}. In fact, even if we promote AμA^{\mu} to a dynamical field, there are no solutions to the equations of motion in Lorentzian signature that would allow a change in winding number [31]. This also guarantees that if AμA^{\mu} is pure gauge asymptotically at t±t\to\pm\infty, the action is invariant under global aa-shifts: we have

S[a+χ]=S[a]λχd4xtrF~F,S[a+\chi]=S[a]-\lambda\chi\int d^{4}x\text{tr}\tilde{F}F, (17)

but

d4xtrF~F\displaystyle\int d^{4}x\text{tr}\tilde{F}F =d4xμCμ=d3xC0|t=t==\displaystyle=\int d^{4}x\partial_{\mu}C^{\mu}=\left.\int d^{3}xC^{0}\right|_{t=-\infty}^{t=\infty}=
=32π2[n(t=+)n(t=)]=0,\displaystyle=32\pi^{2}\left[n(t=+\infty)-n(t=-\infty)\right]=0, (18)

since the homotopy class of the gauge field does not change in time.

However, if AμaA_{\mu}^{a} is promoted to a dynamical field that should also be quantized, the argument in the paragraph above cannot be applied. If we also quantize the gauge field, we cannot exclude the possibility of quantum tunneling between gauge configurations of different homotopy classes. This is precisely what happens in quantum chromodynamics (QCD) and general non-Abelian gauge theories, as manifested by instanton solutions. This is the reason why the vacuum structure of gauge theories is non-trivial [37].

With a quantized gauge field, the system cannot be fixed in given homotopy class because there is always a probability for the gauge field to tunnel to configurations of different classes. In other words, the partition function includes a sum over the instantons connecting different homotopy classes,

Z=ν\displaystyle Z=\sum_{\nu}\int [dμ(ν)dA(ν)da]exp(iSgauge)×exp(iθν)\displaystyle[d\mu_{(\nu)}dA_{(\nu)}da]\exp\left({iS_{\text{gauge}}}\right)\times\exp\left({i\theta\nu}\right) (19)
×\displaystyle\times exp[i(12trF212(a)2λaF~F)],\displaystyle\exp\left[{i\int\left(-\frac{1}{2}\text{tr}F^{2}-\frac{1}{2}(\partial a)^{2}-\lambda a\tilde{F}F\right)}\right],

where we wrote the (homotopically trivial) gauge fixing and Fadeev-Popov contributions in SgaugeS_{\text{gauge}} and [dμ(ν)][d\mu_{(\nu)}]. Then, the vacuum state of the whole system depends on θ\theta, but all physical observables would be the same regardless of the vacuum θ\theta we use to compute them. In the quantum theory defined above, the specific value of theta is unphysical because of the shift symmetry at the quantum level: since

ν=132π2d4xF~F,\nu=\frac{1}{32\pi^{2}}\int d^{4}x\tilde{F}F, (20)

a shift in aa, which should not change any physical observable, induces a shift in θ\theta. Another way of seeing this is by looking at how QaQ_{a} acts on a given θ\theta vacuum of the pure gauge sector (in the absence aa, there is no shift symmetry and θ\theta is observable). Under a gauge transformation in the homotopy class n=1n=1, we have

|θeiθ|θ.|\theta\rangle\to e^{-i\theta}|\theta\rangle. (21)

So, since QaQ_{a} transforms as Eq. (16), we have [32]

exp(iθ16π2λQa)|θ=|θ+θ.\exp\left({\frac{i\theta^{\prime}}{16\pi^{2}\lambda}Q_{a}}\right)|\theta\rangle=|\theta+\theta^{\prime}\rangle. (22)

As QaQ_{a} is conserved (commutes with the Hamiltonian), different theta vacua are degenerated in energy and thus physically indistinguishable.

This degeneracy between the choice of θ\theta and the zero mode (constant part) of aa has a different nature than the usual quantum symmetries. Consider, for instance, the symmetry breaking potential V(Φ)=λΦ(|Φ|2μΦ2/2λΦ)2V(\Phi)=\lambda_{\Phi}\left(|\Phi|^{2}-\mu^{2}_{\Phi}/2\lambda_{\Phi}\right)^{2} for the complex scalar field Φ(x)=r(x)eiα(x)\Phi(x)=r(x)e^{i\alpha(x)}. In this case, a vacuum expectation value for r(x)r(x) spontaneously breaks the symmetry while α\alpha parametrizes the vacuum manifold, with shifts in α\alpha being the non-linear realizations of U(1)(1) phase rotation of the complex scalar field. If we think of aa as a Goldstone mode after some phase transition (as is the case in the Peccei-Quinn (PQ) proposal), then shifts around the vacuum manifold at constant θ\theta are generated by the CS-independent part of QaQ_{a}, that we will call Q~a\tilde{Q}_{a}:

Qa=Q~a+QCS,Q~a=𝑑Sa˙.Q_{a}=\tilde{Q}_{a}+Q_{\rm CS},\quad\tilde{Q}_{a}=\int dS\dot{a}. (23)

Due to Eq. (22), Q~a\tilde{Q}_{a} and QCSQ_{\rm CS} are not individually conserved, and so the conservation of QaQ_{a} implies that Q~a\tilde{Q}_{a} is not conserved: shifts of the vacuum manifold for a given value of θ\theta are not symmetries of the θ\theta invariant quantum theory. Instanton effects explicitly break the continuous θ\theta symmetry transformations. Therefore, there are no Goldstone bosons associated with this instanton induced explitict symmetry breaking. Note that it is a mistake to claim that the axion is a Golsdone boson associated with the anomalous breaking by these instanton effects.

These arguments are at the core of ’t Hooft’s solution to the U(1)(1) problem in QCD [38, 39]. The discussion about how the U(1)(1) symmetry is affected by non-perturbative effects is very similar to the chiral symmetry of massless fermions coupled to non-Abelian gauge fields. The standard chiral rotation combines with the instanton effects to make physics independent of θ\theta. Equivalently, the chiral rotation alone is explicitly broken by instantons.

Returning to dCS gravity, we can apply the following non-Abelian axion results to the dCS pseudoscalar: (ii) The dCS action has a continuous constant shift (quasi) symmetry with an associated divergent-free current jdCSμj^{\mu}_{\text{dCS}}; (iiii) If gravity is quantized in the path integral for the theory, the CS-independent part of the total UdCS(1)\text{U}_{\text{dCS}}(1) charge is not conserved due to gravitational instanton effects444See [40] for a discussion on the gravitational theta sector.; (iiiiii) The total charge QdCSQ_{\text{dCS}} (7) is conserved, but not invariant under homotopically non-trivial gauge transformations.

So, we see that the U(1)dCS(1)_{\text{dCS}} shift symmetry is always present, while the CS-independent part of the charge Q~dCS\tilde{Q}_{\text{dCS}} is only conserved in flat spacetime in the absence of gravitational instantons (which are absent at the semiclassical level). In that regime, Q~dCS\tilde{Q}_{\text{dCS}} is not conserved because the local coupling of ϕ\phi with the gravitational Pontryagin density changes the dynamical evolution of the metric and the dCS scalar. For instance, in regions where R~R\tilde{R}R is non-trivial, it will source ϕ\phi. This does not mean the symmetry under ϕ\phi translations is broken: the total charge QdCSQ_{\text{dCS}} is conserved (even after quantizing ϕ\phi). Similarly, one does not worry about non-trivial topological effects when using the axion coupling to F~F\tilde{F}F to compute the axion-photon interactions or birefringence effects [41, 42, 14].

In the next discussion, we shall couple a Dirac fermion with dCS and show that gravitational instantons do not affect a combination of the dCS-shift and fermion-axial currents.

III Adding Fermions

Before addressing SM fermions, we first consider a single Dirac fermion ψ\psi in addition to the dCS action,

S[eaμ,ϕ,ψ]\displaystyle S[e^{\mu}_{a},\phi,\psi] =d4xg[12κ2R12μϕμϕ\displaystyle=\int d^{4}x\sqrt{-g}\left[\frac{1}{2\kappa^{2}}R-\frac{1}{2}\partial_{\mu}\phi\partial^{\mu}\phi-\right.
λϕR~μνρσRμνρσψ¯eaμγa(μ+14ωμabγab)ψ],\displaystyle\left.-\lambda\phi\tilde{R}_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma}-\bar{\psi}e^{\mu}_{a}\gamma^{a}\left(\partial_{\mu}+\frac{1}{4}\omega^{ab}_{\mu}\gamma_{ab}\right)\psi\right], (24)

where eaμe^{\mu}_{a} is the vierbein field and γab=γ[aγb]\gamma_{ab}=\gamma_{[a}\gamma_{b]} are the Lorentz group generators for Dirac fermions. This classical theory has the global symmetry group U(1)dCS×U(1)V×U(1)A\text{U}(1)_{\text{dCS}}\times\text{U}(1)_{\text{V}}\times\text{U}(1)_{\text{A}}. The new vector and axial U(1)\text{U}(1)’s are phase and chiral phase transformations of the fermion, ψeiαψ\psi\to e^{i\alpha}\psi and ψeiβγ5ψ\psi\to e^{i\beta\gamma_{5}}\psi, respectively. Note that the fermion does not transform under ϕ\phi shifts. The new Noether currents are the ordinary vector and axial currents,

jVμ=iψ¯γμψ,jAμ=iψ¯γμγ5ψ,j^{\mu}_{\text{V}}=i\bar{\psi}\gamma^{\mu}\psi,\quad j^{\mu}_{\text{A}}=i\bar{\psi}\gamma^{\mu}\gamma_{5}\psi, (25)

which are classically conserved upon using the equations of motion of the fermion.

If we consider this theory at the quantum level, the currents will satisfy

μ(jAμ1192π2Kμ)\displaystyle\nabla_{\mu}\left(\langle j^{\mu}_{\text{A}}\rangle-\frac{1}{192\pi^{2}}K^{\mu}\right) =0,\displaystyle=0, (26a)
μ(j~dCSμ+λKμ)\displaystyle\nabla_{\mu}\left(\langle\tilde{j}^{\mu}_{\text{dCS}}\rangle+\lambda K^{\mu}\right) =0,\displaystyle=0, (26b)
μjVμ\displaystyle\nabla_{\mu}\langle j_{\text{V}}^{\mu}\rangle =0.\displaystyle=0. (26c)

The appearance of KμK^{\mu} in the divergence of jAμj^{\mu}_{\text{A}} is due to a non-trivial path-integral Jacobian under axial transformations. Given its dependence on the curvature tensor, which has the same structure as the ϕ\phi-curvature coupling in dCS, one might wonder if a chiral rotation can modify the value of λ\lambda. In other words, there is a linear combination of the scalar shift and chiral rotation that is not affected by instanton effects at the quantum level, and one might wonder whether we can use such a combined U(1)(1) to hide the dCS scalar coupling inside the fermion chiral phase. If so, the pseudoscalar-Pontryagin coupling in dCS gravity would not be observable. This is not the case, as we shall see below, since a constant chiral rotation cannot make a whole dynamic coupling disappear.

We write the path integral for the quantum theory as:

Z[e,ω]=eiSEH[dϕdψ¯dψ]exp{id4xg[12μϕμϕλϕR~μνρσRμνρσψ¯eaμγa(μ+14ωμabγab)ψ]},Z[e,\omega]=e^{iS_{\text{EH}}}\int[d\phi d\bar{\psi}d\psi]\exp\left\{i\int d^{4}x\sqrt{-g}\left[-\frac{1}{2}\partial_{\mu}\phi\partial^{\mu}\phi-\lambda\phi\tilde{R}_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma}-\bar{\psi}e^{\mu}_{a}\gamma^{a}\left(\partial_{\mu}+\frac{1}{4}\omega^{ab}_{\mu}\gamma_{ab}\right)\psi\right]\right\}, (27)

where we brought the Einstein-Hilbert action SEHS_{\text{EH}} contribution to outside the path integral. Under the field redefinition

ψ(x)ψ~(x)=eiβ(x)γ5ψ(x),ψ¯(x)ψ¯~=ψ¯(x)eiβ(x)γ5,\psi(x)\to\tilde{\psi}(x)=e^{i\beta(x)\gamma_{5}}\psi(x),\;\bar{\psi}(x)\to\tilde{\bar{\psi}}=\bar{\psi}(x)e^{i\beta(x)\gamma_{5}}, (28)

we have [dψ¯dψ][dψ¯~dψ~]=[dψ¯dψ]J[β][d\bar{\psi}d\psi]\to[d\tilde{\bar{\psi}}d\tilde{\psi}]=[d\bar{\psi}d\psi]J[\beta], where [43, 44]

J[β]=exp[id4xgβ(x)192π2R~μνρσRμνρσ].J[\beta]=\exp\left[-i\int d^{4}x\sqrt{-g}\frac{\beta(x)}{192\pi^{2}}\tilde{R}_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma}\right]. (29)

Together with the fact that Z[e,ω]Z[e,\omega] does not change by a change of integration variables, this result implies Eq. (26a). That Z[e,ω]Z[e,\omega] is invariant also explains why there is no choice of β(x)\beta(x) in the Jacobian that kills the ϕR~R\phi\tilde{R}R term in the action, as can be seen from

Z[e,ω]=eiSEH[dϕdψ¯~dψ~]\displaystyle Z[e,\omega]=e^{iS_{\text{EH}}}\int[d\phi d\tilde{\bar{\psi}}d\tilde{\psi}] exp{id4xg[12μϕμϕλϕR~μνρσRμνρσψ¯~eaμγa(μ+14ωμabγab)ψ~]}\displaystyle\exp\left\{i\int d^{4}x\sqrt{-g}\left[-\frac{1}{2}\partial_{\mu}\phi\partial^{\mu}\phi-\lambda\phi\tilde{R}_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma}-\tilde{\bar{\psi}}e^{\mu}_{a}\gamma^{a}\left(\partial_{\mu}+\frac{1}{4}\omega^{ab}_{\mu}\gamma_{ab}\right)\tilde{\psi}\right]\right\}
=eiSEH[dϕdψ¯dψ]\displaystyle=e^{iS_{\text{EH}}}\int[d\phi d\bar{\psi}d\psi] exp{id4xg[12μϕμϕ(λϕ+β192π2)R~μνρσRμνρσψ¯eaμγa(μ+14ωμabγab)ψ\displaystyle\exp\left\{i\int d^{4}x\sqrt{-g}\left[-\frac{1}{2}\partial_{\mu}\phi\partial^{\mu}\phi-\left(\lambda\phi+\frac{\beta}{192\pi^{2}}\right)\tilde{R}_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma}-\bar{\psi}e^{\mu}_{a}\gamma^{a}\left(\partial_{\mu}+\frac{1}{4}\omega^{ab}_{\mu}\gamma_{ab}\right)\psi-\right.\right.
μβjAμ]},\displaystyle-\left.\left.\partial_{\mu}\beta j^{\mu}_{\text{A}}\right]\right\}, (30)

so, even if we choose β(x)=192π2λϕ\beta(x)=-192\pi^{2}\lambda\phi, we will find

Z[e,ω]\displaystyle Z[e,\omega] =eiSEH[dϕdψ¯dψ]exp{id4xg[12μϕμϕψ¯eaμγa(μ+14ωμabγab)ψ192π2λϕμjAμ]}\displaystyle=e^{iS_{\text{EH}}}\int[d\phi d\bar{\psi}d\psi]\exp\left\{i\int d^{4}x\sqrt{-g}\left[-\frac{1}{2}\partial_{\mu}\phi\partial^{\mu}\phi-\bar{\psi}e^{\mu}_{a}\gamma^{a}\left(\partial_{\mu}+\frac{1}{4}\omega^{ab}_{\mu}\gamma_{ab}\right)\psi-192\pi^{2}\lambda\phi\nabla_{\mu}j^{\mu}_{\text{A}}\right]\right\}
=eiSEH[dϕdψ¯dψ]exp{id4xg[12μϕμϕψ¯eaμγa(μ+14ωμabγab)ψλϕR~μνρσRμνρσ]},\displaystyle=e^{iS_{\text{EH}}}\int[d\phi d\bar{\psi}d\psi]\exp\left\{i\int d^{4}x\sqrt{-g}\left[-\frac{1}{2}\partial_{\mu}\phi\partial^{\mu}\phi-\bar{\psi}e^{\mu}_{a}\gamma^{a}\left(\partial_{\mu}+\frac{1}{4}\omega^{ab}_{\mu}\gamma_{ab}\right)\psi-\lambda\phi\tilde{R}_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma}\right]\right\}, (31)

where to get the second equality we used the fact that, inside the path integral, the axial current satisfies Eq. (26a). So, we cannot simply set λ\lambda to zero, and the ϕR~R\phi\tilde{R}R coupling of dCS is quantum mechanically non-trivial even after the coupling with fermions. Equivalently, we could have redefined ϕϕ(192π2λ)1β\phi\to\phi-(192\pi^{2}\lambda)^{-1}\beta in the last line of Eq. (III), keeping the R~R\tilde{R}R coupling.

One redefinition that can kill some part of ϕ\phi, however, is to assume a constant β\beta. In this case, the last term inside the exponent of the last line in Eq. (III) will not contribute, and thus, we cannot recover a R~R\tilde{R}R term from it. This is because a field redefinition of the form Eq. (28) with a constant β\beta is a genuine global transformation. Thus, quantum mechanically, the (constant) zero mode ϕ0\phi_{0} of ϕ(xμ)\phi(x^{\mu}) can be absorbed into the chiral phase of the fermions. Although such a phase is not observable, phase differences are, so λ\lambda might still be measurable across regions where ϕ0\phi_{0} changes abruptly (e.g., domain wall solutions). More interestingly, if ψ\psi has a mass, only a combination of ϕ0\phi_{0} and the mass phase will be observable, and any measurement or observation sensible to the phase of the mass and ϕ0\phi_{0} will also constrain the value of λ\lambda. Note that the mass term’s explicit chiral symmetry breaking will not affect these conclusions. On top of that, the results above are naturally extended for couplings with multiple fermions.

The results in this section are a direct consequence of the fact that the path integral is invariant under redefinitions of the integrating fields. Our goal here was to explicitly check that if one wants to use the nontrivial Jacobian (anomaly) to remove the CS coupling with ϕ\phi then, by consistency, one cannot neglect that μjAμ0\nabla_{\mu}j^{\mu}_{\text{A}}\neq 0, which ends up reintroducing the CS coupling. The same calculations also allow us to conclude that, when there is a mass term for the fermions, ϕ0\phi_{0} alone is not physical since only a combination of zero mode of ϕ\phi and the mass matrix phase is observable.

IV A UV completion of a dCS sector in the SM

In this section, we will propose a renormalizable four-dimensional model in which a dCS sector appears after integrating right-handed neutrinos Ψ\Psi in an extension of the SM. We start by identifying the dCS scalar as the phase of a complex scalar whose Yukawa coupling with Ψ\Psi gives them a Majorana mass after spontaneous symmetry breaking. We then couple the sterile right-handed neutrinos with the left-handed neutrinos via another Yukawa coupling and deduce the EFT for the latter and the dCS pseudoscalar below the symmetry-breaking scale.

IV.1 DCS gravity from integrating out heavy fermions

Similar to axions arising from a PQ symmetry breaking, we can realize the dCS pseudoscalar as a Goldstone boson associated with the spontaneous breaking of a U(1)(1) symmetry. The remaining constant scalar shift is then a non-linear realization of the global symmetry below the breaking scale. We review this picture in the following model, recently discussed in Ref. [45] (see also Ref. [46] for the non-Abelian axion case),

S=d4xg\displaystyle S=\int d^{4}x\sqrt{-g} [Ψ¯eaμγa(μ+14ωμcdγcd)Ψ|Φ|2\displaystyle\left[-\bar{\Psi}e^{\mu}_{a}\gamma^{a}\left(\partial_{\mu}+\frac{1}{4}\omega^{cd}_{\mu}\gamma_{cd}\right)\Psi-|\partial\Phi|^{2}\right.
yΨ¯(ΦPR+Φ¯PL)ΨV(|Φ|)],\displaystyle\left.-y\bar{\Psi}\left(\Phi P_{R}+\bar{\Phi}P_{L}\right)\Psi-V(|\Phi|)\right], (32)

where PL,R=(1±γ5)/2P_{L,R}=(1\pm\gamma_{5})/2 are the chiral projectors and V(|Φ|)V(|\Phi|) is a symmetry-breaking potential for the complex field Φ\Phi whose explicit form is not important in the following. The action in Eq. (IV.1) is invariant under the global transformation

Ψeiβγ5Ψ,Ψ¯Ψ¯eiβγ5,Φe2iβΦ.\displaystyle\Psi\to e^{i\beta\gamma_{5}}\Psi,\quad\bar{\Psi}\to\bar{\Psi}e^{i\beta\gamma_{5}},\quad\Phi\to e^{-2i\beta}\Phi. (33)

However, the associated path integral is not invariant as the symmetry is anomalous at the quantum level due to the gravitational contribution to the chiral anomaly. This symmetry is spontaneously broken by the vacuum expectation value of Φ\Phi. Writing Φ(xμ)=(1/2)[v+ρ(xμ)]exp[iϕ(xμ)/v]\Phi(x^{\mu})=(1/\sqrt{2})[v+\rho(x^{\mu})]\exp\left[{i\phi(x^{\mu})/v}\right], with vv a constant, we obtain

S=d4xg\displaystyle S=\int d^{4}x\sqrt{-g} [Ψ¯eaμγa(μ+14ωμcdγcd)Ψ\displaystyle\left[-\bar{\Psi}e^{\mu}_{a}\gamma^{a}\left(\partial_{\mu}+\frac{1}{4}\omega^{cd}_{\mu}\gamma_{cd}\right)\Psi\right. (34)
12(ρ)212(1+ρ/v)2(ϕ)2\displaystyle\left.-\frac{1}{2}(\partial\rho)^{2}-\frac{1}{2}(1+\rho/v)^{2}(\partial\phi)^{2}\right.
y2(v+ρ)Ψ¯(eiϕ/vPR+eiϕ/vPL)ΨV(ρ)],\displaystyle\left.-\frac{y}{\sqrt{2}}(v+\rho)\bar{\Psi}\left(e^{i\phi/v}P_{R}+e^{-i\phi/v}P_{L}\right)\Psi-V(\rho)\right],

where ϕ\phi is massless as expected from a Goldstone mode. In order to obtain dCS gravity, we now integrate out the fermion Ψ\Psi and massive scalar ρ\rho. In doing so, we obtain an effective field theory (EFT) for the pseudoscalar ϕ\phi below the symmetry-breaking scale,

Z[ϕ]=[dΨ¯dΨdρ]exp{id4xg[Ψ¯(eaμγaμ+y2(v+ρ)ei(ϕ/v)γ5)Ψ12(ρ)212(1+ρ/v)2(ϕ)2V(ρ)]},Z[\phi]=\int[d\bar{\Psi}d\Psi d\rho]\exp\left\{i\int d^{4}x\sqrt{-g}\left[-\bar{\Psi}\left(e^{\mu}_{a}\gamma^{a}\nabla_{\mu}+\frac{y}{\sqrt{2}}(v+\rho)e^{i(\phi/v)\gamma_{5}}\right)\Psi-\frac{1}{2}(\partial\rho)^{2}-\frac{1}{2}(1+\rho/v)^{2}(\partial\phi)^{2}-V(\rho)\right]\right\}, (35)

where we simplified the Yukawa coupling by rearranging the phases times chiral projectors into exp[i(ϕ/v)γ5]\exp\left[{i(\phi/v)\gamma_{5}}\right]. Using the field redefinition

Ψeiϕ2vγ5Ψ,Ψ¯Ψ¯eiϕ2vγ5,\Psi\to e^{-i\frac{\phi}{2v}\gamma_{5}}\Psi,\quad\bar{\Psi}\to\bar{\Psi}e^{-i\frac{\phi}{2v}\gamma_{5}}, (36)

we get

Z[ϕ]=[dΨ¯dΨdρ]\displaystyle Z[\phi]=\int[d\bar{\Psi}d\Psi d\rho] exp{id4xg[Ψ¯(eaμγaμ+y2(v+ρ)i2vμϕγμγ5)Ψ12(ρ)212(1+ρ/v)2(ϕ)2\displaystyle\exp\left\{i\int d^{4}x\sqrt{-g}\left[-\bar{\Psi}\left(e^{\mu}_{a}\gamma^{a}\nabla_{\mu}+\frac{y}{\sqrt{2}}(v+\rho)-\frac{i}{2v}\partial_{\mu}\phi\gamma^{\mu}\gamma_{5}\right)\Psi-\frac{1}{2}(\partial\rho)^{2}-\frac{1}{2}(1+\rho/v)^{2}(\partial\phi)^{2}-\right.\right.
+V(ρ)+1192π2ϕ2vRR~]},\displaystyle+\left.\left.V(\rho)+\frac{1}{192\pi^{2}}\frac{\phi}{2v}R\tilde{R}\right]\right\}, (37)

where the last term in the exponential appears due to the non-trivial Jacobian of the path integral measure under the change of variables above and we identified the spinor-covariant derivative μΨ=(μ+ωμabγab/4)Ψ\nabla_{\mu}\Psi=(\partial_{\mu}+\omega_{\mu}^{ab}\gamma_{ab}/4)\Psi. To leading order, the path integral picks the classical saddle where Ψ=0=ρ\Psi=0=\rho, so the path integral then becomes

Z[ϕ]exp{id4xg[12(ϕ)2+ϕ384π2vRR~+]},Z[\phi]\propto\exp\left\{{i\int d^{4}x\sqrt{-g}\left[-\frac{1}{2}(\partial\phi)^{2}+\frac{\phi}{384\pi^{2}v}R\tilde{R}+\ldots\right]}\right\}, (38)

where dots stand for higher-derivative corrections suppressed by the fermionic and ρ\rho masses.555Since the fermionic path integral was done in curved space, there will be also extra curvature contributions, see Refs. [47, 48] and references therein. We see that the dCS gravity term is obtained from integrating-out heavy fermions. This model is thus analogous to the Kim-Shifman-Vainshtein-Zakharov (KSVZ) model for the QCD axion [49, 50]. Again, the shift symmetry of the dCS pseudoscalar is a non-linear realization of the non-anomalous phase transformation of Φ\Phi in Eq. (33). It is worth noting that the theories defined the action in Eq. (34) and Eq. (IV.1) are dual since their path integrals are the same (they differ by a fermionic Ψ\Psi field redefinition), with the latter having no λϕRR~\lambda\phi R\tilde{R} coupling. However, once we integrate out Ψ\Psi, there is no fermionic field redefinition left, and the λϕRR~\lambda\phi R\tilde{R} coupling is a genuine low-energy operator of the effective action in Eq. (38).

IV.2 Coupling a DCS sector to the SM

In what follows, we wish to make a connection with the Standard Model where the neutrinos get a Majorana mass term at energies below the symmetry-breaking scale, with ϕ\phi appearing as a phase in such a mass term. We do so after assuming that the heavy fermion Ψ\Psi plays the role of a right-handed sterile neutrino that couples with the lepton and Higgs doublets, LL and HH, respectively,

L=(LνLe),H~=(ϕ0ϕ+),L=\begin{pmatrix}L_{\nu}\\ L_{e}\end{pmatrix},\quad\tilde{H}=\begin{pmatrix}{\phi^{0}}^{*}\\ -{\phi^{+}}^{*}\end{pmatrix}, (39)

as

Δ=y~(L¯H~PRΨ+Ψ¯PLH~L),\Delta\mathcal{L}=-\tilde{y}(\bar{L}\tilde{H}P_{R}\Psi+\bar{\Psi}P_{L}\tilde{H}^{\dagger}L), (40)

where LνL_{\nu}, LeL_{e}, and Ψ\Psi are Majorana spinors:

Lf=(fLiσ2fL),Ψ=(iσ2νRνR).L_{f}=\begin{pmatrix}f_{L}\\ i\sigma^{2}f_{L}^{*}\end{pmatrix},\quad\Psi=\begin{pmatrix}-i\sigma^{2}\nu_{R}^{*}\\ \nu_{R}\end{pmatrix}. (41)

Note that we need the conjugated H~\tilde{H} instead of the Higgs doublet HT=(ϕ0ϕ+)H^{T}=(\phi^{0}\;\;\phi^{+}) in Eq. (40) to get a hypercharge neutral coupling. We consider a single lepton flavor for simplicity since a generalization for all three generations is straightforward (see the discussion section).

With the new Yukawa coupling (40) and the kinetic terms for LL and HH, the action in Eq. (IV.1) becomes

S=d4xg\displaystyle S=\int d^{4}x\sqrt{-g} [12Ψ¯∇̸ΨμΦμΦ¯12L¯∇̸L\displaystyle\left[-\frac{1}{2}\bar{\Psi}\not{\nabla}\Psi-\partial_{\mu}\Phi\partial^{\mu}\bar{\Phi}-\frac{1}{2}\bar{L}\not{\nabla}L-\right.
|H|2yΨ¯(ΦPR+Φ¯PL)Ψ\displaystyle\left.-|\partial H|^{2}-y\bar{\Psi}\left(\Phi P_{R}+\bar{\Phi}P_{L}\right)\Psi-\right.
y~(L¯H~PRΨ+Ψ¯PLH~L)V(|Φ|)].\displaystyle\left.-\tilde{y}(\bar{L}\tilde{H}P_{R}\Psi+\bar{\Psi}P_{L}\tilde{H}^{\dagger}L)-V(|\Phi|)\right]. (42)

This action is invariant under Eq. (33) provided LL transforms as

Leiβγ5L.\displaystyle L\to e^{i\beta\gamma_{5}}L. (43)

Note that since both LL and ψ\psi are four-component Majorana spinors, the global chiral transformation acting on those is a phase rotation associated with their fermionic number. We discuss further aspects of this symmetry in the discussion section. After performing the field redefinition in Eq. (36), the path integral is now

Z=[dψ¯dψda]exp\displaystyle Z=\int[d\bar{\psi}d\psi da]\text{exp} {id4xg[12Ψ¯(eaμγaμ+y2(v+ρ)i2vμϕγμγ5)Ψ12(ρ)2\displaystyle\left\{i\int d^{4}x\sqrt{-g}\left[-\frac{1}{2}\bar{\Psi}\left(e^{\mu}_{a}\gamma^{a}\nabla_{\mu}+\frac{y}{\sqrt{2}}(v+\rho)-\frac{i}{2v}\partial_{\mu}\phi\gamma^{\mu}\gamma_{5}\right)\Psi-\frac{1}{2}(\partial\rho)^{2}-\right.\right.
|H|2L¯eaμγaμL12(1+ρ/v)2(ϕ)2V(ρ)+1192π2ϕ2vRR~\displaystyle\left.-|\partial H|^{2}-\bar{L}e^{\mu}_{a}\gamma^{a}\nabla_{\mu}L-\frac{1}{2}(1+\rho/v)^{2}(\partial\phi)^{2}-V(\rho)+\frac{1}{192\pi^{2}}\frac{\phi}{2v}R\tilde{R}-\right.
y~(L¯H~eiϕ2vγ5PRΨ+Ψ¯PLeiϕ2vγ5H~L)]}.\displaystyle\left.\left.-\tilde{y}\left(\bar{L}\tilde{H}e^{-i\frac{\phi}{2v}\gamma_{5}}P_{R}\Psi+\bar{\Psi}P_{L}e^{-i\frac{\phi}{2v}\gamma_{5}}\tilde{H}^{\dagger}L\right)\right]\right\}. (44)

We thus see that χ¯=y~L¯H~eiϕ2vγ5\bar{\chi}=-\tilde{y}\bar{L}\tilde{H}e^{-i\frac{\phi}{2v}\gamma_{5}} acts like a source for the right-hand part of Ψ\Psi. Moreover, using L¯=Lc=LTC\bar{L}=L^{c}=L^{T}C, where CC is the charge conjugation matrix, we have

Ψ¯PLeiϕ2vγ5H~L=ΨTPLeiϕ2vγ5H~L¯T\bar{\Psi}P_{L}e^{-i\frac{\phi}{2v}\gamma_{5}}\tilde{H}^{\dagger}L=-\Psi^{T}P_{L}e^{-i\frac{\phi}{2v}\gamma_{5}}\tilde{H}^{\dagger}\bar{L}^{T} (45)

and so η=eiϕ2vγ5H~L\eta=e^{-i\frac{\phi}{2v}\gamma_{5}}\tilde{H}^{\dagger}L acts like a source for the left-hand part of Ψ\Psi. But since Ψ\Psi is a Majorana spinor, its left and right components are related, and both χ\chi and η\eta contribute as sources to the path integral of Ψ\Psi. Hence, the coupling in Eq. (40) will give the contribution

12χ¯(M)1PRχ12η¯(M)1PLη\displaystyle-\frac{1}{2}\bar{\chi}(\not{D}-M)^{-1}P_{R}\chi-\frac{1}{2}\bar{\eta}(\not{D}-M)^{-1}P_{L}\eta
=12mΨ(χ¯PRχ+η¯PLη)+\displaystyle=-\frac{1}{2m_{\Psi}}\left(\bar{\chi}P_{R}\chi+\bar{\eta}P_{L}\eta\right)+\cdots (46)

to the effective Lagrangian. Note that the charge conjugation of the chiral sources fixes χ\chi and η\eta,

χ\displaystyle\chi Cχ¯T=y~eiϕ2vγ5H~TL,\displaystyle\equiv C\bar{\chi}^{T}=-\tilde{y}e^{-i\frac{\phi}{2v}\gamma_{5}}\tilde{H}^{T}L, (47)
η¯\displaystyle\bar{\eta} ηTC=L¯H~eiϕ2vγ5.\displaystyle\equiv\eta^{T}C=\bar{L}\tilde{H}^{*}e^{-i\frac{\phi}{2v}\gamma_{5}}. (48)

The effective action is given by

S\displaystyle S =d4xg[12(ϕ)2|H|212L¯eaμγaμL\displaystyle=\int d^{4}x\sqrt{-g}\left[-\frac{1}{2}(\partial\phi)^{2}-|\partial H|^{2}-\frac{1}{2}\bar{L}e^{\mu}_{a}\gamma^{a}\nabla_{\mu}L-\right.
y~22mΨ(eiϕvL¯H~PRH~TL+eiϕvL¯H~PLH~L)\displaystyle-\left.\frac{\tilde{y}^{2}}{2m_{\Psi}}\left(e^{i\frac{\phi}{v}}\bar{L}\tilde{H}P_{R}\tilde{H}^{T}L+e^{-i\frac{\phi}{v}}\bar{L}\tilde{H}^{*}P_{L}\tilde{H}^{\dagger}L\right)-\right.
ϕ384π2vRR~+],\displaystyle\left.-\frac{\phi}{384\pi^{2}v}R\tilde{R}+\cdots\right], (49)

where the heavy fermion mass mΨm_{\Psi} can be read from Eq. (35) to be mΨ=yv/2m_{\Psi}=yv/\sqrt{2}. The resulting effective action has a modified Weinberg operator [51] that includes a coupling with ϕ\phi.

Notice that the pseudoscalar arises as a complex phase in the neutrino mass matrix and as a coupling to the gravitational Pontryagin density. A measurement of the mass phase will not constrain its coupling to the gravitational Pontryagin density because strong gravitational effects from compact binary can source a large gravitational Pontryagin density which will source large gradients of ϕ\phi.

This is analog to QCD, as the axion couples to chromo-electric and magnetic field via the Pontryagin density, astrophysical sources may enhance the axion and a similar story applies to dCS. Here, there are two observable effects to consider. First, an Earth-based measurement of the dynamical field can constrain the value of the pseudoscalar, which is predicted to be small due to the fact that the Pontryagin term is zero in Minkowski space-time. However in compact binary systems strong gravity can give a non-vanishing RR~R\tilde{R} which will in turn source gravitational waves and non-vanishing pseudoscalar amplitude. Also, in models of Higgs-Inflation, Electroweak symmetry is broken during inflation, and the phase can also be non-vanishing in the early universe. In these strong gravity, early universe regimes the mass matrix phase will also be non-vanishing but obviously difficult to measure. Still, it may be possible to study the effects of the phase and its impact on neutrino cosmology. We leave this possibility for future investigation.

V Discussion and Conclusion

In this paper, we have discussed many symmetry aspects of dCS gravity and anomalies that are relevant when trying to get the dCS scalar-gravity coupling from a UV theory. Most of the necessary concepts and results have analogs in the physics of axions and their coupling with gauge fields. We reviewed some known results on axions in detail to better understand their relation and applicability to the gravitational case. We also proposed an extension of the SM model that gives origin to a dCS sector and the coupling of the dCS pseudoscalar with neutrinos at low energies. In the rest of this section, we discuss the model’s symmetry and its gauging.

The modified Weinberg operator in Eq. (IV.2) breaks the lepton number symmetry (as the SM lepton number symmetry alone is anomalous, we consider the baryonic minus lepton number symmetry). But there is a combination of the U(1)BL(1)_{\text{B}-\text{L}} and the ϕ\phi shift symmetry that leaves the modified Weinberg term invariant and is a quasi-symmetry of the action Eq. (IV.2). However, such a symmetry is anomalous and the ϕRR~\phi R\tilde{R} coupling cannot be eliminated by redefining LL in the path integral for the same reasons explained in Section III: such a redefinition would bring a term proportional to the divergence of the LL’s chiral current (after integration by parts), giving back the dCS coupling. We comment on how the symmetry combination above can be gauged from the low-energy effective perspective in Appendix A.

On the other hand, the action in Eq. (IV.2) is invariant under the simultaneous rotation of Ψ\Psi, LL, and Φ\Phi (with appropriate charges) Eq. (33) and Eq. (43), and such symmetry is non-anomalous because the gravitational contribution of LL and Ψ\Psi to the anomaly cancel each other. Hence, we can gauge this non-anomalous symmetry combination, which we call U(1)BLϕ(1)_{\text{B}-\text{L}-\phi}. This begs the question of what happens with this gauge symmetry after integrating Ψ\Psi. At energies well below the Ψ\Psi mass, the only fermion left in the spectrum is the chiral LL, which contributes to the anomaly of the gauge symmetry. As discussed in Ref. [52], despite this anomaly, this effective chiral gauge theory is unitary (but not renormalizable) because the gauge boson acquires its mass via a Higgs mechanism: the gauge symmetry is non-linearly realized as a shift symmetry of the gauge boson longitudinal mode. The main consequence of this fact is that ϕ\phi appears in the EFT as the longitudinal mode of the gauge boson in a Stueckelberg (Aμ,ϕ)(A^{\mu},\phi) sector and hence gets eaten by the gauge field in the unitary gauge. Moreover, before integrating Ψ\Psi, the Lagrangian includes the term AμjΨμ=AμiΨ¯γμγ5ΨA_{\mu}j^{\mu}_{\Psi}=A_{\mu}i\bar{\Psi}\gamma^{\mu}\gamma_{5}\Psi and the gauge transformation of the effective action necessary to select the unitary gauge will give rise to a term proportional to ϕμjΨμ\phi\partial_{\mu}j^{\mu}_{\Psi}, which after the Ψ\Psi path-integration will appear in the effective Lagrangian as ϕ(RR~FF~)\phi(R\tilde{R}-F\tilde{F}). The gauge variation of this term cancels the non-trivial Jacobian from the LL’s fermionic measure, making the effective action gauge invariant at the quantum level.

The results in the previous sections can be straightforwardly generalized to an arbitrary number of heavy fermions. If we consider the case of all left-handed SM leptons LiL_{i} and three sterile right-handed neutrinos Ψi\Psi_{i}, we have the simplest Majoron model [53, 54], where ϕ\phi play the role of a Majoron666See Refs. [55, 56] for a proposal of identifying the QCD axion with the Majoron.. The yy and y~\tilde{y} Yukawa couplings will be promoted to matrices yijy_{ij} and y~ij\tilde{y}_{ij} with flavor indices and the Weinberg operator in the effective Lagrangian would read

y~imMmn1y~nj(eiϕvLi¯H~PRH~TLj+eiϕvL¯iH~PLH~Lj),-\tilde{y}_{im}M^{-1}_{mn}\tilde{y}_{nj}\left(e^{i\frac{\phi}{v}}\bar{L_{i}}\tilde{H}P_{R}\tilde{H}^{T}L_{j}+e^{-i\frac{\phi}{v}}\bar{L}_{i}\tilde{H}^{*}P_{L}\tilde{H}^{\dagger}L_{j}\right), (50)

withe Mmn1=ymnv/2M_{mn}^{-1}=y_{mn}v/\sqrt{2}. An interesting question is how ϕ\phi affects the seesaw mechanism, a topic that we leave for future investigations.

The Majoron model above has the global non-anomalous symmetry U(1)BLϕ(1)_{\text{B}-\text{L}-\phi} that, if gauged, would hide the dCS scalar ϕ\phi in the longitudinal mode of a massive gauge boson777For gauged Majoron models, see Refs. [57, 58]. Note that if we start with a different number of heavy fermions, unmatched with the three left-handed neutrinos, we cannot straightforwardly gauge U(1)BLϕ(1)_{\text{B}-\text{L}-\phi} because this symmetry will be anomalous from the on-set. Moreover, even for three right-handed neutrinos, gauging the symmetry is a choice888We thank Michael Peskin for bringing up this point to our attention. and there is nothing inconsistent with having a global, non-anomalous, U(1)BLϕ(1)_{\text{B}-\text{L}-\phi}. In this case, ϕ\phi is a dynamical field that couples with the Pontryagin density. This is analogous to the PQ symmetry and axions: although there are models where U(1)PQ(1)_{\text{PQ}} is gauged and the axion is eaten to give the gauge boson a mass, this has not precluded the search for axions in particle physics, astrophysics and cosmology.

Although there are motivations for gauging global symmetries from the string theory perspective, the same point of view also motivates us to rethink the origin of the dCS: we know that a dCS sector should appear in any heterotic four-dimensional compactification due to the Green-Schwarz mechanism (see discussion in Ref. [1]). This is because ϕ\phi can be identified with the zero mode of the ten-dimensional Kalb-Rammond field on the internal space, and this identification is compactification-independent (ϕ\phi is referred to as the model-independent axion) [59]. Then, in the most basic example of how to get dCS from string theory, the way the global dCS shift symmetry becomes local is by a promotion to the gauge symmetry of a 2-form at energies close to and larger than the compactification scale. Thus, at least from a theoretical perspective, a dCS sector can appear at low energies without any inconsistency between its global symmetry and quantum gravity. More realistically, since ϕ\phi pairs with the dilaton, it gets a potential (generated by non-perturbative effects, see Ref. [60] for a review) that breaks its continuous global symmetry down to a discrete group, again in analogy to the PQ axion, and so it cannot be the longitudinal mode of any gauge field.

Acknowledgements.
We thank Chris Hull for enlightening discussions, Peter Adshead, Robert Brandenberger, Keshav Dasgupta, Jim Gates, Greg Gabadadze, Humberto Gilmer, Michael Peskin, Sanjaye Ramgoolam, and Matt Strassler for discussion and feedback on a draft of this paper.

Appendix A Baryon and Lepton Number Anomalies in the SM

In this appendix, we state some results on global anomalies in the Standard Model and comment on how the anomaly changes after introducing a dCS sector. The electroweak sector is a chiral gauge theory due to the chiral coupling of SU(2)L(2)_{\text{L}} gauge bosons with quarks and leptons. As it is well known, chiral gauge theories have gauge anomalies which have to be cancelled, and such is the case in the standard model after summing over all representations and hypercharges of the fermions. Chiral gauge theories have also a fermion number global anomaly, and in the SM, the global baryonic and leptonic number symmetries are indeed anomalous. This is because under redefinitions of the form

ψL(x)ψ~L(x)\displaystyle\psi_{L}(x)\to\tilde{\psi}_{L}(x) =eiα(x)ψL(x),\displaystyle=e^{i\alpha(x)}\psi_{L}(x), (51a)
ψR(x)ψ~R(x)\displaystyle\psi_{R}(x)\to\tilde{\psi}_{R}(x) =eiα(x)ψR(x),\displaystyle=e^{i\alpha(x)}\psi_{R}(x), (51b)

the left and right contributions to the Jacobian of the path integral measure do not cancel when only one component (left or right) couple with the gauge field. So, the number current jnμ(x)=iψ¯(x)γμψ(x)j_{n}^{\mu}(x)=i\bar{\psi}(x)\gamma^{\mu}\psi(x) will satisfy [44]

μjnμ\displaystyle\partial_{\mu}\langle j_{n}^{\mu}\rangle =μ(iψ¯LγμψL+iψ¯RγμψR)\displaystyle=\partial_{\mu}\left(\langle i\bar{\psi}_{L}\gamma^{\mu}\psi_{L}\rangle+\langle i\bar{\psi}_{R}\gamma^{\mu}\psi_{R}\rangle\right)
=±tr(TaTb)16π2F~μνaFμνb,\displaystyle=\pm\frac{\text{tr}(T^{a}T^{b})}{16\pi^{2}}\tilde{F}^{a}_{\mu\nu}F^{\mu\nu b}, (52)

where the negative (positive) sign case is for when only the left (right) chirality couples with the gauge field. The trace is over the representation of the ψ\psi chirality that contributes to the anomaly. In the SM, only the left-hand part of the quark and lepton field doublets qq and ll couple with SU(2)L(2)_{\text{L}}:

SSMd4x\displaystyle S_{\text{SM}}\supset\int d^{4}x [q¯Liγμ(μ+igAμaσa2)qLi\displaystyle\left[-\bar{q}_{L}^{i}\gamma_{\mu}\left(\partial_{\mu}+igA^{a}_{\mu}\frac{\sigma^{a}}{2}\right)q_{L}^{i}-\right.
l¯Liγμ(μ+igAμaσa2)lLi].\displaystyle\left.-\bar{l}_{L}^{i}\gamma_{\mu}\left(\partial_{\mu}+igA^{a}_{\mu}\frac{\sigma^{a}}{2}\right)l^{i}_{L}\right]. (53)

The indices in the doublets are generation indices and the fundamental-representation indices are suppressed. For instance,

lL1=(νeeL),qL1=(uLdL).l_{L}^{1}=\begin{pmatrix}\nu_{e}\\ e_{L}\end{pmatrix},\quad q_{L}^{1}=\begin{pmatrix}u_{L}\\ d_{L}\end{pmatrix}. (54)

The SM baryon number current satisfies

μjBμ=3Nc32π2F~μνaFμνa,\partial_{\mu}\langle j_{\text{B}}^{\mu}\rangle=\frac{3N_{c}}{32\pi^{2}}\tilde{F}_{\mu\nu}^{a}F^{\mu\nu a}, (55)

where NcN_{c} is the number of quark colors. The lepton number current will also be anomalous

μjLμ=332π2F~μνaFμνa.\partial_{\mu}\langle j_{\text{L}}^{\mu}\rangle=\frac{3}{32\pi^{2}}\tilde{F}_{\mu\nu}^{a}F^{\mu\nu a}. (56)

The factor of 33 in the numerator of the anomalies above comes from the contributions of the 33 quark and lepton generations. This tells us that the baryon and lepton numbers are not conserved and, in particular, their associated U(1)(1) transformations cannot be straightforwardly gauged.

However, the combination jBLμ=jBμ/3jLμj_{\text{B}-\text{L}}^{\mu}=j_{\text{B}}^{\mu}/3-j^{\mu}_{\text{L}} satisfies

μjBLμ=Nc332π2F~μνaFμνa,\partial_{\mu}\langle j^{\mu}_{\text{B}-\text{L}}\rangle=\frac{N_{c}-3}{32\pi^{2}}\tilde{F}_{\mu\nu}^{a}F^{\mu\nu a}, (57)

and so is conserved in the SM, since Nc=3N_{c}=3 (note that the factor of 1/31/3 in the baryonic contribution to jBLμj_{\text{B}-\text{L}}^{\mu} is because quarks have a baryon number B=1/3B=1/3). So, the Standard Model has a potential non-anomalous global U(1)BL(1)_{\text{B}-\text{L}} symmetry. In beyond SM models, one can add extra breaking effects for U(1)B(1)_{\text{B}} and U(1)L(1)_{\text{L}} individually, but the U(1)BL(1)_{\text{B}-\text{L}} should be non-anomalous. One can even contemplate gauging such a symmetry.

The discussion above is modified when we consider the Standard Model in curved spacetimes. Since there are no right-handed neutrinos, there is a mismatch in the number of left and right chiralities in the fermionic sector. Recall that a single chiral fermion contributes to the gravitational anomaly of an otherwise non-anomalous current as [43]

μiψ¯L,RγμψL,R=±1384π2R~μνρσRμνρσ,\nabla_{\mu}\langle i\bar{\psi}_{L,R}\gamma^{\mu}\psi_{L,R}\rangle=\pm\frac{1}{384\pi^{2}}\tilde{R}_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma}, (58)

where the negative (positive) sign appears for the left (right) chirality contribution. So, for a non-chiral fermion with propagating left and right-parts, we have

μiψ¯γμψ\displaystyle\nabla_{\mu}\langle i\bar{\psi}\gamma^{\mu}\psi\rangle =μ[iψ¯LγμψL+iψ¯RγμψR]=0,\displaystyle=\nabla_{\mu}\left[\langle i\bar{\psi}_{L}\gamma^{\mu}\psi_{L}\rangle+\langle i\bar{\psi}_{R}\gamma^{\mu}\psi_{R}\rangle\right]=0, (59a)
μiψ¯γμγ5ψ\displaystyle\nabla_{\mu}\langle i\bar{\psi}\gamma^{\mu}\gamma^{5}\psi\rangle =μ[iψ¯LγμψLiψ¯RγμψR]\displaystyle=\nabla_{\mu}\left[\langle i\bar{\psi}_{L}\gamma^{\mu}\psi_{L}\rangle-\langle i\bar{\psi}_{R}\gamma^{\mu}\psi_{R}\rangle\right]
=1192π2R~μνρσRμνρσ.\displaystyle=\frac{1}{192\pi^{2}}\tilde{R}_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma}. (59b)

Let us apply these results to jBLμj^{\mu}_{\text{B}-\text{L}} in the SM. The left-handed contribution to the lepton number overcomes the right-handed one, such that the total lepton number current satisfies

μjLμ=332π2F~μνaFμνa+3384π2R~μνρσRμνρσ.\nabla_{\mu}\langle j^{\mu}_{L}\rangle=\frac{3}{32\pi^{2}}\tilde{F}_{\mu\nu}^{a}F^{\mu\nu a}+\frac{3}{384\pi^{2}}\tilde{R}_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma}. (60)

where the factor of 33 in the numerator of the gravitational contribution is the excess of left over right neutrinos. However, since the number of left and right-handed quarks is the same, there is no gravitational contribution to the jBμj_{\text{B}}^{\mu} anomaly. Hence

μjBLμ=3384π2R~μνρσRμνρσ.\nabla_{\mu}\langle j^{\mu}_{B-L}\rangle=\frac{3}{384\pi^{2}}\tilde{R}_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma}. (61)

In other words, since the lepton number current anomaly has a gravitational contribution while the baryonic current does not possess one, the gravitational contribution to the anomaly in jBLμj^{\mu}_{\text{B}-\text{L}} cannot be canceled.

Now, after introducing the dCS sector, following the discussion of Section III, there will be a combination of U(1)dCS(1)_{\text{dCS}} and U(1)BL(1)_{\text{B}-\text{L}} that is not affected by instanton effects at the quantum level, i.e., it is non-anomalous. We will call this combination U(1)BLϕ(1)_{\text{B}-\text{L}-\phi}. Then, a constant part of ϕ\phi can always be absorbed in the phase of the fermions, and only a combination of such a zero mode and the phase of the mass matrix determinant is observable. This is very similar to the CP violation in QCD with the zero mode of ϕ\phi playing the role of the QCD θ\theta.

In some extensions of the Standard Model, the U(1)BL(1)_{\text{B}-\text{L}} is ultimately gauged, leading to new effects. While this can be done in flat space, the gravitational contribution to the chiral anomaly must be consistently canceled for the gauging in curved spacetimes. However, with an added dCS sector, the non-anomalous U(1)BLϕ(1)_{\text{B}-\text{L}-\phi} can seemingly be gauged without any constraints on the background, which looks attractive, for instance, if one wants to contemplate quantizing gravity. Such a procedure is subtle because despite U(1)BLϕ(1)_{\text{B}-\text{L}-\phi} being non-anomalous at the quantum level, there is no contribution from the chiral anomaly to the conservation of the classical jBLϕμj^{\mu}_{\text{B}-\text{L}-\phi} current. In other words, the partition function would be gauge invariant but not the classical action. Notwithstanding, global anomalous symmetries can be gauged at the expense of non-renormalizability [52, 61], so the resulting gauge theory should be an EFT valid only below some cut off scale.

References