Testing a Lepto-Hadronic Two-Zone Model with Extreme High-Synchrotron Peaked BL Lacs and Track-like High-Energy Neutrinos

: Numerous studies suggest that high-energy (HE) neutrinos and ultra-high-energy (UHE) cosmic rays could originate from extremely high-synchrotron peaked (EHSP) BL Lacs, which have been identiﬁed as effective particle accelerators. Due to the discovery of HE-neutrinos by the IceCube telescope, these hypotheses may shortly have the opportunity to be tested. In this work, we use a two-zone leptohadronic model to explain the spatial coincidence of three EHSP BL Lac: 1RXS J09462.5+010459, 1ES 1101-232, and 3HSP J095507.9+355101 with the arrival of track-like neutrinos. Our results for 1RXS J09462.5+010459 and 1ES 1101-232 indicate that the model accurately describes the electromagnetic emission and neutrino events without increasing the ﬂuxes in the measured bands. In addition, the X-ray ﬂaring state of 3HSP J095507.9+355101 can be explained by our model, but the measured ultraviolet ﬂux during the neutrino arrival time window cannot be explained. For all cases, the broadband emission and neutrino arrival are better described by hard proton distributions ≈ 1.5. Finally, the proton luminosity required to explain the neutrino ﬂuxes is slightly higher than the Eddington limit with a photopion efﬁciency of ≈ 0.1 for non-ﬂaring state cases. On the other hand, for the ﬂaring state of 3HSP J095507.9+355101, the proton luminosity must be higher than the Eddington limit at least by one order of magnitude, even if the photopion efﬁciency reaches unity.


Introduction
A subset of Active galactic nuclei (AGN) known as blazars have relativistic jets that point to our direction of sight [1].They are divided into the Flat-spectrum radio quasars (FSRQ) and BL Lac objects sub-classes [2][3][4].These sources include those for which neutrinos have been observed ( [5,6]).When one of the two jets is aligned with our line of sight, its emission is strongly beamed and Doppler boosted, making it visible from radio up to very-high-energy (VHE) gamma-ray wavelengths [7,8] and exhibiting a wide range of variability.In both their quiescent and flaring phases, these objects display striking differences in distinct bands across the electromagnetic spectrum.
Two peaks can be seen in the spectral energy distribution (SED) of blazars, each one at a different frequency [9][10][11].The first peak is typically well-fit assuming synchrotron emission, and is used to classify blazars as low-synchrotron peaked (LSP), intermediatesynchrotron peaked (ISP), or high-synchrotron peaked (HSP).Using these criteria, FSRQ objects are categorized as LSP blazars since the broadband SED shows a peak at infrared wavelengths.While BL Lac sources can be classified as LSP, ISP, or HSP [12].The synchrotron peak is detected in LSP blazars at frequencies of log (ν peak ) < 14; in ISP blazars, this occurs at frequencies between 14 < log (ν peak ) < 15; and in HSP blazars, this occurs at frequencies of log (ν peak ) > 15.
Furthermore, BL Lacs with a synchrotron peak at log (ν peak ) > 17, are often referred to as extreme high-synchrotron peaked (EHSP) [13].Moreover, Refs.[14,15] have identified a novel category of BL Lacs, depending on the location of the second peak with a HE peak at E peak > 1 TeV; these objects are known as hard-TeV BL Lacs (TBL) or extreme-TeV BL Lacs.There is evidence to suggest that extreme behaviors are not permanent states and that they do not always coexist.For instance, such characteristics have been observed in Markarian 501, 1ES 1727+502, and 1ES 1741+196 [16][17][18][19].Extreme BL Lacs present a significant challenge to the presently understood acceleration and emission processes in blazar jets [20].The mechanism of VHE emission may be better understood if it is considered in conjunction with HE neutrinos.
To date, only two HE-neutrino sources have been found (i.e, TXS 0506 +056 and NGC 1068; [5,6,31]), however several astronomical objects were suggested as possible sources (for a review see [32,33]).The particular mechanism by which the HE neutrinos are created in such environments is not yet known, however several possibilities have been investigated (e.g., [34][35][36][37][38][39][40][41]).According to some research, BL Lacs are major contributors to the HE-neutrinos observed by IceCube.For instance, Ref. [42] found that BL Lacs can account for roughly 10% of IceCube's neutrino background at energies above 0.5 PeV when using the blazar simplified view approach.At lower energies, they estimated that BL Lacs may generate a level of ≈20 percent.Alternatively, more recently, Buson's research shows that blazars are definitely linked to HE astrophysical neutrinos, with an unprecedented chance probability [43].Furthermore, the results reported by [44] pointed out the evidence for IceCube neutrinos connection with HSP and ISP blazars.However, Ref. [45] point out that Buson's result may potentially be attributed to a statistical fluctuation, as well as the presence of spatial and flux non-uniformities within the blazar sample.Additionally, the utilization of the updated catalog yields alternative outcomes that also indicate reduced statistical significance [46].
In addition, Ref. [47] discovered a correlation between the neutrino catalog and EHSP accounting for roughly (1-20)% of the neutrino signal up to that time, and the most plausible neutrino candidates were also found to be ≈6 EHSP.It is interesting to note that [48] had already suggested H2356-309 as a candidate for the IC10 event's emitter source.This finding led [49] to propose a leptohadronic model with realistic parameter values to characterize the SED of the neutrino candidates, which they found to be consistent with the IceCube data for some of the sources.
Moreover, an EHSP was recently discovered coincident with the identification of the HE-neutrino IceCube 200107A; however, there is insufficient evidence to conclusively link the flaring activity of 3HSPJ095507.1+355100to the neutrino alert [50].Ref. [51] studied the predicted multifrequency emission during the X-ray flare and the neutrino event, focusing on the one-zone leptohadronic model, the proton synchrotron, and the intergalactic cascade scenarios, to explain the neutrino connection with the flaring activity of the EHSP BL Lac 3HSPJ095507.1+355100.On the other hand, this link has also been explained using different models including two zones (e.g., [35]).Moreover, a recent study by [52] reported the temporal and spatial coincidence, within a 2-degree range, between the neutrino event IceCube-211208A and the high-flux phase of the blazar PKS 0735+178.This alignment was seen across many electromagnetic bands, including optical, ultraviolet, X-ray, and GeV gamma-ray.
A new model has been developed by Aguilar-Ruiz (hereafter AR2022 [53]), which explain the multiwavelength observations of the six well-known extreme BL Lacs by a leptohadronic two-zone model.The model includes an inner blob of photo-hadronic activity and an outer blob of SSC processes to describe the VHE gamma-ray observations.When protons in the inner blob are accelerated, they collide with photons from the annihilation line in a sub-relativistic pair plasma, producing photo-hadronic interactions.Furthermore, the authors extrapolate their findings to the remaining TBLs and determine the neutrino flux generated by the present TBL population under the Aguilar-Ruiz scenario (herafter AR2023 [54]).Taking into account the gamma-ray luminosity function, they found that if protons are accelerated to energies less than ε p 1 PeV, the estimated neutrino flux of the TBL population is consistent with the estimated neutrino diffuse flux.However, they noted that the inner blob might boost proton velocities by as much as ∼300 PeV.
The findings from a recent study on X-ray polarization (e.g., [55,56]), indicate a discrepancy between the degree value of X-ray polarization and the corresponding radio-to-optical value for Mrk 501 and Mrk 421.These results provide further evidence for the validity of blazar's jet models that extend beyond the one-zone framework.Therefore, by considering those track-like neutrino events with temporal and spatially association with EHSP BL Lacs, we apply this novel scenario (AR2022 and AR2023) to EHSP BL Lac type, i.e., those where the hard-TeV-gamma-rays spectrum is missing, to investigate the link with HE neutrinos.
The structure of this paper is as follows: in Section 2, we provide a comprehensive summary of our model.In Section 3, we present and model our sample, and finally, Section 4 discusses the results and conclusions are presented.The Λ-CDM cosmology was used throughout this entire study, where H 0 = (67.4± 0.5) km s −1 Mpc −1 , matter energy density Ω m = 0.315 ± 0.007, and dark energy density

