Draft version August 30, 2021 Typeset using LATEX twocolumn style in AASTeX62
High-Energy Neutrinos from Magnetized Coronae of Active Galactic Nuclei and Prospects for Identification of Seyfert Galaxies and Quasars in Neutrino Telescopes
Ali Kheirandish,1, ∗ Kohta Murase,2, † and Shigeo S. Kimura3, ‡
1Department of Physics; Department of Astronomy & Astrophysics; Center for Multimessenger Astrophysics, Institute for Gravitation and the Cosmos, The Pennsylvania State University, University Park, PA 16802, USA 2Department of Physics; Department of Astronomy & Astrophysics; Center for Multimessenger Astrophysics, Institute for Gravitation and the Cosmos, The Pennsylvania State University, University Park, PA 16802, USA & Center for Gravitational Physics, Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto, Kyoto 606-8502, Japan 3Frontier Research Institute for Interdisciplinary Sciences, Tohoku University, Sendai 980-8578, Japan & Astronomical Institute, Tohoku University, Sendai 980-8578, Japan
ABSTRACT Particles may be accelerated in magnetized coronae via magnetic reconnections and/or plasma turbu- lence, leading to high-energy neutrinos and soft gamma rays. We evaluate the detectability of neutrinos from nearby bright Seyfert galaxies identified by X-ray measurements. In the disk-corona model, we find that NGC 1068 is the most promising Seyfert galaxy in the Northern sky, where IceCube is the most sensitive, and show prospects for the identification of aggregated neutrino signals from Seyfert galaxies bright in X-rays. Moreover, we demonstrate that nearby Seyfert galaxies are promising tar- gets for the next generation of neutrino telescopes such as KM3NeT and IceCube-Gen2. For KM3NeT, Cen A can be the most promising source in the Southern sky if a significant fraction of the observed X-rays come from the corona, and it could be identified in few years of KM3NeT operation. Our results reinforce the idea that hidden cores of supermassive black holes are the dominant sources of the high-energy neutrino emission and underlines the necessity of better sensitivity to medium-energy ranges in future neutrino detectors for identifying the origin of high-energy cosmic neutrinos.
Keywords: astroparticle physics – galaxies: active – galaxies: jets – gamma rays: galaxies – neutrinos – radiation mechanisms: non-thermal
1. INTRODUCTION identified, the multimessenger data indicate that their The observation of high-energy cosmic neutrinos in sources should be opaque to high-energy γ-rays in the IceCube (Aartsen et al. 2013a,b) has revealed that GeV-TeV range, i.e., hidden to the γ-ray telescopes op- hadronic interactions play a major role in the non- erating at these energies (Murase et al. 2016; Capanema thermal emission in the high-energy universe. While the et al. 2020a,b). Otherwise, the contribution of the γ-ray magnitude of the measured spectrum at high energies, flux leaving transparent sources to the diffuse isotropic γ-ray background (IGRB) would create an excess that & 100 TeV, was found to be consistent with some the- oretical predictions for astrophysical neutrinos (Loeb & would overshoot the measured diffuse flux by Fermi- arXiv:2102.04475v2 [astro-ph.HE] 26 Aug 2021 Waxman 2006; Murase et al. 2008; Kotera et al. 2009; LAT (Ackermann et al. 2015). Consequently, scenarios Murase et al. 2013), the emergence of an order of magni- such as galaxy clusters/groups and starburst galaxies tude higher flux at medium energies in the 10 TeV range as the dominant contributors to the high-energy cos- (Aartsen et al. 2015a,b, 2020a) pointed towards sources mic neutrino flux below 100 TeV are excluded because with efficient neutrino emission. While the dominant they will overwhelm the observed IGRB flux (Sen 2015; sources of high-energy cosmic neutrinos are yet to be Murase et al. 2016; Bechtol et al. 2017). On the other hand, the flux above 100 TeV can still be explained by such cosmic-ray (CR) reservoir scenarios (Murase et al. ∗ [email protected] 2013; Xiao et al. 2016; Murase & Waxman 2016; Fang † [email protected] & Murase 2018; Liu et al. 2018; Peretti et al. 2020). ‡ [email protected] 2
Cores of the active galactic nuclei (AGN), which are emission in Fermi cannot accommodate the reported optically thick for GeV-TeV γ-rays are one of the best neutrino flux. In the meantime, neutrino flux predic- candidates as the source of the high-energy neutrino flux tions from starburst activity are directly correlated with at medium energies. The isotropic distribution of the the Fermi γ-rays and face a similar obstacle. arrival direction of high-energy neutrinos, together with In this paper, we present the expected neutrino emis- the constraints on the Galactic component of the neu- sion from the brightest Seyfert galaxies in the disk- trino flux, points at extragalactic sources hidden in γ- corona model by incorporating the measured intrinsic rays as the origin of this flux. Moreover, in order to X-ray luminosities. We estimate the parameters relevant produce this level of flux, efficient neutrino production for emission from the magnetized coronae and establish is required. Such conditions can be met with large col- the parameter space compatible with the reported neu- umn densities of target material and radiation, which trino spectrum of NGC 1068 and the medium-energy form dense environments suppressing the γ-ray flux at excess in the cascade data. We further employ the pa- the site of production. rameters that can accommodate the emission from NGC In Seyfert galaxies, accretion dynamics and magnetic 1068 for the brightest Seyfert galaxies and study their dissipation will form a magnetized corona above the likelihood of observation, individually and when stacked disk, see e.g., Miller & Stone(2000); Merloni & Fabian together. (2001); Liu et al.(2002a); Blackman & Pessah(2009); Io While providing a consistent level of neutrino flux, our & Suzuki(2014); Jiang et al.(2014); Jiang et al.(2019). modeling finds NGC 1068 as the most promising source Thanks to the dense environments near the supermas- in the Northern hemisphere, where IceCube is most sen- sive black holes, the accelerated CRs in the corona in- sitive. While the majority of the bright sources reside teract efficiently with the gas and target photons, pro- in the Southern sky, we show that the stacking analy- ducing charged and neutral pions that promptly decay sis has a good chance of identifying neutrino emission to neutrinos and γ-rays. The neutrinos will escape the from this class of sources. In addition, the significance environment while GeV-TeV γ-rays are cascaded down of the search for neutrino emission from such analysis inside the source. can inform the principal mechanism for the high-energy The disk-corona model for high-energy neutrino emis- emission. sion from the core of AGNs can successfully accommo- We use the predicted neutrino emission to examine the date the flux of cosmic neutrinos at medium energies in prospects for observation of the nearby, bright Seyfert the 10 100 TeV range (Murase et al. 2020). In ad- galaxies in future neutrino telescopes, KM3NeT and the − dition, modeling of the high-energy emission from the next generation of IceCube at the South Pole: IceCube- hidden cores of AGNs presents a robust connection to Gen2. γ-ray emission at MeV energies. While indicating that Seyfert galaxies are the dom- The dominant role of AGN cores to the flux of high- inant contributors to the high magnitude flux of cos- energy neutrinos at medium energies has been reinforced mic neutrinos at energies below 100 TeV, our prediction by the IceCube’s most recent search for the sources of provides a testable scenario in current and future neu- high-energy neutrinos using 10 years of data (Aartsen trino telescopes. The subset of bright Seyfert galaxies is et al. 2020b). The all-sky (untriggered) and the source likely to be observed with the continuous operation of list (triggered) time-integrated analysis of the data finds the IceCube detector and commissioning of the neutrino the direction of NGC 1068 as the most significant point telescopes in the Northern hemisphere. in the sky. The excess in the direction of NGC 1068 is In the next section, we discuss details of acceleration found local pre-trial p-value of 1.8 10−5, yielding a mechanism and neutrino production in the magnetized × 2.9σ after correcting for trials. corona model. Then, in Sec.3 we present the neu- NGC 1068, aka Messier 77, is a nearby, bright Seyfert trino emission under 3 distinct scenarios from NGC 2 galaxy. The observed infrared luminosity of NGC 1068 1068 and the brightest Seyfert galaxies identified by the is similar to starburst galaxies. It has been considered X-ray surveys. Subsequently, in Sec.4, we will study as a potential CR accelerator and motivated study of the prospects for identification of the neutrino emis- neutrino and γ-ray emission, see , e.g., Yoast-Hull et al. sion in current and future neutrino telescopes. Finally, (2014); Murase & Waxman(2016); Liu et al.(2018); we will discuss the implications of observing neutrinos Lamastra et al.(2019). The best-fit neutrino flux mea- from these sources and the relation to the total flux in sured by IceCube, however, exceeds the observed γ-ray Secs. 5.2. emission by Fermi-LAT (Ackermann et al. 2012). There- fore, the models that are built upon the measured γ-ray High-Energy Neutrinos from Seyfert Galaxies 3
2. MAGNETIZED CORONA MODEL ε , L is the X-ray luminosity, is the volume of the inj X V Progress in X-ray observations of AGN led to the corona, and δ(x) is the delta function. Note that in establishment of the magnetized coronal paradigm this scenario, the injection fraction is proportional to (Haardt & Maraschi 1991; Liu et al. 2002b; Ricci et al. the dissipation rate in the corona, and therefore, is pro- 2018). This is partly confirmed by recent magneto- portional to LX . The diffusion coefficient is assumed to scale as D εq , where q is the power-law index of hydrodynamic (MHD) simulations (e.g., Io & Suzuki εp ∝ p 2014; Jiang et al. 2019). The magnetorotational insta- the diffusion coefficient in momentum space. If the gyro- bility (MRI; Balbus & Hawley 1991) supports that the resonant scattering is dominant, which is not necessarily coronae above the optically thick accretion disk are nat- the case, it corresponds to the spectral index of the tur- urally expected to be hot, magnetized, and turbulent. bulence power spectrum. For demonstrative purposes, Such a plasma can be collisionless in the sense that the we use q = 5/3 motivated by the Kolmogorov turbu- accretion onto a central black hole is faster than their lence, but q = 3/2 (that can be expected for isotropic Coulomb collisions time, and ions may be accelerated via MHD turbulence in the magnetically dominated regime) stochastic acceleration (e.g., Lynn et al. 2014; Kimura or q = 2 (that is a hard-sphere type; see, e.g., Kimura et al. 2016; Comisso & Sironi 2018; Zhdankin et al. et al. 2019b) may be possible. Note that in general the 2019; Kimura et al. 2019b; Wong et al. 2020; Lemoine & CR spectrum can be a power law if the trappiing is Malkov 2020) and/or magnetic reconnection processes significant (Lemoine & Malkov 2020). As cooling pro- (e.g., Hoshino 2013, 2015; Guo et al. 2016; Werner et al. cesses, we take into account Bethe-Heitler, pp inelastic 2018; Ball et al. 2018b). Stcochastic acceleration and collisions, and photomeson (pγ) production processes. magnetic reconnection may naturally coexist in highly We consider advective (infall to a central black hole) magnetized plasmas. Although either may become a and diffusive escapes. The maximum energy of the CRs dominant channel for particle acceleration (Comisso & is set by the balance between the acceleration and these Sironi 2019; Zhdankin et al. 2019; Li et al. 2021), they losses. For the latter, the principal component is the cannot easily be separated. In this study, for demon- Bethe-Heitler process due to abundant disk photons, as stration, we consider particle acceleration and neutrino emphasized in Murase et al.(2020). This results in a production under each scenario separately. cutoff feature in the CR spectrum below PeV energies. In this scenario, we introduce two parameters, the CR 2.1. Stochastic acceleration scenario pressure to the thermal pressure PCR/Pth and turbulent strength η−1. The former determines the normalization It is widely believed that turbulent magnetic fields tur of CRs, whereas the latter is related to the maximum generated by the MRI is responsible for the angular mo- energy of CRs. Note that we use the thermal pressure at mentum transport in accretion flows, which has been the virial temperature, i.e., P = n kT , where n is confirmed by MHD simulations (Stone & Pringle 2001; th p vir p the thermal proton density and T is the virial temper- Machida & Matsumoto 2003; Ohsuga & Mineshige 2011; vir ature, and the realistic thermal pressure can be reduced Narayan et al. 2012; Porth et al. 2019). Turbulence can if the ion temperature is lower than T . See Supple- also be generated through magnetic reconnections re- vir mental Material of Murase et al.(2020) for details (see sulting from magnetic fields amplified by the MRI. CRs also Section III-B of Kimura et al. 2019a, in the context in the MHD turbulence randomly change their energies of radiatively inefficient accretion flows). through interactions with MHD waves. If the acceler- ation in momentum space homogeneously occurs with- out significant trapping (Lemoine & Malkov 2020), the 2.2. Magnetic reconnection (fast acceleration) scenario phenomenon can be described by the following diffusion Recent Particle-In-Cell (PIC) Simulations revealed equation in energy space: that relativistic magnetic reconnections in the ion- electron plasma, where the ion magnetization param- 3 ! ∂ 1 ∂ ∂ ε 2 2 p 2 p p p ˙ eter σi = B /(8πnpmpc ) & 1, can accelerate non- F = 2 εpDεp F + p F + p,inj, ∂t εp ∂εp ∂εp tcool F − tesc F thermal relativistic particles very efficiently, given that (1) the plasma β, which is defined by β 8πn kT /B2, ≡ p p where is the CR distribution function (dN/dε = is smaller than 0.01 (Guo et al. 2016; Werner Fp p ∼ 4πp2 /c), D is the diffusion coefficient in energy et al. 2018; Ball et al. 2018b). For AGN coronae, Fp εp space, t is the cooling timescale, and t is the β can be sufficiently small (e.g., β 1 3) but we cool esc . − escape timescale. The injection function is given by typically expect that the plasma is semi-relativistic, ˙ = f L δ(ε ε )/[4π(ε /c)3 ], where ε is σ = 2kT /(m c2β) < 1 (Werner et al. 2018). It is still Fp,inj inj X p − inj inj V inj i p p the injection energy and finj is the injection fraction at uncertain whether protons can be accelerated above 4
m c2, and further acceleration via turbulence would where η is the acceleration efficiency parameter and ∼ p acc be necessary (Comisso & Sironi 2018; Zhdankin et al. rL is the Larmor radius. At sufficient high energies 2019; Wong et al. 2020; Lemoine & Malkov 2020). Not the Fermi-like acceleration, and ηacc can be as small only the stochastic acceleration but also the first-order as 3(c/V )(Giannios 2010; Zhang et al. 2021), al- ∼ rec Fermi acceleration may operate due to scatterings with though details are uncertain in non-relativistic recon- reconnection outflows (e.g., Pisokas et al. 2018; Hoshino nections (c.f. Hoshino 2012; del Valle et al. 2016). In 2012). In this work, for the purpose of phenomeno- this work, we assume ηacc = 300, but we will see that max rec logical studies, we simply assume that the injection Ep = Ep to better explain the observation of NGC spectrum of protons is given by a single power-law with 1068 (see Section 3.1). The normalization of CRs are R 2 an exponential cutoff, providing the power-law index, provided such that N˙ injεpdεp = CRMc˙ , where CR is s, as a parameter. We assume that the corona can the energy fraction carried by CRs and M˙ is the mass be regarded as electron-proton plasma. X-ray obser- accretion rate in the accretion disk. Only the frac- vational data seem consistent with the electron-proton tion of the accretion energy is spent to accelerate CRs. plasma (Ricci et al. 2018), although better quality data Therefore, we use the total mass accretion rate onto the are necessary to clarify it. For electron-positron-proton BH to deduce the CR normalization. The total mass plasma, the non-thermal proton production efficiency accretion rate is proportional to the bolometric luminos- via magnetic reconnections is expected to be lower than ity, which is estimated from intrinsic X-ray luminosity. that for electron-proton plasma (e.g., Petropoulou et al. Note that the corona is connected to the accretion disk 2019). through magnetic fields, and a significant fraction of the To calculate neutrino emission, we solve the transport accretion energy may be dissipated in the disk-corona equation for the non-thermal protons with single-zone interface. and steady-state assumptions: The previous works focused on the stochastic acceler- ation (Murase et al. 2020) or shock acceleration mech- ε N N d p εp ˙ εp anism (see below), and the magnetic reconnection sce- = Nεp,inj , (2) −dεp tcool − tesc nario was not considered in detail. In this work, we show that the diffuse neutrino flux can be explained by where N = dN/dε , N˙ ε−s is the injection εp p εp,inj ∝ p the magnetic reconnection scenario in the AGN corona term. model. In Sec.6 we discuss the possibility of explaining We consider the same cooling and escape processes the diffuse neutrino flux at medium energies with this as those in the stochastic acceleration scenario. The scenario and comment on the feasibility of identifying maximum proton energy of the injection spectrum is this scenario in the near future. determined by either the size of the acceleration re- gion or the balance between acceleration and the to- max 2.3. Remarks on other acceleration scenarios tal losses via cooling and escape processes, as Ep = rec cool min[Ep ,Ep ]. In general, particle acceleration con- This work focuses on the standard magnetized corona sists of multi-phase processes including the electric field scenario based on recent global MHD simulations, which acceleration and Fermi acceleration (Li et al. 2021). In does not involve accretion shocks. Alternatively, accre- the absence of losses, the former energy can be writ- tion shocks have been discussed as a CR production site ten as Erec eBl , where l is the reconnection in the accretion flows (Begelman et al. 1990; Stecker p ∼ rec rec region thickness that can be smaller than the system et al. 1991; Inoue et al. 2019). Originally, this scenario size (Hoshino 2012; del Valle et al. 2016). The latter is was proposed to explain X-ray observations of AGNs determined by tacc = tcool, where tacc is the proton ac- by CR-induced electromagnetic cascades (Berezinskii & celeration time and tcool is the proton cooling time. The Ginzburg 1981; Zdziarski 1986). However, cutoff and acceleration time depends on details of the magnetic softening features were discovered in bright AGNs in reconnection process. For example, particle bounces 1990s (Maisack et al. 1993; Zdziarski et al. 2000), which may occur either in converging magnetic fluxes with re- ruled out the hadron-induced electromagnetic cascade connection velocity (de Gouveia dal Pino & Lazarian scenario for the observed X-rays. The accretion shock 2005; del Valle et al. 2016; Medina-Torrejon et al. 2021) could still be viable to be a neutrino production site and/or by outflows with Alfv´envelocity (Drake et al. as long as the CR-induced cascade flux is below the 2006; Hoshino 2012). However, recent simulations have X-ray and gamma-ray data (Stecker 2005, 2013; Inoue shown that the acceleration is efficient especially for et al. 2019). However, there are several problems in relativistic reconnections (Hoshino & Lyubarsky 2012). this scenario. First, the accretion shocks are not seen For demonstrative purposes, we assume tacc = ηaccrL/c, in any global MHD simulations dedicated for accretion High-Energy Neutrinos from Seyfert Galaxies 5
flows (e.g., Narayan et al. 2012; Kimura et al. 2014; The high-energy neutrino flux modeled in this study Jiang et al. 2019). Although there are solutions with incorporates both hadronuclear (pp) and photohadronic an accretion shock in the one-dimensional hydrodynamic (pγ) interactions. The inclusion of both the processes in equation system with a steady state assumption (Becker the production of high-energy neutrinos, together with et al. 2011; Chattopadhyay & Kumar 2016), such so- the cooling processes in the source, results in features lutions are not realized in these dedicated simulations. in the neutrino flux from Seyfert galaxies. Thus, the The second point is the angular momentum of the ac- spectrum of high-energy neutrinos would deviate from creting matter. For example, Inoue et al.(2019) explic- simplistic single power laws. itly assumed the free-fall spherical accretion (see their The stochastic acceleration scenario for the produc- Equation 2). However,. we expect a large specific an- tion of high-energy neutrinos may yield a spectrum that gular momentum at an outer scale, which prevents the peaks at 10 TeV energies and steeply falls around ∼ ∼ matter from freely falling to a vicinity of the super- 100 TeV. This feature deviates sharply from the single massive black hole. An anomalously efficient angular power-law assumption for the spectrum incorporated in momentum transport mechanism is necessary to form the search for neutrino sources. Meanwhile, the other an accretion shock, which contradicts with the free-fall scenario considered here, the magnetic reconnection sce- assumption made by Inoue et al.(2019); Inoue et al. nario, can yield a neutrino spectrum that mimics a (2020). Lastly, the accretion rate through the accretion power-law with a cutoff at high energies. Neutrino pro- shocks is assumed to be extremely high to explain the duction in both scenarios may also receive a significant IceCube data (Inoue et al. 2019). It is comparable to contribution from the photomeson production process. or even higher than the accretion rate expected in the This often leads to a bumpy feature before the cutoff in geometrically thin, optically thick disk, which requires the spectrum. two accretion components with totally different angular In this section, we first estimate the neutrino flux momentum distributions. Note that PCR < 0.5Pth is from NGC 1068 and adopt parameters in each accel- required in this scenario, otherwise hot coronae cannot eration scenario that can describe the reported neutrino be maintained. spectrum from the source while maintaining other con- We also note that the Bethe-Heitler process, which straints on the parameters in a physical range. The main was often ignored, is relevant. The 10 100 TeV neu- parameter in our modeling that constrains the neutrino − trino flux from the photomeson production with coro- spectrum is the ratio of the CR to thermal pressures. Af- nal X-rays is suppressed by Bethe-Heitler interactions ter finding the parameters for the neutrino flux that can with disk photons (Murase et al. 2020). On the con- accommodate the neutrino flux from NGC 1068, we es- trary, given that shock acceleration leads to higher max- timate the neutrino flux for the brightest Seyfert galax- imum energies, PeV neutrinos are mainly produced via ies, similar to NGC 1068. Using the BAT AGN Spectro- the photomeson production with disk photons (Stecker scopic Survey (BASS), we select the ten sources that are et al. 1991). classified as Seyfert galaxies and pose the highest X-ray Also, the realistic coronal structure is inhomogeneous. flux. We incorporate the X-ray luminosity measured in As shown in Guti´errezet al.(2021), radio emission is 2-10 keV band. This range suits our selection criteria as more likely to come from outer coronal regions. Such we are considering the brightest X-ray sources whose lu- 42 a multi-zone situation is commonly invoked to explain minosities are LX & 10 erg/s. Less bright sources such 39 radio emission from blazars and radio galaxies because as Sgr A*, with LX < 10 erg/s are better characterized it is known that radio emission cannot be explained by by measurements at lower-energy bands. Through- single-zone models. out this work, we estimate the distance of the Seyfert Note that Kalashev et al.(2015) assumed electric ac- galaxies based on their redshift provided by BASS. To celeration in a gap formed in the magnetosphere. The evaluate the luminosity distance, we employ the follow- −1 −1 gap formation is possible only in low-luminosity AGNs ing cosmological parameters: H0 = 70 km s Mpc , such as M87, but it is likely to be screened in the case Ωm = 0.3, ΩΛ = 0.7. of Seyferts and quasars because the accretion rate is so high that the plasma satisfies the quasi-neutrality con- 3.1. Neutrino emission from NGC 1068 dition (Levinson & Rieger 2011). NGC 1068 is a Seyfert 2 galaxy at z = 0.00303 (Tully 1988), with a heavily obscured nucleus. It is considered as one of the best studied AGNs, thanks to its role in the 3. NEUTRINO EMISSION FROM BRIGHT AGN unification scenarios, which was initially proposed SEYFERT GALAXIES for explaining the broad optical lines in its polarized 6
25 −2 light. The column density of NGC 1068 is & 10 cm , implying a Compton thick environment. The high den- 43 1 6 NGC 1068 (LX = 10 erg s− ) sity environment at the core of NGC 1068 and the high 10− level of mass accretion rate provide favorable environ-
] IceCube 10 yr
1 7
ments for efficient production of high-energy neutrinos. − 10− s In order to estimate the neutrino spectrum in the 2 − Stochastic 8 (high CR pressure) AGN corona model for NGC 1068, we use the intrin- 10− sic X-ray luminosity, LX , suggested by Swift-BAT AGN
Spectroscopic Survey (BASS) (Ricci et al. 2017). The [GeVcm 9 Magnetic µ 10− ¯ BASS catalog provides X-ray fluxes and luminosities in ν Reconnection + µ a wide range, covering both soft (below 2 keV) and hard ν
Φ 10 (above 10 keV) X-rays. Here, we use measured spectrum 2 10− and luminosity in the 2-10 keV band, which provides a E 11 Stochastic medium range. This choice suits the disk-corona model- 10− (modest CR pressure) ing, as the abundance of data in this energy band enables us to use empirical correlation. 2 3 4 5 6 The intrinsic X-ray luminosity, L measured from the 10 10 10 10 10 X E [GeV] direction of NGC 1068 by BAT is 1043 erg/s for the ν ∼ 2-10 keV band, although it could be larger depending on estimates of the column density. We use this luminosity Figure 1. Modeled neutrino spectrum for NGC 1068 com- to model the neutrino spectrum from NGC 1068 in three pared to the best fit flux (yellow band) reported by the Ice- different scenarios. Cube Collaboration 10 yr point source study (Aartsen et al. Assuming a single power-law spectrum, the Ice- 2020b). The red line shows the expected flux in the stochas- Cube Collaboration has reported a best-fit flux of tic acceleration scenario matching IceCube’s best fit at TeVs. The purple line depicts the flux that would give the medium- 3 10−11 TeV−1cm−2s−1 at 1 TeV with an index of × energy neutrino flux, compatible with the total neutrino flux 3.2 (Aartsen et al. 2020b). As a criterion, we adopt ' reported in the cascade analysis (Aartsen et al. 2020a). The parameters that maintain the flux within the best-fit blue line presents the flux expected for the magnetic recon- reported flux and its uncertainties reported by the Ice- nection scenario. Cube Collaboration. First, we model neutrino production assuming the stochastically accelerated particles, where instead of stochastic acceleration scenario. As mentioned earlier, matching the flux at TeVs, we match the diffuse neu- in this scenario, the neutrino spectrum has a more com- trino flux at tens of TeV, motivated by the medium- plicated shape than a single power-law. Accommodat- energy excess in the neutrino spectrum. In this case, as ing the IceCube flux at TeV energies requires a relatively shown previously (Murase et al. 2020) we adopt param- high normalization while the spectrum has to cutoff fast eters that can explain the high-energy neutrino flux ex- enough that the spectrum drops around 100 TeV. Such cess observed at medium energies (Aartsen et al. 2020a). conditions would result in a high level of CR pressure in In this case, the P /P is set to 0.01. Here, the neu- CR th ' the corona model. trino spectrum peaks at 40 TeV, which corresponds to ∼ In order to maintain realistic scenarios, we restrict a lower level of neutrino flux compared to the previous ourselves to the range of parameters for which the ra- scenario. We refer to this case as “Modest CR pressure” tio of the CR pressure (PCR) to the thermal pressure hereafter. (Pth) is bound to less than 0.5. In this limit, the non- These results are compatible with the spectra pre- thermal energy is equal to half of the gravitational bind- sented previously by Murase et al.(2020) where the ing energy at the coronal radius without leaving room CR pressure considered to explain the medium-energy for thermal particles. Although the coronal plasma may neutrino flux and NGC 1068 are found at the level be heated more through magnetic fields connected to the of 1 and 30 percent of the thermal pressure, re- ∼ ∼ inner disk, we assume 0.5 as the maximal case in this spectively. Here, we allow the pressure ratio to be as work, and the neutrino spectrum peaks at 5 TeV and high as 50% to explain the soft spectrum reported for ∼ falls sharply around 20 TeV. We refer to this scenario as NGC 1068 by the IceCube Collaboration (Aartsen et al. “High CR pressure” 2020b). Note that, in principle, both the High CR pres- We consider the second scenario for neutrino emis- sure and Modest CR pressure cases can be viable within sion from NGC 1068 assuming coronal emission from the same stochastic acceleration scenario. For exam- High-Energy Neutrinos from Seyfert Galaxies 7
NGC 1068 NGC 1068 Stochastic (High CR pressure) 1 10− Stochastic (Modest CR pressure) 2 Magnetic Reconnection 10 − 3σ 0 Atmospheric 10 3 10−
4 10−
1 5 10− 10− p-value
events/bin/yr 6 10− 5σ
7 2 10− 10− Stochastic (High CR pressure) 8 Stochastic (Modest CR pressure) 10− Magnetic Reconnection 10 9 103 104 105 106 − 5 6 7 8 9 10 11 12 13 14 15 Eν [GeV] t [yr]
Figure 2. Expected event distribution per energy and year Figure 3. Prospects for identification of neutrino emission in IceCube is shown for the stochastic and reconnection sce- from NGC 1068. Local p-values are shown for stochastic with narios presented for NGC 1068 in Figure1. The atmo- high CR pressure (red), stochastic with modest CR pressure spheric neutrino background event distribution is also shown (violet), and reconnection (blue) scenarios considered in this for comparison (grey). study. ple, the Modest CR pressure may be realized in average trinos found from the direction of NGC 1068. The AGN, whereas some sources such as NGC 1068 may have shaded area shows the uncertainty on the fitted spec- a large CR pressure. trum as reported by IceCube. As shown, all modeled Finally, we consider the magnetic reconnection sce- neutrino spectra are within the 68% uncertainty of the nario for particle acceleration. In this case, the neu- measured spectrum. The parameters that we adapt in trino flux approximately follows mainly the initial CR each scenario for particle acceleration and interaction ef- spectrum until the pγ process becomes the dominant ficiency are presented in Table1. The common param- channel for the production of pions. Therefore, this sce- eters among different scenarios are the same as Murase nario leads to the spectrum close to a power-law spec- et al.(2020). The injected CR, i.e., proton, differen- trum with a cutoff at high energies. For the injected CR tial luminosity for the three scenarios shown in Fig.1 is spectrum, we assume a spectral index of 2. The nor- presented in the Appendix (see Fig. 14). malization and CR maximum energy are set such that We should note that a single power-law spectrum is the modeled flux is constrained to the IceCube steep not a realistic spectral energy distribution for neutrino spectrum reported for NGC 1068 while the PCR/Pth is emission from individual astrophysical objects. While rec bound to be smaller than 0.5. We find Ep 5 PeV neutrino and γ-ray spectra may, in general, reflect the rec ≈ for this purpose. Smaller values of Ep cannot accom- initial CR spectrum, the shape of neutrino and γ-ray modate the IceCube flux without violating the CR to fluxes depends on the nature of the interaction, thresh- thermal pressure maximum band. Larger values, on the olds, and the opacity of the source. The neutrino spectra other hand, would create an excess at high energies that provided in this study take all these into account. On are disfavored by the steep spectra reported for NGC the other hand, the diffuse flux of high-energy neutri- 1068. As described in Sec.2, we set ηacc = 300 for mag- nos (or γ-rays) over a specific range of energies may be cool netic reconnection acceleration. For NGC 1068, Ep is explained by a power-law since the superposition of the too high to match the IceCube data. individual sources would wash out the features. Figure1 shows the three modeled neutrino fluxes from We use the modeled neutrino spectra for NGC 1068 NGC 1068. We also projected the best-fit spectrum to compare with the findings of the IceCube 10 yr point reported by the IceCube Collaboration. The best-fit source study. In addition, we investigate the prospects power-law spectrum corresponds to the 51 excess neu- ∼ 8
of events and creates an excess beyond 100 TeV where NGC 1068 atmospheric neutrinos are scarce. Stochastic (High CR pressure) We present the p-value1 for the observation of each 6 IceCube 10 yr 10− scenario in 5 to 15 years of operation of IceCube in ]
1 Fig.3. We use the method described in ATLAS(2011) − s 2 7 to estimate the statistical significance for identification − 10− of the neutrino emission; for details, see Appendix. The stochastic acceleration scenario with the high CR pres- 8 [GeVcm 10− sure, compatible with the IceCube 10 year flux measure- µ ¯ ν +
µ ments, provides the most likely scenario for observation ν Φ
2 9 of NGC 1068 in IceCube, reaching to 5σ level in less 10− E than 12 years of data. This is in accordance with the reported local p-value in IceCube 10 year point source 10 10− study. The magnetic reconnection scenario is expected to be identified at better than 3σ within the 10 years, 1 2 3 4 5 10 10 10 10 10 while near a decade more observation might be required E [GeV] ν to establish a significance at the discovery level. Finally, the stochastic acceleration scenario with the modest CR Figure 4. Neutrino spectrum for NGC 1068, assuming dif- pressure seems difficult to identify with enough signifi- ferent values of the X-ray luminosity. Here, we study how the cance. However, the prospects for this scenario may be neutrino spectrum depends on the uncertainties associated conservative given uncertainties on the X-ray luminos- with X-ray observations. The line shows expected neutrino ity. As demonstrated in Fig.4, the neutrino flux would emission for the stochastic acceleration scenario, for L be- X be higher for larger intrinsic X-ray luminosities. tween 1042.8 erg s−1 and 1043.8 erg s−1. The key parameter in modeling the neutrino emis- sion from NGC 1068 is the intrinsic X-ray luminosity. for identification of each neutrino emission scenario in The column density (NH) for NGC 1068 is very high. the next decade of IceCube operation. Therefore, it is difficult to estimate the intrinsic X-ray In order to find the p-value for observation of neutri- luminosity, and measurements often carry large uncer- nos from NGC 1068 over the background of atmospheric tainties. The dominant source of uncertainty in mea- neutrinos, we calculate the number of signal neutrinos suring the X-ray luminosity from NGC 1068 is the frac- using the publicly available effective area for the Ice- tion of X-rays that is scattered into our line of sight Cube point source selection (Aartsen et al. 2017). We (Janssen et al. 2015). X-ray luminosities as high as also estimate the expected number of background at- 1044 erg/s have been considered for NGC 1068. NuS- mospheric neutrinos using the zenith dependent atmo- TAR and XMM?Newton monitoring campaigns (Mar- spheric neutrino flux reported by Honda et al.(2007), inucci et al. 2016) found L 7+7 1043 erg/s with ◦ X ' −4 × assuming a resolution of 0.7 . This is indeed larger than N 1025cm−2. Here, we used a central value of 1043 ◦ H ' IceCube’s nominal angular uncertainty of 0.5 at high erg/s, which is a conservative choice. We further in- energies (& 100 TeV). We chose this value because most vestigate the variation of the measured intrinsic X-ray neutrinos in our predicted spectra are found at the 1-10 luminosity for the first scenario. Figure4 shows how the TeV range. The energy distribution of the events per spectrum changes when the X-ray luminosity is varied year corresponding to the three neutrino scenarios we between 1042.8 to 1043.8 erg/s. Here, we set the CRs to consider here is shown in Fig.2. thermal pressure to the maximal value of 0.5. The larger The energy distribution of the events is a convolu- values of X-ray luminosity will increase the flux by al- tion of the spectrum and the IceCube effective area in most an order of magnitude, while the peak of the spec- the direction of NGC 1068. The stochastic acceleration trum is slightly shifted towards higher energies. With a scenario with the high CR pressure case results in an larger value of LX , smaller values of ηtur and PCR/Pth event distribution that peaks around 10 TeV. The mod- are demanded to match the IceCube data, which may est CR pressure case for the stochastic acceleration yield a lower rate, and the majority of the events are found with energies around 30 TeV. On the other hand, the 1 Note that the p-values presented here are pre-trial local p- harder spectrum at higher energies in the magnetic re- values and do not account for the trials associated with point source studies. connection scenario leads to a relatively flat distribution High-Energy Neutrinos from Seyfert Galaxies 9
rec Model PCR/Pth (CR) q ηtur s ηacc Ep Stochastic acceleration with high CR pressure 0.5 (0.03) 5/3 50 --- Stochastic acceleration with modest CR pressure 0.008 (0.0009) 5/3 10 --- Magnetic Reconnection 0.5 (0.01) -- 2.0 300 5 PeV Table 1. Model-dependent parameters and quantities adapted for estimating the neutrino flux from NGC 1068 in each scenario of particle acceleration shown in Fig.1. help relaxing the extreme CR production efficiency for 10 ] the high CR pressure case. 1
s 9
As we discussed earlier, in the magnetic reconnection 2 case, the neutrino spectrum follows the CRs flat spec- 8 trum with a cutoff at very high energies. IceCube’s 10 cm 3
1 7 year observation disfavors a hard spectrum and suggests 6 2 that the flux is suppressed at energies beyond 100 TeV.
We evaluate the neutrino detection rate using a neutrino GeV 5 9 spectrum of a simple power law with an exponential cut- 4 off to constrain the normalization and cutoff energy for [10 3 the neutrino flux. We present in Fig.5 the 2 and 3 σ ex- clusion regions for 10 years of IceCube observation. As 1 TeV 2 such, magnetic reconnection scenarios with large nor- 1 1 2 3 4 malization and/or large values of maximum energies are 10 10 10 10 not preferred. The constraints on Eν,cut will also be use- EE⌫cut,cut[TeV][TeV] ful for testing the accretion shock scenario. The diffusive shock acceleration theory predicts tacc ηaccrL/c with Figure 5. Exclusion region of the parameter space of muon 2 ≈ neutrino flux when spectrum is a single power-law (E−2 with ηacc (20/3)(c/Vs) , where Vs is the shock velocity. ≈ exponential cutoff for NGC 1068. The 2 and 3σ confidence This leads to Emax a few PeV in NGC 1068 (Inoue p ∼ level exclusion regions are obtained assuming that the hard et al. 2020). Thus, our results imply that fast acceler- spectrum with a cutoff is not found in the direction of NGC ation models including both magnetic reconnection and 1068. accretion shock scenarios can critically be tested with near-future neutrino observations. distance as reported by Harris et al.(2010); Karachent- 3.2. Neutrino emission from bright Seyfert galaxies sev et al.(2013); Tully et al.(2013). We should note that the neutrino flux is insensitive to the distance, be- We now use the parameters of the modeled neutrino cause the neutrino flux is roughly proportional to the flux from NGC 1068 to estimate neutrino spectra and X-ray flux. prospects for observations of other bright Seyfert galax- We implement the parameters that we adopted in ies. We first focus on IceCube as it is currently the the previous section for modeling the neutrino spectrum major operating neutrino telescope. We will extend our from NGC 1068, listed in Table1. Again, we consider 3 study into KM3NeT and the next generation of IceCube, scenarios: the stochastic acceleration scenario compati- IceCube-Gen2, in the next section. ble with the 10-year IceCube observation of NGC 1068 For this purpose, we rank Seyfert galaxies by their (high CR pressure), the stochastic acceleration scenario intrinsic X-ray fluxes measured by BAT AGN Spec- compatible with the medium-energy excess in the diffuse troscopic Survey (Ricci et al. 2017). We select the 10 neutrino flux (modest CR pressure), and the magnetic brightest sources for our study and extract their intrin- reconnection scenario. Similar to NGC 1068, we impose sic X-ray luminosities from the BASS catalog. Table2 the physical constraint on the ratio of CR pressure to presents the list of sources and their declination, dis- the thermal pressure to 0.5. Therefore, the normaliza- tance, redshift, intrinsic X-ray flux, and luminosity in tion of the CR spectrum is adjusted to maintain this 2-10 keV band. In order to estimate the luminosity dis- bound whenever this ratio exceeds 0.5 with NGC 1068 tance for source, we incorporate the redshift reported parameters. Figure6 demonstrates the spectrum for the in BASS together with the cosmological parameters in- source list considered in this study. For the stochastic troduced in the Introduction, with the exception of Cen acceleration scenario with the high CR pressure, we have A, Circinus Galaxy, and NGC 4945, for which we em- adjusted flux to match the ratio of CR pressure to the ployed astronomical measurements for evaluation of the 10
Source Declination z dL Intrinsic flux log(Intrinsic luminosity) [deg] [Mpc] [10−12erg cm−2 s−1] [erg s−1] Circinus Galaxy -65.34 0.0014 4.2? 984.4 42.31 ESO 138-1 -59.23 0.0091 39.2 671.3 44.09 NGC 7582 -42.37 0.0052 22.4 507.6 43.48 Cen A -43.02 0.00136 3.8? 347.3 42.39 NGC 1068 -0.013 0.00303 13.0 268.3 42.93 NGC 424 -38.08 0.0118 51.0 188.1 43.77 CGCG 164-019 27.03 0.0299 131.0 179.5 44.57 UGC 11910 10.23 0.0267 116.7 157.5 44.41 NGC 4945 -49.47 0.0019 3.6? 149.4 41.36 NGC 1275 41.51 0.0176 76.4 132.8 43.98 Table 2. Seyfert galaxies with the highest intrinsic X-ray fluxes in the range of 2-10 keV in BAT AGN Spectroscopic Survey (BASS) (Ricci et al. 2017). Here, we show the declination, redshift, and the intrinsic X-ray fluxes and luminosities reported by BASS in 2-10 keV band. For sources closer than 10 Mpc, specified by ?, we employ distances reported by astronomical measurements. For other sources, the luminosity distances are evaluated from the reported redshift.
thermal pressure to 0.5 by rescaling finj, which are used the Northern hemisphere, the event rates for the rest of in our analyses hereafter. The measured intrinsic X- the sources are low, weakening the likelihood of identi- ray luminosity and the distance of the source define the fying individual sources in IceCube. shape and magnitude of the spectra. In the stochastic Using the expected signal and background rates, we 43.5 −1 acceleration scenario, sources with LX & 10 erg s estimate the likelihood for observations of these sources peak at relatively higher energies with a narrow width in IceCube. Table3 summarizes the expected p-values while sources with smaller intrinsic X-ray luminosity under each emission scenario for 10 years of IceCube demonstrate a broader spectrum peaking at lower en- operation. The listed p-values show that for all three ergies. In this scenario, Cen A has the highest level acceleration scenarios, NGC 1068 is the brightest source of flux in the list of bright Seyfert galaxies. For the in IceCube. While the prospects for the identification of stochastic acceleration with modest CR pressure, we most sources are not much promising due to the suppres- employ finj value that yields the CR to thermal pres- sion of events in the Southern hemisphere, with contin- sure of 0.008 for NGC 1068. That corresponds to ued collection data CGCG 164-019 and NGC 1275 are ∼ f 10−7(5 10−5) at 2m c2 (1000m c2). Finally, in likely to be observed at 3σ level in 20 years of IceCube inj ' × p p the magnetic reconnection scenario, the larger X-ray lu- operation. However, we should note that the likelihood minosity leads to a larger contribution of photohadronic of observations depends on the neutrino emission sce- processes. For each source, as the maximum energy in nario: the stochastic acceleration scenario with the high this scenario depends on the size of the system as well CR pressure, compatible with NGC 1068 parameters, as the cooling effects. We evaluate Erec by rescaling would yield 3σ. p ∼ with M as Erec Erec−NGC 1068(M/M )1/2. Another source in the list worth discussing is NGC BH p ≈ p NGC 1068 This scaling is motivated by the model assumptions with 4945. The IceCube 10 year analysis found an excess of l /(2GMc−2) and β constant. 1 event in its direction, corresponding to a p-value of rec ∼ The majority of bright nearby Seyfert galaxies in Ta- 0.48, which is consistent to the expectations found in our ble2 are located in the Southern hemisphere, as pointed study. NGC 4945 is a starburst galaxy. We should note out by Murase et al.(2020). IceCube’s event selection that similar to NGC 1068, the neutrino flux from this is optimal for the Northern sky, where the Earth acts source cannot be explained by the starburst scenarios as a shield for the atmospheric muons. In the Southern such as Eichmann & Becker Tjus(2016). hemisphere, the event selection imposes a higher-energy Other than NGC 1068, our predictions indicate that threshold on the energy of the neutrinos to suppress the even optimistic scenarios are not strong enough to yield atmospheric muon background. This feature suppresses a statistically significant measurement of the neutrino the event rate for the majority of the luminous Seyfert emission from the rest of the bright nearby sources. As galaxies. We show the expected events from the sources such, a stacking search for neutrino emission from the in this list in Fig. 15 in Appendix. Except for the sources bright Seyfert galaxies is going to offer the best chance NGC 1275, UGC 11910, and CGCG 164-019 that are in for identifying these sources in IceCube. Stacking anal- High-Energy Neutrinos from Seyfert Galaxies 11
Stochastic (High CR pressure) Stochastic (Modest CR pressure) Magnetic Reconnection
NGC 424 NGC 1275 NGC 4945 ] 1 −
s 8 8 8 2 10− 10− 10− −
9 9 9 10− 10− 10− [GeVcm µ ¯ ν +
µ 10 10 10 ν 10− 10− 10− Φ 2 E 11 11 11 10− 10− 10− 103 104 105 106 103 104 105 106 103 104 105 106 Cen A Circinus Galaxy CGCG 164-019 ] 1 −
s 8 8 8 2 10− 10− 10− −
9 9 9 10− 10− 10− [GeVcm µ ¯ ν +
µ 10 10 10 ν 10− 10− 10− Φ 2 E 11 11 11 10− 10− 10− 103 104 105 106 103 104 105 106 103 104 105 106 ESO 138-1 UGC 11910 NGC 7582 ] 1 −
s 8 8 8 2 10− 10− 10− −
9 9 9 10− 10− 10− [GeVcm µ ¯ ν +
µ 10 10 10 ν 10− 10− 10− Φ 2 E 11 11 11 10− 10− 10− 103 104 105 106 103 104 105 106 103 104 105 106 Eν [GeV] Eν [GeV] Eν [GeV]
Figure 6. Neutrino flux for bright Seyfert galaxies considered in this study. Here, we show the high (red) and modest (purple) CR pressure stochastic acceleration scenarios as well as the magnetic recoonection scenario (blue), which provide compatible fluxes with the best-fit flux for NGC 1068 or the total neutrino spectrum measurement for parameters presented in Table1. yses are widely used to study the correlation of the ar- years of IceCube. The significance can reach a discov- rival direction of high-energy neutrinos and a catalog of ery level of 5σ within the 20 years of operation for this sources. The cumulative signal from all the sources im- scenario. proves the sensitivity and makes it possible to examine The most conservative scenario here is the stochastic a class of sources with lower fluxes. acceleration scenario with the modest CR pressure. This We project the prospects for identifying each scenario scenario would not yield a signal at a significant level of neutrino emission in a stacking search for 5 to 20 years with the current generation of IceCube. Mainly due to of IceCube in Fig.7. The solid lines show the p-value for the fact that most sources are in the Southern hemi- the duration of IceCube observation for each scenario. sphere. We investigate the likelihood of identifying this The most promising scenario is the stochastic accelera- scenario in the next generation of neutrino telescopes in tion scenario with the high CR pressure case. The high- the next section. level neutrino flux in this scenario would result in an We should note that the significance for the stacking excess of signal that can reach a discovery level of 5σ analysis for stochastic acceleration with high CR pres- in less than 10 years of IceCube operation. The hard sure scenario is dominated by the neutrino flux in the spectrum in the magnetic reconnection scenario would direction of NGC 1068. However, it is worth mention- result in an excess distinguishable from the background ing that for this scenario, a stacking analysis of nearby at high energies that can yield a 3σ evidence in about 7 bright Seyfert galaxies when NGC 1068 is excluded is 12
p-value Source Stochastic (High CR pressure) Stochastic (Modest CR pressure) Magnetic reconnection NGC 1068 10−6 0.09 1.8×10−4 NGC 1275 0.03 0.3 0.1 CGCG 164-019 0.04 0.3 0.1 UGC 11910 0.1 0.4 0.09 Cen A 0.5 0.2 0.2 Circinus Galaxy 0.5 0.3 0.3 NGC 7582 0.5 0.5 0.1 ESO 138-1 0.5 0.5 0.09 NGC 424 0.5 0.5 0.5 NGC 4945 0.5 0.5 0.5 Table 3. Prospects for observations of bright nearby Seyfert galaxies in 10 years of IceCube operations. expected to identify neutrino emission at the level of 3σ with 15 years of IceCube data. Such analysis can test IceCube Stacking 10 1 this scenario independently, without relying on the neu- − Stochastic trino emission reported from the direction NGC 1068, 2 (modest CR pressure) 10− which this scenario is based upon. 3σ Our predicted neutrino emission was built upon the 10 3 disk-corona model of AGN. Two of the sources in the − Magnetic list of bright Seyfert galaxies, Cen A and NGC 1275, are 4 10− Reconnection seen with a high jet activity and it is likely that the X- p-value 5 ray emission arises from the jet rather than the disk. We 10− will further investigate the prospects for observations of 6 bright Seyfert galaxies if the disk-corona emission would 10− 5σ not be dominant. The dashed lines in Fig.7 show the 7 Stochastic p-value when Cen A and NGC 1275 are not considered 10− (high CR pressure) in the source list. While the likelihood of observation is decreased, the cumulative neutrino emission is strong 5 10 15 20 t [yr] and the stacking analysis of the rest of the bright Seyfert galaxies can still reach a significant level in the lifetime of the IceCube detector. Figure 7. Projected p-values for stacking of 10 brightest For Cen A, we should note that although one zone Seyfert galaxies in 5 to 20 years of IceCube operation, for the three acceleration scenario considered in this study. The models typically attribute the high-energy to the solid lines show the p-values for the 10 bright sources listed jet Abdo et al.(2010), it may be difficult to explain some in Table2. The dashed line presents the prospects when the of the X-ray properties with such scenario (Tachibana emission from Cen A and NGC 1275 is excluded. et al. 2015). The observed soft lags in X-rays from Cen A may indicate their coronal origin. We should note with IceCube, which could reveal the dominant sources that despite the large magnitude of the flux predicted responsible for the medium-energy excess in the spec- under both stochastic acceleration and magnetic recon- trum of high-energy cosmic neutrinos. The likelihood nection scenarios for Cen A in our study, our prediction for observations of these sources will be enhanced by is compatible with current upper limits imposed by the the improvements in the event selection in the Southern ANTARES and IceCube Collaboration. In Fig.8 we hemisphere. We will address this in Sec.6. compare the predicted neutrino flux with the current limits from the IceCube cascade source search in the 4. FUTURE NEUTRINO TELESCOPES Southern sky (Albert et al. 2020). The upper limit is In this section, we explore the prospects for identify- obtained for a power-law with exponential cutoff, which ing bright nearby Seyfert galaxies in the next generation is the closest shape to the predicted spectra in our study. of neutrino telescopes. The major development in the In summary, observations of stacked neutrino emis- high-energy neutrino astrophysics would be driven by sion from nearby, bright Seyfert galaxies is promising KM3NeT, currently under construction in the Mediter- High-Energy Neutrinos from Seyfert Galaxies 13
Cen A NGC 1068 – KM3NeT 5 10− Stochastic (High CR pressure) Stochastic (High CR pressure) 6 Stochastic (Modest CR pressure) 10− Stochastic (Modest CR pressure)
] Magnetic Reconnection 1 0 Magnetic Reconnection −
s 7 10 2 10− PCR/Pth 0.5 Atmospheric − ' IceCube cascade U.L. 8 10−
[GeVcm 1 µ
¯ 9 ν 10− 10− + µ ν Φ events/bin/yr 2 10 10− PCR/Pth 0.2 E '
11 2 10− 10−
12 10− 103 104 105 106 103 104 105 106 Eν [GeV] Eν [GeV]
Figure 8. Neutrino spectra from Cen A compared with the discovery potential of a joint IceCube cascade source search Figure 9. Expected signal and background events in for the sources in the Southern sky (Albert et al. 2020). KM3NeT neutrino spectra presented in Fig.1. ranean, and the next generation of IceCube: IceCube- take into account the visibility of the source when cal- Gen2. culating the event rate and evaluating the prospects for KM3NeT is a cubic km-scale water Cherenkov detec- their identification. tor, which is designed to enhance the exposure in the We first examine the likelihood for observation of Southern hemisphere. The larger effective area and im- NGC 1068 in KM3NeT. Figure9 shows the expected proved angular resolution, compared to the current neu- signal and background rate. Given that the source vis- trino detector in the Mediterranean, ANTARES, are ex- ibility is 50%, it would not be identified in KM3NeT, pected to provide better sensitivity to the sources in even for the most optimistic scenario, i.e, the stochas- the Southern sky (Adrian-Martinez et al. 2016). This tic acceleration with the high CR pressure which would advancement would enhance the potential of identify- yield a p-value of few percent after 5 years. ing neutrino emission from bright Seyfert galaxies as We present the p-values for observations of all sources the majority of them are located in the Southern hemi- in 1 and 3 years of operation of KM3NeT in Table4. sphere. The brightest source in KM3NeT in the list of bright In order to estimate the likelihood of identifying neu- Seyfert galaxies is found to be Cen A given that its X- trino emission, we incorporate the effective area for up- rays come from the coronal region (that may not be going muon neutrinos in KM3NeT. We use the effective the case). Thanks to the higher level of neutrino flux, area reported in (Adrian-Martinez et al. 2016) and in- the stochastic acceleration scenario with the high CR clude the trigger efficiency for background rejection in pressure is going to be found at 3σ in the first year point source analyses. The publicly available effective of observation, and the significance can reach 5σ with area for the two blocks of ARCA is averaged over the additional two years of observation. zenith covering the upgoing events. Therefore, in our es- Circinus Galaxy is found to be the second brightest timation of the atmospheric background events, we use source in the Southern sky. In addition to its high level the average atmospheric neutrino flux instead of zenith of neutrino flux at TeV energies, the source benefits from dependent ones that we used for IceCube. We estimate a 100% visibility in KM3NeT. Therefore, the likelihood the signal and background events from the sources listed for its observation is high, which can exceed 3σ in 3 years in Table2. Unlike IceCube, KM3NeT is not located of operation for the stochastic acceleration scenario with at the geographic pole and sources’ zenith angle vary. the high CR pressure. Thus, we need to take into account the duration that As the signal events from the rest of the sources in the source is positioned below the horizon. As such, we the list fall short of yielding a statistical significance in 14
p-value 1 yr (3 yr) Source Visibility Stochastic (high CR pressure) Stochastic (Modest CR pressure) Magnetic Reconnection Cen A 0.7 0.001 (9.3×10−8) 0.2 (0.07) 0.2 (0.05) Circinus Galaxy 1.0 0.008 (1.9×10−5) 0.2 (0.09) 0.2 (0.07) ESO 138-1 1 0.1 (0.02) 0.4 (0.3) 0.3 (0.08) NGC 7582 0.7 0.2 (0.04) 0.4 (0.3) 0.4 (0.2) NGC 1068 0.5 0.2 (0.05) 0.4 (0.4) 0.4 (0.2) NGC 4945 0.8 0.5 (0.2) 0.5 (0.4) 0.5 (0.4) NGC 424 0.7 0.4 (0.2) 0.5 (0.4) 0.5 (0.4) UGC 11910 0.5 0.4 (0.4) 0.5 (0.5) 0.5 (0.5) CGCG 164-019 0.4 0.4 (0.3) 0.5 (0.5) 0.5 (0.5) NGC 1275 0.3 0.4 (0.4) 0.5 (0.5) 0.5 (0.5) Table 4. Prospects for observation of nearby bright Seyfert galaxies in one years of KM3NeT observations.
tive area for IceCube-Gen2 is 5 times larger than the 0 ∼ 10 current IceCube detector. IceCube-Gen2 KM3NeT We present the p-values expected for modest CR pres- 10 1 − sure stochastic acceleration scenario for KM3NeT to- 2 gether with the ones for IceCube-Gen2 in Fig. 10. We 10− 3σ project the prospects for identification of neutrino emis- 3 sion from the bright sources assuming an angular res- 10− olution of 0.3 (solid) and 0.7 degrees (dashed) for each 4 detector. We should note that our estimation of the 10−
p-value prospects for identifying Seyfert galaxies are quite con- 5 10− servative, given that an angular resolution of 0.3 degrees or better is not that far-fetched for KM3NeT. The ex- 6 10− 5σ pected improvements in the angular reconstruction in IceCube-Gen2 will also make it easier to identify these δψ = 0.7 10 7 sources. In fact, our estimates indicate that achieving − δψ = 0.3 finer angular resolutions at 10 30 TeV is crucial for ∼ − 2 4 6 8 10 the identification of neutrino emission from these sources t [yr] especially in the modest CR pressure case. We further show the growth of significance for a given resolution in Sec.6. Figure 10. Prospects for observation of the bright Seyfert galaxies in the next-generation neutrino telescopes: KM3NeT and IceCube-Gen2. The solid (dashed) lines show 5. DISCUSSION expectations for 0.3◦ (0.7◦) angular resolution for the mod- est CR pressure scenario. The tick lines show the prospects 5.1. Aggregated Fluxes for identification of the 10 nearby bright sources in Table2 Highly magnetized and turbulent coronae can be pos- in a stacking analysis. The thin lines show the prospects for sible sites of particle acceleration. The system is calori- identification of the sources in the absence of a signal from metric in the sense that sufficiently high-energy CRs are disk-corona model from Cen A and NGC 1275. depleted via hadronuclear and photohadronic interac- tions. The high magnitude of the neutrino flux at 10– 3 years, we now turn into the prospects for observation 100 TeV makes this scenario a primary candidate for the of neutrino emission in a stacking analysis. We only medium-energy neutrino flux observed in IceCube at the consider the modest CR pressure scenario in stochastic level of E2Φ 10−7 GeV cm−2 s−1 sr−1 (Murase et al. ν ν ∼ acceleration since emission under either of the other two 2020). The diffuse flux mainly originates from AGN at scenarios should be identified by IceCube. In addition high redshifts (with z 1 2), which are too far to de- ∼ − to KM3NeT, we consider IceCube-Gen2 for the stacking tect as individual sources. The contribution from local search in this scenario. Here, we assume that the effec- sources is small, but it is still of interest to ask their aggregated flux. High-Energy Neutrinos from Seyfert Galaxies 15
Stochastic (high CR pressure) Stochastic (modest CR pressure) Magnetic Reconnection ] ] ] 1 1 7
7 7 1