Journal of Cosmology and Astroparticle Physics

The capture of dark matter, and its subsequent annihilation, can heat old, isolated neutron stars. In order for kinetic heating to be achieved, the captured dark matter must undergo sufficient scattering to deposit its kinetic energy in the star. We find that this energy deposit typically occurs quickly, for most of the relevant parameter space. In order for appreciable annihilation heating to also be achieved, the dark matter must reach a state of capture-annihilation equilibrium in the star. We show that this can be fulfilled for all types of dark matter-baryon interactions. This includes cases where the scattering or annihilation cross sections are momentum or velocity suppressed in the non-relativistic limit. Importantly, we find that capture-annihilation equilibrium, and hence maximal annihilation heating, can be achieved without complete thermalization of the captured dark matter. For scattering cross sections that saturate the capture rate, we find that capture-annihilation equilibrium is typically reached on a timescale of less than 1 year for vector interactions and 10⁴ years for scalar interactions.

We provide ingredients and recipes for computing signals of TeV-scale Dark Matter annihilations and decays in the Galaxy and beyond. For each DM channel, we present the energy spectra of at production, computed by high-statistics simulations. We estimate the Monte Carlo uncertainty by comparing the results yielded by the Pythia and Herwig event generators. We then provide the propagation functions for charged particles in the Galaxy, for several DM distribution profiles and sets of propagation parameters. Propagation of e^± is performed with an improved semi-analytic method that takes into account position-dependent energy losses in the Milky Way. Using such propagation functions, we compute the energy spectra of e^±, and at the location of the Earth. We then present the gamma ray fluxes, both from prompt emission and from Inverse Compton scattering in the galactic halo. Finally, we provide the spectra of extragalactic gamma rays. All results areavailable in numerical form and ready to be consumed.

We show that by allowing our Universe to merge with other universes one is lead to modified Friedmann equations that explain the present accelerated expansion of our Universe without the need of a cosmological constant.

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands centered at: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes and one large-aperture 6-m telescope, with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The small aperture telescopes will target the largest angular scales observable from Chile, mapping ≈ 10% of the sky to a white noise level of 2 μK-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, r, at a target level of σ(r)=0.003. The large aperture telescope will map ≈ 40% of the sky at arcminute angular resolution to an expected white noise level of 6 μK-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the Large Synoptic Survey Telescope sky region and partially with the Dark Energy Spectroscopic Instrument. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.

One of the most exciting and pressing issues in cosmology today is the discrepancy between some measurements of the local Hubble constant and other values of the expansion rate inferred from the observed temperature and polarization fluctuations in the cosmic microwave background (CMB) radiation. Resolving these differences holds the potential for the discovery of new physics beyond the standard model of cosmology: Lambda Cold Dark Matter (ΛCDM), a successful model that has been in place for more than 20 years. Given both the fundamental significance of this outstanding discrepancy, and the many-decades-long effort to increase the accuracy of the extragalactic distance scale, it is critical to demonstrate that the local measurements are convincingly free from residual systematic errors. We review the progress over the past quarter century in measurements of the local value of the Hubble constant, and discuss remaining challenges. Particularly exciting are new data from the James Webb Space Telescope (JWST), for which we present an overview of our program and first results. We focus in particular on Cepheids and the Tip of the Red Giant Branch (TRGB) stars, as well as a relatively new method, the JAGB (J-Region Asymptotic Giant Branch) method, all methods that currently exhibit the demonstrably smallest statistical and systematic uncertainties. JWST is delivering high-resolution near-infrared imaging data to both test for and to address directly several of the systematic uncertainties that have historically limited the accuracy of extragalactic distance scale measurements (e.g., the dimming effects of interstellar dust, chemical composition differences in the atmospheres of stars, and the crowding and blending of Cepheids contaminated by nearby previously unresolved stars). For the first galaxy in our program, NGC 7250, the high-resolution JWST images demonstrate that many of the Cepheids observed with the Hubble Space Telescope (HST) are significantly crowded by nearby neighbors. Avoiding the more significantly crowded variables, the scatter in the JWST near-infrared (NIR) Cepheid PL relation is decreased by a factor of two compared to those from HST, illustrating the power of JWST for improvements to local measurements of H₀. Ultimately, these data will either confirm the standard model, or provide robust evidence for the inclusion of additional new physics.