The Model
AR2022 proposes a model in which two emission zones are used to represent the complete SED of six TBLs in their quiescent state.In order to reduce the tension in the parameter required by the one-zone SSC model, the authors posited the need for two dissipation regions, which they dubbed "inner" and "outer".The inner blob is located close to the jet's base, while the outer blob is located distant from the engine.The formation of a pair plasma which arises and is launched above the accretion disc was also utilized by the author.This pair plasma produces a spectrum with a sharp peak at an energy of around ε pl = 511 keV.
Electrons and protons moving at relativistic speeds in the inner blob scatter photons from the pair plasma.Electrons interact via the Compton scattering mechanism, whereas protons engage in photohadronic processes (photopion and photopair).The inner blob's maximal proton energy in AR2022 is chosen to be ε p,max < 100 TeV, hence, the photopion process yields neutral pions, which decay into gamma rays with energy approximately ∼1-30 TeV.Electrons in the inner blob may also have a MeV signature caused by external Compton scattering, while secondary electrons cool down via synchrotron radiation with a signal at the radio-to-optical and MeV bands.The pair-plasma radiation greatly attenuates contributions below the GeV band (see Figure 2b in [53] ), hence MeV-band fluxes are not very prominent.In addition, the outer blob formed by relativistic electrons is explained by a SSC model; proton emission is unimportant in this blob.
Nonetheless, AR2023 notes that maximal proton energies inside the inner blob might reach values of ε p,max ∼ 100 PeV, if this is the case, we should see a peak in gammaray and neutrino emission at energies over ∼5 PeV.The authors note, however, that in order for the estimates of neutrino diffuse flux and the total flux predicted by the existing population of TBLs to be compatible, the maximum proton energy cannot be greater than ε p,max ≈ 1 PeV (see Figure 3 in AR2023 [54] ).Based on this finding, TBL is incapable of producing neutrinos with energy greater than E ν,max ≈ 150 TeV (D i /3).It is worth noting that this criterion does not apply if an EHSP does not have a VHE hard-spectrum peaking at 1 TeV; hence, in certain EHSP, the neutrino spectrum may extend even to E ν 5 PeV.
The observed, pair plasma, and blob (inner or outer) frames used here are the same as those used in AR2022.For observed quantities, we use Latin capital letters followed by the superscript "ob" but the AGN frame will use only the letters themselves.In Greek, the plasma, inner, and outer blobs are represented by the lowercase letters (unprimed), (prime), and (two-prime), respectively.We also take into account the on-axis case for relativistic blobs, where the viewing angle is θ obs 1/Γ and the Doppler factor defined by D = [Γ(1 − β cos θ obs )] −1 , with Γ the Lorentz factor in the blob.

The Inner Blob
VHE gamma rays are thought to be generated in an inner blob by a photohadronic process in the AR2022 model.Here, we suppose that there is one electron for every proton, n e = n p and the spectral indexes of electrons and protons in their distributions are identical α e = α p .
The radiation field of the pair-plasma only affects locations within a few to dozens of Schwarzschild radii, R g , from the center of the engine, with R g = GM • /c 2 ∼ 10 14 cm the Schwarzschild radius for a SMBH mass of M • = 10 9 M ; hence the location r i must be closer to the central engine than the outer blob.The photopion efficiency and the peak of the blob's VHE band emission are both determined by its Lorentz factor, which is favourable to mildly relativistic speeds.If two regions are inside the acceleration and collimation zone, then our scenario is possible, as the Lorentz factor in the inner region is smaller than in the outer region (Γ i < Γ o ).Lorentz factors of Γ i ≈ 1.5 − 3 are preferred for producing VHE-band gamma rays with high photopion efficiency f π ∼ 0.1, as shown in works AR2022 and AR2023.
The magnetic field strength is a crucial factor; we made an estimate of it by assuming that the magnetic luminosity is constant along the jet, L B,i ∼ L B,o .As a result, there is a relationship between the values in the inner and outer blobs, as follows: then, our model's typical values allow us to write it as the magnetic field strength, however, is greater when the equipartition requirement, L B = L p , is applied, thus, we choose these two numbers to be the limits.

