{"id":"9580af19-dae8-4b10-8a11-edb282690617","arxiv_id":"1909.02239","paper_version":4,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"NGC 1068's IceCube neutrinos are attributed to non-thermal particle acceleration in the black hole corona, with coronal parameters derived from a millimeter excess and a testable MeV gamma-ray prediction.","lead":"The paper argues that the high-energy neutrinos IceCube attributes to the Seyfert galaxy NGC 1068 are produced by accelerated particles in the hot corona around its central black hole. The model is consistent with millimeter observations and predicts a detectable MeV gamma-ray glow that future telescopes can test.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Screened geometry is doing all the work: the model's gamma-ray consistency fails if neutrino emission is co-spatial with the synchrotron-emitting corona that sets Rc.","rationale":"The paper builds a plausible case that the high neutrino flux relative to GeV gamma-rays forces a compact, photon-dense production site, and the coronal parameters derived from the mm excess are broadly consistent with those inferred for other Seyferts. The manuscript is honest about the tentative nature of the neutrino signal (2.9 sigma) and the mm excess. However, the decisive test of the model is its consistency with the GeV gamma-ray data. Section 4 explicitly shows the uniform-emissivity case violates those data; only the screened case works. The screened case is therefore not a harmless detail but the load-bearing element of the central claim. My concern sharpens the reader's: the screened geometry is not merely unconstrained but sits uneasily with the same model's parameters. The mm synchrotron excess is modelled as coming from non-thermal electrons distributed over a corona of size Rc=10 Rs; the protons responsible for the pion-production neutrinos are assumed to be injected with equal power, presumably in the same region. Co-spatial emission would require the uniform attenuation factor. To avoid this, the gamma-ray production region must be much smaller than the synchrotron region, an ad hoc separation not justified by the acceleration mechanism. A concrete re-calculation with an extended emitter and uniform attenuation, scanning the allowed eta_g range, would reveal whether any parameter set can satisfy the gamma-ray limits; if not, the model's success hinges entirely on an unsupported geometric assumption. This does not warrant rejection, because MeV observations and improved neutrino spectra can test the scenario, but it reinforces the conditional verdict.","tokens_in":9899,"tokens_out":7449,"duration_ms":80370,"concrete_test":"Using the same particle distribution that reproduces the mm excess (volume radius Rc=10 Rs), recompute the gamma-ray spectrum of Section 4 with the attenuation factor for a uniform, spatially extended emitter embedded in the coronal photon field (e.g., Monte Carlo radiative transfer or the analytic 3u(tau)/tau factor), and scan the gyrofactor over 30 <= eta_g <= 3e4. If no choice of eta_g yields a gamma-ray flux below the Fermi-LAT 4FGL/3FHL data points while the neutrino flux remains consistent with IceCube, then the screened geometry is a necessary and unsupported assumption; if some parameter set satisfies the gamma-ray limits, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central consistency claim in Section 4 depends entirely on the 'screened' attenuation case, exp(-tau), because the 'uniform emissivity' case is stated to violate the low-energy gamma-ray data. However, the coronal parameters used for the hadronic emission — Rc = 10 Rs, B = 100 G, p = 2.7, fnth = 0.03 — are fixed by modelling the mm synchrotron excess with non-thermal electrons filling that same corona (Section 3). The natural reading is that protons and electrons are accelerated throughout the same coronal volume; if so, the gamma-ray production region is co-spatial with the synchrotron-emitting region, and the uniform-mixing attenuation factor (approximately 3/tau for large tau) should apply. The screened case implicitly shrinks the gamma-ray/neutrino production region to a small inner sub-volume surrounded by the X-ray/UV photon field, but no physical mechanism or independent constraint supports this separation. The appeal to disk temperature gradients motivates a surrounding photon field, not a compact hadronic emission zone. Thus the agreement with GeV data is not a prediction of the model but an assumption chosen to make the model work.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that the 2.9-sigma IceCube neutrino excess from NGC 1068 is produced by non-thermal protons accelerated in the corona of the central supermassive black hole. It argues that the reported neutrino flux exceeding the GeV gamma-ray flux forces the emission region to be embedded in a very dense X-ray photon field, and that X-ray binaries cannot provide enough sources, leaving the SMBH corona as the only viable candidate. Using ALMA mm-band measurements, the authors fit a coronal synchrotron component with Rc = 10 Rs, B = 100 G, p = 2.7, and fnth = 0.03, importing the remaining coronal parameters from their earlier work. They then compute hadronic gamma-ray and neutrino emission and show that, in a 'screened' geometry where the gamma-ray production region is surrounded by the attenuating photon field, the model can reproduce the IceCube flux for gyrofactor 30 <= eta_g <= 3e4 without violating GeV data. The uniform-emissivity version of the model is stated to violate low-energy gamma-ray data, so the screened geometry is essential to the claimed consistency. The paper explicitly acknowledges that the mm excess is not firmly detected and that the coronal geometry requires further study.","tokens_in":10183,"tokens_out":5729,"duration_ms":63810,"significance":"If correct, the paper would identify the non-thermal corona of Seyfert galaxies as a high-energy neutrino production site, connecting the ALMA-observed mm excess with IceCube astrophysical neutrinos and motivating MeV gamma-ray observations. The manuscript is careful to state its own limitations: the mm excess is not a firm detection, the uniform-emissivity model violates GeV data, and the screened geometry is not independently constrained. However, the central consistency claim rests on this screened geometry, and the neutrino flux match spans a very wide range of the gyrofactor. The paper is therefore best read as a plausible scenario rather than a definitive identification; its main value is in sharpening the observational tests that could confirm or rule out the coronal neutrino hypothesis.","major_comments":[{"comment":"The central consistency with the GeV gamma-ray data depends entirely on the 'screened' attenuation case, exp(-tau), because the paper states that the uniform-emissivity model violates the low-energy gamma-ray data. However, the coronal parameters used for the hadronic emission (Rc = 10 Rs, B = 100 G, p = 2.7, fnth = 0.03) are fixed in Sec. 3 by fitting the mm synchrotron excess with non-thermal electrons filling that same corona. The natural reading is that the proton population responsible for neutrinos and gamma-rays occupies the same volume as the synchrotron-emitting electrons, in which case the uniform-mixing attenuation factor should apply. The screened case implicitly postulates a compact inner hadronic emission region surrounded by the attenuating photon field, but no physical mechanism or independent constraint is given for this separation; the disk temperature gradient argument motivates a surrounding photon field, not a compact hadronic zone. Since the model's gamma-ray consistency fails under the uniform assumption, the agreement with GeV data is an input assumption rather than a prediction. Please provide a concrete physical model for the radial stratification of the acceleration/emission region, or an independent observable that can test the screened geometry.","section":"Sec. 4, Fig. 2"},{"comment":"The mm excess that determines the coronal parameters is not a firm detection; the authors state that 'we cannot claim a firm detection of this component in NGC 1068, because of a paucity of flux measurements, mixture of beam sizes, and the complex source structure.' The neutrino and gamma-ray flux normalizations are tied to the non-thermal electron energy fraction and coronal parameters inferred from this excess. If the mm excess is not coronal synchrotron emission, the non-thermal electron population may have very different parameters, and the hadronic neutrino flux normalization changes accordingly. The paper should quantify how the predicted neutrino flux scales with the uncertain mm measurements and specify what remains of the conclusion if the excess is due to a different component.","section":"Sec. 3 and Sec. 4"},{"comment":"The model is compared with the IceCube data over a very wide parameter range, 30 <= eta_g <= 3e4, with the lower end eta_g = 30 imported from the authors' previous fit to the diffuse neutrino background (Inoue et al. 2019). As a result, the agreement demonstrates broad consistency rather than a sharp test of the model, and it does not discriminate between the coronal scenario and other possible origins of the 2.9-sigma neutrino excess. The paper should either narrow the allowed eta_g range using independent physical constraints or explicitly present the match as a consistency check rather than evidence that the corona is the unique source.","section":"Sec. 4, Fig. 2"}],"minor_comments":[{"comment":"The text says 'micorolensing observation' in the discussion of coronal size constraints; this should read 'microlensing observation.'","section":"Sec. 3"},{"comment":"The caption contains 'for for comparison'; the duplicate 'for' should be removed.","section":"Figure 2 caption"},{"comment":"The statement that the reported neutrino flux is higher than the GeV gamma-ray flux should specify whether this refers to integrated energy flux or E^2 dN/dE at a particular energy, since the two can behave differently across the spectrum.","section":"Sec. 2"},{"comment":"The definition of the optical depth tau in Eq. (1) is described as computed from the center of the corona, but the uniform-emissivity attenuation factor 3u(tau)/tau is stated without deriving the relationship; a brief definition or reference for u(tau) would improve clarity.","section":"Sec. 4"},{"comment":"The IceCube Collaboration et al. (2019) paper is cited as arXiv:1910.08488; the published version may have different page/table numbers, so the reference should be updated to the final journal version if available.","section":"References"},{"comment":"The abstract and conclusion state that the coronal parameters are 'consistent with the spectral excess registered in the millimeter range,' but Sec. 3 explicitly says the excess is not a firm detection; the wording should be softened to 'possibly associated with the reported mm excess' to avoid overstating the observational support.","section":"Abstract and Sec. 5"}],"recommendation":"major_revision","confidential_remarks":"The paper is a timely scenario piece, but the load-bearing screened geometry is not independently supported, and the mm excess used to fix the coronal parameters is explicitly not a firm detection. I do not think these issues are outright fatal, because the authors are transparent about the limitations and the model yields falsifiable MeV predictions. However, the manuscript in its current form presents as a definite identification of the coronal origin while the actual evidence supports only a conditional consistency. A major revision that either supplies a physical model for the screened geometry or reframes the central claim as a scenario with clearly stated assumptions would bring the paper in line with its own caveats."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper applies the authors' earlier coronal neutrino model to NGC 1068 and claims the IceCube hotspot can be explained by non-thermal protons in the SMBH corona, with parameters (Rc=10 Rs, B=100 G, p=2.7) taken from fitting the millimeter excess. The new contribution is source-specific: using ALMA data to pin down the corona parameters and showing that a screened geometry can match both the IceCube flux and the Fermi GeV upper limits, with a concrete MeV gamma-ray prediction that future missions can test. The qualitative argument that neutrino flux exceeding GeV flux demands attenuation in a dense X-ray photon field, and that the only viable site is the AGN corona, is sound and well presented.\n\nThat said, the central consistency claim is fragile. The screened geometry is doing all the work. The paper admits that the uniform emissivity case violates the low-energy gamma-ray data; only the screened case, where the production region sits inside a surrounding photon field and attenuation goes as exp(-tau), agrees with observations. But the coronal parameters used for hadronic emission are fixed by synchrotron radiation from non-thermal electrons that fill the same corona. If protons and electrons are co-spatial, the uniform mixing attenuation factor should apply. The paper's appeal to disk temperature gradients can motivate a photon field gradient but not a compact hadronic emission zone. Without an independent constraint on this geometry, the consistency with GeV data is effectively assumed, not predicted. The authors are transparent about this—they call for further study—but it remains a load-bearing weakness.\n\nThe other soft spots are minor by comparison: the neutrino signal is 2.9 sigma, the mm excess is explicitly not firm, and several inputs (fnth, temperature, optical depth, equal proton/electron power) are taken from previous work without independent verification for this source. The broad eta_g range (30 to 3e4) means the \"match\" is not sharp.\n\nVerdict: the paper is worth serious review. It is honest, well-written, and has a falsifiable prediction. A referee should push on the screened geometry and ask for a physical mechanism that would confine the hadronic production region while the synchrotron electrons occupy the full corona. That would likely soften the claim but not kill the idea. I would send it to review.","headline":"Plausible coronal neutrino model for NGC 1068, but the GeV consistency depends on an unconstrained screened geometry.","tokens_in":10723,"tokens_out":2639,"would_cite":true,"duration_ms":27859,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper argues that NGC 1068's IceCube neutrinos come from a non-thermal corona around its central black hole, with the accompanying gamma rays hidden by the corona's dense X-ray field.","keywords":["NGC 1068","Seyfert galaxies","active galactic nuclei","corona","high-energy neutrinos","IceCube","gamma-ray attenuation","millimeter excess"],"falsifier":"Observe the 0.1-10 MeV band of NGC 1068 with a future MeV telescope. The screened corona model predicts a measurable MeV gamma-ray component, because the internal optical depth is negligible at those energies; a null detection well below the predicted flux, or a neutrino spectral shape demanding $\\eta_g \\gg 3\\times10^4$, would refute the coronal origin.","tokens_in":9680,"feed_emoji":"🔭","tokens_out":8615,"duration_ms":86613,"temperature":0.7,"pith_summary":"The paper sets out to prove that the neutrino signal IceCube sees from the Seyfert galaxy NGC 1068 is made in the non-thermal corona of its central supermassive black hole. The core clue is that the reported neutrino flux exceeds the GeV gamma-ray flux, which cannot happen unless the neutrino source sits inside a dense keV-scale X-ray photon field that absorbs the accompanying gamma rays; only the region within tens of Schwarzschild radii of a compact object offers such a field. After counting X-ray binaries and finding them orders of magnitude too few, the paper identifies the corona as the only viable site. With a corona radius of 10 Schwarzschild radii, a 100 G magnetic field, and an electron index of 2.7, the model matches both the millimeter excess and the IceCube neutrino flux, provided the gamma-ray source is screened by the surrounding photon field. If correct, it explains why this heavily obscured Seyfert is the hottest neutrino spot and predicts a detectable MeV gamma-ray glow.","feed_headline":"Neutrinos from NGC 1068 trace back to the black hole's corona","feed_subtitle":"A shielded corona matching the millimeter excess reproduces the IceCube flux without breaking gamma-ray limits.","key_machinery":"The central object is the non-thermal corona of the supermassive black hole, treated as a compact source of synchrotron, Compton, and hadronic emission whose parameters are fixed by the millimeter excess. The structural load is carried by the gamma-ray attenuation formula $\\tau \\simeq 10^5\\,(\\epsilon_X/1\\,\\mathrm{keV})^{-1}(L_X/L_{\\mathrm{Edd}})(R_s/R)$ and by the choice between two attenuation geometries: the uniform-emissivity case, which attenuates gamma rays by the much larger factor $3u(\\tau)/\\tau$ and violates the low-energy gamma-ray data, and the screened case, which attenuates by $\\exp(-\\tau)$ and permits the neutrino flux. The screened geometry, with the source inside and the absorbing X-ray/UV photon field outside, lets the same corona produce neutrinos while hiding its GeV-TeV gamma-ray counterparts.","core_discovery":"The central claim is that the IceCube neutrinos from NGC 1068 are produced in the non-thermal corona of its supermassive black hole. Because hadronic neutrino production is accompanied by gamma-ray emission, the observed neutrino flux being higher than the GeV gamma-ray flux forces the source to sit inside a dense keV-scale X-ray photon field that absorbs the gamma rays: $\\tau \\simeq 10^5\\,(\\epsilon_X/1\\,\\mathrm{keV})^{-1}(L_X/L_{\\mathrm{Edd}})(R_s/R)$. Such a field can exist only within tens of Schwarzschild radii of a compact object, and the paper rules out X-ray binaries by number, leaving the corona. With $R_c = 10\\,R_s$, $B = 100$ G, electron index $p = 2.7$, and a screened geometry in which the attenuating X-ray/UV field surrounds the emission region, the model matches the IceCube neutrino flux for gyro factors $30 \\le \\eta_g \\le 3\\times10^4$ and satisfies the gamma-ray upper limits; the uniform-emissivity variant does not.","pith_inferences":["A testable ranking follows from the model: among nearby Seyferts, the neutrino flux should track the absorption-corrected intrinsic X-ray luminosity rather than the observed hard X-ray flux, because obscuration hides X-rays but not neutrinos.","If the screened geometry is generic, many neutrino-emitting coronae could be missing from GeV-TeV gamma-ray catalogs, making the diffuse neutrino background harder to pinpoint with gamma-ray surveys alone.","Extending the same logic to other type-2 Seyferts predicts a population of compact MeV gamma-ray sources that are heavily absorbed at higher energies and spatially coincident with obscured active nuclei; stacking their MeV flux should match the IceCube hotspot map."],"forward_implications":["If the coronal model is correct, NGC 1068's neutrinos are a byproduct of the same non-thermal acceleration that produces the coronal synchrotron millimeter excess, so millimeter variability should accompany changes in the neutrino output.","The model requires the neutrino and gamma-ray production region to be hidden behind a screen of X-ray and UV photons, so the GeV-TeV emission observed from NGC 1068 must have a separate origin, such as star formation, a jet, or a disk wind.","Because only the screened geometry survives the gamma-ray constraint, determining the coronal geometry through future spectral or variability observations decides whether the model lives or dies.","If the gyro factor is near 30, Seyfert coronae contribute substantially to the diffuse neutrino background up to a few hundred TeV; if future neutrino spectra require $\\eta_g \\gg 3\\times10^4$, Seyferts are only subdominant contributors.","The paper's logic implies that heavily obscured, intrinsically bright Seyferts are the best neutrino candidates, since obscuration hides X-rays but not neutrinos."],"supporting_citations":[{"why":"Reports NGC 1068 as the hottest spot in the 10-year survey and supplies the neutrino flux that the coronal model must reproduce.","marker":"IceCube Collaboration et al. 2019"},{"why":"Provides the X-ray binary luminosity function and the number of bright binaries in NGC 1068, used to rule out stellar-mass sources.","marker":"Swartz et al. 2011"},{"why":"Supplies the intrinsic 2-10 keV luminosity LX = 7e43 erg/s that sets the X-ray photon field and Eddington ratio.","marker":"Marinucci et al. 2016"},{"why":"Identifies the pc-scale free-free emission at cm bands, anchoring the subtraction of non-coronal radio components.","marker":"Gallimore et al. 2004"},{"why":"Provides ALMA 256 GHz core flux measurements that contribute to the millimeter excess.","marker":"García-Burillo et al. 2016"},{"why":"Provides ALMA 694 GHz core flux that constrains the synchrotron self-absorption break and corona parameters.","marker":"Impellizzeri et al. 2019"},{"why":"Establishes coronal synchrotron emission as a detected phenomenon in other Seyferts and gives prior determinations of corona size and magnetic field.","marker":"Inoue & Doi 2018"},{"why":"Supplies the coronal gamma-ray and neutrino emission model, the uniform-emissivity treatment, and the fiducial thermal corona assumptions used here.","marker":"Inoue et al. 2019"},{"why":"Provides the GeV gamma-ray flux that the screened model must stay below and the uniform model overproduces.","marker":"The Fermi-LAT collaboration 2019"},{"why":"Provides TeV gamma-ray upper limits that constrain the high-energy extension of the coronal model.","marker":"MAGIC Collaboration et al. 2019"}],"fun_headline_variants":["Neutrino burst from NGC 1068 traced to black hole's corona","Corona of black hole powers NGC 1068's neutrino glow","Neutrinos from NGC 1068: born in the black hole's corona","Black hole corona makes NGC 1068 shine in neutrinos"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the corona is screened: the gamma-ray and neutrino production region sits inside a surrounding X-ray and UV photon field, so attenuation goes as $\\exp(-\\tau)$ rather than the much less forgiving uniform-mixing factor; the paper motivates this with a disk temperature gradient but offers no independent constraint on the geometry.","fun_headline_variants_meta":{"raw":{"variants":["Neutrino burst from NGC 1068 traced to black hole's corona","Corona of black hole powers NGC 1068's neutrino glow","Neutrinos from NGC 1068: born in the black hole's corona","Black hole corona makes NGC 1068 shine in neutrinos"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00034,"raw_usage":{"total_tokens":1893,"prompt_tokens":981,"completion_tokens":912,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":597,"completion_tokens_details":{"reasoning_tokens":833}},"tokens_in":597,"tokens_out":912,"duration_ms":8154,"temperature":1.0,"reasoning_tokens":833,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:56:57.132709+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe the 0.1-10 MeV band of NGC 1068 with a future MeV telescope. The screened corona model predicts a measurable MeV gamma-ray component, because the internal optical depth is negligible at those energies; a null detection well below the predicted flux, or a neutrino spectral shape demanding $\\eta_g \\gg 3\\times10^4$, would refute the coronal origin.","supporting_citations":[{"cited_title":"A., Soria, R., Tennant, A","cited_arxiv_id":null,"evidence_quote":"Provides the X-ray binary luminosity function and the number of bright binaries in NGC 1068, used to rule out stellar-mass sources."},{"cited_title":"F., Baum, S","cited_arxiv_id":null,"evidence_quote":"Identifies the pc-scale free-free emission at cm bands, anchoring the subtraction of non-coronal radio components."},{"cited_title":"2019, ApJ, 880, 40","cited_arxiv_id":null,"evidence_quote":"Supplies the coronal gamma-ray and neutrino emission model, the uniform-emissivity treatment, and the fiducial thermal corona assumptions used here."}],"review_version":1}