Search NASA⌕ Search

SEARCH · Search NASA

Results for “Astronomy data modeling”

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

Non-propagating structures and propagating waves in solar wind turbulence revealed by simulations and observations

Structures and waves are common features of solar wind turbulence at various scales. The interplay between structures and waves is important for processes such as the turbulent energy cascade, plasma heating, and particle scattering. Our understanding of turbulence has been advanced by not only new space missions and numerical simulations, but also techniques that have been developed to interpret the rapidly growing turbulence data. We review basic models of turbulence with a specific focus on the analysis methods for understanding magnetic structures and waves. MHD and kinetic waves in single-spacecraft time series measurements can be identified through mode decomposition or their characteristic polarization signatures. The structures in this paper are considered as zero-frequency, non-propagating or convected modes embedded in the solar wind. The synergy between observations and simulations is most evident in the application of spatial-temporal analysis to multi-spacecraft observation and turbulence simulations. The spatial-temporal analysis has greatly improved our understanding of structures and waves in turbulence. We conclude by discussing prospects for future research.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Cosmological implications of DESI DR2 BAO measurements in light of the latest ACT DR6 CMB data

We report cosmological results from the Dark Energy Spectroscopic Instrument (DESI) measurements of baryon acoustic oscillations (BAO) when combined with recent data from the Atacama Cosmology Telescope (ACT). By jointly analyzing ACT and Planck data and applying conservative cuts to overlapping multipole ranges, we assess how different 𝑃⁢𝑙⁢𝑎⁢𝑛⁢𝑐⁢𝑘 + ACT dataset combinations affect consistency with DESI. While ACT alone exhibits a tension with DESI exceeding 3⁢𝜎 within the Λ⁢ CDM model, this discrepancy is reduced when ACT is analyzed in combination with Planck. For our baseline DESI DR2 BAO + 𝑃⁢𝑙⁢𝑎⁢𝑛⁢𝑐⁢𝑘 P⁢R⁢4 + ACT likelihood combination, the preference for evolving dark energy over a cosmological constant is about 3⁢𝜎, increasing to over 4⁢𝜎 with the inclusion of type Ia supernova data. While the dark energy results remain quite consistent across various combinations of Planck and ACT likelihoods with those obtained by the DESI collaboration, the constraints on neutrino mass are more sensitive, ranging from ∑𝑚 𝜈 < 0.061 eV in our baseline analysis, to ∑𝑚 𝜈 < 0.077 eV (95% confidence level) in the CMB likelihood combination chosen by ACT when imposing the physical prior ∑𝑚 𝜈 > 0 eV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Model-independent measurement of the matter-radiation equality scale in DESI 2024

The peak of the matter power spectrum, known as the turnover (TO) scale, is determined by the horizon size at the time of matter-radiation equality. This scale can serve as a standard ruler, independent of other features in the matter power spectrum, such as baryon acoustic oscillations (BAO). Here, in this work, we present the first detection of the turnover in the galaxy autopower spectrum, utilizing the distribution of quasars and luminous red galaxies (LRG) measured by the Dark Energy Spectroscopic Instrument (DESI) during its first year of survey operations in a model-independent manner. To avoid confirmation bias, we first analyze the data using data blinding methods designed for the DESI baryon acoustic oscillation, redshift space distortion and scale-dependent bias signals. We measure the angle-averaged dilation distance 𝐷 V ⁡(𝑧 = 1.651) = (38.1 ± 2.5)⁢𝑟 H from the quasars and 𝐷 V ⁡(𝑧 = 0.733) = (21.8 ± 1.0)⁢𝑟 H from the LRG sample in units of the horizon 𝑟 H at the matter-radiation-equality epoch. Combining these two constraints and assuming a flat Λ⁢CDM model with three standard neutrino species, we can translate this into a constraint of Ω m ⁢ℎ 2 = 0.13⁢9$^{+0.036}_{−0.046}$. We can break the Ω m −𝐻 0 degeneracy with low-redshift distance measurements from type-Ia supernova (SN) data from Pantheon+, we obtain a sound-horizon free estimate of the Hubble-Lemaître parameter of 𝐻 0 = 65.2$^{+4.9}_{−6.2}$ km/s/Mpc, consistent with sound-horizon dependent DESI measurements. On the other hand, combining the DESI BAO and TO, we find a truly DESI-only measurement of 𝐻 0 = 74.0$^{+7.2}_{−3.5}$ km/s/Mpc, in line with DESI-only full-shape results where the sound-horizon scale is marginalized out. This discrepancy in 𝐻 0 can be reconciled in a 𝑤 0 ⁢𝑤 𝑎 ⁢ CDM cosmology, where the combination of DESI BAO and TO data yields 𝐻 0 = 66.5 ± 7.2 km/s/Mpc.

