Search NASA⌕ Search

SEARCH · Search NASA

Results for “High-Energy Astrophysics”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Towards modelling AR Sco: calibration – reproducing high-energy pulsar emission and testing convergence to Aristotelian electrodynamics

In recent years, kinetic simulations have been crucial to further our understanding of pulsar electrodynamics. Yet, due to the large-scale separation between the gyro-period and the stellar rotation period, resolving the particle gyration has been computationally unfeasible for realistic pulsar parameters. The main aim of this work is comparing our gyro-phase-resolved model with a gyro-centric pulsar model, where our model solves the general equations of motion with included radiation reaction using a higher order numerical solver with adaptive time-steps. Specifically, we aim to (i) reproduce a pulsar’s high-energy emission maps, namely one with 10 per cent of the surface B-field strength of Vela, and the spectra produced by an independent gyro-centric pulsar emission model; and (ii) test convergence of these results to the radiation-reaction limit of Aristotelian electrodynamics. (iii) Additionally, we identify the effect that a large $E_{\parallel }$-field has on the trajectories and radiation calculations. We find that we can reproduce the curvature radiation emission maps and spectra well, using 10 per cent field strengths of the Vela pulsar and injecting our particles at a higher altitude in the magnetosphere. Using sufficiently large $E_{\parallel }$-fields, our numeric results converge to the analytic radiation-reaction limit trajectories. Additionally, we illustrate the importance of accounting for the $\mathbf {E}\times \mathbf {B}$-drift in the particle trajectories and radiation calculations, validating the Harding and collaborators’ model approach. Lastly, we found that our model deals very well with the high-radiation-reaction and high-field regimes present in pulsars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Cosmic-ray cooling in active galactic nuclei as a new probe of inelastic dark matter

We present a novel way to probe inelastic dark matter using cosmic-ray (CR) cooling in active galactic nuclei (AGNs). Dark matter (DM) in the vicinity of supermassive black holes may scatter off CRs, resulting in the rapid cooling of CRs for sufficiently large cross sections. This in turn can alter the high-energy neutrino and gamma-ray fluxes detected from these sources. We show that AGN cooling bounds obtained through the multimessenger data of NGC 1068 and TXS 0506 + 056 allows us to reach unprecedently large mass splittings for inelastic DM (≳ TeV), orders of magnitude larger than those probed by direct detection experiments and DM capture in neutron stars. Furthermore, we demonstrate that cooling bounds from AGNs can probe thermal light DM with small mass splittings. This provides novel and complementary constraints in parts of a parameter space accessible solely by colliders and beam-dump experiments.

astronomical black holes↗
DOE OSTI2024

High-energy neutrino deep inelastic scattering cross sections

We present a state-of-the-art prediction for cross sections of neutrino deep inelastic scattering (DIS) from nucleon at high neutrino energies, E v , up to 1000 EeV (10 12 GeV). Our calculations are based on the latest CT18 NNLO parton distribution functions (PDFs) and their associated uncertainties. To make predictions for the highest energies, we extrapolate the PDFs to small x according to several procedures and assumptions, thus affecting the uncertainties at ultrahigh E v ; we quantify the uncertainties corresponding to these choices. Similarly, we quantify the uncertainties introduced by the nuclear corrections that are required to evaluate neutrino-nuclear cross sections for the neutrino observatories. These results can be applied to currently running astrophysical neutrino observatories, such as IceCube and KM3NeT, as well as various future experiments that have been proposed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Optimizing fictitious states for Bell inequality violation in bipartite qubit systems with applications to the t t ¯ system

There is a significant interest in testing quantum entanglement and Bell inequality violation in high-energy experiments. Since the analyses in high-energy experiments are performed with events statistically averaged over phase space, the states used to determine observables depend on the choice of coordinates through an event-dependent basis and are thus not genuine quantum states, but rather “fictitious states.” We find that the basis which diagonalizes the spin-spin correlations is optimal for constructing fictitious states to test the violation of Bell’s inequality. This result is applied directly to the bipartite qubit system of a top and antitop produced at a hadron collider. We show that the beam axis is the optimal basis choice near the t t ¯ threshold production for measuring Bell inequality violation, while at high transverse momentum the basis that aligns along the momentum direction of the top is optimal. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

High-energy neutrinos from millicharged dark matter annihilation in the Sun