Proton Distribution
It is assumed that the Fermi mechanism is responsible for the acceleration of protons in this blob, leading to a spectrum with a spectral index of α p ≥ 1.5 [58].We assume the proton distribution is a single PL distribution, with the form given by where the K p is the normalization constant in units of [GeV −1 cm −3 ], α p is the spectral index and ε p,min = m p c 2 GeV and ε p,max represent the lowest and highest possible proton energies, respectively.The normalization constant may be determined by integrating the proton distribution to produce the apparent bolometric proton luminosity, L p = 4πR 2 i D 4 i ε p ε p N p (ε p ), then for α p = 2, we obtain the expression The competition between the acceleration and loss processes is used to estimate the maximal energy that protons achieve via the Fermi acceleration process.Therefore, by comparing t acc = ε p /(ηeB c) with the dynamical timescale t ad = R /c, we obtain that the maximum proton energy reached by Fermi mechanism is where e is electron's charge and η is the acceleration efficiency, here we take a value of η ∼ 0.1 [59].This value is three orders of magnitude larger than the one needed to describe the TBL spectrum in AR2022.However, in this study we suppose the maximum proton energy is dictated by where χ p→ν ≈ 0.05 is mean the fraction of energy transfer from parent proton to the neutrino, and z is the redshift of the source.Note that the above expression indicates that the value of ε p,max may be lower than the aforementioned value given by Equation (5).
It is important to mention that, due to the neutrino's energies for the events used in this work only represent the reconstructed energy, then Equation ( 6) only provides the minimum value that ε p,max can take in order to explain that events (The reconstructed energy for through-going muon neutrinos events represents the lower limit of the neutrino energy.This is because there are events initiated outside the detector where neutrino signals cannot be recollected [29]).

Seed Photons by the Pair-Plasma
An essential component of the AR2022 model is the existence of a pair plasma that produces an annihilation line with a peak at ε pl = 511 keV (e.g., [60][61][62][63]).It is assumed that the shape of the photon distribution is extremely narrow, to the point that it may be expressed in the delta approximation as [37] where ) is the photon energy density of the pair-plasma.At the photosphere zone, photons escape from the pair-plasma, which moves with velocity β pl ≈ 0.3 − 0.5 [60], and radius of the R ph ∼ R g .
For the pair plasma, the most crucial assumptions are: (i) at the photosphere, it moves at a mildly relativistic speed β pl = 0.5 (Γ pl = 1.15), (ii) it is generated when the luminosity of the disc above 511 keV exceeds where L E ≈ 4 × 10 46 erg s −1 is the Eddington luminosity equivalent to the typical SMBH's mass of BL Lacs M • ∼ 3 × 10 8 M (i.e., [64,65]), hereafter we assume this value for our EHSP sample, this agrees with the value reported for 3HSP J095507.9+355101[66].
When seen in the inner blob, photons from the pair-plasma undergo a redshift.Consequently, the inner blob's measured energy and energy density are [67] ε pl ε pl /(2Γ rel ) and respectively, where Γ rel is the Lorentz factor between the pair plasma and the inner blob, denoted by Γ rel , may be written as where β i is the velocity of the inner blob and Γ i is the Lorentz factor for that blob.In this case, we may express the relationship between the Doppler effect and the Lorentz factor as D i 2Γ i .

Photopion Process
The photopion process results in the production of nonstable secondary products, i.e., π 0 and π ± mesons [68]: where n 0 , n − , n + are the multiplicities of neutral, negative and positive charged pions, respectively.Pions decays into final stable particles as If we consider the photons of the pair-plasma in the inner blob frame, we can calculate the minimum proton energy required to produce a single rest mass pion by photopion interaction, Moreover, efficiency of photopion process is determined by the expression [69] where ε th ≈ 150 MeV is the threshold energy, ε r is the photon energy in the proton restframe, σ pπ is the photopion cross-section and K pπ is the inelasticity coefficient.We use the two-step approximation, given by σ pπ K pπ = 68 ¯barn for 199 MeV ≤ ε r ≤ 500 MeV and σ pπ K pπ = 72 ¯barn for ε r > 500 MeV.

High-Energy Neutrinos
Neutrino generation rates at HE energies are estimated using the approach of [70], then the observed spectrum is given by the proton source is related with proton distribution as Q p 4πR 2 i cN p , and the observed neutrino energy is related with the proton energy as ε p E ν (1 + z)/( χ p→ν D i ).

Gamma-Rays
Relativistic kinematics is used to obtain a rough approximation of the energy of gamma-rays for the decay of neutral pions.Then, gamma-rays are produced with energies above whereas the maximum energy of gamma-rays is determined by the greatest proton energy and hence by the observed neutrino energy, E ob γ 2 E obs ν .This last assumptions differ with the AR2022 work because they were assume ε p,max ≤ 100 TeV and therefore gamma-rays with energies higher than E ob γ ∼ 10 TeVD i hardly can be produced.Finally, the following equation is used to derive the gamma-ray spectra [70] where C abs is the absorption coefficient from either pair-plasma photons or the light from the cosmos in general (EBL).

Secondary Electrons/Positrons
The electron-positron production rate by photopion process is given by where χ pπ p→e 0.05 is mean the fraction of energy transferred from parent proton to positrons/electrons.
The photopion proton threshold may be used as a rough approximation of the lowest Lorentz factor of the electrons As the photopion efficiency is flat above the threshold (see Figure 1 in [54],) the form of the electron spectrum must be adequately tracked by the proton spectrum, and the maximum value of photopion electrons is therefore dictated by the maximum proton energy γ max e,pπ χ p→e ε p,max /(m e c 2 ), therefore this value must be linked to the observed neutrino energy.

Photopair Process
Similarly, the efficiency due to photopair (pe) or Bethe-Heitler (BH) process, p + γ → p + e − + e + , is calculated as [71,72] where ε th 2m e c 2 /γ p is the threshold energy, α f the fine structure constant, the function ϕ is given and parametrized by [72], and σ T is the Thompson cross-section.Therefore, in our scenario electrons are created with Lorentz factor however, the maximum value is much larger by an order of magnitude.
The electron-positron production rate by Bethe-Heitler processes is given by where χ BH p→e m e /m p is mean the fraction of energy transfer from parent proton to positrons/electrons [73].

