Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fermi acceleration”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

AI-Enabled Operations at Fermi Complex: Multivariate Time Series Prediction for Outage Prediction and Diagnosis

The Main Control Room of the Fermilab accelerator complex continuously gathers extensive time-series data from thousands of sensors monitoring the beam. However, unplanned events such as trips or voltage fluctuations often result in beam outages, causing operational downtime. This downtime not only consumes operator effort in diagnosing and addressing the issue but also leads to unnecessary energy consumption by idle machines awaiting beam restoration. The current threshold-based alarm system is reactive and faces challenges including frequent false alarms and inconsistent outage-cause labeling. To address these limitations, we propose an AI-enabled framework that leverages predictive analytics and automated labeling. Using data from $2,703$ Linac devices and $80$ operator-labeled outages, we evaluate state-of-the-art deep learning architectures, including recurrent, attention-based, and linear models, for beam outage prediction. Additionally, we assess a Random Forest-based labeling system for providing consistent, confidence-scored outage annotations. Our findings highlight the strengths and weaknesses of these architectures for beam outage prediction and identify critical gaps that must be addressed to fully harness AI for transitioning downtime handling from reactive to predictive, ultimately reducing downtime and improving decision-making in accelerator management.

Jain, Milan [PNL, Richland] (ORCID:000000021676111↗
DOE OSTI2021

Suppression of the vacuum space-charge effect in fs-photoemission by a retarding electrostatic front lens

The performance of time-resolved photoemission experiments at fs-pulsed photon sources is ultimately limited by the e–e Coulomb interaction, downgrading energy and momentum resolution. Here, we present an approach to effectively suppress space-charge artifacts in momentum microscopes and photoemission microscopes. A retarding electrostatic field generated by a special objective lens repels slow electrons, retaining the k-image of the fast photoelectrons. The suppression of space-charge effects scales with the ratio of the photoelectron velocities of fast and slow electrons. Fields in the range from -20 to 01100 V/mm for E kin = 100 eV to 4 keV direct secondaries and pump-induced slow electrons back to the sample surface. Ray tracing simulations reveal that this happens within the first 40 to 3 μm above the sample surface for E kin = 100 eV to 4 keV. An optimized front-lens design allows switching between the conventional accelerating and the new retarding mode. Time-resolved experiments at E kin = 107 eV using fs extreme ultraviolet probe pulses from the free-electron laser FLASH reveal that the width of the Fermi edge increases by just 30 meV at an incident pump fluence of 22 mJ/cm2 (retarding field -21 V/mm). For an accelerating field of +2 kV/mm and a pump fluence of only 5 mJ/cm 2 , it increases by 0.5 eV (pump wavelength 1030 nm). At the given conditions, the suppression mode permits increasing the slow-electron yield by three to four orders of magnitude. The feasibility of the method at high energies is demonstrated without a pump beam at E kin = 3830 eV using hard x rays from the storage ring PETRA III. The approach opens up a previously inaccessible regime of pump fluences for photoemission experiments.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Fermi-LAT Detection of the Supernova Remnant G312.4-0.4 in the Vicinity of 4FGL J1409.1-6121e

Abstract Gamma-ray emission provides constraints on the nonthermal radiation processes at play in astrophysical particle accelerators. This allows both the nature of accelerated particles and the maximum energy that they can reach to be determined. Notably, it remains an open question to what extent supernova remnants contribute to the sea of Galactic cosmic rays. In the Galactic plane, at around 312° of Galactic longitude, Fermi-LAT observations show an extended source (4FGL J1409.1−6121e) around five powerful pulsars. This source is described by one large disk of 0.°7 radius with a high significance of 45σin the 4FGL-DR3 catalog. Using 14 yr of Fermi-LAT observations, we revisited this region with a detailed spectro-morphological analysis in order to disentangle its underlying structure. Three sources have been distinguished, including the supernova remnant G312.4−0.4 whose gamma-ray emission correlates well with the shell observed at radio energies. The hard spectrum detected by the LAT, extending up to 100 GeV without any sign of a cutoff, is well reproduced by a purely hadronic model.

Astronomy & Astrophysics↗
DOE OSTI2024

Thermomagnetic responses of semimetals

Solid-state thermomagnetic modules operating based on the Nernst–Ettingshausen effects are an alternative to conventional solid-state thermoelectric modules. These modules are appropriate for low-temperature applications where the thermoelectric modules are not efficient. Here, we briefly discuss the application, performance, similarities, and differences of thermoelectric and thermomagnetic materials and modules. We review thermomagnetic module design, Nernst coefficient measurement techniques, and theoretical advances, emphasizing the Nernst effect and factors influencing its response in semimetals such as carrier compensation, Fermi surface, mobility, phonon drag, and Berry curvature. The main objective is to summarize the materials design criteria to achieve high thermomagnetic performance to accelerate thermomagnetic materials discovery.

42 ENGINEERING↗
DOE OSTI2020

Discovery of a Gamma-Ray Black Widow Pulsar by GPU-accelerated Einstein@Home

We report the discovery of 1.97 ms period gamma-ray pulsations from the 75 minute orbital-period binary pulsar now named PSR J1653-0158. The associated Fermi Large Area Telescope gamma-ray source 4FGL J1653.6-0158 has long been expected to harbor a binary millisecond pulsar. Despite the pulsar-like gamma-ray spectrum and candidate optical/X-ray associations—whose periodic brightness modulations suggested an orbit—no radio pulsations had been found in many searches. The pulsar was discovered by directly searching the gamma-ray data using the GPU-accelerated Einstein@Home distributed volunteer computing system. The multidimensional parameter space was bounded by positional and orbital constraints obtained from the optical counterpart. More sensitive analyses of archival and new radio data using knowledge of the pulsar timing solution yield very stringent upper limits on radio emission. Any radio emission is thus either exceptionally weak, or eclipsed for a large fraction of the time. The pulsar has one of the three lowest inferred surface magnetic-field strengths of any known pulsar with B surf ≈ 4 × 10 7 G. The resulting mass function, combined with models of the companion star's optical light curve and spectra, suggests a pulsar mass gsim 2 M ⊙ . The companion is lightweight with mass ~0.01 M ⊙ , and the orbital period is the shortest known for any rotation-powered binary pulsar. This discovery demonstrates the Fermi Large Area Telescope's potential to discover extreme pulsars that would otherwise remain undetected.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Oscillations of sterile neutrinos from dark matter decay eliminates the IceCube-Fermi tension

IceCube has observed a flux of cosmic neutrinos, with a “bump” in the energy range 10 ≲ E/TeV≲100 that creates a 3σ tension with γ-ray data from the Fermi satellite. This has been interpreted as evidence for a population of hidden cosmic-ray accelerators. Here, we propose an alternative explanation of this conundrum on the basis of cold dark matter which decays into sterile neutrinos after oscillations produce the bump in the cosmic neutrino spectrum.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Mixed Precision Fermi-Operator Expansion on Tensor Cores from a Machine Learning Perspective

Here we present a second-order recursive Fermi-operator expansion scheme using mixed precision floating point operations to perform electronic structure calculations using tensor core units. A performance of over 100 teraFLOPs is achieved for half-precision floating point operations on Nvidia’s A100 tensor core units. The second-order recursive Fermi-operator scheme is formulated in terms of a generalized, differentiable deep neural network structure, which solves the quantum mechanical electronic structure problem. We demonstrate how this network can be accelerated by optimizing the weight and bias values to substantially reduce the number of layers required for convergence. We also show how this machine learning approach can be used to optimize the coefficients of the recursive Fermi-operator expansion to accurately represent the fractional occupation numbers of the electronic states at finite temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Nonresonant particle acceleration in strong turbulence: Comparison to kinetic and MHD simulations

Collisionless, magnetized turbulence offers a promising framework for the generation of nonthermal high-energy particles in various astrophysical sites. Yet, the detailed mechanism that governs particle acceleration has remained subject to debate. By means of 2D and 3D particle-in-cell, as well as 3D (incompressible) magnetohydrodynamic (MHD) simulations, we test here a recent model of nonresonant particle acceleration in strongly magnetized turbulence, which ascribes the energization of particles to their continuous interaction with the random velocity flow of the turbulence, in the spirit of the original Fermi model. To do so, we compare, for a large number of particles that were tracked in the simulations, the predicted and the observed histories of particles momenta. The predicted history is that derived from the model, after extracting from the simulations, at each point along the particle trajectory, the three force terms that control acceleration: the acceleration of the field line velocity projected along the field line direction, its shear projected along the same direction, and its transverse compressive part. Overall, we find a clear correlation between the model predictions and the numerical experiments, indicating that this nonresonant model can successfully account for the bulk of particle energization through Fermi-type processes in strongly magnetized turbulence. Additionally we also observe that the parallel shear contribution tends to dominate the physics of energization in the particle-in-cell simulations, while in the magnetohydrodynamic incompressible simulation, both the parallel shear and the transverse compressive term provide about equal contributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Nested Pebble Bed Blanket (NesPeB)

Recent advances in magnetic confinement fusion technology have attracted billions of dollars of investments in startups from venture capitals and corporations, resulting in the development of devices aiming to demonstrate net energy gain in a self-heated burning plasma, such as SPARC (under construction) and others. However, future fusion power plants must operate in regimes that will require technologies far beyond current experience. According to a National Academies of Science, Engineering, and Medicine report, to have nuclear fusion power plants contributing in a timely manner to the planned reduction of atmospheric carbon dioxide, a pilot plant should be built by 2035, and it should demonstrate fusion power production and the performance of the tritium fuel system (requiring a high enough tritium breeding) by 2040. A recognized key technology gap by [26] is the fusion first wall and blanket since no current blanket concept is considered satisfactory or has been built and proven. The first wall and blanket in magnetic fusion reactors form a vital and complex system, as it must satisfy different functions such as power extraction, tritium breeding, plasma containment, radiation shielding, and safety. The list of design requirements is even longer: high enough tritium production for fusion self-sufficiency, low material activation, decay heat and shutdown dose rates, high thermal efficiency, high-capacity factor, high magnets-divertor-vacuum vessel-first wall life, low corrosion, low cost, and intrinsically safe (requiring minimal licensing). Despite fifty-plus years of research, the first wall and blanket concepts proposed suffer from fundamental technical problems and immaturity (TRL=2-3) that jeopardize the timely delivery of a commercial fusion power plant. A fusion first-wall blanket has never been built nor tested, and a "winning", practical functioning design requires enough engineering margins (high enough tritium breeding considering the uncertainty, etc.), manufacturing simplicity, ease of continuous operation, maintenance, and low cost. A new, groundbreaking blanket concept called "Nested Pebble Bed Blanket" (NesPeB) was developed at ORNL under the successful ARPA-E GAMOW FERMI project (patent application allowed by the USPTO). The NesPeB blanket concept addresses current blanket concepts' shortcomings and technical immaturity, paving the way for accelerated delivery of fusion power plants. NesPeB is based on nested pebbles, which are binary-sized lithium-ceramic pebbles enclosed in "Beryllide" perforated and coated spherical shells, which are also binary-sized, stacked on top of each other, forming a "bed" and cooled by Nitrogen gas also "sweeping" the Helium and Tritium generated by the neutron irradiation of Lithium; the vacuum vessel plasma facing material is Molybdenum-96 and -97 with the first wall cooled by Helium while the divertor armor is made of Tungsten. The simulations of the NesPeB blanket using Fusion Reactors Models Integrator (FERMI) are encouraging as they estimate a tritium breeding ratio (TBR) greater than 1.2 using natural Lithium, acceptable pressure drop, and excellent heat transfer properties. Furthermore, the NesPeB blanket is not limited by magneto-hydro-dynamics (MHD) effects, is designed for online refueling, relies on existing tritium extraction technologies, has a simple construction, and limits the corrosion and chemical reactivity problems. NesPeB has the potential to be transformational and disruptive since it can solve all the main, challenging technical problems of fusion device blankets and accelerate a pilot plant delivery for 10 or more years.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Evidence for flare-accelerated particles in large scale loops in the behind-the-limb gamma-ray solar flare of September 29, 2022

We report on the detection of the gamma-ray emission above 100 MeV from the solar flare of September 29, 2022, by Fermi LAT with simultaneous coverage in HXR by Solar Orbiter STIX. The Solar Orbiter-Earth separation was 178° at the time of the flare as seen from Earth, with Solar Orbiter observing the east limb. Based on STIX imaging, the flare was located 16° behind the eastern limb as seen from Earth. The STIX and GBM non-thermal emission and the LAT emission above 100 MeV all show similarly shaped time profiles, and the Fermi profiles peaked only 20 s after the STIX signal from the main flare site, setting this flare apart from all the other occulted flares observed by Fermi LAT. The radio spectral imaging based on the Nançay Radioheliograph and ORFEES spectrograph reveal geometries consistent with a magnetic structure that connects the parent active region behind the limb to the visible disk. We studied the basic characteristics of the gamma-ray time profile, in particular, the rise and decay times and the time delay between the gamma-ray and HXR peak fluxes. We compared the characteristics of this event with those of four Fermi LAT behind-the-limb flares and with an on-disk event and found that this event is strikingly similar to the impulsive on-disk flare. Based on multiwavelength observations, we find that the gamma-ray emission above 100 MeV originated from ions accelerated in the parent active region behind the limb and was transported to the visible disk via a large magnetic structure connected to the parent active region behind the limb. Our results strongly suggest that the source of the emission above 100 MeV from the September 29, 2022 flare cannot be the CME-driven shock.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Causes of fuel–ablator mix inferred from modeling of monochromatic time-gated radiography of OMEGA cryogenic implosions

Here, we present evidence, in the context of OMEGA cryogenic target implosions, that laser imprint, known to be capable of degrading laser-direct-drive target performance, plays a major role in generating fuel–ablator mix. OMEGA cryogenic target implosions show a performance boundary correlated with acceleration-phase shell stability; for sufficiently low adiabats (where the adiabat is the ratio of the pressure to the Fermi pressure) and high in-flight aspect ratios (IFAR's), the neutron-weighted shell areal density and neutron yield relative to the clean simulated values sharply decline. Direct evidence of Rayleigh–Taylor fuel–ablator mixing was previously obtained using a Si Heα backlighter driven by an ~20-ps short pulse generated by OMEGA EP. The shadow cast by the shell shortly prior to stagnation, as diagnosed using backlit radiographs, shows a softening near the limb, which is evidence of an ablator–fuel mix region for a low-adiabat implosion (α ~ 1.9, IFAR = 14) but not for a moderate adiabat implosion (α ~ 2.5, IFAR = 10). We find good agreement between experimental and synthetic radiographs in simulations that model laser imprint and account for uncertainty in the initial ablator thickness. Furthermore, we further explore the role of other mechanisms such as classical instability growth at the fuel–ablator interface, species concentration diffusion, and long-wavelength drive and target asymmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Timing of pulsars in the globular cluster omega centauri

We present the timing of the first five millisecond pulsars discovered in the globular cluster Omega Centauri and the discovery of a pulsar with a spin period of 3.68 ms. With a timing baseline of ~3.5 yr we are able to measure the derivative of the spin frequency ($\dot{\nu }$) for the first five pulsars. Upper limits on the pulsar line-of-sight acceleration are estimated and compared with predictions based on analytical models of the cluster. We find that PSRs J1326-4728B and D show large negative accelerations, which are in tension with the minimum acceleration predicted by analytical models. We searched for pulsed γ-ray signals using 14.3 yr of data from the Fermi Large Area Telescope. Although we found no evidence for γ-ray pulsations, PSRs J1326-4728A, B, C, and E are associated with X-ray sources. This suggests that the observed γ-ray emission from Omega Centauri is likely caused by the emission of the ensemble of MSPs. Finally, the linearly polarized emission from PSR J1326-4728A yields a rotation measure of -18 ± 8 rad m -2 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Development of High Performance Medium Velocity Superconducting Elliptical Cavities for Hadron Linacs

There is a solid scientific need for continuous wave (CW) proton or heavy ion linear accelerators (linacs) in the energy range up to about 1 GeV per nucleon (GeV/u). Only superconducting (SC) cavities can provide cost-efficient acceleration in CW mode. While low velocity, bc < 0.5 (b is the relative velocity of a charged particle and c is the speed of light), and high velocity, bc > 0.9, SC accelerating structures are established, there is a need in high-performance medium velocity SC accelerating structures. Currently, no superconducting (SC) CW proton or ion accelerators operate in the medium velocity range. We have developed a robust technology for producing high-performance medium-velocity elliptical cavities operating in the frequency range between 600 and 700 MHz. The project was performed in collaboration with Fermi and Argonne National Laboratories. The 5-cell, βopt = 0.65, 644 MHz elliptical cavities built at MSU have been used for this research. An advanced surface treatment method based on nitrogen doping to achieve a high-quality factor and a high accelerating gradient has been demonstrated. The surface treatment has been performed at ANL, and nitrogen doping at FNAL. A stable resonance control with the slow and fast tuner system has also been verified. FNAL developed a dedicated resonance control system for the stable operation of the combined tuner.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Hunting extreme BL Lacertae blazars with Fermi-Large Area Telescope

ABSTRACT The emission of very-high-energy (VHE) photons ($E\ \gt\ 100\ \mathrm{GeV}$) in blazars is closely connected to the production of ultra-relativistic particles and the role of these γ-ray sources as cosmic particle accelerators. This work focuses on a selection of 22 γ-ray objects from the 2BIGB catalogue of high-synchrotron-peaked sources, which are classified as blazar candidates of uncertain type in the 4FGL-DR2 catalogue. We study these sources by means of a re-analysis of the first 10 yr of γ-ray data taken with the Fermi-Large Area Telescope, including the attenuation by the extragalactic background light. Their broad-band spectral energy distributions are also evaluated, using multiwavelength archival data in the radio, optical, and X-ray bands, in terms of one-zone synchrotron-self-Compton models, adding an external Compton component when needed. Out of this analysis, we identify 17 new extreme high-synchrotron-peaked (EHSP) candidates and compare their physical parameters with those of prototypical EHSP blazars. Finally, the resulting models are used to assess their detectability by the present and future generation of ground-based Imaging Atmospheric Cherenkov Telescopes. We find two VHE candidates within the reach of the current and next generation of Cherenkov telescopes: J0847.0−2336 and J1714.0−2029.

Nievas Rosillo, M. (ORCID:0000000283219168)↗
DOE OSTI2025

Elucidation of Ultrafast Decay, Vibrational Beating, and Slow Decay Processes of Excited Azulene

Azulene’s nonradiative decay dynamics and kinetics from its singlet excited-states were studied using mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) combined with trajectory surface-hopping nonadiabatic molecular dynamics (NAMD) and Fermi’s golden rule (FGR) rate theory. The NAMD dynamics reproduce experimental observations that the S 1 → S 0 decay is accelerated and exhibits a crossover from mono- to biexponential kinetics with increasing excess vibrational energy. Minimum-energy-path analyses reveal a continuous S 1 /S 0 conical-intersection seam slightly above the S 1 minimum, providing readily accessible funnels. Time-resolved normal-mode projections reveal selective energy funneling into C–C stretching modes at 1272, 1528, and 1692 cm –1 . Constructive combinations of the latter two modes appear to promote rapid access to the conical-intersection seam, whereas their beating at ∼ 160 fs imposes a natural limit on the decrease of the decay time upon a further increase of excess energies, suggesting that their interference delays the decay. Here, the FGR rate calculation data for the S 1 → S 0 transition reaffirm that the nonradiative decay proceeds mainly near or through the conical intersection, rather than via simple nonadiabatic derivative coupling around the minimum of S 1 . On the other hand, the FGR rates indicate that S 2 decays predominantly to S 0 , as the S 2 –S 1 vibronic coupling is exceptionally weak, which serves as a primary cause for azulene’s characteristic anti-Kasha emission.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Superfluid condensate fraction and pairing wave function of the unitary Fermi gas

The unitary Fermi gas is a many-body system of two-component fermions with zero-range interactions tuned to infinite scattering length. Despite much activity and interest in unitary Fermi gas and its universal properties, there have been great difficulties in performing accurate calculations of the superfluid condensate fraction and pairing wave function. In this paper, we present auxiliary-field lattice Monte Carlo simulations using a lattice interaction which accelerates the approach to the continuum limit, thereby allowing for robust calculations of these difficult observables. As a benchmark test, we compute the ground-state energy of 33 spin-up and 33 spin-down particles. As a fraction of the free Fermi gas energy $E_{\text{FG}}$, we find $E_0/E_{\text{FG}}$ = 0.369(2), 0.372(2), using two different definitions of the finite-system energy ratio, in agreement with the latest theoretical and experimental results. We then determine the condensate fraction by measuring off-diagonal long-range order in the two-body density matrix. We find that the fraction of condensed pairs is α = 0.43(2). Further, we also extract the pairing wave function and find the pair correlation length to be $ζ_pk_F$ = 1.8(3)ℏ, where $k_F$ is the Fermi momentum. Provided that the simulations can be performed without severe sign oscillations, the methods we present here can be applied to superfluid neutron matter as well as more exotic $\textit{P}$-wave and $\textit{D}$-wave superfluids.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Very-high-energy Emission from Pulsars

Abstract Air-Cherenkov telescopes have detected pulsations at energies above 50 GeV from a growing number of Fermi pulsars. These include the Crab, Vela, PSR B1706−44, and Geminga, with the first two having pulsed detections above 1 TeV. In some cases, there appears to be very-high-energy (VHE) emission that is an extension of the Fermi spectra to high energies, while in other cases, additional higher-energy spectral components that require a separate emission mechanism may be present. We present results of broadband spectral modeling using global magnetospheric fields and multiple emission mechanisms that include synchro-curvature (SC) and inverse Compton scattered (ICS) radiation from accelerated particles (primaries) and synchrotron self-Compton (SSC) emission from lower-energy pairs. Our models predict three distinct VHE components: SC from primaries whose high-energy tail can extend to 100 GeV, SSC from pairs that can extend to several TeV, and ICS from primary particles accelerated in the current sheet that scatter pair synchrotron radiation, which appears beyond 10 TeV. Our models suggest that H.E.S.S.-II and MAGIC have detected the high-energy tail of the primary SC component that produces the Fermi spectrum in Vela, Geminga, and PSR B1706−44. We argue that the ICS component peaking above 10 TeV from Vela has been seen by H.E.S.S. Detection of this emission component from the Crab and other pulsars is possible with the High Altitude Water Cherenkov Observatory and Cherenkov Telescope Array, and will directly measure the maximum particle energy in pulsars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Disentangling the Hadronic Components in NGC 1068

The recent detection of high-energy neutrinos by IceCube in the direction of the nearby Seyfert/starburst galaxy NGC 1068 implies that radio-quiet active galactic nuclei can accelerate cosmic-ray ions. Dedicated multimessenger analyses suggest that the interaction of these high-energy ions with ambient gas or photons happens in a region of the galaxy that is highly opaque for GeV–TeV gamma rays. Otherwise, the GeV–TeV emission would violate existing constraints provided by the Fermi Large Area Telescope (LAT) and the Major Atmospheric Gamma Imaging Cherenkov. The conditions of high optical depth are realized near the central supermassive black hole (SMBH). At the same time, the GeV emission detected by the Fermi LAT is likely related to the galaxy’s sustained star formation activity. In this work, we derive a 20 MeV–1 TeV spectrum of NGC 1068 using 14 yr of Fermi LAT observations. We find that the starburst hadronic component is responsible for NGC 1068's emission above ~500 MeV. However, below this energy, an additional component is required. In the 20–500 MeV range, the Fermi LAT data are consistent with hadronic emission initiated by non-thermal ions interacting with gas or photons in the vicinity of the central SMBH. This highlights the importance of the MeV band to discover hidden cosmic-ray accelerators.

79 ASTRONOMY AND ASTROPHYSICS↗