{'fraction_non_alphanumeric': 0.039265440657127164, 'fraction_numerical': 0.037224998038035945, 'mean_word_length': 4.277194919933739, 'pattern_counts': {'":': 0, '<': 9, '<?xml version=': 0, '>': 4, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
A two-zone blazar radiation model for"orphan"neutrino flares

In this work, we investigate the 2014-2015 neutrino flare associated with the blazar TXS 0506+056 and a recently discovered muon neutrino event IceCube-200107A in spatial coincidence with the blazar 4FGL J0955.1+3551, under the framework of a two-zone radiation model of blazars where an inner/outer blob close to/far from the supermassive black hole are invoked. An interesting feature that the two sources share in common is that no evidence of GeV gamma-ray activity is found during the neutrino detection period, probably implying a large opacity for GeV gamma rays in the neutrino production region. In our model, continuous particle acceleration/injection takes place in the inner blob at the jet base, where the hot X-ray corona of the supermassive black hole provides target photon fields for efficient neutrino production and strong GeV gamma-ray absorption. We show that this model can self-consistently interpret the neutrino emission from both two blazars in a large parameter space. In the meantime, the dissipation processes in outer blob are responsible for the simultaneous multi-wavelength emission of both sources. In agreement with previous studies of TXS 0506+056 and, an intense MeV emission from the induced electromagnetic cascade in the inner blob is robustly expected to accompany the neutrino flare in our model could be used to test the model with the next-generation MeV gamma-ray detector in the future.


INTRODUCTION
The origin of extragalactic high-energy neutrinos detected by IceCube Neutrino Observatory is still unclear (Aartsen et al. 2013;IceCube Collaboration 2013). As one of the most powerful astrophysical persistent objects, blazars are widely considered as source candidates for the origin of extragalactic high-energy cosmic rays and neutrinos (e.g., Mannheim et al. 1992;Atoyan, & Dermer 2001;Murase et al. 2014;Padovani et al. 2015;Petropoulou et al. 2015;Padovani et al. 2016).
Although stacking analysis done by Aartsen et al. (2017) suggests that blazars cannot account for the dominant diffuse neutrino background (see also Luo & Zhang 2020, but see Paladino et al. (2019), some individual blazars are still likely to be powerful high-energy cosmic-ray accelerators and emit high-energy neutrinos (Kadler et al. 2016;IceCube Collaboration et al. 2018b;Garrappa et al. 2019;Oikonomou et al. 2019;Giommi et al. 2020a).
A breakthrough in neutrino astronomy was made by IceCube in 2017. For the first time, a highenergy (290 TeV) muon neutrino event, IceCube-170922A, was detected in both spatial and temporal coincidence with the γ-ray flare of the known blazar TXS 0506+056 (IceCube Collaboration et al. 2018b) at the significance level of ∼ 3σ, triggering extensive studies on the neutrino-blazar association (Ansoldi et al. 2018;Keivani et al. 2018;Padovani et al. 2018;Sahakyan 2018;Zhang et al. 2018;Banik & Bhadra 2019;Cerruti et al. 2019;Gao et al. 2019;Laha 2019;Liu et al. 2019;Padovani et al. 2019;Righi et al. 2019;Xue et al. 2019a;Cao et al. 2020). More interestingly, IceCube Collaboration et al. (2018a) investigated the archive data of IceCube in the direction of the same blazar, and discovered a 3.5σ excess of 13±5 high-energy neutrino events, also known as 2014-2015 (hereafter 14-15) neutrino flare, in the period between 2014 September and 2015 March. It is worth noting that no evidence of multi-wavelength activity is found during that time (i.e., "orphan" neutrino flare), and the neu-trino flux is about 5 times higher than the average γray flux, probably implying a strong absorption of GeV photons by intense X-ray radiation field (Wang et al. 2018a;Reimer et al. 2019). Such a phenomenon poses a great challenge to the conventional one-zone models. The main difficulty lies in reconciling the comparatively high flux of neutrinos and comparatively low flux of multiwavelength electromagnetic (EM) emissions, because the emission of the EM cascade which accompanies neutrino emission, would easily overproduce the X-ray flux or the GeV gamma-ray flux of the blazar during the neutrino flare. Indeed, Rodrigues et al. (2019) study external photon fields in the broad-line region (including broad emission lines and isotropized accretion disk radiation) as the possible targets, and they found that the derived neutrino flux is not sufficient to explain the neutrino flare on the premise that the model does not overpredict the multiwavelength flux; Petropoulou et al. (2020a) speculated that the non-thermal photons from the layer of a structured jet could be possible targets although a detailed modelling is yet to be further explored.
Given that the conventional one-zone model is only an approximate description to the radiation processes of a parsec-scale (or even larger) jet residing in the complex environment of the very central region of the galaxy hosting an active supermassive black hole (SMBH), it may not contain all the physical processes relevant with the blazar's emission. If the association between the neutrino flare and the blazar is true, additional physical processes and/or new physical pictures have to be introduced. In fact, efforts have already been made to explore possible explanations to the neutrino flare beyond the framework of the conventional one-zone model. The jet-cloud/star model (Wang et al. 2018a) and the neutral beam model ) provide possible solutions but special or extreme conditions need to be satisfied. An additional compact core with enhanced activity inside a conventional emitting blob is considered ) but the resulting neutrino flux is still insufficient. On the other hand, Liu et al. (2019) and Xue et al. (2019a) considered that more than one dissipation (particle acceleration) zones could occur in the blazar's jet and form multiple radiation zones. By invoking two physically distinct emission zones in the jet, with an inner blob close to or inside of the broad-line region (BLR) and an outer blob far away from the BLR (hereafter, we denote the model by the "inner-outer blob model"), they show that the different physical conditions of the two blobs can help suppress the EM cascade emission in the X-ray band. Consequently the theoretical neutrino flux in the model can be enhanced by at least one order of magnitude compared to that in one-zone models, when explaining the simultaneously multiwavelength emission of TXS 0506+056. A detailed compar-ison between the inner-outer blob model and one-zone models is given in Xue et al. (2019a). In the meantime, a rough analytical estimation given in Xue et al. (2019a) suggests that the "inner-outer" blob model may have the potential to explain the 14-15 neutrino flare as well. If the inner blob happens to occur at the jet base, the intense X-ray photons from hot corona of the accretion disk of the SMBH could be the targets for the neutrino production and also absorb the GeV photons so that no gamma-ray flare would be accompanied with the neutrino flare.
More recently, the second possible association between a blazar, 4FGL J0955.1+3551, and a muon neutrino event, IceCube-200107A, is reported (Giommi et al. 2020a;Paliya et al. 2020). 4FGL J0955.1+3551 is a highsynchrotron peaked (HSP; Abdo et al. 2010) BL Lac object at redshift z = 0.557 (Paiano et al. 2020). Its Xray flux was found in a high state with a factor of 2.5 larger than the average flux in 2012-2013 when Swift started the follow-up observation the day after the arrival of IceCube-200107A (Giommi et al. 2020b). Assuming the spectral index of neutrino spectrum to be −2 (N ν (E) ∝ E −2 ) and the effective area for HESE starting tracks (Blaufuss et al. 2019), Giommi et al. (2020a) calculate that the all-flavour neutrino energy flux between 65 TeV and 2.6 PeV for 30 days, 250 days and 10 years are 3 × 10 −9 erg cm −2 s −1 , 4 × 10 −10 erg cm −2 s −1 and 3 × 10 −11 erg cm −2 s −1 , respectively. The GeV gamma-ray flux is at the level of 10 −12 erg cm −2 s −1 which is more than one order of magnitude lower than the neutrino flux even assuming a 10-yr emission period. Giommi et al. (2020a) argue that IceCube-200107A can only be interpreted as an upward fluctuation (see also Petropoulou et al. 2020b) or a random event detected from a numerous population of faint sources (i.e., the Eddington bias, Strotjohann et al. 2019) in the conventional one-zone models, given the huge contrast between the neutrino flux and the GeV gamma-ray flux. Moreover, with analytical estimations, Paliya et al. (2020) conclude that any theoretical model on explaining IceCube-200107A from 4FGL J0955.1+3551 must involve an intense external photon field.
In this work, we aim to explore that whether the 14-15 neutrino flare and IceCube-200107A could physically originate, respectively, from TXS 0506+056 and 4FGL J0955.1+3551 under the inner-outer blob model, without ascribing the detection to the Eddington bias. The rest of this paper is organized as follows. In Section 2, we give a brief description of the model and the setup in this work; in Section 3, we apply the model to explain the multi-messenger emission of 14-15 neutrino flare from TXS 0506+056 and IceCube-200107A from 4FGL J0955.1+3551 and perform an extensive investigation on the parameter space of the inner blob based on the constraints from the EM observations; in Section 4, we discuss our results and give a conclusion. Throughout the paper, the ΛCDM cosmology with H 0 = 70km s −1 Mpc −1 , Ω m = 0.3, Ω Λ = 0.7 is adopted.

MODEL DESCRIPTION
Let us begin with a brief description about the physical picture of the inner-outer blob model and its basic assumptions. Henceforth, physical quantities with the superscript "AGN" are measured in the AGN frame, whereas quantities without the superscript are measured in the jet's comoving frame, unless specified otherwise.
Following the setup in Xue et al. (2019a), two emitting blobs in the jet, i.e., inner blob and outer blob, are considered in the modeling. We assume that these two blobs are spherical plasmoid of different radii (R in < R out ), moving with the same bulk Lorentz factor Γ j along the jet axis. Assuming the observer view the blob at an angle θ with respect to the jet axis, the Doppler factor of the two blobs are δ D = 1/Γ j (1 − β j cos θ) 1 . The two blobs are filled with uniformly entangled magnetic fields (B in > B out ), relativistic electrons and protons. Electrons and protons are assumed to be injected into both blobs with a broken power-law distribution for γ e,min < γ e < γ e,max , and a power-law distribution respectively. The free parameters for the spectral shape of electrons are the minimum, break, and maximum Lorentz factors γ e,min/b/max and the two spectral indices s e,1/2 , respectively. For the spectral shape of protons, the free parameters are the minimum, and maximum Lorentz factors γ p,min/max and the spectral index s p , respectively. The outer blob is assumed to be far away from the BLR, therefore the energy density of external photons in the outer blob is so weak that it can be neglected in the corresponding leptonic and hadronic processes. Therefore, the hadronic emission can be neglected and relativistic electrons radiate mainly through synchrotron radiation and synchrotron self-Compton (SSC) scatterings, giving rise to multi-wavelength emission. The inner blob is closer to the SMBH, so the radiation of the BLR and the dusty torus (DT) could enter the blob with a significant amplification of the energy density due to the Doppler effect. Those external radiation fields provide copious targets for hadronic interactions of protons as well as the external Compton (EC) scattering of electrons (see the upper panel of Fig. 2). The former process gives rise to neutrino emission and the latter process is relevant for gamma-ray emission. The energy density of BLR (u BLR ) and DT (u DT ) in the comoving frame as a function of the distance along the jet, r AGN in , can be approximately written as (Hayashida et al. 2012): and where φ BLR = φ DT = 0.1 are the fractions of the disk luminosity L d reprocessed into the BLR and DT radiation, respectively. It is assumed that the characteristic distance of BLR in the AGN frame is r AGN BLR = 0.1(L d /10 46 erg s −1 ) 1/2 pc and the characteristic distance of the DT is r AGN DT = 2.5(L d /10 46 erg s −1 ) 1/2 pc (e.g., . The radiation from both the BLR and DT is taken as an isotropic blackbody with a peak at 2 × 10 15 Γ j Hz  and 3 × 10 13 Γ j Hz (Cleary et al. 2007) in the comoving frame. In addition, since no gamma-ray activities are detected by Fermi-LAT for both 14-15 neutrino flare and IceCube-200107A, GeV gamma-ray photons must be absorbed, implying a compact size for the inner blob and an environment immersed in the intense X-ray radiation. Here, we consider that the inner blob is near or at the jet base with a distance comparable to only a few times larger than the Schwarzschild radius of the central SMBH (r Sch ∼ (M SMBH /5×10 8 M ⊙ )×10 14 cm). Then the X-ray photons emitted by the hypothetical corona surrounding the accretion disk (Heckman & Best 2014) would enter the inner blob and interact with the gamma-ray photons (Righi et al. 2019). The emission of the hot corona typically has a power-law spectrum (Kamraj et al. 2018;Ricci et al. 2018) and the corresponding energy density in the comoving frame can be estimated as Note that the size of the hot corona is generally only several times Schwarzschild radius, so the time that a relativistically moving blob need to pass through the hot corona, as measured by an observer, would be only ∼ r Sch /cΓ j δ D ∼ 10 2 s which is much shorter than the duration of the 14-15 neutrino flare or the considered emission period of IceCube-200107A in the previous literature. Therefore, multiple blobs need to form in the Figure 1. A schematic illustration (not to scale) of the inner-outer blob model. The outer blob (green circle) is far away from the BLR and the inner blob (yellow star) is near or at the jet base with a distance comparable to only a few times larger than the Schwarzschild radius of the SMBH. The external photons from hot corona would provide copious targets for hadronic interactions of protons as well as the EC scattering of electrons. On the other hand, the intensive X-ray photons from hot corona also absorb the GeV-TeV gamma rays and generate significant cascades emission, with energy around MeV band. The low-energy synchrotron radiation is suppressed due to the dominant Compton cooling of electrons. The observed multi-wavelength radiation arises from the synchrotron emission and SSC emission of electrons (e.syn and e.SSC) from the outer blob where the EC radiation is negligible.
hot corona region consecutively, or alternatively the inner blob could be related to certain dissipation process that occurs at a quasi-stationary feature such as a standing shock at the bottom of the jet (Fromm et al. 2011), so that the neutrino radiation in the inner blob can last a period comparable to the observation. As a consequence, the inner blob would appear as a quasi-stationary emitting zone, and the Doppler boosting is different from that for the outer blob. For the inner blob, we have F obs ν ∝ δ 2 D F in ν (e.g., Marscher 2009;Liu et al. 2020), while for the outer blob, there is another factor of δ D from the time compression since the entire emission region is moving towards us, then we have F obs ν ∝ δ 3 D F out ν . A sketch of the physical picture of our model is shown in Fig. 1.
To avoid introducing too many free parameters, the strategies for reducing the free parameters in Xue et al. (2019a) are applied here.
1. We assume the two blobs have the same spectral shape for the injected particles, the same electron injection luminosity L e,inj , and the same proton injection luminosity L p,inj .
2. We set the BLR luminosities L AGN BLR of TXS 0506+056 and 4FGL J0955.1+3551 are both 5 × 10 43 erg/s, as estimated by Padovani et al. (2019) and Paliya et al. (2020). Then the disk luminosity L d can be estimated through L d = 10 × L AGN BLR as proposed by .
3. Since the X-ray radiation from hot corona of blazars is not directly observed, here we simply as-sume that its spectral index α = 1, which is the observational typical value of AGN (Netzer 2013) and L 1keV is comparable to BLR luminosity. Thus, L 1keV is set as 5 × 10 43 erg/s for the two objects.
4. Since we already assume the external photons from BLR and DT are not relevant for the leptonic and hadronic processes, we do not specify the location of the outer blob r AGN out . For the inner blob, its distance from SMBH r AGN in is set to be a free parameter but restricted to be a few times of the Schwarzschild radius r Sch . 5. Since the inner blob is located at the jet base, the radius of the inner blob should be R in ≤ r AGN in . On the other hand, the external radiation field is the dominant targets for the hadronic interactions and EC scatterings, and its energy density is only related to position of inner blob r AGN in , therefore the specific value of R in will not significantly change the fitting results. For simplicity, we assume R in = r AGN in 2 .
6. For the minimum and maximum electron Lorentz factors, we set γ e,min = 50 and γ e,max = 10 7 , since our modeling results are not sensitive to them.
7. The minimum proton Lorentz factor γ p,min is fixed to be 1 and we set spectral index s p of proton energy distribution at injection to be 2. Upper panel: Timescales of various cooling processes for electrons (green curves) and protons (red curves) in the inner blob of TXS 0506+056 as a function of the particle energy. Both the particle energy and timescale are measured in the jet comoving frame. The parameters are the same as those used in Table 1. The black horizontal lines denote the escape timescale of particles in the ballistic propagation limit. The meaning of all curves is explained in the inset legend. Lower panel: Comparison of the opacity for γγ annihilation as a function of the photon energy in the observer's frame that contributed by the photons from synchrotron emission of primary electrons (dot-dashed green curve), hot corona (solid red curve), broad-line region (dashed purple curve) and dusty torus (dotted yellow curve), respectively.
8. γ p,max is obtained by equating the acceleration and the cooling or escape time-scales t acc = min{t pγ , t BH , t esc , t p,syn }.
If we assume the particle acceleration is dominated by Fermi-type acceleration (e.g., Rieger et al. 2007), the acceleration time-scale t acc can be expressed as t acc = ηγ p m p c/eB in = 2.6 × 10 −3 γ p ( η 50 )( Bin 2 G ) −1 , where η is an efficiency factor characterizing the acceleration rate. t pγ ≃ [cn soft < σ pγ κ pγ >] −1 is the pγ energy loss time-scale, where n soft is the number density of the soft photons and < σ pγ κ pγ >≃ 10 −28 cm 2 is the inelasticity-weighted pγ interaction crosssection. t BH ≃ [cn soft < σ BH κ BH >] −1 is the BH pair-production cooling time-scale, where σ BH and κ BH are the cross section and inelasticity for the BH pair production process (Chodorowski et al. 1992).
t esc = R in /c is the escape timescale. t p,syn = 6πm e c 2 /(cσ T B 2 γ p )(m p /m e ) 3 is the proton-synchrotron cooling timescale. It can be seen from the upper panel of Fig. 2 that due to the very high energy density of the external photons (especially the photons from hot corona), the pγ energy loss above tens TeV is very efficient, dominating over escape. Therefore, under adopted parameters in Fig. 2, γ p,max is determined by the equality t acc = t pγ .
At last, the number of free parameters is reduced to eleven, namely δ D , B in , B out , R in , R out , L e,inj , L p,inj , s e,1 , s e,2 , γ e,b , and η.
With the above setups, the synchrotron, SSC and EC photons from primary electrons, and external photons from hot corona, BLR and DT are all considered as the targets for pγ interactions and γγ pair-production in the inner blob. The photopion production and Bethe-Heitler (BH) pair-production processes are calculated following Kelner & Aharonian (2008). Triggered by the photons from neutral pion decay, secondary electron/positron pairs from charged pion decay, and BH pair-production, the Compton-supported cascades in the inner blob are evaluated using a semi-analytical method developed by Böttcher et al. (2013) and Wang et al. (2018b). Since the very high-energy gamma-ray photons escape from the inner or the outer blobs will be absorbed by the extragalactic background light (EBL), we calculate the attenuation in the GeV-TeV band using the EBL model of Domínguez et al. (2011).

RESULTS
Observationally, during the 14-15 neutrino flare, no EM activity in any wavelength is observed from TXS 0506+056. This feature constrains that the EM emission from neutrino production region can not contribute significantly to the observed energy band. The UV/X-ray emission of 4FGL J0955.1+3551 is found in the high-state at the time of the arrival of IceCube-200107A. However, no activity is found in other wavelengths. Besides, if the neutrino emission period is assumed to be 250 days or 10 years (Giommi et al. 2020a), the high-state X-ray emission may not be causally linked to the neutrino emission. We therefore speculate that IceCube-200107A is produced in a similar process to that for the 14-15 neutrino flare of TXS 0506+056. In the modeling, we find that there are many sets of parameters may satisfactorily reproduce the multi-messenger observations of TXS 0506+056 and 4FGL J0955.1+3551. Here we firstly show the fitting results with a specific set of pa- Note-The kinetic luminosity in relativistic electrons L k e,in and in relativistic protons L k p,in in the inner blob are calculated as L k e,in ≃ πR 2 in δ 2 D mec 3 γeNe(γe)dγe and L k p,in ≃ πR 2 in δ 2 D mpc 3 γpNp(γp)dγp, where Ne(γe) and Np(γp) are the steady-state electrons and protons energy distributions, respectively.  Fig. 3. The corresponding parameters are shown in Table 1 and Table 2. From Fig. 3, it can be seen that, as expected, the multiwavelength SEDs with measurements are dominated by the synchrotron and SSC emission emitted by the primary electrons from the outer blob, which do not emit high-energy neutrinos, while the neutrinos originate from the inner blob, accompanied by a comparable flux in the MeV band. The high MeV flux is a consequence of the EM cascades developed in the intense X-ray radiation field from the hot corona (see also the lower panel of Fig. 2): given that the X-ray spectrum of the hot corona can extend up to ∼ 100 keV, it can absorb γ-ray photons starting from E γ,th ≈ 2(mec 2 ) 2 100 keV(1+z ) 2 ≃ 2 MeV in the observer's frame with an opacity at 2 MeV being τ γγ = σγγ L 100 keV 4πr AGN in c100 keV ≃ 0.5. The high-energy gammaray photons and electron/positron pairs generated in the pγ interactions will initiate EM cascades and generate lower and lower energy photons or electrons/positrons with the development of the cascade, until the energy of the newly generated photons drop close to E γ,th . These photons will be absorbed by the X-ray radiation field and give birth to the last-generation electron/positron pairs of the EM cascade, with energy E γ,th /2 1 MeV. Therefore, most energies of the first-generation gamma rays and electron/positron pairs that co-produced with neutrinos will deposit to the pairs of 1 MeV, which mainly radiate in the MeV band via the IC scattering off the X-ray photons of the hot corona because the energy density of hot corona dominate over the energy densities of magnetic field and other photon fields 3 (see also Figure 10 of Reimer et al. 2019). Therefore, in the absence of measurements at MeV band, it seems as if the blazars underwent "orphan" neutrino flares.
From the upper panel of Fig. 3, we can see that the multi-wavelength SED of TXS 0506+056 is well fitted by the leptonic emission from outer blob except the highest energy GeV data point which shows a spectral hardening, although a hardening feature may not be significant (Garrappa et al. 2019). This is because the inner blob has a large gamma-ray opacity while the outer blob, which is far away from the BLR and DT radiation field, can hardly produce a hard spectrum extending up to 10 GeV (or above) solely with the SSC process. This is the limiation of the present model. However, in the physical picture of the inner-outer model, there could be additional dissipative regions being formed between the inner and outer blobs (i.e., at a distance ∼ r AGN BLR from the SMBH), which can in principle reproduce the highest-energy GeV data point with the EC process. We refrain ourselves from the further discussion on the interpretation of the last GeV data point and focus back on the neutrino flare. We can see that the inferred neutrino flux during the 14-15 flare of TXS 0506+056 can be well reproduced by our model. For 4FGL J0955.1+3551, the predicted neutrino flux is much higher that of onezone model (Paliya et al. 2020). This is because of the differences between the one-zone model and the innerouter blob model. As indicated in Xue et al. (2019a), by ascribing the low-energy emission (radio to UV band) 3 In the modeling, r AGN in /r AGN BLR ≈ 4 × 10 −3 , therefore ucorona/u BLR ≈ (r AGN in /r AGN BLR ) 2 ∼ 10 5 . The ratio of ucorona to u DT can be estimated in the same way, which is ∼ 10 8 . Upper panel: The multi-messenger emissions of TXS 0506+056 predicted by the two-zone model. The quasi-simultaneous data is taken from Padovani et al. (2018) and the violet bow tie shows the muon-neutrino flux as measured during 14-15 flare (IceCube Collaboration et al. 2018a). Lower panel: The multi-messenger emissions of 4FGL J0955.1+3551 predicted by the two-zone model. The quasi-simultaneous data is taken from Paliya et al. (2020), and the violet dotted, dashed and solid horizontal lines show the ν µ +ν µ neutrino flux for 30 days, 250 days and 10 years (Giommi et al. 2020a). The line styles in both panels have the same meaning. The dashed red and dot-dashed red curves represent the synchrotron and SSC emission from primary electrons in the outer blob, respectively. The dotted green curve represents the emission from pair cascades in the inner blob. The double dot-dashed green curve shows the muon neutrino energy spectrum. The solid black curve is the total emission from both blobs.
of the blazar to the outer blob, synchrotron emissivity of primary electrons in the inner blob could be much lower than that in the one-zone models, especially in the presence of dense external photon fields. As a result, the γγ absorption opacity for sub-PeV/PeV photons is significantly reduced and less than unity (see the lower panel of Fig.2) in the inner blob. This sup-presses the cascade emission by reducing the number of secondary electrons generated in γγ absorption. Also, the secondary electrons generated in the cascade preferentially emit MeV-GeV gamma rays through the IC scattering in the dense external photon fields rather than in the X-ray band through the synchrotron radiation, and hence the restrictive constraints by the X-ray observation is avoided. In our modeling, the predicted neutrino flux is well above the one assuming an emission period of 10 years given by Giommi et al. (2020a). As shown in Fig. 2, since the energy density of the hot corona is very high, the pγ interaction is very efficient (f pγ = t esc /t pγ ≈ 3 for the tens TeV to PeV protons). Therefore, the required proton injection luminosity in the inner-outer blob model is moderate and reasonable from the point of view of energy budget. In our modeling, the baryon loading factors, defined as L p,inj /L e,inj , of TXS 0506+056 and 4FGL J0955.1+3551 are 125 and 600, respectively, which are significant lower than that of the conventional one-zone proton synchrotron model (e.g., Böttcher et al. 2013) and one-zone pγ model (e.g., Xue et al. 2019b). Provided the mass of the SMBH being 3 × 10 8 M ⊙  , the jet's kinetic luminosity employed here is also below the Eddington limit.
To further explore the general applicability of the inner-outer blob model in explaining the blazar "orphan" neutrino flare, we perform an extensive scan of the parameter space for the two sources separately. Since both these two emissions are relevant with the inner blob, we only search the parameter space of the inner blob. In addition, leptonic emissions from the inner blob do not contribute to the observed SED, and hence is neglected here. We scan the parameter space of five parameters for the properties of the inner blob and injected relativistic protons energy distribution within physically plausible ranges, which are 10 δ D 30, 3 × 10 14 r AGN in /cm 3 × 10 15 , 1 B in /G 100, 10 η 100, and 10 45 δ 2 D L p,inj /erg s −1 10 47 . The former two parameters are spaced into five bins linearly while the latter three are spaced into five bins logarithmically. Since the observed SEDs are ascribed to the emissions from the outer blob, the flux of the cascade emission should be limited below 1 σ upward error bar of the data points. Generally speaking, a larger δ D and a smaller r AGN in would increase the energy density of the external photon field and improve the pγ interaction efficiency f pγ . A smaller B in would reduce the predicted Xray flux through the synchrotron radiation of secondary pairs co-produced with neutrinos, making the result less constrained by the X-ray data and leaving more room for the neutrino production.
The scanning results are shown in Fig. 4. For the 14-15 neutrino flare of TXS 0506+056, since the observa-  Fig. 3. The yellow/red solid bands represent the muon neutrino energy spectra with different parameter sets in the inner blob, which can/cannot produce a sufficient flux to explain the neutrino measurements, while the cyan/blue solid bands represent the corresponding emission from pair cascades. The parameter resulting in the highest (lowest) neutrino flux for TXS 0506+056 are δ D = 30(10), r AGN in = 3 × 10 14 (10 15 ) cm, B in = 1(100) G, η = 100(10), and δ 2 D L p,inj = 10 47 (10 45 ) erg/s; the parameters resulting in the highest (lowest) neutrino flux for 4FGL J0955.1+3551 are δ D = 30(10), r AGN in = 3 × 10 14 (10 15 ) cm, B in = 1(100) G, η = 100(10), and δ 2 D L p,inj = 10 46.5 (10 45 ) erg/s. tionally inferred neutrino events is 13 ± 5 during about 5 months, we show the neutrino flux under parameters that are able to produce at least 8 events in 10 TeV -10 PeV by yellow band while the red band show that unable to produce more than 8 events. The corresponding EM fluxes are shown in cyan and blue respectively. The result for 4FGL J0955.1+3551 is exhibited in a sim-ilar way except coloring the bands based on whether the neutrino flux can give rise to 1 event detection by Ice-Cube (yellow and cyan) or not (red and blue) in the range of 65 TeV -2.6 PeV for an emission period of 250 days. The IceCube point-source effective area for (anti-)muon neutrinos (Carver 2019) in the declination −5 • − 30 • and 30 • − 90 • are employed for TXS 0506+056 and 4FGL J0955.1+3551 respectively to calculate the expected neutrino event 4 . It can be seen that there are amount of parameter sets able to explain the 14-15 neutrino flare from TXS 0506+056 and IceCube-200107A from 4FGL J0955.1+3551. For the latter, the most optimistic parameters can result in a non-negligible probability of 20% for detection of one event by IceCube even assuming an emission period of 30 days. One robust and remarkable feature in the resulting SED is the MeV flux. As is shown in Fig. 4, when the observed neutrino energy spectrum can be interpreted under certain parameter set (the yellow band), the corresponding cascade emission will contribute significant radiation in the MeV band (the cyan band). It should be also noted that a spectral bump in the MeV band is not the unique feature predicted in our model. In some other lepto-hadronic models of the 14-15 neutrino flare of TXS 0506+056, an intense MeV emission is also expected provided the presence of a dense external radiation field in the UV to X-ray band (e.g. Wang et al. 2018a;Reimer et al. 2019;Rodrigues et al. 2019;Petropoulou et al. 2020a). On the other hand, the flux level of the MeV emission is generally comparable to the neutrino flux and might be used as an indicator of neutrino flux and distinguish different models. Currently, there is no active instruments able to measure the accompanying MeV emission during the neutrino emission period, but the next generation MeV gamma-ray detector (e.g., Wu et al. 2016;Moiseev, & Amego Team 2017;de Angelis et al. 2018) would be sensitive enough to measure this emission and serve as a critical test to our model.

DISCUSSION AND CONCLUSIONS
In Section 3, we have shown that the inner-blob model may in principle explain the 14-15 neutrino flare from TXS 0506+056 and IceCube-200107A from 4FGL J0955.1+3551, provided that the inner blob occurs in the X-ray corona around the SMBH. One may then raise a question: how much could such a phenomenon contribute to the all-sky diffuse neutrino background? The answer highly depends on how frequently such a phenomenon could take place in a blazar, and is actually unknown without further understanding on its formation mechanism. Let us parameterize the occurrence probability of certain dissipative process in the X-ray corona of a blazar's SMBH to be f dis , the diffuse neutrino flux from the entire population of BL Lacs can be estimated as (Tavecchio et al. 2014) where E ν is the observed neutrino energy, L γ is the observed γ-ray luminosity in the range 10 44 < L γ /erg s −1 < 10 46 and j[L γ , E(1 + z), z] = Σ(L γ , z) Lν (Eν ) Eν is the luminosity-dependent comoving volume neutrino emissivity. Σ(L γ , z) is derived using the luminosity function and parameters for its luminosity-dependent evolution for BL Lacs provided by Ajello et al. (2014). We assume the neutrino luminosity L ν and the gamma-ray luminosity L γ scale as L ν = ξL γ , where L ν is the differential luminosity of muon and anti-muon neutrino predicted by our model for the 14-15 neutrino flare of TXS 0506+056 5 as shown in the upper panel of Fig. 3. Comparing the measured gamma-ray flux and the predicted neutrino flux we have ξ = 2.5. We then find that f dis = 2.6% is needed in order to explain the diffuse neutrino flux measured by IceCube (Haack et al. 2017) as shown in Fig. 5. It would immediately rule out the possibility of such kinds of blazar "orphan" neutrino flares as the main source of diffuse neutrino flux, because the inferred f dis would lead to a very low source density of < 10 −9 Mpc −3 and hence violates the constraint from non-detection of neutrino multiplets in the diffuse neutrino events (Ahlers & Halzen 2014;Murase & Waxman 2016). It suggests that this kind of blazar "orphan" neutrino flare must be very rare, implying that dissipations at the jet base cannot occur often or sustain for a very long time. If a dissipative process occurs at an equal probability per unit length along the jet, the probability of forming a blob within the X-ray coronal (r 10 −4 pc) would be f dis ∼ 0.001% for a (the inner, highly relativistic) jet length of ∼10 pc, resulting in a subdominant contribution to the diffuse neutrino background (see Fig. 5). Nevertheless, the X-ray corona of SMBH could be a promising neutrino production sites provided that protons are accelerated in the region. For example, Inoue et al. (2019) and Murase et al. (2020) suggest that protons could be accelerated via plasma turbulence or magnetic reconnection in the corona of AGN and provide a steady neutrino emission. Different from our model, the neutrino emission proposed by these authors is not beamed so the neutrino flux of each source is low but Contribution of the blazar "orphan" neutrino flare to the all-sky diffuse neutrino background. Green circles represent the measured intensity of the high-energy ν µ +ν µ taken from Haack et al. (2017) and green triangles indicate the upper limits. The solid curve represents the case with f dis = 2.6% following the contrast between the gamma-ray flux and neutrino flux of TXS 0506+056. The dashed curve shows the resulting diffuse neutrino flux with f dis = 0.001%. See Section 4 for detailed discussion.
To summarize, we studied the multi-messenger emissions of blazar TXS 0506+056 in the period of the 2014-2015 neutrino flare and blazar 4FGL J0955.1+3551 during the conjectural emission period of the neutrino event IceCube-200107A with the inner-outer blob model. Considering that certain dissipative process occurs at the jet base, an inner blob may form within the hot corona around the SMBH and produce a neutrino flare via the pγ interactions with the X-ray radiation of the corona, while the non-flare-state multi-wavelength flux of the blazar measured in the same period is ascribed to an outer blob which is far from the SMBH. We found that, on the premise that the multi-wavelength SED being generally reproduced, the 14-15 neutrino flare of TXS 0506+056 could be interpreted in a wide range of parameter space, showing the feasibility that the multi-messenger emission of TXS 0506+056 in the period of the 14-15 neutrino flare can be interpreted in the same physical picture of that for IceCube-170922A. The same model could be also applied to the neutrino event IceCube-200107A and explain its association of another blazar 4FGL J0955.1+3551. By searching the parameter space, we found that the model could well reproduce the neutrino emission of these blazars with an amount of parameter sets. A robust fea-ture predicted by our inner-outer blob model is a simultaneous MeV gamma-ray flare generated by the EM cascade emission with the comparable flux level to that of neutrinos, in agreement with some other models for the 14-15 neutrino flare of TXS 0506+056. In the future, sensitive MeV gamma-ray instrument may be able to catch the MeV flare around the arrival time of a neutrino event from a blazar and serve as a critical test to the inner-outer blob model for the blazar's "orphan" neutrino flares with simultaneous multi-wavelength observation.