Synchrotron Emission of Secondary Pairs
We investigate in depth the emission of secondary pairs (from now on referred to as secondary electrons); such emission must span from near the MeV band to the GeV band.At those energy, pair-plasma photons significantly decrease the flux of MeV-to-GeV gamma-rays.Yet, as previously explained, the observed flux below MeV energies is only partly absorbed.On the other hands, gamma-rays with energies greater than hundreds of GeV can escape from the inner blob without attenuation.The emission of secondary pairs is discussed more below.
Since synchrotron and adiabatic losses dominate electron cooling, the distribution of secondary pairs formed by photohadronic processes (photopion and photopair) is computed as where γe,syn and γe,ad are the rates of synchroton and adiabatic losses.Additionally, by taking into account the typical magnetic field in the inner blob, we may calculate the energy of synchrotron photons The maximum energy of synchrotron photons is defined by the maximum proton energy, and applying the condition of χ pπ p→e χ p→ν , we get because this number is outside of the energy range where gamma-rays are optically thick to γγ absorption, GeV-gamma-rays should aid in constraining our model.On the other hand, photopair pairs produce synchrotron photons with energies In a condition of strong flare, when the fluxes of VHE gamma rays and HE neutrinos are both high, it may be possible to measure the flux of synchrotron secondary pairs.Produced fluxes from photopion and photopair might be connected to proton luminosity as L pπ γ,syn ≈ (1/8) f pπ L p and L BH γ,syn ≈ f BH L p , respectively (see [37,74]) in an environment like the inner blob, where electrons are effectively cooled down by the synchrotron process.
It is important to note that the cross-sections for photopion and photopair processes lead to a flux ratio of L fl,pß /L fl,pe ∼ σ pk pß /(8σ pk pe ) ∼ 10 2 .Our model predicts that the photopair process may only lead to a signature during extremely powerful flare states.
Gamma-ray and secondary-pair fluxes are connected in the photopion scenario as L pπ γ ≈ (4/3)L ν and L pπ γ,syn ≈ (1/3)L ν [70].This suggests that gamma rays and the synchrotron flux of secondary pairs may leave a signature when the HE-neutrino flux is high enough to be detected.

Proton Synchrotron
Accelerated protons inside a dissipation region with a very high magnetic field must radiate via synchrotron, using the Equation ( 6) the characteristic photon energy produced by proton-synchrotron is given by furthermore, the loss timescale is calculated by the following equation which is larger than the dynamical timescale, t dyn R i /c ≈ 6.7 × 10 3 s.However, note that proton synchrotron efficiency must be about t p,syn /t dyn ∼ 10 −2 if the magnetic field intensity is higher, for instance, ∼10 4 G (as the predicted by the equipartion condition indicated by Equation ( 2).This efficiency is similar to that for gamma-ray luminosity produced by secondary electrons is

The Outer Blob
Although the outer blob of the AR2022 model is irrelevant to neutrino generation, it does create the remainder of the SED's emission, which we explain below along with the other important aspects considered and applied in this work: (i) As the electron population is larger than the proton population (n e n p ), only leptonic processes are taken into account.
(ii) The blob travels at relativistic speeds, thus we use Γ o = 5 and D o 10 as Lorentz and Doppler boost factors, respectively.Our choice agrees with the average value, D ∼ 13, reported by [75] who conducted a study with a sample of 217 blazars and utilized a one-zone leptonic model.Even though if D o take greater values, our model could perform the same fitting, but at the cost of needs smaller emission regions or stronger magnetic field.And and interesting consequence is the reduce of the value for the equipartition parameter U B /U e .(iii) Using the assumption that the blob's variability timescales is R o D o t var c with t var ∼ one day and D o ∼ 10, we find that the blob's distance from the SMBH is r o = 2Γ o R o ≈ 2.6 × 10 17 cm, this is equivalent to around ∼ 10 3 R g for a SMBH mass of M • = 3 × 10 8 M .This separation is in close proximity to the hypothesized acceleration and collimation zone (e.g., [76,77]) suggesting that Γ o may be close to the terminal Lorentz factor of the jet.(iv) It was possible to determine the magnetic field by looking at the peak of the synchrotron bump, which was around ν syn ∼ 10 17 Hz.Now when the timescales of the synchrotron and the adiabatic losses are equivalent, we have the electron energy break, which is written as hence, this blob's magnetic field strength is with γ e,min , γ e,br , γ e,max are the minimum, break and maximum Lorentz factors of ultrarelativisc electrons, respectively.
Using the blob's properties as constraints, we may estimate the parameters that defined the electron distribution.If we use approximation Equation (30) to calculate the break Lorentz factor, we obtain Nevertheless, the spectral indices and minimum and maximum Lorentz factors are calculated by fitting the broadband SEDs.

High-Energy Neutrinos Coincident with Extreme High Synchrotron Peaked BL Lacs
This work focuses on the study of track-like astrophysical neutrinos that have the potential to be linked to an EHSP BL Lac object.Two of the events, namely IC190819A and IC190922A, are included in the latest IceCube's catalog [29].However, the third event, i.e., IC200107A, is not listed, we took it into consideration in view of the sequence of publications published by many authors examining the potential correlation [35,51,78].The Table 1 contains information regarding neutrino events and the EHSP BL Lac associated with them that are considered in this study.
Table 1.This table shows the neutrino events utilized for this study.The neutrino event appears in column (1).The value in column (2) represents the reconstructed neutrino energy.The coordinates for the best fit given for right ascension and declination are displayed in columns ( 3) and (4).Column ( 5) is the effective area of muon neutrinos at the neutrino energy event [29,79].The EHSP BL Lac associated with each neutrino event is shown in column (6).Column (7) is the synchrotron peak frequency.Column (8) is the redshift.Column (9) is the luminosity distance and column (10) is the all-favor neutrino luminosity level at the observed neutrino energy, which is calculated using 10 yrs of window time.16 TeV.This neutrino event was discovered at the same location (within 2.01 degrees) as the object 4FGL J0946.2+0104 in the Large Area Telescope (LAT) (This instrument is onboard Fermi Satellite) catalog, which is linked to the EHSP BL Lac 1RXS J09462.5+010459.Analysis of the Fermi-LAT data on one-day, one-month, and months-to-years timescales revealed that this object did not exhibit an increase in the gamma-rays flux during the neutrino detection [81].In addition, the High-Altitude Water Cerenkov (HAWC) observatory looked for both steady and transient sources near the neutrino alert site.Steady source searches in records from November 2014 to May 2018 yielded a 95 percent confidence level upper limit of E 2 dN/dE = 3.2 × 10 −13 TeV cm −2 s −1 at 1 TeV for a spectral index of −2.3 [82].Additionally, the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL) multi-MESSENGER collaboration reports that this neutrino event has no counterpart [83].Finally, Figure 1 shows the light curve of 1RXS J09462.5+010459for Fermi-LAT observations between October 2019 and February 2023.This figure also shows the time of neutrino IC190819A's arrival.It is essential to observe that there is not notably strong flux variation of gamma rays over around the neutrino arrival.[84], which was chosen by the Gold alert stream, with position at R.A. = 167.30+2.81 −2.72 and Dec = −22.27+3.39 −3.31 and reconstructed energy E ob ν = 3.1 PeV.There are three gamma-ray sources in the 4FGL Fermi-LAT catalog that are within the 90% uncertainty region of the event.The distances from the best-fit location to the sources are 2.4 degrees for 4FGL J1120.0-2204,2.8 degrees for 4FGL J1100.0-2044, and 1.77 degrees for 4FGL J1103.6-2329.The first two are unknown sources, while third is the EHPS BL Lac tagged as 1ES 1101-232 is the only known source among the three.