Millicharged dark matter particles can be efficiently captured by the Sun, where they annihilate into tau leptons, leading to the production of high-energy neutrinos. In contrast to the Earth, the high temperature of the Sun suppresses the fraction of millicharged particles that are bound to nuclei, allowing for potentially high annihilation rates. We recast existing constraints from the IceCube Neutrino Observatory and use this information to place new limits on the fraction of the dark matter that is millicharged. This analysis excludes previously unexplored parameter space for masses of m χ ∼ ( 5 – 100 ) GeV , charges of q χ ∼ 10 − 3 – 10 − 2 , and fractional abundances as small as f DM ∼ 10 − 5 . Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Final state radiation from high and ultrahigh energy neutrino interactions

Charged leptons produced by high-energy and ultrahigh-energy neutrinos have a substantial probability of emitting prompt internal bremsstrahlung ν ℓ + N → ℓ + X + γ . This can have important consequences for neutrino detection. We discuss observable consequences at high- and ultrahigh-energy neutrino telescopes and the Large Hadron Collider’s (LHC’s) Forward Physics Facility. Logarithmic enhancements can be substantial (e.g., ∼ 20 % ) when either the charged lepton’s energy or the rest of the cascade is measured. We comment on final state radiation’s impacts on measuring the inelasticity distribution, ν / ν ¯ flux ratio, throughgoing muons, and double-bang signatures for high-energy neutrino observation. Furthermore, for ultrahigh-energy neutrino observation, we find that final state radiation increases the overall detectable energy by as much as 20%, affects flavor measurements, and decreases the energy of both Earth-emergent tau leptons and regenerated tau neutrinos. Many of these have significant impacts on measuring neutrino fluxes and spectra. Finally, for the LHC’s Forward Physics Facility, we find that final state radiation will impact future extractions of strange quark parton distribution functions. Final state radiation should be included in future analyses at neutrino telescopes and the Forward Physics Facility. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Energy gain scale calibration of the XRISM Resolve microcalorimeter spectrometer: ground calibration results and on-orbit comparison

The Resolve instrument aboard the X-ray Imaging and Spectroscopy Mission (XRISM) is a 36-pixel microcalorimeter spectrometer that provides nondispersive spectroscopy with ∼ 5 eV spectral resolution in the soft X-ray waveband. Resolve has a requirement to provide an absolute energy-scale calibration of ± 2 eV from 0.3 to 12 keV. We describe our ground calibration strategy and results of a subset of the ground calibration campaigns, including a discussion of improvements in the energy scale ground calibration compared with Hitomi’s. These improvements include calibration of the low-energy band below 4 keV with the instrument in the flight dewar and the dewar aperture door open, which was not performed for Hitomi, and thorough measurements over an extended high-energy waveband to 22 keV. We also developed an improved technique for gain calibration of “mid-res” secondary events, which have suppressed gain due to proximity to a preceding X-ray event (18 to 70 ms) on the same pixel. We provide a discussion of the on-orbit energy scale monitoring campaigns and an assessment of the Resolve energy scale uncertainties, a key parameter for astrophysics analysis. Energy-scale calibration approaches for future space-based instruments, including the X-ray Integral Field Unit on Athena and microcalorimeter spectrometers proposed or under discussion for future X-ray observatory concepts, have heritage in the calibration of XRISM. We briefly comment on lessons learned from Resolve calibration that are relevant for these future instruments.

Aluminum↗
DOE OSTI2024

Cherenkov photon background for low-noise silicon detectors in space

Future space observatories dedicated to direct imaging and spectroscopy of extra-solar planets will require ultra-low-noise detectors that are sensitive over a broad range of wavelengths. Silicon charge-coupled devices (CCDs), such as EMCCDs, Skipper CCDs, and Multi-Amplifier Sensing CCDs, have demonstrated the ability to detect and measure single photons from ultra-violet to near-infrared wavelengths, making them candidate technologies for this application. In this context, we study a relatively unexplored source of low-energy background coming from Cherenkov radiation produced by energetic charged particles traversing a silicon detector. In the intense radiation environment of space, energetic cosmic rays produce high-energy tracks and more extended halos of low-energy Cherenkov photons, which are detectable with ultra-low-noise detectors. We present a model of this effect that is calibrated to laboratory data, and we use this model to characterize the residual background rate for ultra-low noise silicon detectors in space. We find that the rate of cosmic-ray-induced Cherenkov photon production is comparable to other detector and astrophysical backgrounds that have previously been considered.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Searching for synchrotron emission from the geminga TeV halo using the planck satellite

