Resonant W and Z Boson Production in FSRQ Jets: Implications for Diffuse Neutrino Fluxes

Blazars, particularly Flat Spectrum Radio Quasars (FSRQs), are well-known for their ability to accelerate a substantial population of electrons and positrons, as inferred from multiwavelength radiation observations. Therefore, these astrophysical objects are promising candidates for studying high-energy electron–positron interactions, such as the production of $W^{\pm}$ and $Z$ bosons. In this work, we explore the implications of electron–positron annihilation processes in the jet environments of FSRQs, focusing on the resonant production of electroweak bosons and their potential contribution to the diffuse neutrino flux. By modeling the electron distribution in the jet of the FSRQ 3C~279 during a flaring state, we calculate the reaction rates for $W^{\pm}$ and $Z$ bosons and estimate the resulting diffuse fluxes from the cosmological population of FSRQs. We incorporate the FSRQ luminosity function and its redshift evolution to account for the population distribution across cosmic time, finding that the differential flux contribution exhibits a pronounced peak at redshift $z \sim 1$. While the expected fluxes remain well below the detection thresholds of current neutrino observatories such as IceCube, KM3NeT, or Baikal-GVD, the flux from $Z$ boson production within the jet blob is many orders of magnitude smaller than the total diffuse astrophysical neutrino flux. These results provide a theoretical benchmark for the role of Standard Model electroweak processes in extreme astrophysical environments, highlighting the interplay between particle physics and astrophysics, and illustrating that even extremely rare high-energy interactions can leave a subtle, theoretically meaningful imprint on the diffuse astrophysical neutrino background.

💡 Research Summary

**
The paper investigates the possibility that high‑energy electron–positron pairs in the jets of flat‑spectrum radio quasars (FSRQs) can resonantly produce electroweak gauge bosons, namely W± and Z, and that the subsequent decays of these bosons could contribute to the diffuse astrophysical neutrino background observed on Earth. The authors focus on the well‑studied FSRQ 3C 279 during a bright flare on 20 December 2013, using its multi‑wavelength spectral energy distribution (SED) to constrain the underlying electron population.

A one‑zone leptonic model is adopted, in which a compact, homogeneous “blob” moves relativistically with a bulk Lorentz factor Γ≈δ_D≈50–70. Within the blob, electrons are accelerated by a combination of diffusive shock acceleration and stochastic (turbulent) acceleration, while losing energy through synchrotron radiation, synchrotron‑self‑Compton (SSC), and external Compton (EC) scattering on photons from the dusty torus and broad‑line region (BLR). The steady‑state electron distribution N_e(γ) is obtained by solving a Fokker‑Planck equation that includes energy diffusion (coefficient D₀), systematic gains/losses (parameter a), and escape (characterized by τ). Two representative parameter sets (Model 1 and Model 2) are taken from previous fits to the 3C 279 flare, yielding electron Lorentz factors up to γ_max≈10⁴–10⁵ and a total electron power of order 10⁴⁵ erg s⁻¹.

With this electron distribution, the authors compute the rate of resonant e⁺e⁻ → W± and e⁺e⁻ → Z reactions. The cross‑section is modeled by a Breit‑Wigner formula centered at the boson masses (M_W≈80.3 GeV, M_Z≈91.2 GeV) with their Standard Model widths (Γ_W≈2.1 GeV, Γ_Z≈2.5 GeV). Because the jet plasma is highly Doppler‑boosted, the center‑of‑mass energies of the colliding pairs are shifted, allowing a fraction of the pairs to sit within the narrow resonance windows. The resulting reaction rates scale as n_{e⁺e⁻}²⟨σv⟩, where n_{e⁺e⁻} is the pair density inside the blob (∼10⁻⁴–10⁻³ cm⁻³). The Z channel enjoys a larger peak cross‑section (∼10⁻³⁴ cm²) than the W channel, but only about 20 % of Z decays produce neutrinos (Z→νν̄), whereas each W decays to a charged lepton and a neutrino with ∼10 % branching per lepton flavor. Consequently, the total neutrino yield from both channels is comparable, with Z contributing slightly more.

To assess the cosmological impact, the single‑source neutrino emissivity is folded with an empirical FSRQ luminosity function Φ(L,z) and its redshift evolution. The integration over luminosity and redshift yields a diffuse neutrino flux that peaks at redshift z≈1, reflecting the epoch of maximal FSRQ activity. However, the predicted flux is many orders of magnitude below the measured astrophysical neutrino flux by IceCube. Quantitatively, the Z‑induced neutrino flux is ≲10⁻⁶ of the total diffuse flux, and the W‑induced component is even smaller. The authors explicitly compare these numbers with the current sensitivities of IceCube, KM3NeT, and Baikal‑GVD, concluding that detection is infeasible with present instruments.

The paper discusses several sources of uncertainty. First, the assumption of an isotropic, homogeneous electron–positron distribution may be oversimplified; real jets likely exhibit anisotropies and temporal variability that could locally enhance the resonant pair density. Second, the use of a hard‑sphere diffusion coefficient (D₀∝γ²) is justified by turbulent simulations, but alternative turbulence spectra (e.g., Kolmogorov) change the steady‑state electron spectrum only modestly (factor ≲2 in the neutrino flux). Third, the Standard Model boson widths are taken as vacuum values; in the dense plasma of a jet, collective effects could modify the effective width or branching ratios, though such corrections are expected to be small. Fourth, the luminosity function and its evolution are subject to observational biases, which could shift the redshift of the flux peak or alter the overall normalization. Finally, the SED fitting itself carries systematic uncertainties (magnetic field strength, external photon field energy densities) that propagate into the inferred pair density.

Despite these caveats, the work provides a valuable benchmark for the role of rare electroweak processes in extreme astrophysical environments. It demonstrates that even when the underlying particle physics is well understood, the astrophysical conditions (pair density, Doppler boosting, redshift distribution) render the contribution to the observable neutrino sky negligible. The study opens a pathway for future investigations that could incorporate multi‑zone jet structures, detailed particle‑in‑cell simulations of e⁺e⁻ collisions in magnetized plasma, and the next generation of neutrino telescopes (e.g., IceCube‑Gen2, KM3NeT‑ARCA) whose improved sensitivity might eventually probe such faint signals. In summary, resonant W and Z production in FSRQ jets is a theoretically sound but practically undetectable source of high‑energy neutrinos, serving as a stringent test of our understanding of both jet physics and electroweak interactions in the cosmos.