HE-Neutrino
During the neutrino detection, these object did not show an increase in the flux of gamma rays, as shown by an analysis of the Fermi-LAT data on 1-day, 1-month, and months-to-years timescales [85].In a similar vein, the HAWC observatory looked for both steady and transient sources around the neutrino alert location.Even though this neutrino was not inside HAWC's field of view, the collaboration nonetheless searched for it in archival data from November 2014 to May 2018 and established a 95 confidence level of E 2 dN/dE = 1.2 × 10 −12 TeV cm −2 s −1 at 1 TeV for a spectral index of −2.3 [86].It has also been reported by the INTEGRAL multi-MESSENGER collaboration that this neutrino event has no counterpart.[87].
Finally, Figure 2 shows the light curve of 1ES 1101-232 for Fermi-LAT observatory observations between October 2010 and February 2023.This figure also shows the time of neutrino IC190922A's arrival.It is essential to observe the flux variations of gamma rays over a period of ten years not notably strong.ν , E ob ν = 333 TeV is the expected energy.There are two Fermi 4FGL sources inside the 90% C.L. reconstructed arrival localization region: The nearest gamma-ray source is 4FGL J0955.1+3551(3HSP J095507.9+355101), which is 0.63 deg from the best-fit event position.4FGL J0957.8+3423 is a second source that is only 1.50 degrees from the best location of the event.Neither object showed any significant flux increases in subsequent observations with HAWC, Fermi-LAT, or INTE-GRAL.Between 2020 January 07 06:02:39 UTC and 2020 January 07 12:24:01 UTC (∼ 6 h), HAWC reported an upper limit of 4.108 × 10 −12 TeV cm −2 s −1 at 1 TeV [89].Despite this, Swift observations show that 4FGL J0957.8+3423 had higher X-rays emission (observed 2−10 keV flux) during the neutrino observation than the archival data from 2012-2013 by a factor of ≈ 2.6, indicating that the source was behaving in flaring state during the neutrino detection [90].Nevertheless, the second source could not be identified, and only an upper limit was set.
Moreover, Ref. [91] highlights potential evidence for hour-scale flux variability in the X-ray band using NuSTAR and NICER light curves.This limits the size of the emission zone to less than R 3.5 × 10 14 cm D. Furthermore, as shown by the soft X-ray data the maximal intensity value is obtained after the neutrino detection and remains constant for nearly 30 days [78].

Modeling the SED of EHSP
In this section, for each EHSP, we use the model outlined in Section 2. The SSC model is used to describe the outer blob's emission [92].For the emission of the inner blob (where neutrinos are produced in our model), the following considerations are taken into account as a consequence of the earlier research works performed by [53,54]: the Doppler boost for a jet aligned with the line of sight is D i = 2.6 for Γ i = 1.5, and the blob's size is R i = 2 × 10 14 cm. Figure 3 shows the proton timescales for photopion, photopair, synchrotron, adiabatic, escape and acceleration processes.We calculate the photopion and photopair timescales by considering an annihilation line luminosity L keV = 3 × 10 44 erg s −1 and Γ i = 1.5.The escape timescale was calculated considering a diffusion process, i.e., t esc R i 2 /(2D diff ) where D diff 1 3 cηr L is the diffusion coefficient in the Bohm regime and r L = ε p /eB is the Larmor radius.As Figure 3 indicates photopion process is the most relevant below ∼1 PeV and ∼100 PeV when we consider a value of magnetic field B = 10 4 G and B = 10 3 G, respectively, while for higher energies synchrotron losses become the dominant process and suppressing photopion process.Using Equation ( 6), the maximal proton energy is calculated as a function of the observed neutrino's energy, By considering the monochromatic spectrum of the 511 keV annihilation line (see Equation ( 7)) and the photopion processes efficiency equation given in Equation ( 15), the photopion efficiency is determined to be The proton luminosity is estimated using the relationship ε p L ε p (8/3) f −1 pπ ε ν L ν , where luminosity of all-flavor neutrino for each event is estimated using the relationship , where ∆T ob is the observation time window, and A µ,eff is the effective area for muon neutrinos dependent of the neutrino energy and the neutrino declination and is given in [31].Furthermore, the bolometric proton luminosity is estimated using the expresion, We consider two different scenarios to model the SED, and we describe each one in the following.