Pulsars convert a significant fraction of their total spin-down power into very high-energy electrons, leading to the formation of TeV halos. While these halos are well characterized at TeV energies, it remains unclear whether pulsars also accelerate electrons efficiently at lower energies and how these particles propagate through their surrounding environments. We aim to test whether synchrotron emission from $\sim 50$–$300 \, \textrm{GeV}$ electrons around the Geminga pulsar can be detected in the frequency range observed by the Planck satellite. This would help constrain low-energy particle acceleration and diffusion in the vicinity of pulsars. We model the expected synchrotron emission from Geminga’s TeV halo based on various diffusion and injection spectrum scenarios and compare these predictions to publicly available multi-frequency Planck data. We find no conclusive evidence of spatially extended synchrotron emission associated with Geminga in any of Planck’s frequency bands. Our calculations show that even under favorable diffusion and injection conditions, the predicted synchrotron flux lies well below Planck’s measured background levels.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Beyond $\rho ^{2/3}$ scaling: microscopic origins and multimessengers of high-density nuclear symmetry energy

The nature and Equation of State (EOS) of dense neutron-rich matter are still very poorly known, while they have broad impacts on many interesting issues in both astrophysics and nuclear physics. In particular, the nuclear symmetry energy $E_{\textrm{sym}}(\rho )$ encoding the cost to make nuclear matter more neutron-rich has been the most uncertain component of the EOS of dense neutron-rich nucleonic matter. It significantly affects the radii, tidal deformations, cooling rates, and frequencies of various oscillation modes of isolated neutron stars as well as the strain amplitude and frequencies of gravitational waves from their mergers, besides its many effects on the structures of nuclei as well as the dynamics and observables of their collisions. Siemens (1970s) observed that $E_{\textrm{sym}}(\rho )$ scales as $(\rho /\rho _0)^{2/3}$ near the saturation density $\rho _0$ of nuclear matter, since both the kinetic part and the potential contribution (quadratic in momentum) exhibit this dependence. The scaling holds if: (1) the nucleon isoscalar potential is quadratic in momentum, and (2) the isovector interaction is weakly density-dependent. After examining many empirical evidences and understanding theoretical findings in the literature, we conclude that: (1) Siemens’ $\rho ^{2/3}$ scaling is robust and serves as a valuable benchmark for both nuclear theories and experiments up to $2\rho _0$ but breaks down at higher densities, (2) Experimental and theoretical findings about $E_{\textrm{sym}}(\rho )$ up to $2\rho _0$ are broadly consistent, but uncertainties remain large for its curvature $K_{\textrm{sym}}(\rho )$ and higher-order parameters, (3) Above $2\rho _0$ , uncertainties grow due to poorly constrained spin-isospin-dependent tensor and three-body forces as well as the resulting nucleon short-range correlations. Looking forward, combining signals from both observations of neutron stars and terrestrial heavy-ion reaction experiments is the most promising path to finally constraining the high-density $E_{\textrm{sym}}(\rho )$ and the EOS of supradense neutron-rich matter. Multiple examples of community efforts to further constrain the high-density $E_{\textrm{sym}}(\rho )$ using both real and mocked data of present and future high-precision observations of neutron stars, as well as heavy-ion collisions involving high-energy rare isotopes, are briefly reviewed.

