⚛️ phenomenology

Multimessenger Constraints on Nuclear Acceleration in Neutrino-Emitting Seyfert Coronae

This paper investigates how stochastic acceleration and radiative cascades constrain the nuclear composition and power fractions in neutrino-emitting Seyfert coronae, demonstrating that limits on nuclear contributions depend critically on the specific nuclear mixture, acceleration conditions, and available power, as illustrated by the distinct behaviors of NGC 1068 and NGC 7469.

Original authors: Chengchao Yuan

Published 2026-10-14

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Chengchao Yuan

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Deep in the centers of distant galaxies, supermassive black holes act as cosmic engines, surrounded by superheated clouds of gas and magnetic fields known as coronae. These regions are violent places where particles are accelerated to energies far beyond what any machine on Earth could produce. When these high-speed particles collide with light or other matter, they can shatter into smaller pieces and create ghostly particles called neutrinos. Neutrinos are unique because they pass through almost everything without stopping, carrying a direct message from the heart of these extreme environments to telescopes on Earth. For years, astronomers have detected a steady stream of these neutrinos coming from a specific galaxy called NGC 1068, but a crucial question has remained unanswered: what kind of particles are being accelerated in the first place? Are they mostly simple hydrogen nuclei, or are they heavier, more complex atoms like helium or iron? The answer matters because the type of particle determines how much energy is available to create the neutrinos we see, and it changes the entire story of how these cosmic accelerators work.

A researcher has now tackled this mystery by building a detailed computer model of the environment inside the corona of NGC 1068 and a second galaxy, NGC 7469. Instead of guessing the composition of the particles, they simulated how different mixtures of atoms would behave under the intense conditions of a black hole's corona. They tested three distinct scenarios: one dominated by helium, a second with a mix of carbon, nitrogen, and oxygen, and a third made almost entirely of iron. The researcher tracked how these particles would be accelerated by magnetic turbulence, how they would lose energy by colliding with light or gas, and how they would break apart into lighter fragments. Crucially, they followed the chain reaction that follows these collisions, calculating how the resulting debris would produce both neutrinos and gamma rays, the high-energy cousins of visible light. By comparing their simulated results with actual observations from the IceCube neutrino detector and gamma-ray telescopes, they could see which mixtures of particles were allowed by the laws of physics and which were ruled out by the data.

The study reveals that the composition of the accelerated particles is not just a detail; it fundamentally changes what the universe looks like to our instruments. When the researcher assumed the corona was filled mostly with helium, the model allowed for a much larger fraction of the total energy to be carried by heavy nuclei without violating the observed limits on gamma rays. However, when they tested heavier mixtures like carbon, nitrogen, oxygen, or pure iron, the constraints became much tighter. In these heavier scenarios, the particles break apart more easily and lose energy faster, meaning they cannot reach the high energies needed to produce the neutrinos we detect unless the total power of the accelerator is significantly different. The researcher found that for NGC 1068, the data favors a scenario where lighter elements like helium play a major role, as this mixture fits the observed neutrino and gamma-ray signals more comfortably than the heavier options.

The situation changes slightly when looking at NGC 7469, a galaxy that appears to emit neutrinos with even higher energies than NGC 1068. Because the particles in this galaxy are being pushed to more extreme speeds, the model suggests that heavier nuclei can survive longer and contribute more significantly to the neutrino signal. In this case, the data allows for a broader range of compositions, including larger fractions of the heavier elements, because the higher energy environment supports a more vigorous cascade of particle interactions. The researcher emphasizes that their findings are based on simulations of how particles move and interact under specific physical laws, rather than a direct measurement of the particles themselves. They show that the answer depends on a delicate balance between how fast the particles are accelerated, how much energy is available, and how quickly they are destroyed by their surroundings.

Ultimately, this work demonstrates that we cannot understand the high-energy universe by looking at neutrinos alone. The type of particle being accelerated leaves a specific fingerprint on the radiation it produces. If the corona is rich in helium, the rules for how much energy can be assigned to heavy nuclei are different than if it is rich in iron. The study provides a new way to test these ideas, showing that by combining neutrino data with gamma-ray limits, astronomers can narrow down the possible ingredients of these cosmic accelerators. While the exact recipe for the particles in NGC 1068 and NGC 7469 is still being refined, this research establishes that the composition of the accelerated matter is a key variable that shapes the entire high-energy emission from these galactic cores. The results suggest that future observations will need to account for these nuclear mixtures to fully decode the messages arriving from the edge of black holes.

Technical Summary: Multimessenger Constraints on Nuclear Acceleration in Neutrino-Emitting Seyfert Coronae