Scenario A:
We use it to describe the steady state of EHSP BL Lacs where no significative electromagnetic signature was observed during the detection of the neutrino events (i.e., 1RXS J09462.5+010459 and 1ES 1101-232).For this scenario we consider the minimum requirement necessary to produce the annihilation line by the pair-plasma, L keV = 3 × 10 −3 L Edd ≈ 1 × 10 44 erg s −1 [60]; It is important to note that in these circumstances, the photopion process has a value of f pπ ≈ 0.1.For these EHSP BL Lacs we take into account the available archival multifrequecy data in taken from https://tools.ssdc.asi.it/SED/accessed on 1 May 2023.

Scenario B:
We use it to describe the multifrequency observations of 3HSP J095507.9+355101 during the detection of the neutrino event.We take into account the X-ray flaring state during the arrival of IC200107A to constrain the parameters of the inner blob; meanwhile the archival data is described by SSC emission in the outer blob.In this scenario, we assume the pair-plasma luminosity grows, L keV > 3 × 10 −3 L Edd , and therefore the photopion efficiency improves to the point where it equals the unity, which happens for L keV 3 × 10 −2 L Edd .
Finally we calculate the SED applying our model, and taking into account the magnetic field strength in the range of values suggested by the Equations ( 1) and (2).Furthermore, the distribution of stable particles will be broken if B i has a very large value, as secondary unstable charged particles created by photopion process may be cooled down by the synchrotron process before they decay.However, if the value of B i is very strong, secondary unstable charged particles produced by photopion processes, might cool down by synchrotron process before they decay into other particles; therefore, this must be reflected in a break in the stable particle distribution.The cooling timescale for muons and charged pions through synchrotron is comparable, however muons decay slower than charged pions, with a mean lifetime of τ µ = 2.197 × 10 −6 s and a lifetime of τ π ± = 2.6 × 10 −8 s, respectively [93].As a result, we focus solely on muons, and the energy break can be roughly estimated by comparing the synchrotron loss timescale, t , with the muon lifetime in the inner blob frame t lt,µ ∼ ε µ /(m µ c 2 )τ µ .Taking into account the fact that neutrinos only carry around a third of the energy of the parent muon [94], we predict that an energy break must occur in the observed all-flavor neutrino spectrum at as a result, this effect is only significant for PeV neutrinos and extremely strong magnetic fields.Furthermore, Figure 3 indicates that for such intense magnetic fields synchrotron losses becomes more relevant than photopion above ε p ∼ 1 PeV making hard the production of neutrinos with energies greater than therefore, PeV neutrinos (as the case of IC190922A) are produced more efficiently in environments with magnetic fields that no exceed B i ∼ 10 3 G.

1RXS J09462.5+010459
1RXS J09462.5+010459 is located at a redshift z = 0.577 which corresponds to a luminosity distance of d L ≈ 3.48 Gpc.Using Equation ( 6), the maximum proton energy must be ε p,max 1.2 PeV 1+z 1.577 , then we use 10 PeV.We model this EHSP BL Lac using historical data and use the scenario A, as it did not show any significant flux increases in any electromagnetic band.Using the effective area at the neutrino event's energy, we calculate the all-flavor neutrino flux to be E ν F ν ≈ 5.9 × 10 −12 erg cm −2 s −1 for ∆T ob = 10 yrs.This is equivalent to a neutrino luminosity E ν L E ν ≈ 8.52 × 10 45 erg s −1 .Then the proton luminosity at ε p,max ≈ 1.2 PeV is ε p L ε p | ε p,max ≈ 2.37 × 10 47 erg s −1 D −4 i with photopion efficiency equal to f pπ ≈ 0.1.Likewise, using Equation ( 36) and taking into account proton spectrum with α p ≤ 2, the bolometric proton luminosities for α p = 1.5 and α p = 2 are L p ≈ 5 × 10 46 erg s −1 and L p ≈ 3.5 × 10 47 erg s −1 , respectively.
Using the parameters that describe the emission of the outer blob and Equation ( 1), the strength of magnetic field is B i ≈ 683 G; But, if we use the equipartition condition given by Equation ( 2) and the proton luminosities above calculated, we obtain B i,eq ≈ 1.1 × 10 4 G and B i,eq ≈ 2.8 × 10 4 G for α p = 1.5 and α p = 2, respectively.Observe that these values are one or two orders of magnitude greater than those derived from Equation (1).It is essential to note that with the larger value of B i predicted by the equipartition condition, the neutrino break must appear at ∼100 TeV energies (see Figure 3) then we adopt B i,max ≈ 1 × 10 4 G to explain the IC190819A.In this work, we use as limit values of our model in the SED of 1RXS J09462.5+010459:α p = 1.5 − 2. and B i = 683 G to B i,max = 1 × 10 4 G for each case.Our results are plotted in Figure 4 and the parameters used to described the outer blob and the inner blob are listed in Table 2 and Table 3, respectively.Figure 4 shows our resulting SED for each case considered for this EHSP BL Lac.Our results indicate differences between each other especially in hard X-rays to MeV gamma-rays, where the model with α p = 2 and B i = 1 × 10 4 G produces the hardest X-ray spectrum, making it the most distinctive signature.The neutrino fluxes for spectral indexes of 1.5 and 2 is displayed in web-green color with solid and dotted line, respectively.The results of modeling the SED of 1RXS 094620.5+010459for various spectral indexes and magnetic fields are presented with black lines (solid, dashed, dotted, and dotted-dashed).Moreover, an orange line designates the annihilation line.In addition, we illustrate the contribution of SSC in the outer blob (gray line), the gamma-rays after correction for the EBL attenuation effects (orange dotted line), the synchrotron of photohadronic pairs (salmon dotted line), and the synchrotron of protons (violet dotted line) for the model with α p = 1.5 and B i = 1 × 10 4 G.