The Einstein Telescope (ET), the European project for a third-generation gravitational-wave detector, has a reference configuration based on a triangular shape consisting of three nested detectors with 10 km arms, where each detector has a 'xylophone' configuration made of an interferometer tuned toward high frequencies, and an interferometer tuned toward low frequencies and working at cryogenic temperature. Here, we examine the scientific perspectives under possible variations of this reference design. We perform a detailed evaluation of the science case for a single triangular geometry observatory, and we compare it with the results obtained for a network of two L-shaped detectors (either parallel or misaligned) located in Europe, considering different choices of arm-length for both the triangle and the 2L geometries. We also study how the science output changes in the absence of the low-frequency instrument, both for the triangle and the 2L configurations. We examine a broad class of simple 'metrics' that quantify the science output, related to compact binary coalescences, multi-messenger astronomy and stochastic backgrounds, and we then examine the impact of different detector designs on a more specific set of scientific objectives.

We present new arcminute-resolution maps of the Cosmic Microwave Background temperature and polarization anisotropy from the Atacama Cosmology Telescope, using data taken from 2013–2016 at 98 and 150 GHz. The maps cover more than 17,000 deg², the deepest 600 deg² with noise levels below 10μK-arcmin. We use the power spectrum derived from almost 6,000 deg² of these maps to constrain cosmology. The ACT data enable a measurement of the angular scale of features in both the divergence-like polarization and the temperature anisotropy, tracing both the velocity and density at last-scattering. From these one can derive the distance to the last-scattering surface and thus infer the local expansion rate, H₀. By combining ACT data with large-scale information from WMAP we measure H₀=67.6± 1.1 km/s/Mpc, at 68% confidence, in excellent agreement with the independently-measured Planck satellite estimate (from ACT alone we find H₀=67.9± 1.5 km/s/Mpc). The ΛCDM model provides a good fit to the ACT data, and we find no evidence for deviations: both the spatial curvature, and the departure from the standard lensing signal in the spectrum, are zero to within 1σ; the number of relativistic species, the primordial Helium fraction, and the running of the spectral index are consistent with ΛCDM predictions to within 1.5–2.2σ. We compare ACT, WMAP, and Planck at the parameter level and find good consistency; we investigate how the constraints on the correlated spectral index and baryon density parameters readjust when adding CMB large-scale information that ACT does not measure. The DR4 products presented here will be publicly released on the NASA Legacy Archive for Microwave Background Data Analysis.

We formulate f(Q,C) gravity and cosmology. Such a construction is based on the symmetric teleparallel geometry, but apart form the non-metricity scalar Q we incorporate in the Lagrangian the boundary term C of its difference from the standard Levi-Civita Ricci scalar  R̊. We extract the general metric and affine connection field equations, we apply them at a cosmological framework, and adopting three different types of symmetric teleparallel affine connections we obtain the modified Friedmann equations. As we show, we acquire an effective dark-energy sector of geometrical origin, which can lead to interesting cosmological phenomenology. Additionally, we may obtain an effective interaction between matter and dark energy. Finally, examining a specific model, we show that we can obtain the usual thermal history of the universe, with the sequence of matter and dark-energy epochs, while the effective dark-energy equation-of-state parameter can be quintessence-like, phantom-like, or cross the phantom-divide during evolution.

The combined fit of the measured energy spectrum and shower maximum depth distributions of ultra-high-energy cosmic rays is known to constrain the parameters of astrophysical models with homogeneous source distributions. Studies of the distribution of the cosmic-ray arrival directions show a better agreement with models in which a fraction of the flux is non-isotropic and associated with the nearby radio galaxy Centaurus A or with catalogs such as that of starburst galaxies. Here, we present a novel combination of both analyses by a simultaneous fit of arrival directions, energy spectrum, and composition data measured at the Pierre Auger Observatory. The model takes into account a rigidity-dependent magnetic field blurring and an energy-dependent evolution of the catalog contribution shaped by interactions during propagation. We find that a model containing a flux contribution from the starburst galaxy catalog of around 20% at 40 EeV with a magnetic field blurring of around 20° for a rigidity of 10 EV provides a fair simultaneous description of all three observables. The starburst galaxy model is favored with a significance of 4.5σ (considering experimental systematic effects) compared to a reference model with only homogeneously distributed background sources. By investigating a scenario with Centaurus A as a single source in combination with the homogeneous background, we confirm that this region of the sky provides the dominant contribution to the observed anisotropy signal. Models containing a catalog of jetted active galactic nuclei whose flux scales with the γ-ray emission are, however, disfavored as they cannot adequately describe the measured arrival directions.

