---
title: 'Leptophilic ALPs: Overview & Experimental Prospects'
url: https://www.emergentmind.com/topics/leptophilic-axion-like-particles-alps
type: topic
---

# Leptophilic ALPs: Overview & Experimental Prospects

Leptophilic axion-like particles (ALPs) are pseudoscalar states associated with a broken global symmetry whose effective interactions are concentrated in the lepton sector. In contemporary effective descriptions they are written either with derivative couplings to lepton currents or, after integration by parts and use of equations of motion, with pseudoscalar couplings proportional to lepton masses, accompanied in many treatments by anomaly-induced couplings to electroweak gauge bosons. Their phenomenology is correspondingly shaped by charged-lepton masses, flavor structure, and the presence or absence of electroweak contact terms, yielding distinct signatures in exotic \(\tau\) decays, charged-lepton-flavor-violation observables, beam dumps, forward detectors, reactor experiments, and anomalous magnetic moments [2509.23949, 2504.05873, 2507.15271].

## 1. Effective-field-theory formulation

A standard starting point is a derivative coupling between the ALP field \(a\) and a leptonic Peccei–Quinn current,
\[
{\cal L}\supset \partial_\mu a\,J^\mu_{\rm PQ,\ell},
\]
with
\[
J^{\mu}_{\rm PQ,\ell}=\frac{c^V_{\ell}}{2\Lambda}\,\bar\ell\gamma^{\mu}\ell
+\frac{c^A_{\ell}}{2\Lambda}\,\bar\ell\gamma^{\mu}\gamma^5\ell
+\frac{c_\nu}{2\Lambda}\,\bar\nu_\ell\gamma^{\mu}P_L\nu_\ell.
\]
After integration by parts and use of the equations of motion, this produces a pseudoscalar lepton coupling \(i\,c^A_{\ell}\,\frac{m_\ell}{\Lambda}\,a\,\bar\ell\gamma^5\ell\), a four-point \(a\bar\ell\gamma^\mu P_L\nu_\ell W^-_\mu\) interaction proportional to \(c^A_{\ell}-c^V_{\ell}+c_\nu\), and anomaly-induced couplings to \(F\tilde F\), \(F\tilde Z\), \(Z\tilde Z\), and \(W^+\tilde W^-\) [2509.23949].

A complementary organization uses two EFT bases for charged leptons. In the derivative basis,
\[
{\cal L}_{\rm D}
=
\tfrac12\,\partial_\mu a\,\partial^\mu a
-\tfrac12\,m_a^2\,a^2
+\frac{\partial_\mu a}{2\Lambda}\sum_{\ell}C_\ell\;\bar\psi_\ell\gamma^\mu\gamma^5\psi_\ell,
\]
whereas in the pseudoscalar basis,
\[
{\cal L}_{\rm P}
=
\tfrac12\,(\partial_\mu a)^2
-\tfrac12\,m_a^2\,a^2
-i \,\sum_\ell g_{a\ell}\,a\,\bar\psi_\ell\gamma^5\psi_\ell,
\qquad
g_{a\ell}\equiv \frac{m_\ell C_\ell}{\Lambda}.
\]
Up to field redefinitions both Lagrangians yield identical tree-level ALP–lepton amplitudes; they differ when ALP–photon interactions are probed. In the derivative basis a chiral field redefinition generates a tree-level \(aF^{\mu\nu}\tilde F_{\mu\nu}\) term, while in the pseudoscalar basis there is no tree-level ALP–photon coupling [2504.05873].

This EFT structure fixes two recurrent features of leptophilic ALP phenomenology. First, direct couplings to charged leptons scale with \(m_\ell\), so rates and lifetimes depend strongly on flavor thresholds. Second, the four-point \(a\ell\nu W\) vertex is not a minor detail: in several benchmark scenarios it controls production, especially in exotic meson and \(\tau\) decays [2412.19195].

## 2. Coupling benchmarks and flavor structure

