---
title: Muon PDF Probes in Collider Physics
url: https://www.emergentmind.com/topics/muon-pdf-probes
type: topic
---

# Muon PDF Probes in Collider Physics

Muon PDF probes are analyses in which the muon is treated as a source of perturbatively calculable partons, or, in a distinct analogy, as a parent particle whose decay spectra are parameterized and constrained through secondary scattering observables. In the collider usage, the relevant objects are electroweak parton distribution functions \(f_{i/\mu}(x,Q)\) for muons, neutrinos, gauge bosons, quarks, gluons, and interference states such as the mixed \(Z/\gamma\) PDF; in the decay-spectrum usage, the Michel parameters \(P_x\) and \(w_x\) compactly encode how the muon’s energy is distributed among decay products and can be probed through coherent elastic neutrino–nucleus scattering (CE\(\nu\)NS) [2410.21383, 2408.13191, 2602.16771, 2502.18175].

## 1. Partonic description of the high-energy muon

At energies far above the muon mass, collinear radiation of gauge bosons and fermion pairs from the incoming muon becomes enhanced by logarithms of the ratio of the hard scale to the particle masses. Emissions such as
\[
\mu \to \mu \gamma,\quad \mu\to\mu Z,\quad \mu\to \nu_\mu W,\quad \gamma\to f\bar f
\]
justify a PDF description of the muon at multi-TeV colliders. The inclusive cross section for a final state \(F\) is written as
\[
\sigma_{\mu^+\mu^-\to F} = \sum_{i,j}\int_0^1 dx_1\,dx_2\; f_{i/\mu^+}(x_1,Q)\, f_{j/\mu^-}(x_2,Q)\, \hat\sigma_{ij\to F}(\hat s, Q,\mu_R),
\]
with \(\hat s=x_1x_2s\). In this framework, the “initial state” is no longer a pure on-shell muon, but a muon plus a cloud of collinear electroweak and QCD partons [2602.16771].

The evolution of these PDFs is governed by generalized DGLAP equations. In the notation used for muon-collider PDF evolution,
\[
Q^2\frac{d f_B(x,Q)}{d Q^2} = P_B^v(Q)\, f_B(x,Q) + \sum_{A,C} \frac{\alpha_{ABC}}{2\pi}\, \widetilde{P}_{BA}^C \otimes f_A + \frac{v^2}{16\pi^2Q^2} \sum_{A,C} \widetilde{U}^C_{BA}\otimes f_A.
\]
The perturbative boundary condition for a muon beam is
\[
f_\mu(x,Q=m_\mu) = \delta(1-x),\qquad f_{i\neq\mu}(x,Q<m_i)=0.
\]
This is a central distinction from proton PDFs: the muon beam has no nonperturbative input comparable to hadronic structure functions [2602.16771].

The partonic content generated by the evolution includes the valence muon, radiatively induced charged leptons and neutrinos, the electroweak gauge bosons \(\gamma\), \(W^\pm\), and \(Z\), quarks and gluons generated via \(\gamma/Z \to q\bar q\), and interference PDFs. In the Standard Model there are two mixed channels singled out by quantum numbers: transverse photon–\(Z\) interference and longitudinal \(Z\)–Higgs interference. In the broken phase, the transverse \((\gamma,Z)\) sector must therefore be treated as a matrix of PDFs,
\[
\rho^{\rm split}_{V_1 V_2} =
\begin{pmatrix}
f_\gamma & f_{Z/\gamma}\\
f_{Z/\gamma} & f_Z
\end{pmatrix},
\qquad V_1,V_2\in\{\gamma,Z\},
\]
rather than as independent scalar distributions [2408.13191].

## 2. Neutrino and gauge-boson content as collider probes

The most distinctive electroweak constituent of the muon is the muon-neutrino PDF generated by the charged-current splitting
\[
\mu_L^- \to \nu_\mu + W^-.
\]
This induces a distribution \(f_{\nu_\mu/\mu}(x,Q^2)\) that is strongly biased toward large \(x\sim 1\), because the emission is dominantly soft in the gauge boson and the neutrino tends to take most of the muon’s energy. At both 3 TeV and 10 TeV muon colliders, the \(\mu^-\mu^+\) and \(\mu_L^-\bar\nu_\mu\) luminosities are the largest at large \(\sqrt{\hat s}\), while gauge-boson–gauge-boson luminosities are subdominant in that regime [2410.21383].

