There have been several reports of a detection of an unexplained excess of X-ray emission at $\simeq$ 3.5 keV in astrophysical systems. One interpretation of this excess is the decay of sterile neutrino dark matter. The most influential study to date analysed 73 clusters observed by the XMM-Newton satellite. We explore evidence for a $\simeq$ 3.5 keV excess in the XMM-PN spectra of 117 redMaPPer galaxy clusters ($0.1 < z < 0.6$). In our analysis of individual spectra, we identify three systems with an excess of flux at $\simeq$ 3.5 keV. In one case (XCS J0003.3+0204) this excess may result from a discrete emission line. None of these systems are the most dark matter dominated in our sample. We group the remaining 114 clusters into four temperature ($T_{\rm X}$) bins to search for an increase in $\simeq$ 3.5 keV flux excess with $T_{\rm X}$ - a reliable tracer of halo mass. However, we do not find evidence of a significant excess in flux at $\simeq$ 3.5 keV in any $T_{\rm X}$ bins. To maximise sensitivity to a potentially weak dark matter decay feature at $\simeq$ 3.5 keV, we jointly fit 114 clusters. Again, no significant excess is found at $\simeq$ 3.5 keV. We estimate the upper limit of an undetected emission line at $\simeq$ 3.5 keV to be $2.41 \times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$, corresponding to a mixing angle of $\sin^2(2\theta)=4.4 \times 10^{-11}$, lower than previous estimates from cluster studies. We conclude that a flux excess at $\simeq$ 3.5 keV is not a ubiquitous feature in clusters and therefore unlikely to originate from sterile neutrino dark matter decay.

The size of the dust torus in Active Galactic Nuclei (AGN) and their high-luminosity counterparts, quasars, can be inferred from the time delay between UV/optical accretion disk continuum variability and the response in the mid-infrared (MIR) torus emission. This dust reverberation mapping (RM) technique has been successfully applied to $\sim 70$ $z\lesssim 0.3$ AGN and quasars. Here we present first results of our dust RM program for distant quasars covered in the SDSS Stripe 82 region combining $\sim 20$-yr ground-based optical light curves with 10-yr MIR light curves from the WISE satellite. We measure a high-fidelity lag between W1-band (3.4 $\mu$m) and $g$ band for 587 quasars over $0.3\lesssim z\lesssim 2$ ($\left<z\right>\sim 0.8$) and two orders of magnitude in quasar luminosity. They tightly follow (intrinsic scatter $\sim 0.17$ dex in lag) the IR lag-luminosity relation observed for $z<0.3$ AGN, revealing a remarkable size-luminosity relation for the dust torus over more than four decades in AGN luminosity, with little dependence on additional quasar properties such as Eddington ratio and variability amplitude. This study motivates further investigations in the utility of dust RM for cosmology, and strongly endorses a compelling science case for the combined 10-yr Vera C. Rubin Observatory Legacy Survey of Space and Time (optical) and 5-yr Nancy Grace Roman Space Telescope 2$\mu$m light curves in a deep survey for low-redshift AGN dust RM with much lower luminosities and shorter, measurable IR lags. The compiled optical and MIR light curves for 7,384 quasars in our parent sample are made public with this work.

We present a measurement of the Hubble constant $H_0$ using the gravitational wave (GW) event GW190814, which resulted from the coalescence of a 23 $M_\odot$ black hole with a 2.6 $M_\odot$ compact object, as a standard siren. No compelling electromagnetic counterpart with associated host galaxy has been identified for this event, thus our analysis accounts for $\sim$ 2,700 potential host galaxies within a statistical framework. The redshift information is obtained from the photometric redshift (photo-$z$) catalog from the Dark Energy Survey. The luminosity distance is provided by the gravitational wave sky map published by the LIGO/Virgo Collaboration. Since this GW event has the second-smallest sky localization area after GW170817, GW190814 is likely to provide the best constraint on cosmology from a single standard siren without identifying an electromagnetic counterpart. Our analysis uses photo-$z$ probability distribution functions and corrects for photo-$z$ biases. We also reanalyze the binary-black hole GW170814 within this updated framework. We explore how our findings impact the $H_0$ constraints from GW170817, the only GW merger associated with a unique host galaxy, and therefore the most powerful standard siren to date. From a combination of GW190814, GW170814 and GW170817, our analysis yields $H_0 = 69.0^{+ 14}_{- 7.5 }~{\rm km~s^{-1}~Mpc^{-1}}$ (68% Highest Density Interval, HDI) for a prior in $H_0$ uniform between $[20,140]~{\rm km~s^{-1}~Mpc^{-1}}$. The addition of GW190814 and GW170814 to GW170817 improves the 68% HDI from GW170817 alone by $\sim 12\%$, showing how well-localized mergers without counterparts can provide a marginal contribution to standard siren measurements, provided that a complete galaxy catalog is available at the location of the event.