We study the full-sky distribution of the radio emission from the stimulated decay of axions which are assumed to compose the dark matter in the Galaxy. Besides the constant extragalactic and CMB components, the decays are stimulated by a Galactic radio emission with a spatial distribution that we empirically determine from observations. We compare the diffuse emission to the counterimages of the brightest supernovæ remnants, and take into account the effects of free-free absorption. We show that, if the dark matter halo is described by a cuspy NFW profile, the expected signal from the Galactic center is the strongest. Interestingly, the emission from the Galactic anti-center provides competitive constraints that do not depend on assumptions on the uncertain dark matter density in the inner region. Furthermore, the anti-center of the Galaxy is the brightest spot if the Galactic dark matter density follows a cored profile. The expected signal from stimulated decays of axions of mass m_a ∼ 10^-6 eV is within reach of the Square Kilometer Array for an axion-photon coupling  g_aγ ≳ (2-3) × 10^-11 GeV^-1.

Ricci scalar being zero is equivalent to the vacuum field equation in Finsler space-time. The Schwarzschild metric can be concluded from the field equation's solution if the space-time conserves spherical symmetry. This research aims to investigate Finslerian Schwarzschild-de Sitter space-time. Recent studies based on Finslerian space-time geometric models are becoming more prevalent because the local anisotropic structure of space-time influences the gravitational field and gives rise to modified cosmological relations. We suggest a gravitational field equation with a non-zero cosmological constant in Finslerian geometry and apprehend that the presented Finslerian gravitational field equation corresponds to the non-zero Ricci scalar. In Finsler geometry, the peer of spherical symmetry is the Finslerian sphere. Assuming space-time to conserve the "Finslerian sphere" symmetry, the counterpart of the Riemannian sphere (Finslerian sphere) must have a constant flag curvature (λ). It is demonstrated that the Finslerian covariant derivative of the geometric part of the gravitational field equation is preserved under a condition using the Chern connection. According to the string theory, string clouds can be defined as a pool of strings made due to symmetry breaking in the universe's early stages. We find that for λ ≠ 1, this solution resembles a black hole surrounded by a cloud of strings. Furthermore, we investigate null and time-like geodesics for λ = 1. In this regard, the photon geodesics are obtained that are the closest paths to the photon sphere of the first photons visible at the black hole shadow limit. Also, circular orbit conditions are obtained for the effective potential.

The thermal plasma filling the early universe generated a stochastic gravitational wave background that peaks in the microwave frequency range today. If the graviton production rate is expressed as a series in a fine-structure constant, α, and the temperature over the Planck mass,  T²/m_pl², then the lowest-order contributions come from single (∼αT²/m_pl²) and double (∼T⁴/m_pl⁴) graviton production via 2 → 2 scatterings. We show that in the Standard Model, single-graviton production dominates if the maximal temperature is smaller than 4 × 10¹⁸ GeV. This justifies previous calculations which relied solely on single-graviton production. We mention Beyond the Standard Model scenarios in which the single and double-graviton contributions could be of comparable magnitudes. Finally, we elaborate on what these results imply for the range of applicability of General Relativity as an effective theory.