Problem Statement
High-energy (HE) neutrinos from Seyfert galaxies, particularly NGC 1068 and the harder-spectrum candidate NGC 7469, provide a direct probe of particle acceleration near supermassive black holes. While the existence of HE neutrinos implies hadronic processes, the composition of the accelerated particles (protons vs. nuclei) remains uncertain. Previous models often assumed fixed injection spectra or specific compositions without fully accounting for the competition between stochastic acceleration, energy losses, escape, and nuclear destruction (fragmentation). A critical gap exists in understanding how the nuclear mixture affects the resulting neutrino and gamma-ray fluxes, specifically whether fragmentation can overcome energy limitations imposed during the acceleration phase without secondary reacceleration.

Methodology
The author constructs a steady-state model of a compact, turbulent Seyfert corona surrounding a supermassive black hole ( $M_{BH} \sim 10^7 M_\odot$ ). The framework integrates three main components:

Stochastic Acceleration and Transport: The model employs a Fokker–Planck (FP) equation to describe the evolution of primary particles (protons and seven nuclear species: $^4$ He, $^{12}$ C, $^{14}$ N, $^{16}$ O, $^{23}$ Na, $^{28}$ Si, $^{56}$ Fe). The FP stage includes acceleration via magnetic turbulence, escape via rigidity-dependent diffusion, and continuous losses (synchrotron, Bethe–Heitler pair production) and catastrophic losses (photomeson production, photodisintegration, gas collisions).
Power Normalization: The model assigns a fraction $f_{nuc}$ of the total cosmic-ray power ( $L_{CR}$ ) to nuclei ( $A>1$ ) and the remainder to protons. The injection rates are normalized such that the total kinetic power injected into the cascade matches the assigned acceleration power, independent of the specific escape rate.
Nuclear and Radiative Cascades: The stationary FP distributions serve as effective injection rates for a nuclear transport code (NuCascade). This code tracks the evolution of primary nuclei, daughter nuclei, and secondary nucleons through photonuclear interactions and decays. The resulting pion injection rates and electron–positron pairs (from Bethe–Heitler interactions) are fed into the time-dependent lepto-hadronic code am3 to calculate neutrino emission and electromagnetic (EM) cascades.

The study evaluates three composition templates: a Solar-like mixture (helium-dominated), a Carbon-Nitrogen-Oxygen (CNO) mixture, and an Iron (Fe) mixture. Constraints are applied using joint neutrino (IceCube) and gamma-ray (Fermi-LAT, MAGIC) flux selections for NGC 1068 and NGC 7469.

Key Results

Energy Limitation by Acceleration: The characteristic energy per nucleon available to nuclear cascades is determined by the balance between acceleration and losses. The author finds that without secondary reacceleration, fragmentation redistributes energy among lighter nuclei and nucleons but cannot increase the characteristic energy per nucleon beyond the scale set during the primary acceleration phase.
Composition-Dependent Constraints (NGC 1068): For NGC 1068, joint neutrino and gamma-ray constraints are highly sensitive to the nuclear mixture. The helium-dominated Solar template permits larger nuclear power fractions ( $f_{nuc}$ ) than the heavier CNO or Fe mixtures. This is because the Solar template yields a higher neutrino flux ratio ( $r_\nu$ ) at 3 TeV for purely nuclear injection compared to heavier elements, better compensating for the reduction in primary proton power.
Harder Spectra and Higher Energies (NGC 7469): For NGC 7469, which has a harder candidate neutrino spectrum, the model yields more energetic nuclear cascades due to higher parent nucleon energies. Consequently, the source-specific flux selections allow for larger nuclear fractions for the heavier (CNO and Fe) mixtures compared to the constraints found for NGC 1068.
Spectral Evolution: Increasing $f_{nuc}$ generally weakens the high-energy neutrino tail and shifts the spectral peak to lower energies, while making lower-energy emission relatively more prominent. However, substantial neutrino emission persists even at $f_{nuc} = 1$ .

Significance and Claims
The paper demonstrates that constraints on the nuclear composition of Seyfert coronae are not solely determined by the total power or acceleration efficiency but are intrinsically linked to the specific nuclear yields and the energy per nucleon supplied by the acceleration mechanism. The author claims that their results show:

Fragmentation cannot overcome energy limitations imposed during acceleration without secondary reacceleration.
Multimessenger constraints (neutrino and gamma-ray) are template-dependent; specifically, helium-dominated injections allow for higher nuclear power fractions in NGC 1068 than heavier mixtures under the adopted prescriptions.
The interplay between acceleration conditions, available power, and nuclear yields must be considered simultaneously to interpret neutrino data from AGNs.

The work provides a framework for assessing the contribution of Seyfert galaxies to the diffuse astrophysical neutrino background and suggests that future population studies should combine escaping nuclear spectra with propagation models to constrain the cosmic-ray composition. The author notes that magnetic deflections hinder direct association of TeV–PeV cosmic rays with individual Seyferts, reinforcing the value of neutrino and EM observations as direct probes of coronal acceleration.

More like this