We present a description of the Australian Dark Energy Survey (OzDES) and summarise the results from its six years of operations. Using the 2dF fibre positioner and AAOmega spectrograph on the 3.9-metre Anglo-Australian Telescope, OzDES has monitored 771 AGN, classified hundreds of supernovae, and obtained redshifts for thousands of galaxies that hosted a transient within the 10 deep fields of the Dark Energy Survey. We also present the second OzDES data release, containing the redshifts of almost 30,000 sources, some as faint as $r_{\mathrm AB}=24$ mag, and 375,000 individual spectra. These data, in combination with the time-series photometry from the Dark Energy Survey, will be used to measure the expansion history of the Universe out to $z\sim1.2$ and the masses of hundreds of black holes out to $z\sim4$. OzDES is a template for future surveys that combine simultaneous monitoring of targets with wide-field imaging cameras and wide-field multi-object spectrographs.

We present DES16C3cje, a low-luminosity, long-lived type II supernova (SN II) at redshift 0.0618, detected by the Dark Energy Survey (DES). DES16C3cje is a unique SN. The spectra are characterized by extremely narrow photospheric lines corresponding to very low expansion velocities of $\lesssim1500$ km s$^{-1}$, and the light curve shows an initial peak that fades after 50 days before slowly rebrightening over a further 100 days to reach an absolute brightness of M$_r\sim -15.5$ mag. The decline rate of the late-time light curve is then slower than that expected from the powering by radioactive decay of $^{56}$Co but is comparable to that expected from accretion power. Comparing the bolometric light curve with hydrodynamical models, we find that DES16C3cje can be explained by either i) a low explosion energy (0.11 foe) and relatively large $^{56}$Ni production of 0.075 M$_{\odot}$ from a $\sim15$ M$_{\odot}$ red supergiant progenitor typical of other SNe II, or ii) a relatively compact $\sim40$ M$_{\odot}$ star, explosion energy of 1 foe, and 0.08 M$_{\odot}$ of $^{56}$Ni. Both scenarios require additional energy input to explain the late-time light curve, which is consistent with fallback accretion at a rate of $\sim0.5\times{10^{-8}}$ M$_{\odot}$ s$^{-1}$.

We present a detailed abundance analysis of the three brightest member stars at the top of the giant branch of the ultra-faint dwarf galaxy Grus~II. All stars exhibit a higher than expected $\mathrm{[Mg/Ca]}$ ratio compared to metal-poor stars in other ultra-faint dwarf galaxies and in the Milky Way halo. Nucleosynthesis in high mass ($\geqslant 20$M$_\odot$) core-collapse supernovae has been shown to create this signature. The abundances of this small sample (3) stars suggest the chemical enrichment of Grus~II could have occurred through substantial high-mass stellar evolution and is consistent with the framework of a top-heavy initial mass function. However, with only three stars it can not be ruled out that the abundance pattern is the result of a stochastic chemical enrichment at early times in the galaxy. The most metal-rich of the three stars also possesses a small enhancement in rapid neutron-capture ($r$-process) elements. The abundance pattern of the $r$-process elements in this star matches the scaled $r$-process pattern of the solar system and $r$-process enhanced stars in other dwarf galaxies and in the Milky Way halo, hinting at a common origin for these elements across a range of environments. All current proposed astrophysical sites of $r$-process element production are associated with high-mass stars, thus the possible top-heavy initial mass function of Grus~II would increase the likelihood of any of these events occurring. The time delay between the $\alpha$ and $r$-process element enrichment of the galaxy favors a neutron star merger as the origin of the $r$-process elements in Grus~II.

