Ao Zhang, Wenlei Chen, Manel Errando

The origin of intergalactic magnetic fields (IGMFs) remains one of the key open questions in cosmology. Gamma-ray pair halos produced by electromagnetic cascades from TeV-emitting blazars provide a powerful indirect probe of these fields. In this work, we present a novel search for pair halos that explicitly exploits their expected anisotropic morphology, aligning with the projected orientation of blazar jets on the sky. Using a Monte Carlo framework to model the spatial distribution of cascade emission, we identify an optimal sample of 21 high-synchrotron-peaked BL Lac objects with well-constrained jet position angles from radio interferometry. By rotating and stacking \textit{Fermi}-LAT observations of these sources along their jet directions, we enhance sensitivity to anisotropic extended emission that would be diluted in traditional orientation-agnostic analyses. Applying a likelihood analysis to the combined dataset, we find evidence for a non-zero IGMF, excluding the null hypothesis at $3.8σ$ level and obtaining a best-fit field strength of $B_0 = 2.8 \times 10^{-16}\,\mathrm{G}$, with a $99\%$ confidence interval of $0.9 \times 10^{-16}\,\mathrm{G} < B_0 < 8.9 \times 10^{-16}\,\mathrm{G}$. Our result is consistent with previous constraints from spectral, spatial, and temporal studies, while demonstrating that incorporating anisotropic information provides a significant gain in sensitivity. This approach opens a new avenue for probing intergalactic magnetism and highlights the potential of future high-angular-resolution gamma-ray observations to directly image pair halos and map magnetic fields in cosmic voids.