Observations of the Lyman-α forest in distant quasar spectra with upcoming surveys are expected to provide significantly larger and higher-quality datasets. To interpret these datasets, it is imperative to develop efficient simulations. One such approach is based on the assumption that baryonic densities in the intergalactic medium (IGM) follow a lognormal distribution. We extend our earlier work to assess the robustness of the lognormal model of the Lyman-α forest in recovering the parameters characterizing IGM state, namely, the mean-density IGM temperature (T₀), the slope of the temperature-density relation (γ), and the hydrogen photoionization rate (Γ₁₂), by comparing with high-resolution Sherwood SPH simulations across the redshift range 2 ≤ z ≤ 2.7. These parameters are estimated through a Markov Chain Monte Carlo (MCMC) technique, using the mean and power spectrum of the transmitted flux. We find that the usual lognormal distribution of IGM densities cannot recover the parameters of the SPH simulations. This limitation arises from the fact that the SPH baryonic density distribution cannot be described by a simple lognormal form. To address this, we extend the model by scaling the linear density contrast by a parameter ν. While the resulting baryonic density is still lognormal, the additional parameter gives us extra freedom in setting the variance of density fluctuations. With this extension, values of T₀ and γ implied in the SPH simulations are recovered at ∼ 1 - σ (≲ 10%) of the median (best-fit) values for most redshifts bins. However, this extended lognormal model cannot recover Γ₁₂ reliably, with the best-fit value discrepant by ≳ 3 - σ for z > 2.2. Despite this limitation in the recovery of Γ₁₂, whose origins we explain, we argue that the model remains useful for constraining cosmological parameters.

We study conformal transformations in the most general parity-preserving models of the New General Relativity type. Then we apply them to analysis of cosmological perturbations in the (simplest) spatially flat cosmologies. Strong coupling issues around Minkowski spacetime are seen for many special cases of these models. At the same time, the behaviour of the most general three-parameter case seems to be very robust, presumably always with only the eight first-class constraints coming from diffeomorphisms. Also the case of the so-called 1-parameter New GR doesn't show any discontinuity between Minkowski and the cosmology, though without showing any deviations from GR which would be observable at this level either.

Non-radial oscillations of Neutron Stars (NSs) provide a means to learn important details regarding their interior composition and equation of state. We consider the effects of Δ-baryons on non-radial f-mode oscillations and other NS properties within the Density-Dependent Relativistic Mean Field formalism. Calculations are performed for Δ-admixed NS matter with and without hyperons. Our study of the non-radial f-mode oscillations revealed a distinct increase in frequency due to the addition of the Δ-baryons with upto 20% increase in frequency being seen for canonical NSs. Other bulk properties of NSs, including mass, radii, and dimensionless tidal deformability (Λ) were also affected by these additional baryons. Comparing our results with available observational data from pulsars (NICER) and gravitational waves (LIGO-VIRGO collaboration), we found strong agreement, particularly concerning Λ.

The thermal plasma filling the early universe generated a stochastic gravitational wave background that peaks in the microwave frequency range today. If the graviton production rate is expressed as a series in a fine-structure constant, α, and the temperature over the Planck mass,  T²/m_pl², then the lowest-order contributions come from single (∼αT²/m_pl²) and double (∼T⁴/m_pl⁴) graviton production via 2 → 2 scatterings. We show that in the Standard Model, single-graviton production dominates if the maximal temperature is smaller than 4 × 10¹⁸ GeV. This justifies previous calculations which relied solely on single-graviton production. We mention Beyond the Standard Model scenarios in which the single and double-graviton contributions could be of comparable magnitudes. Finally, we elaborate on what these results imply for the range of applicability of General Relativity as an effective theory.

It is a fundamental unsolved question in general relativity how to unambiguously characterize the effective collective dynamics of an ensemble of fluid elements sourcing the local geometry, in the absence of exact symmetries. In a cosmological context this is sometimes referred to as the averaging problem. At the heart of this problem in relativity is the non-uniqueness of the choice of foliation within which the statistical properties of the local spacetime are quantified, which can lead to ambiguity in the formulated average theory. This has led to debate in the literature on how to best construct and view such a coarse-grained hydrodynamic theory. Here, we address this ambiguity by performing the first quantitative investigation of foliation dependence in cosmological spatial averaging. Starting from the aim of constructing slicing-independent integral functionals (volume, mass, entropy, etc.) as well as average functionals (mean density, average curvature, etc.) defined on spatial volume sections, we investigate infinitesimal foliation variations and derive results on the foliation dependence of functionals and on extremal leaves. Our results show that one may only identify fully foliation-independent integral functionals in special scenarios, requiring the existence of associated conserved currents. We then derive bounds on the foliation dependence of integral functionals for general scalar quantities under finite variations within physically motivated classes of foliations. Our findings provide tools that are useful for quantifying, eliminating or constraining the foliation dependence in cosmological averaging.