1ES 1101-232
The redshift of 1ES 1101-232, which is spatially coincident with the neutrino event IC190922A, is z = 0.186, which corresponds to a luminosity distance of d L ≈ 879 Mpc.Using Equation ( 6), the maximum proton energy is ε p,max 38.61 PeV 1+z 1.186 , then we use 100 PeV.We model 1ES 1101-232 using the archival data as it did not show significant flux increases in any electromagnetic band.The associated all-flavor neutrino flux is E ν F ν ≈ 5.9 × 10 −11 erg cm −2 s −1 and the corresponding neutrino luminosity is E ν L E ν ≈ 5.5 × 10 45 erg s −1 taking into account for a window time of 10 yrs.Applying the scenario A, the photopion efficiency is f pπ ≈ 0.1, hence the proton luminosity at Therefore the bolometric proton luminosities for α p = 1.5 and α p = 2. are L p ≈ 2 × 10 46 erg s −1 and L p ≈ 1.6 × 10 47 erg s −1 , respectively.
Using the parameters that describe the emission of the outer blob and Equation ( 1), the strength of magnetic field is B i ≈ 986 G; But, if we use the equipartition condition given by Equation ( 2) and the proton luminosities above calculated, we obtain B i,eq ≈ 8.6 × 10 3 G and B i,eq ≈ 2.5 × 10 4 G for α p = 1.5 and α p = 2, respectively.It is important to be noted that the magnetic field cannot exceed B i, ≈ 1 × 10 3 G in order to successfully explain the energy of the neutrino event IC190922A and prevent the neutrino break (see Equation (37) and Figure 3).
Finally, similar to the case of 1RXS J09462.5+010459,we model the SED of 1ES 1101-232 using spectral indexes of α p = 1.5 − 2. and magnetic field value of B i = 1 × 10 3 G. Figure 5 shows our results, and Tables 2 and 3 enumerate the parameters used to describe the outer blob and the inner blob, respectively.
Our results are plotted in Figure 5 and the parameter used to described the outer blob and the inner blob are listed in Table 2 and Table 3, respectively.The results shown in Figure 5 indicate that for the cases considered, there are electromagnetic differences, particularly in the MeV and GeV energy bands, where cases with α p = 2 produce higher fluxes than those with α p = 1.5.Hence, models with α p = 2 may be ruled out because the resulting flux exceeds of HAWC's upper limit (UL) set between November 2014 and May 2018.
The maximum proton energy assumed for the rationale of IC200107A is ε p,max 3.4 PeV 1+z 1.557 , then we use 10 PeV.In contrast to the two other EHSP, this one exhibits an increase in X-ray flux during neutrino detection, and the light curves indicate that it remains constant for at least thirty days after neutrino arrival [78].The fact that neutrinos were produced during flaring implies that the luminosity of neutrinos must be significantly greater than that required for a window time of 10 yrs.In our analysis, we therefore consider a time window of ∆T = 250 days, hence the all-flavor neutrino luminosity is E ν L E ν ≈ 4.7 × 10 47 erg s −1 (Our approximate values are comparable to those reported by [78], they calculate the bolometric neutrino luminosity for 30 days, 250 days and 10 yrs obtaining L ν ≈ 4 × 10 48 erg s −1 , 5 × 10 47 erg s −1 , 3 × 10 46 erg s −1 , respectively).Hence, in the context of scenario B and the previously mentioned neutrino luminosity, the bolometric proton luminosities for α p = 1.5 and α p = 2 are L p ≈ 2.7 × 10 47 erg s −1 and L p ≈ 2.1 × 10 48 erg s −1 , respectively.
During this flaring state, the equipartion magnetic field is in the range B i,eq ≈ (2.5 − 7.0) × 10 4 G.However, such an intense magnetic field prohibits the production of neutrinos as energetic as the IC200107A event.To avoid muons cooled by synchrotron losses, the magnetic field's maximal value must be less than B i,max ≈ 1 × 10 4 G (see Equation ( 37)).On the other hand, using the Equation ( 1), the strength of magnetic field is B i ≈ 655 G.
We use the following values to model the SED of 3HSP J095507.9+355101:α p = 1.5 − 2 and B i = (0.66 − 10) × 10 3 G.Finally, our results are shown in Figure 6, and Tables 2 and 3 list the parameter used to describe the outer blob and inner blob, respectively.Figure 6 demonstrates that our model produces different results depending on the selection of the set of parameters, with the most significant differences occurring in the radio/IR, MeV, and VHE gamma-rays bands, where models with α p = 2 predict higher fluxes than those with α p = 1.5.Therefore, during future observations of blazar flares and neutrino events, the radio/infrared signature produced by the photopair process could be utilized to test our model.It is important to note that our results indicate that the evolution of X-rays can be explained by varying the magnetic field strength from B i ≈ 655 G to ≈2 × 10 4 G while keeping α p = 1.5 constant.In addition, our results indicate that none of our hypotheses can account for the ultraviolet (UV) emission observed during neutrino detection, indicating that this phenomenon may result from a different mechanism.

