Dwarf Spheroidal Galaxies Research Articles

Semiannihilation of fermionic dark matter

The continued nonobservation of events emanating from dark matter annihilations in various direct and indirect detection experiments calls into question the mechanism for determining the relic density of a weakly interacting massive particle. However, if the relic density is determined primarily by a semiannihilation process, as opposed to the usual annihilation, this tension can be ameliorated. Here, we investigate a ${Z}_{3}$ symmetric effective field theory incorporating a fermionic dark matter that semiannihilates to right-handed neutrinos. The dynamics of the right-handed neutrinos and the impact of its late decays are also scrutinized while obtaining the correct dark matter relic. Finally, indirect detection bounds on the semiannihilation cross sections are drawn from the gamma-ray observations in the direction of dwarf spheroidal galaxies (Fermi-LAT), and including the projections obtained for the H.E.S.S. and the CTA detectors.

VTSCat: The VERITAS Catalog of Gamma-Ray Observations

The ground-based gamma-ray observatory Very Energetic Radiation Imaging Telescope Array System (VERITAS, https://veritas.sao.arizona.edu/) is sensitive to photons of astrophysical origin with energies in the range between ≈85 GeV and ≈30 TeV. The instrument consists of four 12 m diameter imaging Cherenkov telescopes operating at the Fred Lawrence Whipple Observatory in southern Arizona. VERITAS started four-telescope operations in 2007 and collects about 1100 hr of good-weather data per year. The VERITAS collaboration has published over 100 journal articles since 2008 reporting on gamma-ray observations of a large variety of objects: Galactic sources like supernova remnants, pulsar wind nebulae, and binary systems; extragalactic sources like star-forming galaxies, dwarf-spheroidal galaxies, and highly variable active galactic nuclei. This note presents VTSCat: the catalog of high-level data products from all VERITAS publications.

Dynamics of dwarf galaxies in f(R) gravity

ABSTRACT We use the kinematic data of the stars in eight dwarf spheroidal galaxies to assess whether f(R) gravity can fit the observed profiles of the line-of-sight velocity dispersion of these systems without resorting to dark matter. Our model assumes that each galaxy is spherically symmetric and has a constant velocity anisotropy parameter β and constant mass-to-light ratio consistent with stellar population synthesis models. We solve the spherical Jeans equation that includes the Yukawa-like gravitational potential appearing in the weak field limit of f(R) gravity, and a Plummer density profile for the stellar distribution. The f(R) velocity dispersion profiles depends on two parameters: the scale length ξ−1, below which the Yukawa term is negligible, and the boost of the gravitational field δ > −1. δ and ξ are not universal parameters, but their variation within the same class of objects is expected to be limited. The f(R) velocity dispersion profiles fit the data with a value $\xi ^{-1}= 1.2^{+18.6}_{-0.9}$ Mpc for the entire galaxy sample. On the contrary, the values of δ show a bimodal distribution that picks at $\overline{\delta }=-0.986\pm 0.002$ and $\overline{\delta }=-0.92\pm 0.01$. These two values disagree at 6σ and suggest a severe tension for f(R) gravity. It remains to be seen whether an improved model of the dwarf galaxies or additional constraints provided by the proper motions of stars measured by future astrometric space missions can return consistent δ’s for the entire sample and remove this tension.

Indirect signals of dark matter can change depending on where you look

We propose that the nature of indirect signals of dark matter (DM) can depend on the Galactic environment they originate from. We demonstrate this possibility in models where DM annihilates into light mediators whose branching fractions depend on a long-range force sourced by ordinary matter. In particular, electromagnetic signals of DM may only arise near the centers of galaxies where the ordinary matter densities, and hence astrophysical background levels, are high. We briefly discuss how our model could explain the Galactic Center gamma ray excess, without leaving much of a trace in baryon-poor environments, like dwarf spheroidal galaxies. Similar spatial dependence of indirect signals can also apply to models featuring metastable DM decay into light mediators.

The metallicity distribution in the core of the Sagittarius dwarf spheroidal: Minimising the metallicity biases

We present the metallicity and radial velocity for 450 bona fide members of the Sagittarius dwarf spheroidal (Sgr dSph) galaxy, measured from high-resolution spectra (R ≃ 18 000) obtained with FLAMES at the VLT. The targets were carefully selected (a) to sample the core of the main body of Sgr dSph while avoiding contamination from the central stellar nucleus, and (b) to prevent any bias on the metallicity distribution by selecting targets based on theirGaiaparallax and proper motions. All the targets selected in this way were confirmed as radial velocity members. We used this sample to derive the first metallicity distribution of the core of Sgr dSph, which is virtually unaffected by metallicity biases. The observed distribution ranges from [Fe/H] ≃ −2.3 to [Fe/H] ≃ 0.0, with a strong, symmetric, and relatively narrow peak around [Fe/H] ≃ −0.5 and a weak and extended metal-poor tail, in which only 13.8 ± 1.9% of the stars have [Fe/H] < −1.0. We confirm previous evidence of correlations between chemical and kinematical properties of stars in the core of Sgr. In our sample, stars with [Fe/H] ≥ −0.6 display a lower velocity dispersion and a higher rotation amplitude than those with [Fe/H] < −0.6, confirming previous suggestions of disk/halo structure for the progenitor of the system.

Investigating the effect of Milky Way dwarf spheroidal galaxies extension on dark matter searches with Fermi-LAT data

Satellite galaxies of the Milky Way with high mass-to-light ratios and little baryon content, i.e., dwarf spheroidal galaxies (dSphs), are among the most promising targets to detect or constrain the nature of dark matter (DM) through its final annihilation products into high-energy photons. Previously, the assumption that DM emission from dSphs is pointlike has been used to set strong constraints on DM candidates using data from the Fermi Large Area Telescope (LAT). However, due to their high DM densities and proximity, dSphs actually have sufficient angular extension to be detected by the Fermi-LAT. Here, we perform a comprehensive analysis about the impact of accounting for angular extension in the search for gamma-ray DM signals toward known dSphs with Fermi-LAT. We show that, depending on the dSph under consideration, limits on the DM cross section can be weakened by up to a factor of 2--2.5, while the impact on the stacked, i.e., combined, limits is at most 1.5--1.8 depending on the annihilation channel. This result is of relevance when comparing dSph limits to other multimessenger DM constraints and for testing the DM interpretation of anomalous ``excesses.''

A probabilistic deep learning model to distinguish cusps and cores in dwarf galaxies

ABSTRACT Numerical simulations within a cold dark matter (DM) cosmology form haloes whose density profiles have a steep inner slope (‘cusp’), yet observations of galaxies often point towards a flat central ‘core’. We develop a convolutional mixture density neural network model to derive a probability density function (PDF) of the inner density slopes of DM haloes. We train the network on simulated dwarf galaxies from the NIHAO and AURIGA projects, which include both DM cusps and cores: line-of-sight velocities and 2D spatial distributions of their stars are used as inputs to obtain a PDF representing the probability of predicting a specific inner slope. The model recovers accurately the expected DM profiles: $\sim 82{{\ \rm per\ cent}}$ of the galaxies have a derived inner slope within ±0.1 of their true value, while $\sim 98{{\ \rm per\ cent}}$ within ±0.3. We apply our model to four Local Group dwarf spheroidal galaxies and find results consistent with those obtained with the Jeans modelling based code GravSphere: the Fornax dSph has a strong indication of possessing a central DM core, Carina and Sextans have cusps (although the latter with large uncertainties), while Sculptor shows a double peaked PDF indicating that a cusp is preferred, but a core cannot be ruled out. Our results show that simulation-based inference with neural networks provide a innovative and complementary method for the determination of the inner matter density profiles in galaxies, which in turn can help constrain the properties of the elusive DM.

Search for primordial black hole dark matter with x-ray spectroscopic and imaging satellite experiments and prospects for future satellite missions

Ultra-light primordial black holes (PBHs) in the mass range of 10$^{16}$ - 10$^{22}$ g are allowed by current observations to constitute a significant fraction, if not all, of the dark matter in the Universe. In this work, we present limits on ultra-light, non-rotating PBHs which arise from the non-detection of the Hawking radiation signals from such objects in the keV-MeV energy band. Namely, we consider observations from the current-generation missions XMM-Newton and INTEGRAL/SPI and discuss the observational perspectives of the future missions Athena, eXTP, and THESEUS for PBH searches. Based on 3.4 Msec total exposure time XMM-Newton observations of Draco dwarf spheroidal galaxy, we conclude that PBH with masses $\lesssim 10^{16}$ g can not make all dark matter at 95% confidence level. Our ON-OFF-type analysis of $>100$ Msec of INTEGRAL/SPI data on the Milky Way halo puts significantly stronger constraints. Only $\lesssim 10$% dark matter can be presented by PBHs with masses $\lesssim 3\cdot 10^{16}$ g while the majority of dark matter can not be represented by PBHs lighter than $7\cdot 10^{16}$ g at 95% confidence level. We discuss the strong impact of systematic uncertainty related to the variations of instrumental and astrophysical INTEGRAL/SPI background on the derived results and estimate its level. We also show that future large-field-of-view missions such as THESEUS/X-GIS will be able to improve the constraints by a factor of 10-100 depending on the level of control under the systematics of these instruments.

Proper Motions, Orbits, and Tidal Influences of Milky Way Dwarf Spheroidal Galaxies

We combine Gaia early data release 3 astrometry with accurate photometry and utilize a probabilistic mixture model to measure the systemic proper motion of 52 dwarf spheroidal (dSph) satellite galaxies of the Milky Way (MW). For the 46 dSphs with literature line-of-sight velocities we compute orbits in both a MW and a combined MW + Large Magellanic Cloud (LMC) potential and identify Car II, Car III, Hor I, Hyi I, Phx II, and Ret II as likely LMC satellites. 40% of our dSph sample has a >25% change in pericenter and/or apocenter with the MW + LMC potential. For these orbits, we use a Monte Carlo sample for the observational uncertainties for each dSph and the uncertainties in the MW and LMC potentials. We predict that Ant II, Boo III, Cra II, Gru II, and Tuc III should be tidally disrupting by comparing each dSph's average density relative to the MW density at its pericenter. dSphs with large ellipticity (CVn I, Her, Tuc V, UMa I, UMa II, UMi, Wil 1) show a preference for their orbital direction to align with their major axis even for dSphs with large pericenters. We compare the dSph radial orbital phase to subhalos in MW-like N-body simulations and infer that there is not an excess of satellites near their pericenter. With projections of future Gaia data releases, we find that dSph's orbital precision will be limited by uncertainties in the distance and/or MW potential rather than in proper motion precision. Finally, we provide our membership catalogs to enable community follow-up.

Indirect detection of eV dark matter via infrared spectroscopy

Infrared spectroscopy has been developed significantly. In particular, infrared photons can be measured with high spectral and angular resolution in state-of-art spectrographs. They are sensitive to monochromatic photons due to the decay and annihilation of particles beyond the Standard Model, such as dark matter (DM), while suppressing background photons that form a continuous spectrum. In this paper, we study the indirect detection of the DM decaying into infrared light using infrared spectrographs. In particular, we show that serious thermal and astrophysical noises can be overcome. As concrete examples, the Warm INfrared Echelle spectrograph to Realize Extreme Dispersion and sensitivity (WINERED) installed at the Magellan Clay 6.5m telescope and Near-Infrared Spectrograph (NIRSpec) at the James Webb Space Telescope (JWST) are discussed. We show that a few hours of measurements of a faint dwarf spheroidal galaxy with WINERED (NIRSpec-like spectrograph) in the Magellan telescope (JWST) can probe an axion-like particle DM in the mass range $m_\phi=1.8 - 2.7\,$eV ($0.5-4\,$eV) with a photon coupling $g_{\phi\gamma\gamma}\gtrsim 10^{-11}{\rm GeV}^{-1}$. Complemental approaches, taking advantage of the high resolutions, such as the measurement of the Doppler shift of the signal photon lines and the possible search of the DM decay around the Milky Way galaxy center with Infrared Camera and Spectrograph (IRCS) at 8.2m Subaru telescope, are also presented.

Constraining the primordial curvature perturbation using dark matter substructure

We investigate the primordial curvature perturbation by the observation of dark matter substructure. Assuming a bump in the spectrum of the curvature perturbation in the wavenumber of k>1 Mpc^{-1}, we track the evolution of the host halo and subhalos in a semi-analytic way. Taking into account possible uncertainties in the evaluation of the tidal stripping effect on the subhalo growth, we find a new robust bound on the curvature perturbation with a bump from the number of observed dwarf spheroidal galaxies in our Galaxy and the observations of the stellar stream. The upper limit on the amplitude of the bump is O(10^{-7}) for k~10^3 Mpc^{-1}. Furthermore we find the boost factor, which is crucial for the indirect detection of dark matter signals, is up to O(10^4) due to the bump that is allowed in the current observational bounds.

Tidally induced velocity gradients in the Milky Way dwarf spheroidal satellites

ABSTRACT We present a kinematic study of six dwarf spheroidal galaxies (dSph) satellites of the Milky Way (MW), namely Carina, Draco, Fornax, Sculptor, Sextans, and Ursa Minor. We combine proper motions (PMs) from the Gaia Data Release 3 (DR3) and line-of-sight velocities (vlos) from the literature to derive their 3D internal kinematics and to study the presence of internal velocity gradients. We find velocity gradients along the line-of-sight for Carina, Draco, Fornax, and Ursa Minor, at ≥1σ level of significance. The value of such gradients appears to be related to the orbital history of the dwarfs, indicating that the interaction with the MW is causing them. Dwarfs that are close to the MW and moving towards their orbits pericentres show, on average, larger velocity gradients. On the other hand, dwarfs that have recently left their orbits pericentres show no significant gradients. Lastly, dwarfs located at large Galactocentric distances show gradients with an intermediate intensity. Our results would indicate that the torque caused by the strong tidal forces exerted by the MW induces a strong velocity gradient when the dwarfs approach their orbits pericentres. During the pericentre passage, the rapid change in the forces direction would disrupt such gradient, which may steadily recover as the galaxies recede. We assess our findings by analysing dwarfs satellites from the TNG50 simulation. We find a significant increase in the intensity of the detected gradients as the satellites approach their pericentre, followed by a sharp drop as they abandon it, supporting our results for the dSphs of the MW.

Motivations for a large self-interacting dark matter cross-section from Milky Way satellites

ABSTRACT We explore the properties of Milky Way (MW) subhaloes in self-interacting dark matter models for moderate cross-sections of 1–5 cm2 g−1 using high-resolution zoom-in N-body simulations. We include the gravitational potential of a baryonic disc and bulge matched to the MW, which is critical for getting accurate predictions. The predicted number and distribution of subhaloes within the host halo are similar for 1 and 5 cm2 g−1 models, and they agree with observations of MW satellite galaxies only if subhaloes with peak circular velocity over all time >7 km s−1 are able to form galaxies. We do not find distinctive signatures in the pericentre distribution of the subhaloes that could help distinguish the models. Using an analytical model to extend the simulation results, we are able to show that subhaloes in models with cross-sections between 1 and 5 cm2 g−1 are not dense enough to match the densest ultrafaint and classical dwarf spheroidal galaxies in the MW. This motivates exploring velocity-dependent cross-sections with values larger than 5 cm2 g−1 at the velocities relevant for the satellites such that core collapse would occur in some of the ultrafaint and classical dwarf spheroidals.

Dark matter indirect detection limits from complete annihilation patterns

While cosmological and astrophysical probes suggest that dark matter would make up for 85% of the total matter content of the Universe, the determination of its nature remains one of the greatest challenges of fundamental physics. Assuming the ΛCDM cosmological model, Weakly Interacting Massive Particles would annihilate into Standard Model particles, yielding γ-rays, which could be detected by ground-based telescopes. Dwarf spheroidal galaxies represent promising targets for such indirect searches as they are assumed to be highly dark matter dominated with the absence of astrophysical sources nearby. Previous studies have led to upper limits on the annihilation cross-section assuming single exclusive annihilation channels. In this work, we consider a more realistic situation and take into account the complete annihilation pattern within a given particle physics model. This allows us to study the impact on the derived upper limits on the dark matter annihilation cross-section from a full annihilation pattern compared to the case of a single annihilation channel. We use mock data for the Cherenkov Telescope Array simulating the observations of the promising dwarf spheroidal galaxy Sculptor. We show the impact of considering the full annihilation pattern within a simple framework where the Standard Model of particle physics is extended by a singlet scalar. Such a model shows new features in the shape of the predicted upper limit which reaches a value of 〈σv〉 = 3.8 × 10-24 cm-3s-1 for a dark matter mass of 1 TeV at 95% confidence level. We suggest considering the complete particle physics information in order to derive more realistic limits.

The Dark Matter Tidal Stripping History of the Sagittarius Core with N-body Simulations

The infall of the Sagittarius (Sgr) dwarf spheroidal galaxy in the Milky Way halo is an unique opportunity to understand how the different components of a dwarf galaxy could be tidally removed. In this work, we reconstruct the Sgr core morphology and kinematics on the basis of a model that has already successfully reproduced the Sgr stream. Here we use a very high resolution model that almost resolves individual stars in the Sgr core. It reproduces most of the observed morphology and kinematic properties, without specific fine tuning. We also show that the dark matter may have been almost entirely stripped by Milky Way tides after two passages at the pericenter. Finally the model predicts that the Sgr core will be fully disrupted within the next 2 Gyr.

A Deep View into the Nucleus of the Sagittarius Dwarf Spheroidal Galaxy with MUSE. III. Discrete Multicomponent Population-dynamical Models Based on the Jeans Equations

We present comprehensive multicomponent dynamical models of M54 (NGC 6715), the nuclear star cluster of the Sagittarius (Sgr) dwarf galaxy, which is undergoing a tidal disruption in the Milky Way halo. Previous papers in this series used a large MUSE mosaic data set to identify multiple stellar populations in the system and study their kinematic differences. Here, we use Jeans-based dynamical models that fit the population properties (mean age and metallicity), spatial distributions, and kinematics simultaneously. They provide a solid physical explanation for our previous findings. Population-dynamical models deliver a comprehensive view of the whole system, and allow us to disentangle the different stellar populations. We explore their dynamical interplay and confirm our previous findings about the build-up of Sgr’s nuclear cluster via contributions from globular cluster stars, Sgr inner field stars, and in situ star formation. We explore various parameterizations of the gravitational potential and show the importance of a radially varying mass-to-light ratio for the proper treatment of the mass profile. We find a total dynamical mass within M54's tidal radius (∼75 pc) of 1.60 ± 0.07 × 106 M ⊙ in excellent agreement with N-body simulations. Metal-poor globular cluster stars contribute about 65% of the total mass or 1.04 ± 0.05 × 106 M ⊙. Metal-rich stars can be further divided into young and intermediate-age populations, which contribute 0.32 ± 0.02 × 106 M ⊙ (20%) and 0.24 ± 0.02 × 106 M ⊙ (15%), respectively. Our population-dynamical models successfully distinguish the different stellar populations in Sgr’s nucleus because of their different spatial distributions, ages, metallicities, and kinematic features.

The Impact of Binaries on the Dynamical Mass Estimate of Dwarf Galaxies

Binary stars are recognized to be important in driving the dynamical evolution of stellar systems and also in determining some of their observational features. In this study, we explore the role that binary stars have in modulating the estimates of the velocity dispersion of stellar systems. To this aim, we developed a tool that allows investigation of the dependence of synthetic velocity dispersion on a number of crucial quantities characterizing the binary content: binary fraction and the distributions of their mass ratio, eccentricity, and semimajor axis. As an application, we evaluate the impact that binary stars have on the estimation of the dynamical mass of dwarf spheroidal and ultrafaint dwarf galaxies, finding that it can be particularly relevant, especially for low-mass and low-density systems. These results bear profound implications for the interpretation of the measured velocity dispersion in such systems, since it weakens or relieves the claim for the need of large amounts of dark matter.

On the use of field RR Lyrae as galactic probes – VI. Mixed mode RR Lyrae variables in Fornax and in nearby dwarf galaxies

ABSTRACT We investigate the properties of the mixed-mode (RRd) RR Lyrae (RRL) variables in the Fornax dwarf spheroidal (dSph) galaxy by using B- and V-band time series collected over 24 yr. We compare the properties of the RRds in Fornax with those in the Magellanic Clouds and in nearby dSphs, with special focus on Sculptor. We found that the ratio of RRds over the total number of RRLs decreases with metallicity. Typically, dSphs have very few RRds with 0.49≲ P0 ≲0.53 d, but Fornax fills this period gap in the Petersen diagram (ratio between first overtone over fundamental period versus fundamental period). We also found that the distribution in the Petersen diagram of Fornax RRds is similar to Small Magellanic Cloud (SMC) RRds, thus suggesting that their old stars have a similar metallicity distribution. We introduce the Period–Amplitude RatioS diagram, a new pulsation diagnostics independent of distance and reddening. We found that Large Magellanic Cloud (LMC) RRds in this plane are distributed along a short- and a long-period sequence that we identified as the metal-rich and the metal-poor component. These two groups are also clearly separated in the Petersen and Bailey (luminosity amplitude versus logarithmic period) diagrams. This circumstantial evidence indicates that the two groups have different evolutionary properties. All the pulsation diagnostics adopted in this investigation suggest that old stellar populations in Fornax and Sculptor dSphs underwent different chemical enrichment histories. Fornax RRds are similar to SMC RRds, while Sculptor RRds are more similar to the metal-rich component of the LMC RRds.

The Pristine Inner Galaxy Survey (PIGS) – IV. A photometric metallicity analysis of the Sagittarius dwarf spheroidal galaxy

ABSTRACT We present a comprehensive metallicity analysis of the Sagittarius dwarf spheroidal galaxy (Sgr dSph) using Pristine CaHK photometry. We base our member selection on Gaia EDR3 astrometry, applying a magnitude limit at G0 = 17.3, and our population study on the metallicity-sensitive photometry from the Pristine Inner Galaxy Survey (PIGS). Working with photometric metallicities instead of spectroscopic metallicities allows us to cover an unprecedented large area (∼100 square degrees) of the dwarf galaxy, and to study the spatial distribution of its members as function of metallicity with little selection effects. Our study compares the spatial distributions of a metal-poor population of 9719 stars with [Fe/H] <−1.3 and a metal-rich one of 30 115 stars with [Fe/H] >−1.0. The photometric Sgr sample also allows us to assemble the largest sample of 1150 very metal-poor Sgr candidates ([Fe/H] <−2.0). By investigating and fitting the spatial properties of the metal-rich and metal-poor population, we find a negative metallicity gradient which extends up to 12 degrees from the Sgr centre (or ∼5.5 kpc at the distance of Sgr), the limit of our footprint. We conclude that the relative number of metal-poor stars increases in the outer areas of the galaxy, while the central region is dominated by metal-rich stars. These findings suggest an outside-in formation process and are an indication of the extended formation history of Sgr, which has been affected by the tidal interaction between Sgr and the Milky Way.

Constraints on dark matter self-interaction from the internal density profiles of X-COP galaxy clusters

The fundamental properties of the postulated dark matter (DM) affect the internal structure of gravitationally bound structures. In the cold dark matter paradigm, DM particles interact only via gravity. Their distribution is well represented by an Einasto profile with shape parameter α ≈ 0.18 in the smallest dwarf galaxies or the most massive galaxy clusters alike. Conversely, if DM particles self-interact via additional forces, we expect the mass density profiles of DM halos to flatten in their central regions, thereby increasing the Einasto shape parameter. We measured the structural properties of 12 massive galaxy clusters from observations of their hot gaseous atmosphere, using the X-ray observatory XMM-Newton, and of the Sunyaev-Zel’dovich effect using the Planck all-sky survey. After removing morphologically disturbed systems, we measured Einasto shape parameters with mean ⟨α⟩=0.19 ± 0.03 and intrinsic scatter σα = 0.06, which is in close agreement with the prediction of the cold dark matter paradigm. We used cosmological hydrodynamical simulations of cluster formation with self-interacting DM (BAHAMAS-SIDM) to determine how the Einasto shape parameter depends on the self-interaction cross section. We used the fitted relation to turn our measurements of α into constraints on the self-interaction cross section, which imply σ/m < 0.19 cm2 g−1 (95% confidence level) at collision velocity vDM − DM ∼ 1000 km s−1. This is lower than the interaction cross section required for DM self-interactions to solve the core-cusp problem in dwarf spheroidal galaxies, unless the cross section is a strong function of velocity.