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

A Beamdump facility at Jefferson Lab

The potential of the intense secondary muon, neutrino, and (hypothetical) light dark matter beams at the Thomas Jefferson National Accelerator Facility (Jefferson Lab) is explored. These are produced in the high-power dumps with high-current electron beams. Light dark matter searches with the approved Beam Dump eXperiment (BDX) are driving the realization of a new underground vault behind Hall A that could be extended to a Beamdump Facility with little additional installations. High-energy muons created via the Bethe–Heitler process uniquely do not proceed through the more common pion production and decay channels. Several possible muon physics applications are highlighted. Neutrino detector technologies and experiments suitable for a beamdump facility are outlined.

Accelerator Physics↗
DOE OSTI2025

X-Ray Polarization of the Magnetar 1E 1841−045

We report on IXPE and NuSTAR observations beginning 40 days after the 2024 outburst onset of magnetar 1E 1841−045, marking the first IXPE observation of a magnetar in an enhanced state. Our spectropolarimetric analysis indicates that both a blackbody (BB) plus double power-law (PL) and a double blackbody plus power-law spectral model fit the phase-averaged intensity data well, with a hard PL tail (Γ = 1.19 and 1.35, respectively) dominating above ≈5 keV. For the former model, we find the soft PL (the dominant component at soft energies) exhibits a polarization degree (PD) of ≈30% while the hard PL displays a PD of ≈40%. Similarly, the cool BB of the 2BB+PL model possesses a PD of ≈15% and a hard PL PD of ≈57%. For both models, each component has a polarization angle (PA) compatible with celestial north. Model-independent polarization analysis supports these results, wherein the PD increases from ≈15% to ≈70% in the 2–3 keV and 6–8 keV ranges, respectively, while the PA remains nearly constant. We find marginal evidence for phase-dependent variability of the polarization properties, namely a higher PD at phases coinciding with the hard X-ray pulse peak. We compare the hard X-ray PL to the expectation from resonant inverse Compton scattering (RICS) and secondary pair cascade synchrotron radiation from primary high-energy RICS photons; both present reasonable spectropolarimetric agreement with the data, albeit the latter does so more naturally. We suggest that the soft PL X-ray component may originate from a Comptonized corona in the inner magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

High Optical-to-X-Ray Polarization Ratio Reveals Compton Scattering in BL Lacertae’s Jet

Blazars, supermassive black hole systems with highly relativistic jets aligned with the line of sight, are the most powerful long-lived emitters of electromagnetic emission in the Universe. We report here on a radio-to-gamma-ray multiwavelength campaign on the blazar BL Lacertae with unprecedented polarimetric coverage from radio to X-ray wavelengths. The observations caught an extraordinary event on 2023 November 10–18, when the degree of linear polarization of optical synchrotron radiation reached a record value of 47.5%. In stark contrast, the Imaging X-ray Polarimetry Explorer found that the X-ray (Compton scattering or hadron-induced) emission was polarized at less than 7.4% (3 σ confidence level). We argue here that this observational result rules out a hadronic origin of the high-energy emission and strongly favors a leptonic (Compton scattering) origin, thereby breaking the degeneracy between hadronic and leptonic emission models for BL Lacertae and demonstrating the power of multiwavelength polarimetry to address this question. Furthermore, the multiwavelength flux and polarization variability, featuring an extremely prominent rise and decay of the optical polarization degree, is interpreted for the first time by the relaxation of a magnetic “spring” embedded in the newly injected plasma. This suggests that the plasma jet can maintain a predominant toroidal magnetic field component parsecs away from the central engine.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Global tuning of hadronic interaction models with accelerator-based and astroparticle data

In high-energy and astroparticle physics, event generators play an essential role, even in the simplest data analyses. As analysis techniques become more sophisticated, e.g. based on deep neural networks, their correct description of the observed event characteristics becomes even more important. Physical processes occurring in hadronic collisions are simulated within a Monte Carlo framework. A major challenge is the modeling of hadron dynamics at low momentum transfer, which includes the initial and final phases of every hadronic collision. QCD-inspired phenomenological models used for these phases cannot guarantee completeness or correctness over the full phase space. These models usually include parameters which must be tuned to suitable experimental data. Until now, event generators have been developed and tuned mainly on the basis of data from high-energy physics experiments at accelerators. The wealth of data available from the latest generation of astroparticle experiments has not yet been fully exploited, and in many cases is not satisfactorily described. Both kinds of data sets are complementary as astroparticle experiments provide sensitivity especially to hadrons produced nearly parallel to the collision axis and cover center-of-mass energies up to several hundred TeV, well beyond those reached at colliders so far. In this report, we provide an overview of state-of-the-art event generators and their tuning, including the most relevant inputs from high-energy accelerator and astroparticle experiments. We present a road map that shows, for the first time, how the unified tuning of event generators with accelerator-based and astroparticle data can be performed.