This structure makes several Standard Model processes into direct probes of the neutrino content of the muon. In
\[
\mu^+\mu^- \to e^- \bar\nu_e + X,
\]
the dominant signal channel is
\[
\mu^-\bar\nu_\mu\to e^- \bar\nu_e,
\]
while the main background is vector boson fusion through \(W^-\gamma\to e^- \bar\nu_e\) and \(W^-Z\to e^- \bar\nu_e\). For
\[
\mu^+\mu^- \to W^- \gamma + X,
\]
the corresponding neutrino-initiated subprocess is
\[
\mu^- \bar\nu_\mu \to W^- \gamma.
\]
In both classes of observables, high-\(p_T\) regions become increasingly dominated by the \(\mu\bar\nu_\mu\) luminosity, so differential distributions in \(p_T\), rapidity, and \(\hat s\) provide an experimental handle on \(f_{\nu_\mu/\mu}\) itself [2410.21383].

A related use of “muon PDF probes” appears in multi-TeV vector boson fusion. At a muon collider, collinear radiation of \(W\) and \(Z\) bosons off energetic muons makes the beams act like sources of effective electroweak gauge bosons, so that heavy-state production can be organized schematically as
\[
\sigma(\mu^+\mu^- \to \mu^+\mu^- h_2) \sim \int f_{V/\mu}(x_1)\, f_{V/\mu}(x_2)\, \hat\sigma(VV\to h_2)\,dx_1dx_2.
\]
Within a model-independent heavy-scalar framework with couplings \(c_2\) and \(c_{12}\), observing
\[
V V \to h_2 \to Z h_1
\]
implies \(c_2\neq 0\) and \(c_{12}\neq 0\), so that the basis-independent invariant
\[
K_{12}=c_1 c_2 c_{12}\simeq c_2 c_{12}
\]
is nonzero. Under the stated assumptions, this is sufficient to establish CP violation in the scalar sector. In this sense, the effective “PDF of the muon into \(W/Z\)” becomes the production-side lever arm for a bosonic CP-violation test [2511.21874].

## 3. The mixed \(Z/\gamma\) PDF and interference structure

The mixed \(Z/\gamma\) PDF is an interference object rather than a probability density in the ordinary sense. It arises because, after integrating over the azimuth of the collinear emission, interference between different helicity states vanishes, but interference between distinct species survives when they share the same unbroken quantum numbers and can enter the same splitting and hard amplitude. In the Standard Model, this selects exactly the transverse \(\gamma\)–\(Z\) sector and the longitudinal \(Z\)–Higgs sector [2408.13191].

At leading order, the effective-vector-boson approximation for \(f_{Z/\gamma}\) is accidentally suppressed by the small vector coupling combination
\[
Q_{\mu_L}^Z+Q_{\mu_R}^Z = -\frac{1}{2}+2 s_W^2 \approx -0.038.
\]
This suppression is not present in the leading-logarithm resummed result. Once electroweak evolution polarizes the muon PDFs and includes higher-order splittings, the mixed \(Z/\gamma\) PDF becomes comparable in size to other electroweak gauge-boson PDFs. Extending the analytic approximation to \(\mathcal O(\alpha^2)\) lifts the accidental cancellation through double gauge-boson emission and Sudakov-enhanced terms of the form \(\alpha^2\log^3(Q^2/m_Z^2)\), bringing the analytic treatment into good agreement with the full LL-resummed numerical evolution [2408.13191].

This interference PDF has direct phenomenological consequences. In high-energy Compton scattering,
\[
\mu^- (Z,\gamma)\to \mu^- \gamma,
\]
the ratio
\[
\Delta_{Z/\gamma}\equiv \frac{\sigma_{Z/\gamma}}{\sigma_{\rm tot}}
\]
is negative and ranges from \(\sim 1\%\) up to tens of percent at large transverse momentum and backward muon rapidity. For a 10 TeV muon collider with \(10\,\mathrm{ab}^{-1}\), the effect exceeds \(3\sigma\) in several \((p_T,y_\mu)\) bins and exceeds \(5\sigma\) in the bin \(p_T\in[2,5]\) TeV and \(y_\mu\in[-1,0]\). The corresponding \(\nu_\mu W^-\) background is small in the most sensitive region [2408.13191].

The same interference structure affects Higgs observables. In associated \(WH\) production, the inclusion of the \(Z/\gamma\) PDF shifts the 10 TeV Standard Model prediction from \(6.81\) fb to \(6.63\) fb in the \(\kappa_W=\kappa_Z=1\) limit, while at 3 TeV the change is at the \(\sim 2\%\) level. This is already comparable to the \(\sim 1\%\) precision anticipated for a high-luminosity \(WH\) measurement, so the mixed PDF must be included consistently in any percent-level extraction of \(\kappa_W\) and \(\kappa_Z\) [2408.13191].

## 4. Beyond-the-Standard-Model distortions of muon PDFs