It is assumed in standard cosmology that the Universe underwent a period of inflation in its earliest phase, providing the seeds for structure formation through vacuum fluctuations of the inflaton scalar field. These fluctuations get stretched by the quasi-exponential expansion of the Universe and become squeezed. The aim of this paper is to deepen the understanding of the squeezing process, considering the effect of self-interactions. Axion-like particles can provide a useful setup to study this effect. Specifically we focus on the consequences that a non-trivial evolution of the background axion field has on the squeezing of the perturbations. We follow the evolution of the axion's fluctuation modes from the horizon exit during inflation to the radiation-dominated epoch. We compute Bogoliubov coefficients and squeezing parameters, which are linked to the axion particle number and isocurvature perturbation. We find that the quantum mechanical particle production and the squeezing of the perturbations are enhanced, if one accounts for anharmonic effects, i.e., the effect of higher order terms in the potential. This effect becomes particularly strong towards the hilltop of the potential.

The excess radio background detected by ARCADE 2 represents a puzzle within the standard cosmological model. There is no clear viable astrophysical solution, and therefore, it might indicate the presence of new physics. Radiative decays of a relic neutrino ν_i (either i=1, or i=2, or i=3) into a sterile neutrino ν_ s, assumed to be quasi-degenerate, provide a solution that currently evades all constraints posed by different cosmological observations and reproduces very well the ARCADE 2 data. We find a very good fit to the ARCADE 2 data with best fit values τ_i = 1.46 × 10^21 s and Δ m_i = 4.0 × 10^-5 eV, where τ_i is the lifetime and Δ m_i is the mass difference between the decaying active neutrino and the sterile neutrino. On the other hand, if relic neutrino decays do not explain ARCADE 2 data, then these place a stringent constraint Δ m_i^3/2τ_i ≳ 2 × 10^14 eV^3/2 s in the range 1.4 × 10^-5 eV<Δ m_i < 2.5 × 10^-4 eV. The solution also predicts a stronger 21 cm absorption global signal than the predicted one from the ΛCDM model, with a contrast brightness temperature T_21 = -238^+21_-20 mK (99% C.L.) at redshift z≃ 17. This is in mild tension with the even stronger signal found by the EDGES collaboration, T_21 = - 500^+200_-500 mK, suggesting that this might have been overestimated, possibly receiving a contribution from some unidentified foreground source.

A singlet majoron can arise from the seesaw framework as a pseudo-Goldstone boson when the heavy Majorana neutrinos acquire masses via the spontaneous breaking of global U(1)_L symmetry. The resulting cosmological impacts are usually derived from the effective majoron-neutrino interaction, and the majoron abundance is accumulated through the freeze-in neutrino coalescence. However, a primordial majoron abundance can be predicted in a minimal setup and lead to distinctive cosmological effects. In this work, we consider such a primordial majoron abundance from relativistic freeze-out and calculate the modification to the effective neutrino number N_eff. We demonstrate that the measurements of N_eff will constrain the parameter space from a primordial majoron abundance in an opposite direction to that from neutrino coalescence. When the contributions from both the primordial abundance and the freeze-in production coexist, the U(1)_L-breaking scale (seesaw scale) f will be pushed into a "sandwiched window". Remarkably, for majoron masses below 1 MeV and above the eV scale, the future CMB-S4 experiment will completely close such a low-scale seesaw window for f ∈ [1,10⁵] GeV. We highlight that any new light particle with a primordial abundance that couples to Standard Model particles may lead to a similar sandwiched window, and such a general phenomenon deserves careful investigation.

Single-field models of inflation might lead to amplified scalar fluctuations on small scales due, for example, to a transient ultra-slow-roll phase. It was argued by Kristiano & Yokoyama in ref. [1] that the enhanced amplitude of the scalar power spectrum on small scales has the potential to induce a sizeable 1-loop correction to the spectrum at large scales. In this work, we repeat the calculation for the 1-loop correction presented in ref. [1]. We closely follow their assumptions but evaluate the loop numerically. This allows us to consider both instantaneous and smooth transitions between the slow-roll and ultra-slow-roll phases. In particular, we generate models featuring realistic, smooth evolution from an analytic inflationary potential. We find that, upon fixing the amplitude of the peak in the power spectrum at short scales, the resulting 1-loop correction is not significantly reduced by considering a smooth evolution. In particular, for a power spectrum with a tree-level peak amplitude potentially relevant for small-scale phenomenology, e.g. primordial black hole production, the 1-loop correction on large scales is a few percent of the tree-level power spectrum.