Discussion and Conclusions
This work proposes a two-zone model to describe the probable relationship between EHSP BL Lacs and track-neutrino detection.We have demonstrated that a recent model proposed by AR2022 and AR2023 may adequately characterize the SED of these objects.In contrast to AR2022, we have considered the emission of secondary pairs produced by photohadronic processes (photopion and photopair).This emission peaks in the hard X-ray band and extends to GeV gamma rays for ε p,max ∼ 10 − 100 PeV.
Such emission could generate a signature in the MeV band that could be detected by future telescopes such as AMEGO or eASTROGAM, mainly either during electromagnetic flare episodes in BL Lacs or intense HE-neutrino flares.In general, our results indicate that our model provides a reasonable explanation for the SED of EHSP from X-rays to gamma-rays and its connection with HE-neutrino observations under the assumption of a prolonged steady state lasting ∼ 10 yr.However, in the case of 1RXS J09462.5+010459,our model fails to adequately explain the data in the radio to optical wavelength range.Important to mention that in this work, we do not consider the emission produced by the host-galaxy, which can remove part of the discrepancy between our model and optical data [15].
During flare states, such as 3HSP J095507.9+35510,our model can explain the X-ray spectrum's hardening due to proton-synchrotron and synchrotron emission of photopion secondary pairs.However, our model cannot explain the UV-to-soft-X-ray flux observed during the neutrino arrival, which may indicate that it is produced by a mechanism not discussed here.
The bolometric proton luminosities required to explain the broadband emission of 1RXS J09462.5+010459 and 1ES 1101-232 are a factor of ∼2 greater and less than the Eddington limit, respectively (∼0.4 ×10 47 erg s −1 assuming a SMBH mass M • = 3 × 10 8 M ) considering f pπ ≈ 0.1 and α p = 1.5.However, the situation changes when we attempt to correlate the neutrino event IC200107A with 3HSP J095507.9+355101within a 250 days time window, i.e., the required luminosity is L p ≈ (3.7 − 26.6) × 10 47 erg s −1 even with the maximum process efficiency f pπ ≈ 1.In general, a harder proton distribution provides a better explanation for the SED considered in this work demanding lower jet power than softer cases.Hadronic models requiring proton luminosities greater than Eddington limit have been pointed out as unavoidable implications for modeling blazar's SED, especially for proton-synchrotron scenario [96], although leptohadronic models offer moderate jet's power for EHSP (e.g., [97]), TeV emission and neutrinos events association demands higher power than Eddington limit (e.g., [51,52,98]).Despite Eddington luminosity provides an estimation of the maximum jet's power, it may not be a strict limit, and its violation might occur as a result of non-spherical geometry or the presence of various types of instabilities (i.e., [99]).Therefore our results of HE-neutrino production in EHSP's jets under our scenario may suggest another type of accretion in such environments.
A crucial and essential aspect of our model is the existence of a 511 keV annihilation line, which is the source of the seed photons for photohadronic interactions.This emission has been observed in micro-quasars and may also exist in SMBH phenomena.In fact, for two EHSPs modeled as steady state, no signature in any electromagnetic frequency can be distinguished from the SSC contribution.The annihilation line is the only characteristic that can be used to test our model.However, this feature is not always observable, especially during intense flare states when masked by the synchrotron emission of protons and secondary pairs, as in 3HSP J095507.9+35510.
In this study, we have considered the minimum luminosity required to generate an optical depth environment that leads to the formation of the annihilation ∼1 × 10 44 erg s −1 .During the flaring state, however, we have supposed this value should be higher, 3 × 10 45 erg s −1 , to achieve the maximum photopion efficiency, f pπ ≈ 1.
The increase in the luminosity of the annihilation line implies an increase in the optical depth to pair-creation via photon-photon collision in the MeV-to-GeV energy range.Out of this band, gamma-rays led to a signature that is not entirely attenuated and can be used to constrain our model with future observations; see Figures 5 and 6.However, we have not considered the effect of the broad-line region, which could reduce the flux at VHE energies.Unfortunately, observations do not reveal evidence of the existence of such region, and only for 3HSP J095507.9+35510 an upper limit, L BLR < 3 × 10 −3 L Edd , is determined.In 3HSP J095507.9+35510, the tentative evidence of rapid variability timescales suggests an emitting region size of R 3.5 × 10 14 cm lends support to our model because our proposed inner blob has comparable typical sizes.

Figure 3 .
Figure 3. Proton timescales in the comoving frame for photopion (red line), photopair (orange line), synchrotron (pink lines), adiabatic (green line), escape via diffusion (magenta line) and acceleration (blue line) processes.For photopion and photopair processes we consider the seed photons from the pair-plasma with an annihilation line luminosity L keV = 3 × 10 −3 L Edd , and boost Lorentz factor of the blob Γ i = 1.5.The vertical dashed line indicates the proton energy at 1 PeV.

Figure 4 .
Figure 4.The neutrino fluxes for spectral indexes of 1.5 and 2 is displayed in web-green color with solid and dotted line, respectively.The results of modeling the SED of 1RXS 094620.5+010459for various spectral indexes and magnetic fields are presented with black lines (solid, dashed, dotted, and dotted-dashed).Moreover, an orange line designates the annihilation line.In addition, we illustrate the contribution of SSC in the outer blob (gray line), the gamma-rays after correction for the EBL attenuation effects (orange dotted line), the synchrotron of photohadronic pairs (salmon dotted line), and the synchrotron of protons (violet dotted line) for the model with α p = 1.5 and B i = 1 × 10 4 G.

Figure 5 .
Figure5.The black lines (solid, dashed, dotted, and dotted-dashed) represent the results of modeling the SED of 1ES 1101-232 for different spectral index and magnetic fields; the corresponding neutrino fluxes for spectral indices of 1.5 and 2 are shown in web-green color (solid and dotted line).In addition, the annihilation line emission is drawn in orange.The contribution of the SSC in the outer blob (grey line), the resultant gamma rays after EBL attenuation (orange dotted line), the synchrotron of photohadronic pairs (salmon dotted line), and proton synchrotron (violet dotted line) are all illustrated for the model with α p = 1.5.

Figure 6 .
Figure 6.Neutrino fluxes for spectral indices of 1.5 and 2 are shown in web-green color (solid and dotted lines), and the results of modeling the SED of 3HSPJ095507.9+355101assuming a time window for the neutrino event of 250 days are presented with black lines (solid, dashed, dotted, and dotted−dashed).Also, an orange line represents the annihilation line, and a grey solid line represents the SSC's contribution from the outer blob.
[78]is values was calculated assuming a E −2 neutrino spectrum, however if the spectrum is taken as E −1 and E −2.7 then neutrinos have energies of 1.40+5.75−1.22PeVand0.16 +0.83 −0.12 PeV, respectively[78].b The 1ES 1101-232 light curves for the Fermi−LAT, Swift−BAT, and ASAS−SN Telescopes are displayed from October 2010 to February 2023.An illustration of neutrino event IC190922A is shown in red.The yellow shadow zone represents average count rate with a significance level of 1σ over 10 year of observation.detected by IceCube on 2020 January 07 at 09:42:18.36UT, in fact, belong to the starting track event category, despite not matching the Bronze or Gold starting track classification [88].R.A. = 1489.18+2.20 −1.83 and Dec = 35.46+1.10 −1.22 are the coordinates of the reconstructed arrival direction, and for a neutrino spectrum of E −2 3.1.3.IC200107A and 3HSP J095507.9+355101Thenew neural network classifier confirmed that a high-energy starting event

Table 2 .
Parameters used to model the outer blob during the steady state of EHSP BL Lac coincident with track-like neutrino events.

Table 3 .
Parameters used to model the inner blob of EHSP BL Lac coincident with track-like neutrino events.