Several benchmark patterns recur across the recent literature. In the \(\tau\)-factory analysis, the electroweak-violating tauphilic benchmark is defined by \(c^V_\tau=c_\nu=0\), \(c^A_\tau\neq0\), with couplings to \(e\) and \(\mu\) set to zero; only the \(a\tau\tau\) axial-vector coupling and anomaly terms survive. The electroweak-preserving right-handed tauphilic benchmark instead takes \(c_\nu=0\) and \(c^V_\tau=c^A_\tau\neq0\), so the \(aW\tau\nu\) four-point interaction is absent. The lepton-flavor-universal benchmark imposes \(c^{V,A}_{e}=c^{V,A}_{\mu}=c^{V,A}_{\tau}\neq0\) and \(c_\nu=c^A_\ell-c^V_\ell\), so that all three charged leptons couple equally and the low-energy \(aW\ell\nu\) amplitudes cancel [2509.23949].

Forward-detector studies organize the same physics somewhat differently. There the electroweak-violating benchmark is \(c^V_\ell=0\), \(c_\nu=0\), \(c^A_\ell\neq0\), which yields a large four-point \(a\ell\nu W\) coupling unsuppressed by \(m_\ell\). The electroweak-preserving benchmark sets \(c^V_\ell=c^A_\ell\neq0\), \(c_\nu=0\), leaving only the mass-suppressed \(a\bar\ell\ell\) and anomaly terms. A third option, the left-right softly-asymmetric model, takes \(c^V_\ell\neq c^A_\ell\) with \(\Delta\equiv |(c_L-c_R)/(c_L+c_R)|\ll1\); a benchmark value \(\Delta=0.1\) is used to enhance production through the four-point interaction without shortening the ALP lifetime through the pseudoscalar decay coupling [2412.19195].

Flavor structure can also be expressed directly in the charged-lepton indices. A general derivative interaction,
\[
\mathcal{L}
=
\frac12(\partial_\mu a)(\partial^\mu a)-\frac12\,m_a^2\,a^2
+\frac{\partial_\mu a}{2\,f_a}\sum_{i,j=1}^3
\Bigl[g_L^{ij}\,\overline{L}_i\gamma^\mu P_L L_j+g_R^{ij}\,\overline{\ell}_i\gamma^\mu P_R \ell_j\Bigr],
\]
supports anarchic LFV, hierarchical LFV, universal LFC, and non-universal LFC patterns. In much of that phenomenology one takes \(g_L^{ij}=-g_R^{ij}\equiv g_{ij}\), so that the interaction is purely axial in the charged-lepton sector [2507.15271].

In anomaly-free flavored models the coupling matrices are built from symmetry charges and the charged-lepton rotation matrices,
\[
C^V_{ij}=\tfrac12\bigl[(U_R^{e\dagger}x^R U_R^e)_{ij}+(U_L^{e\dagger}x^L U_L^e)_{ij}\bigr],
\qquad
C^A_{ij}=\tfrac12\bigl[(U_R^{e\dagger}x^R U_R^e)_{ij}-(U_L^{e\dagger}x^L U_L^e)_{ij}\bigr].
\]
A sufficient condition for absence of an electromagnetic anomaly with only SM fermions plus right-handed neutrinos is \(\sum_i x^L_i=0\) and \(\sum_i x^R_i=0\) [2007.08834].

## 3. Decays, branching fractions, and lifetime systematics