Despite vast improvements in the measurement of the cosmological parameters, the nature of dark energy and an accurate value of the Hubble constant (H$_0$) in the Hubble-Lema\^itre law remain unknown. To break the current impasse, it is necessary to develop as many independent techniques as possible, such as the use of Type II supernovae (SNe II). The goal of this paper is to demonstrate the utility of SNe II for deriving accurate extragalactic distances, which will be an asset for the next generation of telescopes where more-distant SNe II will be discovered. More specifically, we present a sample from the Dark Energy Survey Supernova Program (DES-SN) consisting of 15 SNe II with photometric and spectroscopic information spanning a redshift range up to 0.35. Combining our DES SNe with publicly available samples, and using the standard candle method (SCM), we construct the largest available Hubble diagram with SNe II in the Hubble flow (70 SNe II) and find an observed dispersion of 0.27 mag. We demonstrate that adding a colour term to the SN II standardisation does not reduce the scatter in the Hubble diagram. Although SNe II are viable as distance indicators, this work points out important issues for improving their utility as independent extragalactic beacons: find new correlations, define a more standard subclass of SNe II, construct new SN II templates, and dedicate more observing time to high-redshift SNe II. Finally, for the first time, we perform simulations to estimate the redshift-dependent distance-modulus bias due to selection effects.

We use imaging from the first three years of the Dark Energy Survey to characterize the dynamical state of 288 galaxy clusters at $0.1 \lesssim z \lesssim 0.9$ detected in the South Pole Telescope (SPT) Sunyaev-Zeldovich (SZ) effect survey (SPT-SZ). We examine spatial offsets between the position of the brightest cluster galaxy (BCG) and the center of the gas distribution as traced by the SPT-SZ centroid and by the X-ray centroid/peak position from Chandra and XMM data. We show that the radial distribution of offsets provides no evidence that SPT SZ-selected cluster samples include a higher fraction of mergers than X-ray-selected cluster samples. We use the offsets to classify the dynamical state of the clusters, selecting the 43 most disturbed clusters, with half of those at $z \gtrsim 0.5$, a region seldom explored previously. We find that Schechter function fits to the galaxy population in disturbed clusters and relaxed clusters differ at $z>0.55$ but not at lower redshifts. Disturbed clusters at $z>0.55$ have steeper faint-end slopes and brighter characteristic magnitudes. Within the same redshift range, we find that the BCGs in relaxed clusters tend to be brighter than the BCGs in disturbed samples, while in agreement in the lower redshift bin. Possible explanations includes a higher merger rate, and a more efficient dynamical friction at high redshift. The red-sequence population is less affected by the cluster dynamical state than the general galaxy population.

For ground-based optical imaging with current CCD technology, the Poisson fluctuations in source and sky background photon arrivals dominate the noise budget and are readily estimated. Another component of noise, however, is the signal from the undetected population of stars and galaxies. Using injection of artificial galaxies into images, we demonstrate that the measured variance of galaxy moments (used for weak gravitational lensing measurements) in Dark Energy Survey (DES) images is significantly in excess of the Poisson predictions, by up to 30\%, and that the background sky levels are overestimated by current software. By cross-correlating distinct images of "empty" sky regions, we establish that there is a significant image noise contribution from undetected static sources (US), which on average are mildly resolved at DES resolution. Treating these US as a stationary noise source, we compute a correction to the moment covariance matrix expected from Poisson noise. The corrected covariance matrix matches the moment variances measured on the injected DES images to within 5\%. Thus we have an empirical method to statistically account for US in weak lensing measurements, rather than requiring extremely deep sky simulations. We also find that local sky determinations can remove the bias in flux measurements, at a small penalty in additional, but quantifiable, noise.

Obtaining accurate distributions of galaxy redshifts is a critical aspect of weak lensing cosmology experiments. One of the methods used to estimate and validate redshift distributions is apply weights to a spectroscopic sample so that their weighted photometry distribution matches the target sample. In this work we estimate the \textit{selection bias} in redshift that is introduced in this procedure. We do so by simulating the process of assembling a spectroscopic sample (including observer-assigned confidence flags) and highlight the impacts of spectroscopic target selection and redshift failures. We use the first year (Y1) weak lensing analysis in DES as an example data set but the implications generalise to all similar weak lensing surveys. We find that using colour cuts that are not available to the weak lensing galaxies can introduce biases of $\Delta~z\sim0.015$ in the weighted mean redshift of different redshift intervals. To assess the impact of incompleteness in spectroscopic samples, we select only objects with high observer-defined confidence flags and compare the weighted mean redshift with the true mean. We find that the mean redshift of the DES Y1 weak lensing sample is typically biased at the $\Delta~z=0.005-0.05$ level after the weighting is applied. The bias we uncover can have either sign, depending on the samples and redshift interval considered. For the highest redshift bin, the bias is larger than the uncertainties in the other DES Y1 redshift calibration methods, justifying the decision of not using this method for the redshift estimations. We discuss several methods to mitigate this bias.