cosmological parameters↗
DOE OSTI2024

Detection of the large-scale tidal field with galaxy multiplet alignment in the DESI Y1 spectroscopic survey

We explore correlations between the orientations of small galaxy groups, or ‘multiplets’, and the large-scale gravitational tidal field. Using data from the Dark Energy Spectroscopic Instrument (DESI) Y1 survey, we detect the intrinsic alignment (IA) of multiplets to the galaxy-traced matter field out to separations of $100\,h^{-1}$ Mpc. Unlike traditional IA measurements of individual galaxies, this estimator is not limited by imaging of galaxy shapes and allows for direct IA detection beyond redshift $z=1$. Multiplet alignment is a form of higher order clustering, for which the scale-dependence traces the underlying tidal field and amplitude is a result of small-scale ($\lt 1h^{-1}$ Mpc) dynamics. Within samples of bright galaxies, luminous red galaxies (LRG) and emission-line galaxies, we find similar scale-dependence regardless of intrinsic luminosity or colour. This is promising for measuring tidal alignment in galaxy samples that typically display no IA. DESI’s LRG mock galaxy catalogues created from the A BACUS S UMMIT N -body simulations produce a similar alignment signal, though with a 33 per cent lower amplitude at all scales. An analytic model using a non-linear power spectrum (NLA) only matches the signal down to 20 $h^{-1}$ Mpc. Our detection demonstrates that galaxy clustering in the non-linear regime of structure formation preserves an interpretable memory of the large-scale tidal field. Multiplet alignment complements traditional two-point measurements by retaining directional information imprinted by tidal forces, and contains additional line-of-sight information compared to weak lensing. This is a more effective estimator than the alignment of individual galaxies in dense, blue, or faint galaxy samples.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

New dark matter analysis of milky way dwarf satellite galaxies with madhatv2

We obtain bounds on dark matter annihilation using 14 years of publicly available Fermi-LAT data from a set of 54 dwarf spheroidal galaxies, using spectral information from 16 energy bins. We perform this analysis using our updated and publicly available code , which can be used to test a variety of models for dark matter particle physics and astrophysics in an accessible manner. In particular, we note that including Carina III in the analysis strengthens constraints on s -wave annihilation into two-body Standard Model final states by a factor of ∼ 3 but broadens the error on the constraint due to the large uncertainty of its J -factor. Our findings illustrate the importance of verifying if Carina III is in fact a dwarf spheroidal galaxy and measuring more precisely its J -factor. More generally, they highlight the significance of forthcoming discoveries of nearby ultrafaint dwarfs for dark matter indirect detection. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2025

Tools for probing new physics with newly discovered gamma-ray targets

Here, we present a computational tool, TweedleDEE, for empirically modeling diffuse gamma-ray background emission in a 1° region of the sky, using publicly available gamma-ray data off-axis from the region of interest. This background model allows a user to perform a purely data-driven search for anomalous localized sources of gamma-ray emission, including new physics. A major application of this tool would be in searching for dark matter annihilation in newly discovered astrophysical targets. For this purpose, we derive a scaling relation for determining velocity-dependent 𝐽-factors using only the stellar parameters, which can be broadly applied to obtain dark matter constraints from new targets. As an application of these tools, we use TweedleDEE and MADHATv2 to perform the first search for dark matter annihilation in the newly discovered Leo VI dwarf spheroidal galaxy, and present model constraints for a variety of choices of the annihilation channel and velocity dependence of the cross section.

dark matter↗
DOE OSTI2026

A New Measurement of the Extragalactic Background Light Using 15 yr of Fermi-large Area Telescope Data