For same-flavor charged-lepton decays, the standard width is
\[
\Gamma(a\to\ell^+\ell^-)
=
\frac{g_{a\ell}^2\,m_a}{8\pi}
\sqrt{1-\frac{4m_\ell^2}{m_a^2}},
\]
or equivalently
\[
\Gamma(a\to\ell^+\ell^-)
=
\frac{g_{\ell\ell}^2}{8\pi f_a^2}\,m_\ell^2\,m_a\,
\sqrt{1-\frac{4m_\ell^2}{m_a^2}}.
\]
For flavor-off-diagonal channels,
\[
\Gamma(a\to\ell_1^\pm\ell_2^\mp)
=
\frac{g_{\ell_1\ell_2}^2}{16\pi f_a^2}\,(m_{\ell_1}+m_{\ell_2})^2\,m_a
\left[1-\frac{(m_{\ell_1}-m_{\ell_2})^2}{m_a^2}\right]^{1/2}
\left[1-\frac{(m_{\ell_1}+m_{\ell_2})^2}{m_a^2}\right]^{1/2}.
\]
The loop-induced diphoton mode is commonly written as
\[
\Gamma(a\to\gamma\gamma)=\frac{|g_{a\gamma\gamma}^{\rm eff}|^2\,m_a^3}{64\pi},
\]
or, in a charged-lepton loop representation,
\[
\Gamma(a\to\gamma\gamma)=
\frac{\alpha^2 m_a^3}{64\pi^3 f_a^2}
\Biggl|\sum_{\ell=e,\mu,\tau}g_{\ell\ell}\,
I_{\gamma\gamma}\!\Bigl(\tfrac{m_a}{2m_\ell}\Bigr)\Biggr|^2.
\]
The total width \(\Gamma_{\rm tot}\) fixes the proper lifetime through \(\tau_a=\hbar/\Gamma_{\rm tot}\), and in the derivative formulations \(\tau_a\propto \Lambda^2/[(c^A_\ell)^2 m_\ell^2 m_a]\) [2504.05873, 2507.15271, 2412.19195].

These formulae encode a strong threshold dependence. Below the corresponding \(2m_\ell\) threshold, the tree-level dilepton channel is closed and the lifetime can become long enough for displaced or far-detector searches. In the tauphilic and LFU benchmarks studied for STCF, the lifetime behavior differs markedly: the figure-level summary states that the tauphilic benchmark remains much longer-lived, while in the LFU case the \(e^+e^-\) and \(\mu^+\mu^-\) channels can shorten the lifetime dramatically [2509.23949].

For universal LFC couplings the low-mass diphoton channel may dominate below \(2m_\mu\), while for a \(\tau\)-philic ALP below \(2m_\tau\) the decay is only \(a\to\gamma\gamma\), with \(\Gamma\propto m_a^7\); the lifetime therefore grows rapidly. This feature is central to the far-detector reach of beam dumps and SHiP-like facilities [2507.15271].

## 4. Production channels

Exotic \(\tau\) decays are among the cleanest direct production mechanisms. In the leptonic four-body mode,
\[
d\Gamma(\tau^-\to\ell^-\bar\nu_\ell\nu_\tau a)
=
\frac1{2m_\tau}\,\overline{|M|^2}\,d\Phi_4,
\]
and the integrated rate is parameterized as
\[
\Gamma_{\tau\to\ell\nu\nu a}
=
c_\tau^{A\,2}\;\frac{m_\tau^5}{512\pi^5\Lambda^2}\;I\bigl(m_a/m_\tau\bigr).
\]
For \(c_\tau^A=1\), the corresponding branching ratios satisfy
\[
{\rm BR}(\tau^-\to e^-\bar\nu_e\nu_\tau a)
\simeq
{\rm BR}(\tau^-\to \mu^-\bar\nu_\mu\nu_\tau a)
\sim 10^{-10}-10^{-8},
\]
nearly flat for \(m_a\lesssim100\) MeV and then falling for heavier \(m_a\). The relevant diagrams are ALP bremsstrahlung off the \(\tau\) leg, a four-point \(aW\tau\nu\) contact term present only in the electroweak-violating case, and an \(aW\ell\nu\) contact term present in LFU but cancelled at low energy against the previous contribution [2509.23949].