We construct and validate the selection function of the MARD-Y3 sample. This sample was selected through optical follow-up of the 2nd ROSAT faint source catalog (2RXS) with Dark Energy Survey year 3 (DES-Y3) data. The selection function is modeled by combining an empirically constructed X-ray selection function with an incompleteness model for the optical follow-up. We validate the joint selection function by testing the consistency of the constraints on the X-ray flux--mass and richness--mass scaling relation parameters derived from different sources of mass information: (1) cross-calibration using SPT-SZ clusters, (2) calibration using number counts in X-ray, in optical and in both X-ray and optical while marginalizing over cosmological parameters, and (3) other published analyses. We find that the constraints on the scaling relation from the number counts and SPT-SZ cross-calibration agree, indicating that our modeling of the selection function is adequate. Furthermore, we apply a largely cosmology independent method to validate selection functions via the computation of the probability of finding each cluster in the SPT-SZ sample in the MARD-Y3 sample and vice-versa. This test reveals no clear evidence for MARD-Y3 contamination, SPT-SZ incompleteness or outlier fraction. Finally, we discuss the prospects of the techniques presented here to limit systematic selection effects in future cluster cosmological studies.

We perform a joint analysis of the counts and weak lensing signal of redMaPPer clusters selected from the Dark Energy Survey (DES) Year 1 dataset. Our analysis uses the same shear and source photometric redshifts estimates as were used in the DES combined probes analysis. Our analysis results in surprisingly low values for $S_8 =\sigma_8(\Omega_{\rm m}/0.3)^{0.5}= 0.65\pm 0.04$, driven by a low matter density parameter, $\Omega_{\rm m}=0.179^{+0.031}_{-0.038}$, with $\sigma_8-\Omega_{\rm m}$ posteriors in $2.4\sigma$ tension with the DES Y1 3x2pt results, and in $5.6\sigma$ with the Planck CMB analysis. These results include the impact of post-unblinding changes to the analysis, which did not improve the level of consistency with other data sets compared to the results obtained at the unblinding. The fact that multiple cosmological probes (supernovae, baryon acoustic oscillations, cosmic shear, galaxy clustering and CMB anisotropies), and other galaxy cluster analyses all favor significantly higher matter densities suggests the presence of systematic errors in the data or an incomplete modeling of the relevant physics. Cross checks with X-ray and microwave data, as well as independent constraints on the observable--mass relation from SZ selected clusters, suggest that the discrepancy resides in our modeling of the weak lensing signal rather than the cluster abundance. Repeating our analysis using a higher richness threshold ($\lambda \ge 30$) significantly reduces the tension with other probes, and points to one or more richness-dependent effects not captured by our model.