We confront massive Proca-Nuevo gravity with cosmological observations. The former is a non-linear theory involving a massive spin-1 field, that can be extended incorporating operators of the Generalized Proca class, and when coupled to gravity it can be covariantized in a way that exhibits consistent and ghost-free cosmological solutions, without experiencing instabilities and superluminalities at the perturbative level. When applied at a cosmological framework it induces extra terms in the Friedmann equations, however due to the special non-linear construction the field is eliminated in favor of the Hubble function. Thus, the resulting effective dark energy sector is dynamical, however it contains the same number of free parameters with the ΛCDM concordance model. We use data from Supernovae Ia (SNIa) and Cosmic Chronometers (CC) observations and we construct the corresponding likelihood-contours for the free parameters. Interestingly enough, application of various information criteria, such as AIC, BIC and DIC, shows that the scenario of massive Proca-Nuevo gravity, although having exactly the same number of free parameters with ΛCDM paradigm, it is more efficient in fitting the data. Finally, the reconstructed dark-energy equation-of-state parameter shows statistical compatibility with the model-independent, data-driven reconstructed one.

We consider black hole formation due to the gravitational collapse produced by large density fluctuations during an epoch of reheating with a stiff equation of state and calculate the induced gravitational wave spectrum. By considering the existing bounds on the total energy density of gravitational waves today, we find constraints on the parameter space of this scenario. We also calculate the lepton asymmetry generated by metric perturbations via the chiral gravitational anomaly present in the Standard Model and find that, once the electroweak sphaleron processes have taken place, the large spectrum of scalar perturbations responsible for black hole formation induces a peak in the baryon asymmetry fluctuations on small scales.

Gravity theories that modify General Relativity in the slow-motion regime can introduce nonperturbative corrections to the stochastic gravitational-wave background (SGWB) from supermassive black-hole binaries in the nano-Hertz band, while not affecting the quadrupolar nature of the gravitational-wave radiation and remaining perturbative in the highly-relativistic regime, as to satisfy current post-Newtonian (PN) constraints. We present a model-agnostic formalism to map such theories into a modified tilt for the SGWB spectrum, showing that negative PN corrections (in particular -2PN) can alleviate the tension in the recent pulsar-timing-array data if the detected SGWB is interpreted as arising from supermassive binaries. Despite being preliminary, current data have already strong constraining power, for example they set a novel (conservative) upper bound on theories with time-varying Newton's constant (a -4PN correction) at least at the level of Ġ/G ≲ 10^-5 yr^-1 for redshift z=[0.1÷1]. We also show that NANOGrav data are best fitted by a broken power-law interpolating between a dominant -2PN or -3PN modification at low frequency, and the standard general-relativity scaling at high frequency. Nonetheless, a modified gravity explanation should be confronted with binary eccentricity, environmental effects, nonastrophysical origins of the signal, and scrutinized against statistical uncertainties. These novel tests of gravity will soon become more stringent when combining all pulsar-timing-array facilities and when collecting more data.

We study the full-sky distribution of the radio emission from the stimulated decay of axions which are assumed to compose the dark matter in the Galaxy. Besides the constant extragalactic and CMB components, the decays are stimulated by a Galactic radio emission with a spatial distribution that we empirically determine from observations. We compare the diffuse emission to the counterimages of the brightest supernovæ remnants, and take into account the effects of free-free absorption. We show that, if the dark matter halo is described by a cuspy NFW profile, the expected signal from the Galactic center is the strongest. Interestingly, the emission from the Galactic anti-center provides competitive constraints that do not depend on assumptions on the uncertain dark matter density in the inner region. Furthermore, the anti-center of the Galaxy is the brightest spot if the Galactic dark matter density follows a cored profile. The expected signal from stimulated decays of axions of mass m_a ∼ 10^-6 eV is within reach of the Square Kilometer Array for an axion-photon coupling  g_aγ ≳ (2-3) × 10^-11 GeV^-1.