If flavor-violating couplings \(g_{\tau\ell}\) are present, two-body \(\tau\to \ell a\) dominates:
\[
\Gamma(\tau\to\ell a)
=
\frac{|g_{\ell\tau}|^2}{64\pi f_a^2}\,m_\tau^3
\Bigl(1+\tfrac{m_\ell}{m_\tau}\Bigr)^2
\Bigl[\Bigl(1-\tfrac{m_\ell}{m_\tau}\Bigr)^2-\tfrac{m_a^2}{m_\tau^2}\Bigr]
\sqrt{\lambda(m_\tau^2,m_\ell^2,m_a^2)}.
\]
For \(m_a\ll m_\tau\),
\[
{\rm Br}(\tau\to\ell a)\sim1.2\times10^{-4}\,
\biggl(\frac{10^7\,\rm GeV}{f_a}\biggr)^2|g_{\ell\tau}|^2.
\]
If \(g_{\tau\ell}\simeq0\), production proceeds instead through \(\tau\)-line bremsstrahlung in ordinary decays such as \(\tau^-\to\pi^-\nu_\tau a\), \(\tau\to\pi\pi^0\nu_\tau a\), and \(\tau\to\ell\nu\nu a\), with branching ratios around \(\sim10^{-12}(10^4/f_a)^2\) in the quoted examples [2507.15271].

At hadron facilities the four-point interaction also drives exotic meson decays. The modes emphasized in forward-detector studies include \(B^0\to K^* a\), \(K^+\to\pi^+ a\), \(K_L\to\pi^0 a\), \(K_S\to\pi^0 a\), and \(P^+\to\ell^+\nu_\ell a\) with \(P=\pi,K,D,D_s\). In electroweak-violating scenarios the three-body charged-meson decays are enhanced because the \((c^A_\ell-c^V_\ell+c_\nu)^2\) term dominates [2412.19195].

For sub-10 MeV masses, reactor environments offer an entirely different production regime. Starting from
\[
\mathcal{L}\supset -\,i\,g_{aee}\,a\,\bar e\,\gamma_5\,e,
\qquad
g_{aee}=\frac{m_e\,C_{ae}}{f_a},
\]
the relevant source process is Compton-like production \(\gamma e\to a e\), with a reactor-core photon flux
\[
\frac{d\Phi_\gamma}{dE_\gamma}
=
5.8\times10^{17}\;\bigl(P_{\rm MW}/{\rm MW}\bigr)\;
e^{-1.1\,E_\gamma/{\rm MeV}}
\quad[{\rm MeV^{-1}\,s^{-1}}].
\]
Detection then proceeds through inverse Compton-like scattering, axio-electric absorption, or \(a\to e^+e^-\) decays if \(m_a>2m_e\) [2010.15712].

## 5. Search strategies and mapped parameter space