Li, Bao-An [East Texas A & M Univ., Commerce, TX (↗
DOE OSTI2026

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Large-scale Cosmic-ray Anisotropies with 19 yr of Data from the Pierre Auger Observatory

Results are presented for the measurement of large-scale anisotropies in the arrival directions of ultra–high-energy cosmic rays detected at the Pierre Auger Observatory during 19 yr of operation, prior to AugerPrime, the upgrade of the observatory. The 3D dipole amplitude and direction are reconstructed above 4 EeV in four energy bins. Besides the established dipolar anisotropy in R.A. above 8 EeV, the Fourier amplitude of the 8–16 EeV energy bin is now also above the 5σ discovery level. No time variation of the dipole moment above 8 EeV is found, setting an upper limit to the rate of change of such variations of 0.3% yr$^{−1}$ at the 95% confidence level. Additionally, the results for the angular power spectrum are shown, demonstrating no other statistically significant multipoles. The results for the equatorial dipole component down to 0.03 EeV are presented, using for the first time a data set obtained with a trigger that has been optimized for lower energies. Finally, model predictions are discussed and compared with observations, based on two source emission scenarios obtained in the combined fit of spectrum and composition above 0.6 EeV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Analysis of Two Models for the Angular Structure of the Outflows Producing the Swift/XRT “Larger-angle Emission” of Gamma-Ray Bursts

The quasi-instantaneous emission from a relativistic surface endowed with a Lorentz factor that decreases away from the outflow symmetry axis can naturally explain the three phases observed by Swift X-Ray Telescope (XRT) in gamma-ray bursts (GRBs) and their afterglows (GRB tail, afterglow plateau, and postplateau) based only on the angular change of the relativistic Doppler boost across the outflow surface. We develop further the analytical formalism of the “larger-angle emission” model for the case of “n-exponential” outflows (where the Lorentz factor Γ dependence of the angular location θ is Γ ~ exp{-(θ/θ c ) n }), and compare its ability to account for the X-ray emission of XRT afterglows relative to that of “power-law” outflows (Γ ∼ θ −g ). Power-law outflows yield longer afterglow plateaus, followed by slower postplateau flux decays than n-exponential outflows, features which may be used in identifying which type of angular structure is at work in a given afterglow. Identifying the Γ(θ) angular structure that accommodates XRT light curves is slightly complicated by the fact that the afterglow X-ray light curve is also determined by how two characteristics of the comoving-frame emission spectrum (peak energy $E'_p$ and peak intensity $i'_p$) change with the angular location or, equivalently, with the Lorentz factor. Here, we assume power-law Γ dependences of those spectral characteristics and find that, unlike power-law outflows, n-exponential outflows cannot account for plateaus with a temporal dynamical range larger than 100 (2 dex in logarithmic space). To capture all the information contained in XRT afterglow measurements (0.3–10 keV unabsorbed flux and effective spectral slope), we calculate 0.3 and 10 keV light curves using a broken-power-law emission spectrum of peak energy and low- and high-energy slopes that are derived from the effective slope measured by XRT. This economical peak energy determination is found to be consistent with the results of more expensive spectral fits. The angular distributions of the Lorentz factor, comoving frame peak energy, and peak intensity (Γ(θ), $E'_p$(θ), $i'_p$(θ)) constrain the (yet-to-be determined) convolution of various features of the production of relativistic jets by solar-mass black holes and of their propagation through the progenitor/circumburst medium, while the $E'_p$(Γ) and $i'_p$(Γ) dependences may constrain the GRB dissipation mechanism and the GRB emission process.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

IceCube Search for Neutrino Emission from X-Ray Bright Seyfert Galaxies

The recent IceCube detection of TeV neutrino emission from the nearby active galaxy NGC 1068 suggests that active galactic nuclei (AGNs) could make a sizable contribution to the diffuse flux of astrophysical neutrinos. The absence of TeV γ-rays from NGC 1068 indicates neutrino production in the vicinity of the supermassive black hole, where the high radiation density leads to γ-ray attenuation. Therefore, any potential neutrino emission from similar sources is not expected to correlate with high-energy γ-rays. Disk-corona models predict neutrino emission from Seyfert galaxies to correlate with keV X-rays because they are tracers of coronal activity. Using through-going track events from the Northern Sky recorded by IceCube between 2011 and 2021, we report results from a search for individual and aggregated neutrino signals from 27 additional Seyfert galaxies that are contained in the Swift's Burst Alert Telescope AGN Spectroscopic Survey. Besides the generic single power law, we evaluate the spectra predicted by the disk-corona model assuming stochastic acceleration parameters that match the measured flux from NGC 1068. Assuming all sources to be intrinsically similar to NGC 1068, our findings constrain the collective neutrino emission from X-ray bright Seyfert galaxies in the northern sky, but, at the same time, show excesses of neutrinos that could be associated with the objects NGC 4151 and CGCG 420-015. These excesses result in a 2.7σ significance with respect to background expectations.

Abbasi, R. [Loyola University, Chicago, IL (United↗
DOE OSTI2025

Energy Spectrum and Mass Composition of Ultra-high-energy Cosmic Rays Originating from Relativistic Jets of Nearby Radio Galaxies

Relativistic jets of radio galaxies (RGs) are possible sources of ultra-high-energy cosmic rays (UHECRs). Recent studies combining relativistic hydrodynamic simulations with Monte Carlo particle transport have demonstrated that UHECRs can be accelerated to energies beyond 10 20 eV through shocks, turbulence, and relativistic shear in jet-induced flows of Fanaroff–Riley type RGs. The resulting time-asymptotic UHECR spectrum is well modeled by a double power law with an “extended” exponential cutoff, primarily shaped by relativistic shear acceleration. In this study, we adopt this novel source spectrum and simulate the propagation of UHECRs from nearby RGs using the CRPropa code. We focus on Virgo A (Vir A), Centaurus A (Cen A), Fornax A (For A), and Cygnus A (Cyg A), expected to be the most prominent UHECR sources among RGs. We then analyze the energy spectrum and mass composition of UHECRs arriving at Earth. We find that, due to the extended high-energy tail in the source spectrum, UHECRs from Vir A, which has a higher Lorentz factor, exhibit a higher flux at the highest energies and a lighter mass composition at Earth compared to those from Cen A and For A with lower Lorentz factors. Despite Cyg A having an even higher Lorentz factor, the large distance limits its contribution. With a small number of nearby prominent RGs, our findings suggest that if RGs are the major sources of UHECRs, the energy spectrum and mass composition of observed UHECRs would exhibit hemispheric differences between the Northern and Southern skies at the highest energies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Nonthermal Ionization of Kilonova Ejecta: Observable Impacts

The characteristic rapid rise and decline at optical wavelengths of a kilonova is the product of the low ejecta mass (≲0.05 M ⊙ ) and high ejecta velocity (≳0.1c). We show that, even at very early times (≲2 days), regions of ejecta fall below critical density and temperature thresholds at which nonlocal thermodynamic equilibrium (NLTE) effects become important. Here, we present an approximate method for calculating the ionization state of the ejecta that accounts for the NLTE impact of high-energy electrons produced in the beta decay of freshly synthesized r-process elements. We find that incorporating ionization from high-energy electrons produces an “inverted” and “blended” ionization structure, where the most highly ionized species are located in the fastest moving homologous ejecta, and multiple ionization states coexist. In radiation transport calculations, the higher degree of ionization reduces line blanketing in optical bands, leading to improved agreement with the light-curve properties of AT 2017gfo, such as the duration, decay rates, brightness, and colors. Our quasi-NLTE implementation helps to alleviate tensions in kilonova modeling: for high-velocity (∼0.3 c) ejecta components, our models require less mass for a given peak brightness in optical bands, by as much as a factor of 3; our models can explain the presence of observed features associated to Sr II, W III, Se III, and Te III under conditions where LTE models would predict only neutral species; and we naturally predict the coexistence of species like Sr II and Ce III without the need for fine-tuning of the ejecta properties.

High energy astrophysics↗
DOE OSTI2026

The Extent of Solar Energetic Particle Irradiation in the Sun’s Protoplanetary Disk

Solar flares emit X-rays and high-energy (MeV to GeV) ions called solar energetic particles (SEPs). Astronomical observations show solar-mass protostellar fluxes are a factor Φ ≈ 3 × 10 2 –3 × 10 3 times higher than the present-day Sun. Constraining Φ in the early solar system is important for modeling ionization in the Sun’s protoplanetary disk, the extent of magnetorotational instability or magnetocentrifugal outflows, or even production of short-lived radionuclides. Recent interpretations of meteoritic data—including cosmogenic Ne in hibonite grains, initial ( 10 Be/ 9 Be) 0 ratios in Ca-rich, Al-rich inclusions (CAIs), or even inferences of live 7 Be in CAIs—suggest values Φ > 10 5 , reaching as high as Φ ≈ 6 × 10 6 , which would make the young Sun extraordinarily active, even for a protostar. We constrain Φ by reexamining these data. We conclude that cosmogenic Ne was produced in hibonite grains as they resided in the disk; 36 Cl was created in Cl-poor grains after the disk dissipated; 10 Be was inherited from the molecular cloud, with almost no (<1%) 10 Be created in the disk; and there is no evidence whatsoever for any live 7 Be in CAIs. We show these data are consistent with a value Φ ≈ 3 × 10 3 for the first >5 Myr of the solar nebula. The early Sun evidently emitted a flux of X-rays and SEPs not atypical for a protostar.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Characterizing the γ-Ray Emission from Low-luminosity AGN

A majority of the active galactic nuclei (AGN) in the local Universe are classified as low-luminosity AGN (LLAGN), having bolometric luminosities ≲10 42 erg s −1 . Although high-energy γ-ray emission is predicted from both the jets and disks of LLAGN, to date only four have been detected by the Fermi Large Area Telescope (Fermi-LAT). In this work, we therefore conduct a comprehensive study of all the LLAGN from the Palomar spectroscopic survey of bright, northern galaxies, including both subthreshold and detected γ-ray sources, using 14.4 yr of LAT data. Our analysis results in a new detection of one LLAGN, as well as a detection of the subthreshold population using a stacking technique. We find that the signal from the subthreshold sample is consistent with being dominated by star formation activity, although a contribution from compact jets or a mixed contribution from jetted and nonjetted systems is also feasible. On the other hand, the individually detected LLAGN are likely dominated by jet emission. We perform detailed spectral modeling for a subset of these sources and find that the γ-ray signal can be explained by synchrotron self-Compton radiation if the inner jet emission region is weakly magnetized, with its total energy density being strongly particle dominated and only slowly moving. With this work we also publicly release our Python-based stacking library for analyzing subthreshold source populations with the LAT based on a proven technique used in numerous studies.

79 ASTRONOMY AND ASTROPHYSICS↗