The extragalactic background light (EBL) from ultraviolet to infrared comprises the emission from all stars, galaxies, and actively accreting black holes in the observable Universe. A precise measurement of the EBL is critically important to probe models of star formation and galaxy evolution. The EBL can be measured via the absorption imprint left on the spectra of γ-ray blazars. In this work, we rely on 15 yr of Fermi-LAT data and 1576 blazars to measure the EBL optical depth in the 0 < z < 4.3 range. We detect the EBL attenuation with ≅23σ significance and measure the optical depth in 19 redshift bins, extending the coverage and improving on our previous results. This allows us to reconstruct the evolution of the EBL and find general consistency with recent EBL models. These results represent the most precise determination of the EBL with GeV γ rays to date.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Cosmology of dark energy radiation

In this work, we quantify the cosmological signatures of dark energy radiation—a novel description of dark energy, which proposes that the dynamical component of dark energy is comprised of a thermal bath of relativistic particles sourced by thermal friction from a slowly rolling scalar field. For a minimal model with particle production emerging from first principles, we find that the abundance of radiation sourced by dark energy can be as large as Ω DER = 0.03, exceeding the bounds on relic dark radiation by three orders of magnitude. Although the background and perturbative evolution of dark energy radiation are distinct from Quintessence, we find that current and near-future cosmic microwave background and supernova data will not distinguish these models of dark energy. We also find that our constraints on all models are dominated by their impact on the expansion rate of the Universe. Considering extensions that allow the dark radiation to populate neutrinos, axions, and dark photons, we evaluate the direct detection prospects of a thermal background comprised of these candidates consistent with cosmological constraints on dark energy radiation. Our study indicates that a resolution of ~6 meV is required to achieve sensitivity to relativistic neutrinos compatible with dark energy radiation in a neutrino capture experiment on tritium. We also find that dark matter axion experiments lack sensitivity to a relativistic thermal axion background, even if enhanced by dark energy radiation, and dedicated search strategies are required to probe new parameter space. We derive constraints arising from a dark photon background from oscillations into visible photons, and find that viable parameter space can be explored with the late dark energy radiation experiment.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Dark Energy Survey Year 6 results: cell-based coadds and METADETECTION weak lensing shape catalogue

We present the metadetection weak lensing galaxy shape catalogue from the 6-yr Dark Energy Survey (DES Y6) imaging data. This data set is the final release from DES, spanning 4422 deg 2 of the southern sky. We describe how the catalogue was constructed, including the two new major processing steps, cell-based image coaddition, and shear measurements with metadetection. The DES Y6 M etadetection weak lensing shape catalogue consists of 151 922 791 galaxies detected over riz bands, with an effective number density of n eff = 8.22 galaxies per arcmin 2 and shape noise of σ e = 0.29. We carry out a suite of validation tests on the catalogue, including testing for point spread function (PSF) leakage, testing for the impact of PSF modelling errors, and testing the correlation of the shear measurements with galaxy, PSF, and survey properties. In addition to demonstrating that our catalogue is robust for weak lensing science, we use the DES Y6 image simulation suite to estimate the overall multiplicative shear bias of our shear measurement pipeline. We find no detectable multiplicative bias at the roughly half-per cent level, with m = (3.4 ± 6.1) x 10 –3 , at 3σ uncertainty. This is the first time both cell-based coaddition and Metadetection algorithms are applied to observational data, paving the way to the Stage-IV weak lensing surveys.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

The DEHVILS in the details: Type Ia supernova Hubble residual comparisons and mass step analysis in the near-infrared