The STCF study combines prompt vertexing and a proposed far hall. The prompt detector is assumed to reconstruct displaced vertices up to \(L_{\max}=1\) m from the interaction point, with ECAL photon detection threshold \(p_T(\gamma)>10\) MeV and vertex resolution of order \(100\,\mu{\rm m}\). The far detector is a cylindrical hall of length \(30\) m, radius \(20\) m, with front face \(10\) m from the interaction point; the decay probability in a distance interval \([l_1,l_2]\) is \(P=\exp(-l_1/\beta\gamma\tau)-\exp(-l_2/\beta\gamma\tau)\). For tauphilic \(a\to\gamma\gamma\), the prompt analysis requires exactly two photons with \(p_T>10\) MeV, a common displaced vertex in \(10^{-4}\,{\rm m}<\ell<1\,{\rm m}\), \(m_{\gamma\gamma}\) within \(\pm10\%\) of \(m_a\), and vetoes around \(\pi^0\), \(\eta\), \(\eta'\), and \(K^0_L\). At \(\sqrt s=3.56\) GeV the machine yields \(\mathcal{O}(10^8)\) \(\tau^+\tau^-\) pairs per year, while at \(\sqrt s=4.2\) GeV it yields \(\mathcal{O}(3.6\times10^9)\) pairs per year. In the tauphilic case the projected sensitivity is centered on an unprobed region around \(m_a\sim1000\) MeV and \(c^A_\tau\sim10^2-10^3\), and the right panel of the projection figure shows a \(15\)–\(75\%\) gain in the electroweak-violating reach relative to the electroweak-preserving reach. In the LFU case, by contrast, the prompt \(\sqrt s=4.2\) GeV trilepton search covers \(100\,{\rm MeV}\lesssim m_a\lesssim200\,{\rm MeV}\) with \(|c^A_\ell|\sim1\)–\(10\), while the far detector offers comparable reach to CHARM for \(m_a\sim1\)–\(10\) MeV and \(|c^A_\ell|\sim1\) [2509.23949].

Forward-detector studies at the LHC emphasize the LLP regime. FASER with \(300\,{\rm fb}^{-1}\) and FASER II with \(3000\,{\rm fb}^{-1}\) assume a background-free \(N\ge3\) signal criterion and require total visible energy \(>500\) GeV; heavy-meson spectra are simulated with Pythia 8, light mesons with EPOS-LHC, and detector acceptance with the foresee package. LHCb uses \(\eta\in[2,5]\), \(p_T(\mu)>0.5\) GeV, a displacement \(\ell_{xy}\in[6,22]\) mm, and a global efficiency of \(0.7\). In the electrophilic case, the electroweak-preserving scenario yields no sensitivity at FASER, while FASER II reaches \(c_e^A\sim10^2\) for \(m_a<5\) MeV. In the electroweak-violating scenario, FASER excludes \(c_e^A\sim10\)–\(10^3\) for \(m_a\lesssim0.1\) GeV, and FASER II reaches \(c_e^A\sim1\) for \(m_a\lesssim1\) GeV. In the muonphilic case above \(2m_\mu\), FASER loses sensitivity because the dimuon mode shortens the lifetime, but the LHCb Upgrade with \(300\,{\rm fb}^{-1}\) in the left-right softly-asymmetric scenario can probe down to \(c_\mu^A\sim10^{-2}\) for \(0.2\lesssim m_a\lesssim0.6\) GeV [2412.19195].

A broader laboratory compilation shows that loop-induced ALP–photon interactions are basis dependent and that experimental exclusions differ sharply between derivative and pseudoscalar formulations. In the derivative basis, LEP monophoton and exotic-\(Z\) bounds constrain \(C_\ell/\Lambda\lesssim10^{-3\text{–}2}\,{\rm GeV^{-1}}\) for \(m_a\lesssim M_Z\), while in the pseudoscalar basis the \(m_\ell^2/M_Z^2\) suppression makes these bounds negligible. E137 sets the strongest limits at small couplings, \(\lesssim10^{-5}\)–\(10^{-4}\,{\rm GeV^{-1}}\), for \(m_a\sim10\)–\(200\) MeV when the ALP decays visibly in the detector; NA64\(_e\) covers larger couplings \(\gtrsim10^{-3}\)–\(10^{-2}\,{\rm GeV^{-1}}\); and NA64\(_\mu\) is the leading probe of muon-only couplings. The schematic exclusion plots indicate that, in the derivative basis, the region \(10^{-6}\lesssim C_\ell/\Lambda\lesssim10^{-2}\,{\rm GeV^{-1}}\) and \(1\lesssim m_a\lesssim200\,{\rm MeV}\) is largely excluded for electron and muon couplings, leaving small windows only in the pseudoscalar basis or for tau-only couplings [2504.05873].

Long-lived ALPs from \(\tau\) decays probe a complementary sector of parameter space. CHARM and BEBC, using \(400\) GeV proton beams and \(\tau\) production from \(D_s\to\tau\nu\), exclude LFV leptophilic ALPs with \(f_a/g_{\ell\tau}\lesssim10^7\)–\(10^8\) GeV for \(m_a\sim0.2\)–\(0.9\) GeV. Projected \(90\%\) C.L. sensitivities extend further: the Fermilab Main Injector neutrino detectors reach \(f_a/g_{\ell\tau}\sim10^7\) GeV at \(m_a\sim0.2\) GeV, ProtoDUNE reaches a few \(\times10^7\) GeV, FASER reaches \(10^6\) GeV at \(m_a\sim0.4\) GeV, FASER-2 extends to \(10^8\) GeV, DUNE Near Detector probes a few \(\times10^7\) GeV, and SHiP reaches \(10^9\) GeV for \(m_a\sim0.2\)–\(0.8\) GeV. In the universal LFC scenario SHiP, FASER-2, ProtoDUNE, and DUNE also cover \(f_a\sim10^2\)–\(10^4\) GeV for \(m_a<2m_\mu\) [2507.15271].

For the MeV-scale electron-coupled regime, reactor searches fill a different gap. The quoted sensitivities show that current or near-future reactor setups with a \(10\) kg Ge detector, \(L\sim10\) m, and backgrounds \(10\)–\(100\) DRU can probe \(g_{aee}\sim10^{-5}\)–\(10^{-6}\) for \(m_a\sim0.5\)–\(1\) MeV. A ton-scale Xe detector with background \(1\)–\(10\) DRU extends to \(g_{aee}\sim10^{-8}\) at \(m_a\sim1\) MeV, entering regions unconstrained by Borexino or cooling arguments [2010.15712].

## 6. Precision anomalies and charged-lepton-flavor violation

The role of leptophilic ALPs in anomalous magnetic moments is highly benchmark dependent. In the universal charged-lepton model analyzed by Ganguly et al., the one-loop contribution to \(\Delta a_\mu\) is always negative because of the \(\gamma_5\) structure,
\[
\Delta a_\mu^{\rm (1\text{-}loop)}
=
-\,\frac{g_{a\mu\mu}^2}{8\pi^2}\,\frac{m_\mu^2}{m_a^2}
\int_0^1dx\;
\frac{x^3}{1-x+\frac{m_\mu^2}{m_a^2}x^2},
\]
whereas the Barr–Zee-type two-loop term is positive. After imposing collider and other laboratory bounds, the only surviving slice able to account for the Fermilab/BNL \(2\sigma\) band is
\[
5\,{\rm GeV}\lesssim m_a\lesssim6\,{\rm GeV},
\qquad
0.08\,{\rm GeV^{-1}}\lesssim c_{\ell\ell}/f_a\lesssim0.15\,{\rm GeV^{-1}}.
\]
The same study notes exclusion from BaBar \(e^+e^-\to\mu^+\mu^-a\to4\mu\) for \(c_{\ell\ell}/f_a\gtrsim0.15\,{\rm GeV^{-1}}\) and \(m_a<4\) GeV, while loop-induced \(a\gamma\gamma\) couplings are tested by mono-\(\gamma\) and tri-\(\gamma\) searches at LEP, CDF, and the LHC [2204.07920].

A different precision target is the anomalous magnetic moment of the electron in the \(m_a>f_a\) regime. The EFT analysis of that region argues that the often-assumed hierarchy \(m_a\ll f_a\) is not a fundamental requirement: the consistency condition is \(m_a<\Lambda\sim g_*f_a\), with the gray exclusion region in the \(m_a\)–\(f_a\) plane given by \(m_a>4\pi f_a\). In that framework a leptophilic ALP may explain the \(-3.8\sigma\) Caesium tension in \(\Delta a_e\) over an orange best-fit band that extends partly into \(m_a>f_a\), while present \(\mu\to e\gamma\) and \(\mu\to e\) conversion bounds imply roughly \(m_a f_a\gtrsim \mathcal{O}(10\,{\rm TeV})^2\) and \(m_a f_a\gtrsim \mathcal{O}(12\,{\rm TeV})^2\), with future Mu2e/COMET pushing to \(\mathcal{O}(25\,{\rm TeV})^2\) [2605.00115].

Charged-lepton-flavor violation provides an additional organizing principle. In anomaly-free flavored constructions, tree-level \(\ell_i\to\ell_j a\) obeys
\[
\Gamma(\ell_i\to\ell_j a)
=
\frac{m_i^3}{16\pi}\,\frac{|C^e_{ij}|^2}{4f_a^2}
\left(1-\frac{m_a^2}{m_i^2}\right)^2,
\]
while one-loop ALP exchange induces \(\ell_i\to\ell_j\gamma\). The quoted phenomenology includes present and projected limits such as \({\rm BR}(\mu\to e a)<2.6\times10^{-6}\), a MEG II forward projection \({\rm BR}(\mu\to e a)<1.3\times10^{-7}\), and \({\rm BR}(\tau\to\mu a)<4.5\times10^{-3}\). In Model II, a keV ALP can fit the XENON1T e-recoil excess while satisfying stellar-cooling and LFV limits through suitable flavor textures. This suggests that the status of leptophilic ALPs in precision observables is controlled as much by flavor alignment as by mass and coupling size [2007.08834].

## 7. Model-building realizations and theoretical issues

Several constructions realize leptophilic ALPs as pseudo–Nambu–Goldstone bosons of broken lepton-sector symmetries. In the anomaly-free flavored framework, the ALP can originate from a flavon field \(\phi\) with
\[
\phi(x)=\frac{v_\phi+s(x)}{\sqrt2}\,e^{i\,a(x)/v_\phi},
\]
so that \(f_a\simeq v_\phi\). Two representative realizations are given. Model I ties the ALP to the flavon that generates charged-lepton masses through Froggatt–Nielsen-like powers \(\epsilon\equiv v_\phi/\Lambda\simeq0.1\), yielding CKM-like charged-lepton mixing. Model II assumes a larger flavor symmetry, allowing \(U_L^e,U_R^e\sim U_{\rm PMNS}\) with large \(\mathcal{O}(1)\) off-diagonals and, for \(x^L=-x^R\), purely axial couplings with \(C^V=0\) [2007.08834].

An inverse-seesaw realization identifies the ALP with the phase of a singlet scalar \(\sigma\),
\[
\sigma(x)=\frac{v_\sigma+\rho(x)}{\sqrt2}\,e^{\,i\,a(x)/v_\sigma},
\]
where \(\langle\sigma\rangle=v_\sigma/\sqrt2\sim10^{9\div14}\,{\rm GeV}\). The same vacuum expectation value generates inverse-seesaw mass scales through Planck-suppressed operators and the low-energy ALP interaction is derivative,
\[
\mathscr{L}_{\rm eff}
=
\tfrac12(\partial_\mu a)^2-\tfrac12m_a^2a^2-\frac{g_{a\gamma}}{4}aF_{\mu\nu}\tilde F^{\mu\nu}
+\frac{\partial_\mu a}{v_\sigma}J_X^\mu.
\]
In this class of models,
\[
g_{a\ell\ell}=\frac{X_\ell m_\ell}{f_a},
\qquad
g_{a\nu_i\nu_i}\simeq \frac{X_L m_{\nu_i}}{f_a},
\]
and a discrete gauge symmetry such as \(Z_{13}\) is used to forbid dangerous low-dimension gravity-induced operators that would otherwise spoil the ALP mass or inverse-seesaw scales [1503.03502].

A persistent conceptual issue is the status of the ALP–photon interaction. Some EFT formulations generate an \(aF\tilde F\) term under field redefinition, while anomaly-free flavored models impose no electromagnetic anomaly at tree level. This suggests that the phenomenological weight of photon-mediated searches is basis- and model-dependent rather than universal [2504.05873, 2007.08834].

Taken together, the recent literature places leptophilic ALPs at the intersection of LLP searches, flavor physics, and precision tests. The parameter space is not organized by a single canonical benchmark: tauphilic, LFU, electrophilic, muonphilic, flavor-violating, anomaly-free, and \(m_a>f_a\) realizations all emphasize different observables. The resulting picture is one of a broad but structured search program in which exotic \(\tau\) decays, far detectors, beam dumps, reactor experiments, and charged-lepton precision measurements probe complementary limits of the same lepton-centered EFT.

Source: https://www.emergentmind.com/topics/leptophilic-axion-like-particles-alps