While many studies have shown a correlation between properties of the light curves of Type Ia SN (SNe Ia) and properties of their host galaxies, it remains unclear what is driving these correlations. We introduce a new direct method to study these correlations by analyzing `parent' galaxies that host multiple SNe Ia 'siblings'. Here, we search the Dark Energy Survey SN sample, one of the largest samples of discovered SNe, and find 8 galaxies that hosted two likely Type Ia SNe. Comparing the light-curve properties of these SNe and recovered distances from the light curves, we find no better agreement between properties of SNe in the same galaxy as any random pair of galaxies, with the exception of the SN light-curve stretch. We show at $2.8\sigma$ significance that at least 1/2 of the intrinsic scatter of SNe Ia distance modulus residuals is not from common host properties. We also discuss the robustness with which we could make this evaluation with LSST, which will find $100\times$ more pairs of galaxies, and pave a new line of study on the consistency of Type Ia supernovae in the same parent galaxies. Finally, we argue that it is unlikely some of these SNe are actually single, lensed SN with multiple images.

The astrophysical $s$-process is one of the two main processes forming elements heavier than iron. A key outstanding uncertainty surrounding $s$-process nucleosynthesis is the neutron flux generated by the ${}^{22}\mathrm{Ne}(\alpha, n){}^{25}\mathrm{Mg}$ reaction during the He-core and C-shell burning phases of massive stars. This reaction, as well as the competing ${}^{22}\mathrm{Ne}(\alpha, \gamma){}^{26}\mathrm{Mg}$ reaction, is not well constrained in the important temperature regime from ${\sim} 0.2$--$0.4$~GK, owing to uncertainties in the nuclear properties of resonances lying within the Gamow window. To address these uncertainties, we have performed a new measurement of the ${}^{22}\mathrm{Ne}({}^{6}\mathrm{Li}, d){}^{26}\mathrm{Mg}$ reaction in inverse kinematics, detecting the outgoing deuterons and ${}^{25,26}\mathrm{Mg}$ recoils in coincidence. We have established a new $n / \gamma$ decay branching ratio of $1.14(26)$ for the key $E_x = 11.32$ MeV resonance in $^{26}\mathrm{Mg}$, which results in a new $(\alpha, n)$ strength for this resonance of $42(11)~\mu$eV when combined with the well-established $(\alpha, \gamma)$ strength of this resonance. We have also determined new upper limits on the $\alpha$ partial widths of neutron-unbound resonances at $E_x = 11.112,$ $11.163$, $11.169$, and $11.171$ MeV. Monte-Carlo calculations of the stellar ${}^{22}\mathrm{Ne}(\alpha, n){}^{25}\mathrm{Mg}$ and ${}^{22}\mathrm{Ne}(\alpha, \gamma){}^{26}\mathrm{Mg}$ rates, which incorporate these results, indicate that both rates are substantially lower than previously thought in the temperature range from ${\sim} 0.2$--$0.4$~GK.

The five-year Dark Energy Survey supernova programme (DES-SN) is one of the largest and deepest transient surveys to date in terms of volume and number of supernovae. Identifying and characterising the host galaxies of transients plays a key role in their classification, the study of their formation mechanisms, and the cosmological analyses. To derive accurate host galaxy properties, we create depth-optimised coadds using single-epoch DES-SN images that are selected based on sky and atmospheric conditions. For each of the five DES-SN seasons, a separate coadd is made from the other 4 seasons such that each SN has a corresponding deep coadd with no contaminating SN emission. The coadds reach limiting magnitudes of order $\sim 27$ in $g$-band, and have a much smaller magnitude uncertainty than the previous DES-SN host templates, particularly for faint objects. We present the resulting multi-band photometry of host galaxies for samples of spectroscopically confirmed type Ia (SNe Ia), core-collapse (CCSNe), and superluminous (SLSNe) as well as rapidly evolving transients (RETs) discovered by DES-SN. We derive host galaxy stellar masses and probabilistically compare stellar-mass distributions to samples from other surveys. We find that the DES spectroscopically confirmed sample of SNe Ia selects preferentially fewer high mass hosts at high redshift compared to other surveys, while at low redshift the distributions are consistent. DES CCSNe and SLSNe hosts are similar to other samples, while RET hosts are unlike the hosts of any other transients, although these differences have not been disentangled from selection effects.

We introduce a galaxy cluster mass observable, $\mu_\star$, based on the stellar masses of cluster members, and we present results for the Dark Energy Survey (DES) Year 1 observations. Stellar masses are computed using a Bayesian Model Averaging method, and are validated for DES data using simulations and COSMOS data. We show that $\mu_\star$ works as a promising mass proxy by comparing our predictions to X-ray measurements. We measure the X-ray temperature-$\mu_\star$ relation for a total of 150 clusters matched between the wide-field DES Year 1 redMaPPer catalogue, and Chandra and XMM archival observations, spanning the redshift range $0.1<z<0.7$. For a scaling relation which is linear in logarithmic space, we find a slope of $\alpha = 0.488\pm0.043$ and a scatter in the X-ray temperature at fixed $\mu_\star$ of $\sigma_{{\rm ln} T_X|\mu_\star}=0.266^{+0.019}_{-0.020}$ for the joint sample. By using the halo mass scaling relations of the X-ray temperature from the Weighing the Giants program, we further derive the $\mu_\star$-conditioned scatter in mass, finding $\sigma_{{\rm ln} M|\mu_\star}=0.26^{+ 0.15}_{- 0.10}$. These results are competitive with well-established cluster mass proxies used for cosmological analyses, showing that $\mu_\star$ can be used as a reliable and physically motivated mass proxy to derive cosmological constraints.

The center determination of a galaxy cluster from an optical cluster finding algorithm can be offset from theoretical prescriptions or $N$-body definitions of its host halo center. These offsets impact the recovered cluster statistics, affecting both richness measurements and the weak lensing shear profile around the clusters. This paper models the centering performance of the \RM~cluster finding algorithm using archival X-ray observations of \RM-selected clusters. Assuming the X-ray emission peaks as the fiducial halo centers, and through analyzing their offsets to the \RM~centers, we find that $\sim 75\pm 8 \%$ of the \RM~clusters are well centered and the mis-centered offset follows a Gamma distribution in normalized, projected distance. These mis-centering offsets cause a systematic underestimation of cluster richness relative to the well-centered clusters, for which we propose a descriptive model. Our results enable the DES Y1 cluster cosmology analysis by characterizing the necessary corrections to both the weak lensing and richness abundance functions of the DES Y1 redMaPPer cluster catalog.

Using archival X-ray observations and a log-normal population model, we estimate constraints on the intrinsic scatter in halo mass at fixed optical richness for a galaxy cluster sample identified in Dark Energy Survey Year-One (DES-Y1) data with the redMaPPer algorithm. We examine the scaling behavior of X-ray temperatures, $T_X$, with optical richness, $\lambda_{RM}$, for clusters in the redshift range $0.2<z<0.7$. X-ray temperatures are obtained from Chandra and XMM observations for 58 and 110 redMaPPer systems, respectively. Despite non-uniform sky coverage, the $T_X$ measurements are $> 50\%$ complete for clusters with $\lambda_{RM} > 130$. Regression analysis on the two samples produces consistent posterior scaling parameters, from which we derive a combined constraint on the residual scatter, $\sigma_{\ln Tx | \lambda} = 0.275 \pm 0.019$. Joined with constraints for $T_X$ scaling with halo mass from the Weighing the Giants program and richness--temperature covariance estimates from the LoCuSS sample, we derive the richness-conditioned scatter in mass, $\sigma_{\ln M | \lambda} = 0.30 \pm 0.04\, _{({\rm stat})} \pm 0.09\, _{({\rm sys})}$, at an optical richness of approximately 70. Uncertainties in external parameters, particularly the slope and variance of the $T_X$--mass relation and the covariance of $T_X$ and $\lambda_{RM}$ at fixed mass, dominate the systematic error. The $95\%$ confidence region from joint sample analysis is relatively broad, $\sigma_{\ln M | \lambda} \in [0.14, \, 0.55]$, or a factor ten in variance.

We use present cosmological observations and forecasts of future experiments to illustrate the power of large-scale structure (LSS) surveys in probing dark matter (DM) microphysics and unveiling potential deviations from the standard $\Lambda$CDM scenario. To quantify this statement, we focus on an extension of $\Lambda$CDM with DM-neutrino scattering, which leaves a distinctive imprint on the angular and matter power spectra. After finding that future CMB experiments (such as COrE+) will not significantly improve the constraints set by the Planck satellite, we show that the next generation of galaxy clustering surveys (such as DESI) could play a leading role in constraining alternative cosmologies and even have the potential to make a discovery. Typically we find that DESI would be an order of magnitude more sensitive to DM interactions than Planck, thus probing effects that until now have only been accessible via $N$-body simulations.

Over the past few decades, an anomalous 511 keV gamma-ray line has been observed from the centre of the Milky Way. Dark matter (DM) in the form of light weakly interacting massive particles (WIMPs) annihilating into electron-positron pairs has been one of the leading hypotheses of the observed emission. Given the small required cross section, a further coupling to lighter particles is required to produce the correct relic density. Here, we derive constraints from the Planck satellite on light WIMPs that were in equilibrium with either the neutrino or electron sector in the early universe. For the neutrino sector, we obtain a lower bound on the WIMP mass of 4 MeV for a real scalar and 10 MeV for a Dirac fermion DM particle, at 95% CL. For the electron sector, we find even stronger bounds of 7 and 11 MeV, respectively. Using these results, we show that, in the absence of additional ingredients such as dark radiation, the light thermally produced WIMP explanation of the 511 keV excess is strongly disfavoured by the latest cosmological data. This suggests an unknown astrophysical or more exotic DM source of the signal.

In the thermal dark matter (DM) paradigm, primordial interactions between DM and Standard Model particles are responsible for the observed DM relic density. In Boehm et al. (2014), we showed that weak-strength interactions between DM and radiation (photons or neutrinos) can erase small-scale density fluctuations, leading to a suppression of the matter power spectrum compared to the collisionless cold DM (CDM) model. This results in fewer DM subhaloes within Milky Way-like DM haloes, implying a reduction in the abundance of satellite galaxies. Here we use very high resolution N-body simulations to measure the dynamics of these subhaloes. We find that when interactions are included, the largest subhaloes are less concentrated than their counterparts in the collisionless CDM model and have rotation curves that match observational data, providing a new solution to the "too big to fail" problem.

Interactions between dark matter (DM) and radiation (photons or neutrinos) in the early Universe suppress density fluctuations on small mass scales. Here we perform a thorough analysis of structure formation in the fully non-linear regime using N-body simulations for models with DM-radiation interactions and compare the results to a traditional calculation in which DM only interacts gravitationally. Significant differences arise due to the presence of interactions, in terms of the number of low-mass DM haloes and their properties, such as their spin and density profile. These differences are clearly seen even for haloes more massive than the scale on which density fluctuations are suppressed. We also show that semi-analytical descriptions of the matter distribution in the non-linear regime fail to reproduce our numerical results, emphasizing the challenge of predicting structure formation in models with physics beyond collisionless DM.

The cold dark matter (CDM) model faces persistent challenges on small scales. In particular, taken at face value, the model significantly overestimates the number of satellite galaxies around the Milky Way. Attempts to solve this problem remain open to debate and have even led some to abandon CDM altogether. However, current simulations are limited by the assumption that dark matter feels only gravity. Here, we show that including interactions between CDM and radiation (photons or neutrinos) leads to a dramatic reduction in the number of satellite galaxies, alleviating the Milky Way satellite problem and indicating that physics beyond gravity may be essential to make accurate predictions of structure formation on small scales. The methodology introduced here gives constraints on dark matter interactions that are significantly improved over those from the cosmic microwave background.

In this paper, we explore the impact of Dark Matter-photon interactions on the CMB angular power spectrum. Using the one-year data release of the Planck satellite, we derive an upper bound on the Dark Matter-photon elastic scattering cross section of sigma_{DM-photon} < 8 x 10^{-31} (m_DM/GeV) cm^2 (68% CL) if the cross section is constant and a present-day value of sigma_{DM-photon} < 6 x 10^{-40} (m_DM/GeV) cm^2 (68% CL) if it scales as the temperature squared. For such a limiting cross section, both the B-modes and the TT angular power spectrum are suppressed with respect to LCDM predictions for l > 500 and l > 3000 respectively, indicating that forthcoming data from CMB polarisation experiments and Planck could help to constrain and characterise the physics of the dark sector. This essentially initiates a new type of dark matter search that is independent of whether dark matter is annihilating, decaying or asymmetric. Thus, any CMB experiment with the ability to measure the temperature and/or polarisation power spectra at high l should be able to investigate the potential interactions of dark matter and contribute to our fundamental understanding of its nature.

We present a new study on the elastic scattering cross section of dark matter (DM) and neutrinos using the latest cosmological data from Planck and large-scale structure experiments. We find that the strongest constraints are set by the Lyman-alpha forest, giving sigma_{DM-neutrino} < 10^{-33} (m_DM/GeV) cm^2 if the cross section is constant and a present-day value of sigma_{DM-neutrino} < 10^{-45} (m_DM/GeV) cm^2 if it scales as the temperature squared. These are the most robust limits on DM-neutrino interactions to date, demonstrating that one can use the distribution of matter in the Universe to probe dark ("invisible") interactions. Additionally, we show that scenarios involving thermal MeV DM and a constant elastic scattering cross section naturally predict (i) a cut-off in the matter power spectrum at the Lyman-alpha scale, (ii) N_eff ~ 3.5 +/- 0.4, (iii) H_0 ~ 71 +/- 3 km/s/Mpc and (iv) the possible generation of neutrino masses.