Abstract

Current Cherenkov telescopes have identified a population of ultra-high-frequency peaked BL Lac objects (UHBLs), also known as extreme blazars, that exhibit exceptionally hard TeV spectra, including 1ES 0229+200, 1ES 0347-121, RGB J0710+591, 1ES 1101-232, and 1ES 1218+304. Although one-zone synchrotron-self-Compton (SSC) models have been generally successful in interpreting the high-energy emission observed in other BL Lac objects, they are problematic for UHBLs, necessitating very large Doppler factors and/or extremely high minimum Lorentz factors of the emitting leptonic population. In this context, we have investigated alternative scenarios where hadronic emission processes are important, using a newly developed (lepto-)hadronic numerical code to systematically explore the physical parameters of the emission region that reproduces the observed spectra while avoiding the extreme values encountered in pure SSC models. Assuming a fixed Doppler factor $\delta=30$, two principal parameter regimes are identified, where the high-energy emission is due to: 1) proton-synchrotron radiation, with magnetic fields $B \sim 1-100\ \textrm{G}$ and maximum proton energies $E_{p;max} \lesssim 10^{19}$ eV; and 2) synchrotron emission from p-$\gamma$-induced cascades as well as SSC emission from primary leptons, with $B \sim 0.1-1\ \textrm{G}$ and $E_{p;max} \lesssim 10^{17}$ eV. This can be realized with plausible, sub-Eddington values for the total (kinetic plus magnetic) power of the emitting plasma, in contrast to hadronic interpretations for other blazar classes that often warrant highly super-Eddington values.