Albrecht, J. [Ruhr U., Bochum, RAPP Ctr.; Ruhr U.,↗
DOE OSTI2025

Delving into the depths of NGC 3783 with XRISM II. Cross-calibration of X-ray instruments used in the large, multi-mission observational campaign

Context. Accurate X-ray spectroscopic measurements are fundamental for deriving basic physical parameters of the most abundant baryon components in the Universe. The plethora of X-ray observatories currently operational enables a panchromatic view of the high-energy emission of celestial sources. However, uncertainties in the energy-dependent calibration of the instrument transfer functions (e.g. the effective area, energy redistribution, or gain) can limit - and historically, did limit - the accuracy of X-ray spectroscopic measurements. Aims. We revised the status of the cross-calibration among the scientific payload on board four operation missions: Chandra, NuSTAR, XMM-Newton, and the recently launched XRISM. XRISM carries the micro-calorimeter Resolve, which yields the best energy resolution at energies ≥2 keV. For this purpose, we used the data from a 10-day-long observational campaign targeting the nearby active galactic nucleus NGC 3783, carried out in July 2024. Methods. We present a novel model-independent method for assessing the cross-calibration status that is based on a multi-node spline of the spectra with the highest-resolving power (XRISM/Resolve in our campaign). We also estimated the impact of the intrinsic variability of NGC 3783 on the cross-calibration status due to the different time coverages of participating observatories and performed an empirical reassessment of the Resolve throughput at low energies. Results. Based on this analysis, we derived a set of energy-dependent correction factors of the observed responses, enabling a statistically robust analysis of the whole spectral dataset. They will be employed in subsequent papers describing the astrophysical results of the campaign. Aims. We revised the status of the cross-calibration among the scientific payload on board four operation missions: Chandra, NuSTAR, XMM-Newton, and the recently launched XRISM. XRISM carries the micro-calorimeter Resolve, which yields the best energy resolution at energies ≥2 keV. For this purpose, we used the data from a 10-day-long observational campaign targeting the nearby active galactic nucleus NGC 3783, carried out in July 2024. Methods. We present a novel model-independent method for assessing the cross-calibration status that is based on a multi-node spline of the spectra with the highest-resolving power (XRISM/Resolve in our campaign). We also estimated the impact of the intrinsic variability of NGC 3783 on the cross-calibration status due to the different time coverages of participating observatories and performed an empirical reassessment of the Resolve throughput at low energies. Results. Based on this analysis, we derived a set of energy-dependent correction factors of the observed responses, enabling a statistically robust analysis of the whole spectral dataset. They will be employed in subsequent papers describing the astrophysical results of the campaign.

Active Galactic Nuclei, individual: NGC 3783↗
DOE OSTI2024

Establishing accretion flares from supermassive black holes as a source of high-energy neutrinos

ABSTRACT The origin of cosmic high-energy neutrinos remains largely unexplained. For high-energy neutrino alerts from IceCube, a coincidence with time-variable emission has been seen for three different types of accreting black holes: (1) a gamma-ray flare from a blazar (TXS 0506+056), (2) an optical transient following a stellar tidal disruption event (TDE; AT2019dsg), and (3) an optical outburst from an active galactic nucleus (AGN; AT2019fdr). For the latter two sources, infrared follow-up observations revealed a powerful reverberation signal due to dust heated by the flare. This discovery motivates a systematic study of neutrino emission from all supermassive black hole with similar dust echoes. Because dust reprocessing is agnostic to the origin of the outburst, our work unifies TDEs and high-amplitude flares from AGN into a population that we dub accretion flares. Besides the two known events, we uncover a third flare that is coincident with a PeV-scale neutrino (AT2019aalc). Based solely on the optical and infrared properties, we estimate a significance of 3.6σ for this association of high-energy neutrinos with three accretion flares. Our results imply that at least ∼10 per cent of the IceCube high-energy neutrino alerts could be due to accretion flares. This is surprising because the sum of the fluence of these flares is at least three orders of magnitude lower compared to the total fluence of normal AGN. It thus appears that the efficiency of high-energy neutrino production in accretion flares is increased compared to non-flaring AGN. We speculate that this can be explained by the high Eddington ratio of the flares.

Astronomy & Astrophysics↗
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 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 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 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↗