Measurements of type Ia supernovae (SNe Ia) in the near-infrared (NIR) have been used both as an alternate path to cosmology compared to optical measurements and as a method of constraining key systematics for the larger optical studies. With the DEHVILS sample, the largest published NIR sample with consistent NIR coverage of maximum light across three NIR bands ( Y, J , and H ), we check three key systematics: (i) the reduction in Hubble residual scatter as compared to the optical, (ii) the measurement of a “mass step” or lack thereof and its implications, and (iii) the ability to distinguish between various dust models by analyzing slopes and correlations between Hubble residuals in the NIR and optical. We produce SN Ia simulations of the DEHVILS sample and find that it is harder to differentiate between various dust models than previously understood. Additionally, we find that fitting with the current SALT3-NIR model does not yield accurate wavelength-dependent stretch-luminosity correlations, and we propose a limited solution for this problem. From the data, we see that (i) the standard deviation of Hubble residual values from NIR bands treated as standard candles are 0.007–0.042 mag smaller than those in the optical, (ii) the NIR mass step is not constrainable with the current sample size of 47 SNe Ia from DEHVILS, and (iii) Hubble residuals in the NIR and optical are correlated in the data. We test a few variations on the number and combinations of filters and data samples, and we observe that none of our findings or conclusions are significantly impacted by these modifications.

Astronomy & Astrophysics↗
DOE OSTI2026

The Thermal and Kinematic Sunyaev–Zeldovich Effect in Galaxy Clusters and Filaments Using Multifrequency Temperature Maps of the Cosmic Microwave Background: A399–A401 Cluster Pair Case Study

We present a multifrequency and multi-instrument methodology to study the physical properties of galaxy clusters and cosmic filaments using cosmic microwave background observations. Our approach enables simultaneous measurement of both the thermal (tSZ) and kinematic Sunyaev–Zeldovich (kSZ) effects, incorporates relativistic corrections, and models astrophysical foregrounds such as thermal dust emission. We do this by jointly fitting a single physical model across multiple maps from multiple instruments at different frequencies, rather than fitting a model to a single Compton-y map. We demonstrate the success of this method by fitting the A399–A401 galaxy cluster pair and filament system using archival data from the Planck satellite and new, targeted deep data from the Atacama Cosmology Telescope, covering 11 different frequencies over 14 maps from 30 GHz to 545 GHz. Our tSZ results are consistent with previous work using Compton-y maps. We measure the line-of-sight peculiar velocities of the cluster–filament system using the kSZ effect and find statistical uncertainties on individual cluster peculiar velocities of ≲600 km s −1 , which are competitive with current state-of-the-art measurements. Additionally, we measure the optical depth of the filament component with a signal-to-noise of 8.5σ and reveal hints of its morphology. This modular approach is well-suited for application to future instruments across a wide range of millimeter and submillimeter wavebands.