Muon PDF probes are also a BSM search strategy. In an \(L_\mu-L_\tau\) model with a new gauge boson \(Z'\), the Lagrangian
\[
\mathcal{L} \supset -\frac{1}{4}F'^{\alpha\beta}F'_{\alpha\beta}+\frac{1}{2}M_{Z'}^2 Z'^\alpha Z'_\alpha
-g' Z'_{\alpha} (\ell_\mu^\dagger\bar{\sigma}^\alpha \ell_\mu + \mu^{c} \sigma^\alpha \mu^{c\,\dagger} - \ell_\tau^\dagger\bar{\sigma}^\alpha \ell_\tau - \tau^{c} \sigma^\alpha \tau^{c\,\dagger})
\]
adds new parton species and mixed states to the DGLAP system. Below the electroweak scale, for \(Q > M_{Z'}\), the evolution includes \(Z'_T\), \(Z'_L\), \(\gamma_T/Z'_T\), \(\nu_\mu\), \(\nu_\tau\), and \(\mu^+\); above the electroweak scale it tracks \(Z'_\pm\), \(Z'_L\), \(\gamma_\pm/Z'_\pm\), \(Z_\pm/Z'_\pm\), \(Z_L/Z'_L\), \(\bar{\nu}_\mu\), and \(\mu^+_{L/R}\). Numerically, the main effects are an enhancement of the muon PDF at \(x<1\) and a suppression of the photon PDF, both at the percent level for the benchmark \(M_{Z'}=50\,\mathrm{GeV}\), \(g'=0.02\) [2602.16771].

The associated collider observable is the invariant-mass fraction
\[
\tau \equiv \frac{\hat s}{s},
\]
through the distributions \(d\sigma_F/d\tau\) in the final states \(F=\mu^+\mu^-, \gamma\gamma, \mu^\pm\gamma, e^+ e^-\). A reference setup uses \(\sqrt{s}=10\) TeV, \(\mathcal L_0=10\,\mathrm{ab}^{-1}\), central cuts \(|\eta|<2.0\) or \(2.5\), and \(m_F\ge 2\,\mathrm{TeV}\). Within that framework, the profile-likelihood analysis is almost insensitive to the luminosity uncertainty because the sensitivity is dominated by the shape of \(d\sigma_F/d\tau\), not by overall normalization. In the mass window \(M_{Z'}\sim 50\)–\(100\) GeV, the PDF-based reach lies significantly below the direct-search curves from \(\tau^+\tau^-\gamma\) and \(\nu\bar\nu\gamma\), so it probes smaller \(g'\) couplings than those associated-production channels at the same collider [2602.16771].

The neutrino PDF also modifies new-physics production channels directly. For a heavy scalar \({\rm SU}(2)_L\) doublet \(S=(S_+,S_-)\), the charged-current mode
\[
\nu_\mu \mu^+ \to W^{+*} \to S_+ S_-^\dagger
\]
is present only because of the muon’s \(\nu_\mu\) content. The inclusive LePDF-based prediction agrees with the fixed-order simulation of
\[
\mu^- \mu^+ \to S_+ S_-^\dagger W^-
\]
within scale uncertainties, and the kinematic distributions show the expected ISR pattern: the radiated \(W\) is mostly forward and low-\(p_T\), while the heavy scalars are central and high-\(p_T\) [2410.21383].

A further example is resonant single production of an axion-like particle \(\phi\) through VBF. In a 10 TeV collider, the \(Z/\gamma\)-induced interference contribution typically modifies the total cross section by \(\sim 10\%\) for \(M_\phi\lesssim 4\) TeV, with the sign and magnitude depending on the combinations \(c_{\phi\gamma\gamma}\), \(c_{\phi ZZ}\), and \(c_{\phi\gamma Z}\). This makes the mixed-PDF sector relevant not only for Standard Model precision studies but also for coupling extraction in BSM resonance searches [2408.13191].

## 5. Muon-decay spectra as “muon PDFs” in the Michel-parameter sense

A distinct use of the phrase “Muon PDF probes” appears in stopped-muon decay. In this context, the energy distributions of \(e^+\), \(\nu_e\), and especially \(\bar\nu_\mu\) are treated as objects analogous to PDFs: they encode how the muon’s energy is distributed among decay products and are parameterized by neutrino Michel parameters \(P_x\) and \(w_x\). These parameters are quadratic functions of the effective muon-decay couplings \(g^X_{\epsilon\eta}\), so deviations from the Standard Model values \(P_x=1\) and \(w_x=0\) diagnose non-\(V\!-\!A\) Lorentz structures, right-handed currents, or right-handed neutrinos in muon decay [2502.18175].

At the Spallation Neutron Source, the delayed \(\nu_e\) and \(\bar\nu_\mu\) fluxes from
\[
\mu^+ \to e^+ + \nu_e + \bar\nu_\mu
\]
feed the CE\(\nu\)NS signal measured by COHERENT. After convolution with the CE\(\nu\)NS cross section, the delayed event rate can be written in terms of effective rescaling factors \(x_{\nu_e}\) and \(x_{\bar\nu_\mu}\), with
\[
\begin{aligned}
x_{\bar\nu_\mu} &= P_{\bar\nu_L} - \frac{4}{3}\,P_{\bar\nu_L} w_{\bar\nu_L} + \frac{4}{3}\,P_{\nu_L} w_{\nu_L},\\
x_{\nu_e} &= P_{\nu_L} - \frac{4}{3}\,P_{\nu_L} w_{\nu_L}  + \frac{4}{3}\,P_{\bar\nu_L} w_{\bar\nu_L}.
\end{aligned}
\]
This makes the recoil-energy and timing distributions of CE\(\nu\)NS into probes of the muon-decay spectrum itself [2502.18175].

Using CsI and LAr data from COHERENT, the experimentally best-constrained linear combination is
\[
0.54\, P_{\nu_L} + 0.84\, P_{\bar\nu_L} + 0.40\, \left(P_{\nu_L}w_{\nu_L} - P_{\bar\nu_L}w_{\bar\nu_L}\right)  = 1.25\pm 0.21.
\]
This is the first direct experimental constraint involving the muon antineutrino Michel parameters \(P_{\bar\nu_L}\) and \(w_{\bar\nu_L}\), and only the second constraint on the electron-neutrino Michel parameters \(P_{\nu_L}\) and \(w_{\nu_L}\). The quoted one-at-a-time \(1\sigma\) ranges are
\[
\begin{aligned}
P_{\nu_L}& = 0.76\left(^{+0.24}_{-0.39}\right), &\qquad w_{\nu_L} &= 0.00\left(^{+0.29}_{-0.00}\right),\\
P_{\bar\nu_L}& = 0.84\left(^{+0.16}_{-0.25}\right), &\qquad w_{\bar\nu_L} &= 0.32 \left(^{+0.43}_{-0.53}\right).
\end{aligned}
\]
Within this usage, CE\(\nu\)NS acts as a precision production-side measurement of the muon’s decay “PDFs,” rather than of electroweak initial-state radiation [2502.18175].

## 6. Experimental foundations, uncertainties, and limitations

Muon PDF probes depend on both beam preparation and theoretical control. A compact front-end muon source optimized for a Neutrino Factory or muon-collider program begins with an 8 GeV proton beam on a liquid mercury jet inside a 20 T solenoid, followed by a 20 T \(\to\) 2 T taper over 6.00 m, a 2 T decay channel, a bent-solenoid chicane, RF bunching and phase rotation, a matching section, and an ionization-cooling channel. The quoted performance is about \(0.12\,\mu^+\) and \(0.12\,\mu^-\) per incident 8 GeV proton within the reference acceptance, a train of \(\sim 21\) well-defined bunches centered around \(p\simeq 245\,\mathrm{MeV}/c\), and a transverse normalized rms emittance reduced by a factor of \(\sim 3\) to \(\varepsilon_{\perp,N}\approx 6\,\mathrm{mm\cdot rad}\). The chicane and absorber reduce downstream energy deposition by more than an order of magnitude compared to a scheme without them. This source is explicitly described as directly relevant to using muons as high-energy probes once further acceleration is provided [1504.00380].

Theoretical systematics remain substantial. For the \(\nu_\mu\) PDF, factorization-scale variations \(Q\to Q/2\) and \(Q\to 2Q\) give relative changes of \(\sim 38\%\) around \(Q=500\) GeV and \(\sim 19\%\) around \(Q=3\) TeV, whereas QED-like pieces are at the \(\sim 8\%\) and \(\sim 7\%\) level, respectively. For mixed \(Z/\gamma\) effects, the scale spread is at the few-percent level in benchmark Compton and \(WH\) studies. In the \(L_\mu-L_\tau\) case study, varying the hard scale by a factor of two changes \(d\sigma_F/d\tau\) at the \(\sim 10\%\) level, larger than the percent-level deviations induced by the benchmark \(Z'\) couplings. These numbers indicate that higher-order electroweak evolution, matched fixed-order matrix elements, and improved treatment of massive and mixed PDFs are required for precision applications [2410.21383, 2408.13191, 2602.16771].

Experimental limitations are equally specific. Muon-collider detector concepts include forward shielding to suppress beam-induced background, so only objects with
\[
10^\circ < \theta < 170^\circ
\]
are reconstructed in some analyses. Many vector-boson-fusion muons and neutrinos therefore escape detection and appear as missing energy. Sensitivity estimates in several studies are purely statistical, with systematics such as background normalization, b-tagging efficiencies, and detector calibrations not yet included. This suggests that muon PDF probes are already structurally well defined, but that their ultimate reach will depend on coordinated advances in beam delivery, electroweak PDF theory, event-generator implementations, and detector realism [2511.21874, 2410.21383, 2602.16771].

Source: https://www.emergentmind.com/topics/muon-pdf-probes