Gill, Ajay S. [National Research Council of Canada↗
DOE OSTI2025

One Million Open-source Cislunar Orbits

Cislunar space, encompassing the region from geosynchronous orbit to beyond the Moon, is poised to become a cornerstone for future exploration, scientific discovery, and national security. Missions in this region, spanning durations from weeks to decades, require robust infrastructure and reliable transit capabilities. The complex gravitational influences of the Moon, Sun, and planets, along with thermal radiation from Earth and the Sun, lead to significant trajectory deviations, resulting in kilometer-scale errors within days. Leveraging the high-performance computing resources at Lawrence Livermore National Laboratory (LLNL), we have simulated one million high-fidelity cislunar trajectories, now publicly available via LLNL’s Green Data Oasis and the Unified Data Library. Generated using the open-source Space Situational Awareness Python package, these trajectories match the precision of commercial tools such as AGI’s Systems Tool Kit and NASA’s General Mission Analysis Tool. This data set is a valuable resource for reference, statistical analysis of cislunar orbit populations, and training machine learning models for rapid orbit classification with minimal observational input. Preliminary analysis reveals stable bands in Keplerian element space, particularly around five geosynchronous radii across a range of inclinations and eccentricities. Beyond this threshold, the Moon’s influence disrupts most unassisted orbits, though co-orbiting L4/L5 Lunar Trojans persist throughout the six-year simulation.

Astronomy and AstroPhysics↗
DOE OSTI2024

A More Precise Measurement of the Radius of PSR J0740+6620 Using Updated NICER Data

PSR J0740+6620 is the neutron star with the highest precisely determined mass, inferred from radio observations to be 2.08 ± 0.07 M ⊙ . Measurements of its radius therefore hold promise to constrain the properties of the cold, catalyzed, high-density matter in neutron star cores. Previously, Miller et al. and Riley et al. reported measurements of the radius of PSR J0740+6620 based on Neutron Star Interior Composition Explorer (NICER) observations accumulated through 2020 April 17, and an exploratory analysis utilizing NICER background estimates and a data set accumulated through 2021 December 28 was presented in Salmi et al. Here we report an updated radius measurement, derived by fitting models of X-ray emission from the neutron star surface to NICER data accumulated through 2022 April 21, totaling ~1.1 Ms additional exposure compared to the data set analyzed in Miller et al. and Riley et al., and to data from XMM-Newton observations. We find that the equatorial circumferential radius of PSR J0740+6620 is ${12.92}_{-1.13}^{+2.09}$ km (68% credibility), a fractional uncertainty ~83% the width of that reported in Miller et al., in line with statistical expectations given the additional data. If we were to require the radius to be less than 16 km, as was done in Salmi et al., then our 68% credible region would become $R={12.76}_{-1.02}^{+1.49}$ km, which is close to the headline result of Salmi et al. Our updated measurements, along with other laboratory and astrophysical constraints, imply a slightly softer equation of state than that inferred from our previous measurements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

SIDDA: SInkhorn Dynamic Domain Adaptation for image classification with equivariant neural networks

Modern neural networks (NNs) often do not generalize well in the presence of a ‘covariate shift’; that is, in situations where the training and test data distributions differ, but the conditional distribution of classification labels given the data remains unchanged. In such cases, NN generalization can be reduced to a problem of learning more robust, domain-invariant features. Domain adaptation (DA) methods include a broad range of techniques aimed at achieving this; however, these methods have struggled with the need for extensive hyperparameter tuning, which then incurs significant computational costs. In this work, we introduce SInkhorn Dynamic Domain Adaptation (SIDDA), an out-of-the-box DA training algorithm built upon the Sinkhorn divergence, that can achieve effective domain alignment with minimal hyperparameter tuning and computational overhead. We demonstrate the efficacy of our method on multiple simulated and real datasets of varying complexity, including simple shapes, handwritten digits, real astronomical observations, and remote sensing data. These datasets exhibit covariate shifts due to noise, blurring, differences between telescopes, and variations in imaging wavelengths. SIDDA is compatible with a variety of NN architectures, and it works particularly well in improving classification accuracy and model calibration when paired with symmetry-aware equivariant NNs (ENNs). We find that SIDDA consistently enhances the generalization capabilities of NNs, achieving up to a ${\approx}40\%$ improvement in classification accuracy on unlabeled target data, while also providing a more modest performance gain of $\lesssim 1\%$ on labeled source data. We also study the efficacy of DA on ENNs with respect to the varying group orders of the dihedral group DN, and find that the model performance improves as the degree of equivariance increases. Finally, if SIDDA achieves proper domain alignment, it also enhances model calibration on both source and target data, with the most significant gains in the unlabeled target domain—achieving over an order of magnitude improvement in the expected calibration error and Brier score. SIDDA’s versatility across various NN models and datasets, combined with its automated approach to domain alignment, has the potential to significantly advance multi-dataset studies by enabling the development of highly generalizable models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Kilonovae and Long-duration Gamma-Ray Bursts

Recent detections of kilonova-like emission following long-duration gamma-ray bursts GRB 211211A and GRB 230307A have been interpreted as originating from the merger of two neutron stars. In this work, we demonstrate that these observations are also consistent with nucleosynthesis originating from a collapsar scenario. Our model is capable of reproducing the observed optical and infrared light curves using a single, weak r-process component. The absence of lanthanide-rich material in our model, consistent with the data, challenges the prevailing interpretation that a red evolution in such transients necessarily indicates the presence of heavy r-process elements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Four Years of Type Ia Supernovae Observed by TESS: Early-time Light-curve Shapes and Constraints on Companion Interaction Models

Abstract We present 307 type Ia supernova (SN) light curves from the first 4 yr of the Transiting Exoplanet Survey Satellite mission. We use this sample to characterize the shapes of the early-time light curves, measure the rise times from first light to peak, and search for companion star interactions. Using simulations, we show that light curves must have noise <10% of the peak flux to avoid biases in the early-time light-curve shape, restricting our quantitative analysis to 74 light curves. We find that the mean power-law index t β 1 of the early-time light curves isβ 1 = 1.93 ± 0.57, and the mean rise time to peak is 15.7 ± 3.5 days. The underlying population distribution forβ 1 may instead consist of a Gaussian component with mean 2.29, width 0.34, and a long tail extending to values less than 1.0. We find that the data can rarely distinguish between models with and without companion interaction models. Nevertheless, we find three high-quality light curves that tentatively prefer the addition of a companion interaction model, but the statistical evidence for the companion interactions is not robust. We also find two SNe that disfavor the addition of a companion interaction model to a curved power-law model. Taking the 74 SNe together, we calculate 3σupper limits on the presence of companion signatures to control for orientation effects that can hide companions in individual light curves. Our results rule out common progenitor systems with companions having Roche lobe radii >31R ⊙ (separations >5.7 × 10 12 cm, 99.9% confidence level) and disfavor companions having Roche lobe radii >10R ⊙ (separations >1.9 × 10 12 cm, 95% confidence level). Lastly, we discuss the implications of our results for the intrinsic fraction of single degenerate progenitor systems.

Astronomy & Astrophysics↗
DOE OSTI2025

An Exceptionally Powerful, Radiatively Driven Ultrafast Outflow in the Rapidly Accreting AGN RE J1034+396

We report the analysis of ∼1 Ms of XMM-Newton observations of the rapidly accreting active galactic nucleus RE J1034+396. The 0.3–9 keV EPIC-pn spectra are well described by a model consisting of steep continuum emission from the corona accompanied by relativistically blurred reflection from a highly ionized accretion disk. The source is known to exhibit strong excess soft X-ray emission, which we show is well represented by thermal disk photons Comptonized by a warm plasma spanning the inner accretion flow. Additionally, the EPIC-pn data provide compelling evidence (ΔC ∼ 60 for four additional parameters) for the presence of an ultrafast outflow (UFO) with a line-of-sight velocity v/c = 0.307$^{+0.001}_{-0.005}$, and an emission signature consistent with reflection of the corona from modestly ionized, outflowing gas. The simultaneous 0.5–2.5 keV RGS spectra show clear absorption lines. Modelling of these data confirms the presence of the UFO and constrains its equivalent hydrogen column density, log N H /(atom cm −2 ) = 21.7$^{+0.1}_{-0.2}$. The RGS data also reveal at least two warm absorber components with a modest outflow velocity (1680$^{+40}_{-50}$ km s -1 ). The measured properties and time evolution of the UFO in RE J1034+396 suggest that it is formed from collisionally ionized plasma, launched from the disk surface and accelerated by radiation pressure. The high terminal velocity and substantial absorbing column density imply that the outflow carries sufficient momentum and energy to transform its environment, being capable of driving out essentially all dust and gas it interacts with along the line of sight, even if the AGN were initially surrounded by a Compton-thick absorber.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The cosmic-ray positron excess and its imprint in the Galactic gamma-ray sky

Abstract We study the origin of the positron excess observed in the local cosmic-ray spectrum at high energies, and relate it to the cosmic rays and gamma-ray emission across the entire Galaxy. In particular, we explore the hypothesis of a single, dominant source accountable for primary electron-positron pairs. Since we are agnostic about the physical nature of the underlying source population, we consider four models that are representative of young pulsars, old stars (as a tracer of millisecond pulsars), and annihilating dark matter particles. In the dark matter hypothesis, we consider both a cored and a cuspy model for the halo in the Milky Way. Then, we compare the associated gamma-ray sky maps withFermi-LAT data. We find that the emission arising from pulsar wind nebulae is fairly concentrated near the mid plane, and therefore additional cosmic-ray sources must be invoked to explain the emission at the center of the Galaxy. If the local positron excess were mainly due to millisecond pulsars, inverse Compton scattering by the particles injected would naturally account for a non-negligible fraction of the central gamma-ray emission. The same process would lead to a tension for a standard NFW dark matter profile, exceeding the gamma-ray data by almost a factor of 2 in some regions of the Galaxy. Although the results for an isothermal, cored profile are in agreement with these data, the cross section needed in both cases is around 2 orders of magnitude above the thermal cross-section, disfavouring the dark matter interpretation.

Astronomy & Astrophysics↗