In this paper we study the density profiles of self-interacting dark matter (SIDM) haloes spanning the full observable mass range, from dwarf galaxies to galaxy clusters. Using realistic simulations that model the baryonic physics relevant for galaxy formation, we compare the density profiles of haloes simulated with either SIDM or cold and collisionless dark matter (CDM) to those inferred from observations of stellar velocity dispersion, gas rotation curves, weak and strong gravitational lensing, and/or X-ray maps. We make our comparison in terms of the maximal surface density of haloes, circumventing the need for semi-analytic or parametric models for dark matter density profiles. We find that the maximal surface density as a function of halo mass is well reproduced by simulations of CDM that include baryons, and that a SIDM cross-section of 1 cm$^2/$g is marginally excluded. This constraint is driven mainly by galaxy clusters, with less massive systems being relatively less affected by DM self-interactions with velocity-independent cross-section.

An unidentified line at energy around 3.5 keV was detected in the spectra of dark matter-dominated objects. Recent work of Dessert et al. [1812.06976] used 30 Msec of XMM-Newton blank-sky observations to constrain the admissible line flux, challenging its dark matter decay origin. We demonstrate that these bounds are overestimated by more than an order of magnitude due to improper background modeling. Therefore the dark matter interpretation of the 3.5 keV signal remains viable.

Traditional direct searches for dark matter, looking for nuclear recoils in deep underground detectors, are challenged by an almost complete loss of sensitivity for light dark matter particles. Consequently, there is a significant effort in the community to devise new methods and experiments to overcome these difficulties, constantly pushing the limits of the lowest dark matter mass that can be probed this way. From a model-building perspective, the scattering of sub-GeV dark matter on nucleons essentially must proceed via new light mediator particles, given that collider searches place extremely stringent bounds on contact-type interactions. Here we present an updated compilation of relevant limits for the case of a scalar mediator, including a new estimate of the near-future sensitivity of the NA62 experiment as well as a detailed evaluation of the model-specific limits from Big Bang nucleosynthesis. We also derive updated and more general limits on DM particles upscattered by cosmic rays, applicable to arbitrary energy- and momentum dependences of the scattering cross section. Finally we stress that dark matter self-interactions, when evaluated beyond the common s-wave approximation, place stringent limits independently of the dark matter production mechanism. These are, for the relevant parameter space, generically comparable to those that apply in the commonly studied freeze-out case. We conclude that the combination of existing (or expected) constraints from accelerators and astrophysics, combined with cosmological requirements, puts robust limits on the maximally possible nuclear scattering rate. In most regions of parameter space these are at least competitive with the best projected limits from currently planned direct detection experiments.

The Lyman-$\alpha$ forest is a powerful tool to constrain warm dark matter models (WDM). Its main observable -- flux power spectrum -- should exhibit a suppression at small scales in WDM models. This suppression, however, can be mimicked by a number of thermal effects related to the instantaneous temperature of the intergalactic medium (IGM), and to the history of reionization and of the IGM heating ("pressure effects"). Therefore, to put robust bounds on WDM one needs to disentangle the effect of free-streaming of dark matter particles from the influence of all astrophysical effects. This task cannot be brute-forced due to the complexity of the IGM modelling. In this work, we model the sample of high-resolution and high-redshift quasar spectra (Boera et al 2018) assuming a thermal history that leads to the smallest pressure effects while still being broadly compatible with observations. We explicitly marginalize over observationally allowed values of IGM temperature and find that (thermal) WDM models with masses above 1.9 keV (at 95% CL) are consistent with the spatial shape of the observed flux power spectrum at $z=4-5$. Even warmer models would produce a suppression at scales that are larger than observed, independently of assumptions about thermal effects. This bound is significantly lower than previously claimed bounds, demonstrating the importance of the knowledge about the reionization history and of the proper marginalization over unknowns.

Observations of the redshifted 21-cm signal (in absorption or emission) allow us to peek into the epoch of "dark ages" and the onset of reionization. These data can provide a novel way to learn about the nature of dark matter, in particular about the formation of small size dark matter halos. However, the connection between the formation of structures and 21-cm signal requires knowledge of stellar to total mass relation, escape fraction of UV photons, and other parameters that describe star formation and radiation at early times. This baryonic physics depends on the properties of dark matter and in particular in warm-dark-matter (WDM) models, star formation may follow a completely different scenario, as compared to the cold-dark-matter case. We use the recent measurements by the EDGES [J. D. Bowman, A. E. E. Rogers, R. A. Monsalve, T. J. Mozdzen, and N. Mahesh, An absorption profile centred at 78 megahertz in thesky-averaged spectrum,Nature (London) 555, 67 (2018).] to demonstrate that when taking the above considerations into account, the robust WDM bounds are in fact weaker than those given by the Lyman-$\alpha$ forest method and other structure formation bounds. In particular, we show that resonantly produced 7 keV sterile neutrino dark matter model is consistent with these data. However, a holistic approach to modelling of the WDM universe holds great potential and may in the future make 21-cm data our main tool to learn about dark matter clustering properties.

Flux ratio anomalies in strong gravitationally lensed quasars constitute a unique way to probe the abundance of non-luminous dark matter haloes, and hence the nature of dark matter. In this paper we identify double imaged quasars as a statistically efficient probe of dark matter, since they are 20 times more abundant than quadruply imaged quasars. Using N-body simulations that include realistic baryonic feedback, we measure the full distribution of flux ratios in doubly imaged quasars for cold (CDM) and warm dark matter (WDM) cosmologies. Through this method, we fold in two key systematics - quasar variability and line-of-sight structures. We find that WDM cosmologies predict a ~6 per cent difference in the cumulative distribution functions of flux ratios relative to CDM, with CDM predicting many more small ratios. Finally, we estimate that ~600 doubly imaged quasars will need to be observed in order to be able to unambiguously discern between CDM and the two WDM models studied here. Such sample sizes will be easily within reach of future large scale surveys such as Euclid. In preparation for this survey data we require discerning the scale of the uncertainties in modelling lens galaxies and their substructure in simulations, plus a strong understanding of the selection function of observed lensed quasars.

We investigate how hypothetical particles - sterile neutrinos - can be produced in the interior of exploding supernovae via the resonant conversion of $\bar\nu_\mu$ and $\bar \nu_\tau$. The novelty of our treatment lies in the proper account of the resulting lepton number diffusion. We compute the yield of sterile neutrinos and find that even after taking into account back reaction, sterile neutrinos can carry out a sizeable fraction of the total energy of the explosion comparable to that of active neutrinos. The production is, however, exponentially sensitive to the temperature in the inner supernovae regions, making robust predictions of challenging. In order to understand whether this production affects supernova evolution and can therefore be constrained, detailed simulations including the effects of sterile neutrinos are needed.

The observed Lyman-$\alpha$ flux power spectrum (FPS) is suppressed on scales below $\sim~ 30~{\rm km~s}^{-1}$. This cutoff could be due to the high temperature, $T_0$, and pressure, $p_0$, of the absorbing gas or, alternatively, it could reflect the free streaming of dark matter particles in the early universe. We perform a set of very high resolution cosmological hydrodynamic simulations in which we vary $T_0$, $p_0$ and the amplitude of the dark matter free streaming, and compare the FPS of mock spectra to the data. We show that the location of the dark matter free-streaming cutoff scales differently with redshift than the cutoff produced by thermal effects and is more pronounced at higher redshift. We, therefore, focus on a comparison to the observed FPS at $z>5$. We demonstrate that the FPS cutoff can be fit assuming cold dark matter, but it can be equally well fit assuming that the dark matter consists of $\sim 7$ keV sterile neutrinos in which case the cutoff is due primarily to the dark matter free streaming.

We report a detection of 3.5 keV line in the Milky Way in 5 regions offset from the Galactic Center by distances from 10' to 35 degrees. We build an angular profile of this line and compare it with profiles of several astrophysical lines detected in the same observations. We compare our results with other detections and bounds previously obtained using observations of the Milky Way.

In this work, we discuss the relation between the strength of the self-interaction of dark matter particles and the predicted properties of the inner density distributions of dark matter haloes. We present the results of $N$-body simulations for 28 haloes performed with the same initial conditions for cold dark matter and for self-interacting dark matter, for a range of cross-sections. We provide a simple phenomenological description of these results and compare with the semi-analytical model typically used in the literature. Using these results, we then predict how the inner dark matter surface density and core radius should depend on the self-interaction cross-section for observed haloes.

We review sterile neutrinos as possible Dark Matter candidates. After a short summary on the role of neutrinos in cosmology and particle physics, we give a comprehensive overview of the current status of the research on sterile neutrino Dark Matter. First we discuss the motivation and limits obtained through astrophysical observations. Second, we review different mechanisms of how sterile neutrino Dark Matter could have been produced in the early universe. Finally, we outline a selection of future laboratory searches for keV-scale sterile neutrinos, highlighting their experimental challenges and discovery potential.

Using direct numerical simulations (DNS), we study laminar and turbulent dynamos in chiral magnetohydrodynamics (MHD) with an extended set of equations that accounts for an additional contribution to the electric current due to the chiral magnetic effect (CME). This quantum phenomenon originates from an asymmetry between left- and right-handed relativistic fermions in the presence of a magnetic field and gives rise to a chiral dynamo. We show that the magnetic field evolution proceeds in three stages: (1) a small-scale chiral dynamo instability; (2) production of chiral magnetically driven turbulence and excitation of a large-scale dynamo instability due to a new chiral effect (alpha_mu effect); and (3) saturation of magnetic helicity and magnetic field growth controlled by a conservation law for the total chirality. The $\alpha_\mu$ effect becomes dominant at large fluid and magnetic Reynolds numbers and is not related to kinetic helicity. The growth rate of the large-scale magnetic field and its characteristic scale measured in the numerical simulations agree well with theoretical predictions based on mean-field theory. The previously discussed two-stage chiral magnetic scenario did not include stage (2) during which the characteristic scale of magnetic field variations can increase by many orders of magnitude. Based on the findings from numerical simulations, the relevance of the CME and the chiral effects revealed in the relativistic plasma of the early universe and of proto-neutron stars are discussed.

The observed surface densities of dark matter halos are known to follow a simple scaling law, ranging from dwarf galaxies to galaxy clusters, with a weak dependence on their virial mass. Here we point out that this can not only be used to provide a method to determine the standard relation between halo mass and concentration, but also to use large samples of objects in order to place constraints on dark matter self-interactions that can be more robust than constraints derived from individual objects. We demonstrate our method by considering a sample of about 50 objects distributed across the whole halo mass range, and by modelling the effect of self-interactions in a way similar to what has been previously done in the literature. Using additional input from simulations then results in a constraint on the self-interaction cross section per unit dark matter mass of about $\sigma/m_\chi\lesssim 0.3$ cm$^2$/g. We expect that these constraints can be significantly improved in the future, and made more robust, by i) an improved modelling of the effect of self-interactions, both theoretical and by comparison with simulations, ii) taking into account a larger sample of objects and iii) by reducing the currently still relatively large uncertainties that we conservatively assign to the surface densities of individual objects. The latter can be achieved in particular by using kinematic observations to directly constrain the average halo mass inside a given radius, rather than fitting the data to a pre-selected profile and then reconstruct the mass. For a velocity-independent cross-section, our current result is formally already somewhat smaller than the range $0.5-5$ cm$^2$/g that has been invoked to explain potential inconsistencies between small-scale observations and expectations in the standard collisionless cold dark matter paradigm.

We use BOSS DR9 quasars to constrain 2 cases of dark matter models: cold-plus-warm (C+WDM) where the warm component is a thermal relic, and sterile neutrinos resonantly produced in the presence of a lepton asymmetry (RPSN). We establish constraints on the relic mass m_x and its relative abundance $F=\Omega_{wdm}/\Omega_{dm}$ using a suite of hydrodynamical simulations in 28 C+WDM configurations. We find that the 3 sigma bounds approximately follow F ~ $0.35 (keV/m_x)^{-1.37}$ from BOSS alone. We also establish constraints on sterile neutrino mass and mixing angle by producing the non-linear flux power spectrum of 8 RPSN models, where the input linear power spectrum is computed directly from the particles distribution functions. We find values of lepton asymmetries for which sterile neutrinos as light as 6.5 keV (resp. 3.5 keV) are consistent with BOSS at the 2 sigma (resp. 3sigma) level. These limits tighten by close to a factor of 2 for lepton asymmetries departing from those yielding the coolest distribution functions. Our Ly-a forest bounds can be strengthened if we include higher-resolution data from XQ-100, HIRES and MIKE. At these smaller scales, the flux power spectrum exhibits a suppression that can be due to Doppler broadening, IGM pressure smoothing or free-streaming of WDM particles. In the current work, we show that if one extrapolates temperatures from lower redshifts via broken power laws in T_0 and gamma, then our 3 sigma C+WDM bounds strengthen to F~ $0.20 (keV/m_x)^{-1.37}$, and the lightest RPSN consistent with our extended data set have masses of 7.0 keV at the 3 sigma level. Using dedicated hydrodynamical simulations, we show that a 7 keV sterile neutrino produced in a lepton asymmetry $L = 8 \times 10^{-6}$ is consistent at 1.9 sigma (resp. 3.1 sigma) with BOSS (resp. BOSS + higher-resolution), for the thermal history models tested in this work.

The magnetohydrodynamic (MHD) description of plasmas with relativistic particles necessarily includes an additional new field, the chiral chemical potential associated with the axial charge (i.e., the number difference between right- and left-handed relativistic fermions). This chiral chemical potential gives rise to a contribution to the electric current density of the plasma (\emph{chiral magnetic effect}). We present a self-consistent treatment of the \emph{chiral MHD equations}, which include the back-reaction of the magnetic field on a chiral chemical potential and its interaction with the plasma velocity field. A number of novel phenomena are exhibited. First, we show that the chiral magnetic effect decreases the frequency of the Alfv\'{e}n wave for incompressible flows, increases the frequencies of the Alfv\'{e}n wave and of the fast magnetosonic wave for compressible flows, and decreases the frequency of the slow magnetosonic wave. Second, we show that, in addition to the well-known laminar chiral dynamo effect, which is not related to fluid motions, there is a dynamo caused by the joint action of velocity shear and chiral magnetic effect. In the presence of turbulence with vanishing mean kinetic helicity, the derived mean-field chiral MHD equations describe turbulent large-scale dynamos caused by the chiral alpha effect, which is dominant for large fluid and magnetic Reynolds numbers. The chiral alpha effect is due to an interaction of the chiral magnetic effect and fluctuations of the small-scale current produced by tangling magnetic fluctuations (which are generated by tangling of the large-scale magnetic field by sheared velocity fluctuations). These dynamo effects may have interesting consequences in the dynamics of the early universe, neutron stars, and the quark--gluon plasma.

We re-analyse high redshift and high resolution Lyman-$\alpha$ forest spectra considered in (Viel et al 2013), seeking to constrain the properties of warm dark matter particles. Compared to this previous work, we consider a wider range of thermal histories of the intergalactic medium. We find that both warm and cold dark matter models can explain the cut-off observed in the flux power spectra of high-resolution observations equally well. This implies, however, very different thermal histories and underlying re-ionisation models. We discuss how to remove this degeneracy.

We study galaxy formation in sterile neutrino dark matter models that differ signifi- cantly from both cold and from `warm thermal relic' models. We use the EAGLE code to carry out hydrodynamic simulations of the evolution of pairs of galaxies chosen to resemble the Local Group, as part of the APOSTLE simulations project. We compare cold dark matter (CDM) with two sterile neutrino models with 7 keV mass: one, the warmest among all models of this mass (LA120) and the other, a relatively cold case (LA10). We show that the lower concentration of sterile neutrino subhalos compared to their CDM counterparts makes the inferred inner dark matter content of galaxies like Fornax (or Magellanic Clouds) less of an outlier in the sterile neutrino cosmologies. In terms of the galaxy number counts the LA10 simulations are indistinguishable from CDM when one takes into account halo-to-halo (or `simulation-to-simulation') scatter. In order for the LA120 model to match the number of Local Group dwarf galaxies, a higher fraction of low mass haloes is required to form galaxies than is predicted by the EAGLE simulations. As the census of the Local Group galaxies nears completion, this population may provide a strong discriminant between cold and warm dark matter models.

N-body dark matter simulations of structure formation in the $\Lambda$CDM model predict a population of subhalos within Galactic halos that have higher central densities than inferred for satellites of the Milky Way, a tension known as the `too big to fail' problem. Proposed solutions include baryonic effects, a smaller mass for the Milky Way halo, and warm dark matter. We test these three possibilities using a semi-analytic model of galaxy formation to generate luminosity functions for Milky Way halo-analogue satellite populations, the results of which are then coupled to the Jiang & van den Bosch model of subhalo stripping to predict the subhalo $V_\mathrm{max}$ functions for the 10 brightest satellites. We find that selecting the brightest satellites (as opposed to the most massive) and modelling the expulsion of gas by supernovae at early times increases the likelihood of generating the observed Milky Way satellite $V_\mathrm{max}$ function. The preferred halo mass is $6\times10^{11}M_{\odot}$, which has a 14 percent probability to host a $V_\mathrm{max}$ function like that of the Milky Way satellites. This probability is reduced to 8 percent for a $1.0\times10^{12}M_{\odot}$ halo and to 3 percent for a $1.4\times10^{12}M_{\odot}$ halo. We conclude that the Milky Way satellite $V_\mathrm{max}$ function is compatible with a CDM cosmology, as previously found by Sawala et al. using hydrodynamic simulations. Sterile neutrino-warm dark matter models achieve a higher degree of agreement with the observations, with a maximum 35 percent chance of generating the observed Milky Way satellite $V_\mathrm{max}$ function. However, more work is required to check that the semi-analytic stripping model is calibrated correctly in the sterile neutrino cosmology, and to check if our sterile neutrino models produce sufficient numbers of faint satellites.

The presence of asymmetry between fermions of opposite handedness in plasmas of relativistic particles can lead to exponential growth of a helical magnetic field via a small-scale chiral dynamo instability known as the chiral magnetic effect. Here, we show, using dimensional arguments and numerical simulations, that this process produces through the Lorentz force chiral magnetically driven turbulence. A k^{-2} magnetic energy spectrum emerges via inverse transfer over a certain range of wavenumbers k. The total chirality (magnetic helicity plus normalized chiral chemical potential) is conserved in this system. Therefore, as the helical magnetic field grows, most of the total chirality gets transferred into magnetic helicity until the chiral magnetic effect terminates. Quantitative results for height, slope, and extent of the spectrum are obtained. Consequences of this effect for cosmic magnetic fields are discussed.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

The sterile neutrino is a viable dark matter candidate that can be produced in the early Universe via non-equilibrium processes, and would therefore possess a highly non-thermal spectrum of primordial velocities. In this paper we analyse the process of structure formation with this class of dark matter particles. To this end we construct primordial dark matter power spectra as a function of the lepton asymmetry, $L_6$, that is present in the primordial plasma and leads to resonant sterile neutrino production. We compare these power spectra with those of thermally produced dark matter particles and show that resonantly produced sterile neutrinos are much colder than their thermal relic counterparts. We also demonstrate that the shape of these power spectra is not determined by the free-streaming scale alone. We then use the power spectra as an input for semi-analytic models of galaxy formation in order to predict the number of luminous satellite galaxies in a Milky Way-like halo. By assuming that the mass of the Milky Way halo must be no more than $2\times10^{12}M_{\odot}$ (the adopted upper bound based on current astronomical observations) we are able to constrain the value of $L_6$ for $M_s\le 8$~keV. We also show that the range of $L_6$ that is in best agreement with the 3.5~keV line (if produced by decays of 7~keV sterile neutrino) requires that the Milky Way halo has a mass no smaller than $1.5\times10^{12}M_{\odot}$. Finally, we compare the power spectra obtained by direct integration of the Boltzmann equations for a non-resonantly produced sterile neutrino with the fitting formula of Viel~et~al. and find that the latter significantly underestimates the power amplitude on scales relevant to satellite galaxies.

The recent discovery of the unidentified emission line at 3.55 keV in galaxies and clusters has attracted great interest from the community. As the origin of the line remains uncertain, we study the surface brightness distribution of the line in the Perseus cluster since that information can be used to identify its origin. We examine the flux distribution of the 3.55 keV line in the deep Suzaku observations of the Perseus cluster in detail. The 3.55 keV line is observed in three concentric annuli in the central observations, although the observations of the outskirts of the cluster did not reveal such a signal. We establish that these detections and the upper limits from the non-detections are consistent with a dark matter decay origin. However, absence of positive detection in the outskirts is also consistent with some unknown astrophysical origin of the line in the dense gas of the Perseus core, as well as with a dark matter origin with a steeper dependence on mass than the dark matter decay. We also comment on several recently published analyses of the 3.55 keV line.

We present results of a search for the 3.5 keV emission line in our recent very long (~ 1.4 Ms) XMM-Newton observation of the Draco dwarf spheroidal galaxy. The astrophysical X-ray emission from such dark matter-dominated galaxies is faint, thus they provide a test for the dark matter origin of the 3.5 keV line previously detected in other massive, but X-ray bright objects, such as galaxies and galaxy clusters. We do not detect a statistically significant emission line from Draco; this constrains the lifetime of a decaying dark matter particle to tau > (7-9) x 10^27 s at 95% CL (combining all three XMM-Newton cameras; the interval corresponds to the uncertainty of the dark matter column density in the direction of Draco). The PN camera, which has the highest sensitivity of the three, does show a positive spectral residual (above the carefully modeled continuum) at E = 3.54 +/- 0.06 keV with a 2.3 sigma significance. The two MOS cameras show less-significant or no positive deviations, consistently within 1 sigma with PN. Our Draco limit on tau is consistent with previous detections in the stacked galaxy clusters, M31 and the Galactic Center within their 1-2 sigma uncertainties, but is inconsistent with the high signal from the core of the Perseus cluster (which has itself been inconsistent with the rest of the detections). We conclude that this Draco observation does not exclude the dark matter interpretation of the 3.5 keV line in those objects.

Recent studies of M31, the Galactic centre, and galaxy clusters have made tentative detections of an X-ray line at ~3.5 keV that could be produced by decaying dark matter. We use high resolution simulations of the Aquarius project to predict the likely amplitude of the X-ray decay flux observed in the GC relative to that observed in M31, and also of the GC relative to other parts of the Milky Way halo and to dwarf spheroidal galaxies. We show that the reported detections from M31 and Andromeda are compatible with each other, and with upper limits arising from high galactic latitude observations, and imply a decay time {\tau} ~10^28 seconds. We argue that this interpretation can be tested with deep observations of dwarf spheroidal galaxies: in 95 per cent of our mock observations, a 1.3 Msec pointed observation of Draco with XMM-Newton will enable us to discover or rule out at the 3{\sigma} level an X-ray feature from dark matter decay at 3.5 keV, for decay times {\tau} < 0.8x10^28 sec.

We report a weak line at 3.52+/-0.02 keV in X-ray spectra of M31 galaxy and the Perseus galaxy cluster observed by MOS and PN cameras of XMM-Newton telescope. This line is not known as an atomic line in the spectra of galaxies or clusters. It becomes stronger towards the centers of the objects; is stronger for Perseus than for M31; is absent in the spectrum of a deep "blank sky'' dataset. Although for each object it is hard to exclude that the feature is due to an instrumental effect or an atomic line, it is consistent with the behavior of a dark matter decay line. Future (non-)detections of this line in multiple objects may help to reveal its nature.

We detect a line at 3.539 +/- 0.011 keV in the deep exposure dataset of the Galactic Center region, observed with the XMM-Newton. The dark matter interpretation of the signal observed in the Perseus galaxy cluster, the Andromeda galaxy [1402.4119] and in the stacked spectra of galaxy clusters [1402.2301], together with non-observation of the line in blank sky data, put both lower and upper limits on the possible intensity of the line in the Galactic Center data. Our result is consistent with these constraints for a class of Milky Way mass models, presented previously by observers, and would correspond to radiative decay dark matter lifetime tau_dm ~ (6-8) x 10^{27} sec. Although it is hard to exclude an astrophysical origin of this line based the Galactic Center data alone, this is an important consistency check of the hypothesis that encourages to check it with more observational data that are expected by the end of 2015.

We study the dynamics of a plasma of charged relativistic fermions at very high temperature $T\gg m$, where $m$ is the fermion mass, coupled to the electromagnetic field. In particular, we derive a magneto-hydrodynamical description of the evolution of such a plasma. We show that, as compared to conventional MHD for a plasma of non-relativistic particles, the hydrodynamical description of the relativistic plasma involves new degrees of freedom described by a pseudo-scalar field originating in a local asymmetry in the densities of left-handed and right-handed fermions. This field can be interpreted as an effective axion field. Taking into account the chiral anomaly we present dynamical equations for the evolution of this field, as well as of other fields appearing in the MHD description of the plasma. Due to its non-linear coupling to helical magnetic fields, the axion field significantly affects the dynamics of a magnetized plasma and can give rise to a novel type of inverse cascade.

We revisit the X-ray spectrum of the central 14' of the Andromeda galaxy, discussed in our previous work [1402.4119]. Recently in [1408.1699] it was claimed that if one limits the analysis of the data to the interval 3-4 keV, the significance of the detection of the line at 3.53 keV drops below 2 sigma. In this note we show that such a restriction is not justified, as the continuum is well-modeled as a power law up to 8 keV, and parameters of the background model are well constrained over this larger interval of energies. This allows for a detection of the line at 3.53 keV with a statistical significance greater than ~3 sigma and for the identification of several known atomic lines in the energy range 3-4 keV. Limiting the analysis to the 3-4 keV interval results in increased uncertainty, thus decreasing the significance of the detection. We also argue that, with the M31 data included, a consistent interpretation of the 3.53 keV line as an atomic line of K XVIII in all studied objects is problematic.

The diversity of structures in the Universe (from the smallest galaxies to the largest superclusters) has formed under the pull of gravity from the tiny primordial perturbations that we see imprinted in the cosmic microwave background. A quantitative description of this process would require description of motion of zillions of dark matter particles. This impossible task is usually circumvented by coarse-graining the problem: one either considers a Newtonian dynamics of "particles" with macroscopically large masses or approximates the dark matter distribution with a continuous density field. There is no closed system of equations for the evolution of the matter density field alone and instead it should still be discretized at each timestep. In this work we describe a method of solving the full 6-dimensional Vlasov-Poisson equation via a system of auxiliary Schroedinger-like equations. The complexity of the problem gets shifted into the choice of the number and shape of the initial wavefunctions that should only be specified at the beginning of the computation (we stress that these wavefunctions have nothing to do with quantum nature of the actual dark matter particles). We discuss different prescriptions to generate the initial wave functions from the initial conditions and demonstrate the validity of the technique on two simple test cases. This new simulation algorithm can in principle be used on an arbitrary distribution function, enabling the simulation of warm and hot dark matter structure formation scenarios.

Large Observatory For X-ray Timing (LOFT) is a next generation X-ray telescope selected by European Space Agency as one of the space mission concepts within the ``Cosmic Vision'' programme. The Large Area Detector on board of LOFT will be a collimator-type telescope with an unprecedentedly large collecting area of about 10 square meters in the energy band between 2 and 100 keV. We demonstrate that LOFT will be a powerful dark matter detector, suitable for the search of the X-ray line emission expected from decays of light dark matter particles in galactic halos. We show that LOFT will have sensitivity for dark matter line search more than an order of magnitude higher than that of all existing X-ray telescopes. In this way, LOFT will be able to provide a new insight into the fundamental problem of the nature of dark matter.

We review the status of sterile neutrino dark matter and discuss astrophysical and cosmological bounds on its properties as well as future prospects for its experimental searches. We argue that if sterile neutrinos are the dominant fraction of dark matter, detecting an astrophysical signal from their decay (the so-called 'indirect detection') may be the only way to identify these particles experimentally. However, it may be possible to check the dark matter origin of the observed signal unambiguously using its characteristic properties and/or using synergy with accelerator experiments, searching for other sterile neutrinos, responsible for neutrino flavor oscillations. We argue that to fully explore this possibility a dedicated cosmic mission - an X-ray spectrometer - is needed.

Cold dark matter models predict the existence of a large number of substructures within dark matter halos. If the cold dark matter consists of weakly interacting massive particles, their annihilation within these substructures could lead to diffuse GeV emission that would dominate over the annihilation signal of the host halo. In this work we search for GeV emission from three nearby galaxy clusters: Coma, Virgo and Fornax. We first remove known extragalactic and galactic diffuse gamma-ray backgrounds and point sources from the Fermi 2-year catalog and find a significant residual diffuse emission in all three clusters. We then investigate whether this emission is due to (i) unresolved point sources; (ii) dark matter annihilation; or (iii) cosmic rays (CR). Using 45 months of Fermi-LAT data we detect several new point sources (not present in the Fermi 2-year point source catalogue) which contaminate the signal previously analyzed by Han et al.(arxiv:1201.1003). Including these and accounting for the effects of undetected point sources, we find no significant detection of extended emission from the three clusters studied. Instead, we determine upper limits on emission due to dark matter annihilation and cosmic rays. For Fornax and Virgo the limits on CR emission are consistent with theoretical models, but for Coma the upper limit is a factor of 2 below the theoretical expectation. Allowing for systematic uncertainties associated with the treatment of CR, the upper limits on the cross section for dark matter annihilation from our clusters are more stringent than those from analyses of dwarf galaxies in the Milky Way. We rule out the thermal cross section for supersymmetric dark matter particles for masses as large as 100 GeV (depending on the annihilation channel).

We study the properties of the diffuse gamma-ray background around the Galactic plane at energies 20 -- 200 GeV. We find that the spectrum of this emission possesses significant spacial variations with respect to the average smooth component. The positions and shapes of these spectral features change with the direction on the sky. We therefore argue, that the spectral feature around 130 GeV, found in several regions around the Galactic Center and in the Galactic plane in [1203.1312, 1204.2797, 1205.1045, 1206.1616], can not be interpreted with confidence as a gamma-ray line, but may be a component of the diffuse background and can be of instrumental or astrophysical origin. Therefore, the dark matter origin of this spectral feature becomes dubious.

We analyze the properties of dark matter halos in the cold-plus-warm dark matter cosmologies (CWDM). We study their dependence on the fraction and velocity dispersion of the warm particle, keeping the free-streaming scale fixed. To this end we consider three models with the same free-streaming: (1) a mixture of 90% of CDM and 10% of WDM with the mass 1 keV; (2) a mixture of 50% of CDM and 50% of WDM with the mass 5 keV; and (3) pure WDM with the mass 10 keV. "Warm" particles have rescaled Fermi-Dirac spectrum of primordial velocities (as non-resonantly produced sterile neutrinos would have). We compare the properties of halos among these models and with a LCDM with the same cosmological parameters. We demonstrate, that although these models have the same free-streaming length and the suppression of matter spectra are similar at scales probed by the Lyman-alpha forest comoving wave-numbers k<3-5 h/Mpc), the resulting properties of halos with masses below 1e11 Msun are different due to the different behaviour of matter power spectra at smaller scales. In particular, we find that while the number of galaxies remains the same as in LCDM case, their density profiles become much less concentrated, and hence in better agreement with current observational constraints. Our results imply that a single parameter (e.g. free streaming length) description of these models is not enough to fully capture their effects on the structure formation process.

This white paper addresses the hypothesis of light sterile neutrinos based on recent anomalies observed in neutrino experiments and the latest astrophysical data.

In thermal equilibrium the ground state of the plasma of Standard Model particles is determined by temperature and exactly conserved combinations of baryon and lepton numbers. We show that at non-zero values of the global charges a translation invariant and homogeneous state of the plasma becomes unstable and the system transits into a new state, containing a large-scale magnetic field. The origin of this effect is the parity-breaking character of weak interactions and chiral anomaly. This situation can occur in the early Universe and may play an important role in its subsequent evolution.

High resolution N-body simulations of galactic cold dark matter haloes indicate that we should expect to find a few satellite galaxies around the Milky Way whose haloes have a maximum circular velocity in excess of 40 kms. Yet, with the exception of the Magellanic Clouds and the Sagittarius dwarf, which likely reside in subhaloes with significantly larger velocities than this, the bright satellites of the Milky Way all appear to reside in subhaloes with maximum circular velocities below 40 kms. As recently highlighted by Boylan-Kolchin et al., this discrepancy implies that the majority of the most massive subhaloes within a cold dark matter galactic halo are much too concentrated to be consistent with the kinematic data for the bright Milky Way satellites. Here we show that no such discrepancy exists if haloes are made of warm, rather than cold dark matter because these haloes are less concentrated on account of their typically later formation epochs. Warm dark matter is one of several possible explanations for the observed kinematics of the satellites.

We show that the evolution of magnetic fields in a primordial plasma, filled with Standard Model particles, at temperatures T > 10 MeV is strongly affected by the quantum chiral anomaly -- an effect that has been neglected previously. Although reactions equilibrating left and right-chiral electrons are in deep thermal equilibrium for T < 80 TeV, an asymmetry between these particle develops in the presence of strong magnetic fields. This results in magnetic helicity transfer from shorter to longer scales. This also leads to an effective generation of lepton asymmetry that may survive in the plasma down to temperatures T ~ 10 MeV, which may strongly affect many processes in the early Universe. Although we report our results for the Standard Model, they are likely to play an important role also in its extensions.

ORIGIN is a proposal for the M3 mission call of ESA aimed at the study of metal creation from the epoch of cosmic dawn. Using high-spectral resolution in the soft X-ray band, ORIGIN will be able to identify the physical conditions of all abundant elements between C and Ni to red-shifts of z=10, and beyond. The mission will answer questions such as: When were the first metals created? How does the cosmic metal content evolve? Where do most of the metals reside in the Universe? What is the role of metals in structure formation and evolution? To reach out to the early Universe ORIGIN will use Gamma-Ray Bursts (GRBs) to study their local environments in their host galaxies. This requires the capability to slew the satellite in less than a minute to the GRB location. By studying the chemical composition and properties of clusters of galaxies we can extend the range of exploration to lower redshifts (z ~ 0.2). For this task we need a high-resolution spectral imaging instrument with a large field of view. Using the same instrument, we can also study the so far only partially detected baryons in the Warm-Hot Intergalactic Medium (WHIM). The less dense part of the WHIM will be studied using absorption lines at low redshift in the spectra for GRBs.

In the recent paper of Hooper & Goodenough (1010.2752) it was reported that gamma-ray emission from the Galactic Center region contains an excess compared to the contributions from the large-scale diffuse emission and known point sources. This excess was argued to be consistent with a signal from annihilation of Dark Matter with a power law density profile. We reanalyze the Fermi data and find instead that it is consistent with the "standard model" of diffuse emission and of known point sources. The main reason for the discrepancy with the interpretation of 1010.2752 is different (as compared to the previous works) spectrum of the point source at the Galactic Center assumed in 1010.2752. We discuss possible reasons for such an interpretation.

We discuss the universal relation between density and size of observed Dark Matter halos that was recently shown to hold on a wide range of scales, from dwarf galaxies to galaxy clusters. Predictions of LambdaCDM N-body simulations are consistent with this relation. We demonstrate that this property of LambdaCDM can be understood analytically in the secondary infall model. Qualitative understanding given by this model provides a new way to predict which deviations from LambdaCDM or large-scale modifications of gravity can affect universal behavior and, therefore, to constrain them observationally.

A signal from decaying dark matter (DM) can be unambiguously distinguished from spectral features of astrophysical or instrumental origin by studying its spatial distribution. We demonstrate this approach by examining the recent claim of 0912.0552 regarding the possible DM origin of the 2.5 keV line in Chandra observations of the Milky Way satellite known as Willman 1. Our conservative strategy is to adopt a relatively large dark mass for Willman 1 and relatively small dark masses for the comparison objects. We analyze archival observations by XMM-Newton of M31 and Fornax dwarf spheroidal galaxy (dSph) and Chandra observations of Sculptor dSph. By performing a conservative analysis of X-ray spectra, we show the absence of a DM decay line with parameters consistent with those of 0912.0552. For M31, the observations of the regions between 10 and 20 kpc from the center, where the uncertainties in the DM distribution are minimal, make a strong exclusion at the level above 10sigma. The minimal estimate for the amount of DM in the central 40 kpc of M31 is provided by the model of 0912.4133, assuming the stellar disk's mass to light ratio ~8 and almost constant DM density within a core of 28 kpc. Even in this case one gets an exclusion at 5.7sigma from central region of M31 whereas modeling all processed data from M31 and Fornax produces more than 14sigma exclusion. Therefore, despite possible systematic uncertainties, we exclude the possibility that the spectral feature at ~2.5 keV found in 0912.0552 is a DM decay line. We conclude, however, that the search for DM decay line, although demanding prolonged observations of well-studied dSphs, M31 outskirts and other similar objects, is rather promising, as the nature of a possible signal can be checked. An (expected) non-observation of a DM decay signal in the planned observations of Willman 1 should not discourage further dedicated observations.

We suggest a new efficient way to constrain a certain class of large scale modifications of gravity. We show that the scale-free relation between density and size of Dark Matter halos, predicted within the LambdaCDM model with Newtonian gravity, gets modified in a wide class of theories of modified gravity.

We present a new universal relation, satisfied by matter distributions at all observed scales, and show its amazingly good and detailed agreement with the predictions of the most up-to-date pure dark matter simulations of structure formation in the Universe. This work extends the previous analysis [0904.4054; 0909.5203] to a larger range of masses, demonstrates a different scaling law, and compares it with numerical simulations. This behaviour seems to be insensitive to the complicated feedback of baryons on dark matter. Therefore, it potentially allows to compare theoretical predictions directly with observations, thus providing a new tool to constrain the properties of dark matter. Such a universal property, observed in structures of all sizes (from dwarf spheroidal galaxies to galaxy clusters), is difficult to explain without dark matter, thus providing new evidence for its existence.

We present a comprehensive overview of an extension of the Standard Model that contains three right-handed (sterile) neutrinos with masses below the electroweak scale [the Neutrino Minimal Standard Model, (nuMSM)]. We consider the history of the Universe from the inflationary era through today and demonstrate that most of the observed phenomena beyond the Standard Model can be explained within the framework of this model. We review the mechanism of baryon asymmetry of the Universe in the nuMSM and discuss a dark matter candidate that can be warm or cold and satisfies all existing constraints. From the viewpoint of particle physics the model provides an explanation for neutrino flavor oscillations. Verification of the nuMSM is possible with existing experimental techniques.

We propose an X-ray mission called Xenia to search for decaying superweakly interacting Dark Matter particles (super-WIMP) with a mass in the keV range. The mission and its observation plan are capable of providing a major break through in our understanding of the nature of Dark Matter (DM). It will confirm, or reject, predictions of a number of particle physics models by increasing the sensitivity of the search for decaying DM by about two orders of magnitude through a wide-field imaging X-ray spectrometer in combination with a dedicated observation program. The proposed mission will provide unique limits on the mixing angle and mass of neutral leptons, right handed partners of neutrinos, which are important Dark Matter candidates. The existence of these particles is strongly motivated by observed neutrino flavor oscillations and the problem of baryon asymmetry of the Universe. In super-WIMP models, the details of the formation of the cosmic web are different from those of LambdaCDM. The proposed mission will, in addition to the search for decaying Dark Matter, provide crucial insight into the nature of DM by studying the structure of the "cosmic web". This will be done by searching for missing baryons in emission, and by using gamma-ray bursts as backlight to observe the warm-hot intergalactic media in absorption.

We revisit Lyman-alpha bounds on the dark matter mass in Lambda Warm Dark Matter (Lambda-WDM) models, and derive new bounds in the case of mixed Cold plus Warm models (Lambda-CWDM), using a set up which is a good approximation for several theoretically well-motivated dark matter models. We combine WMAP5 results with two different Lyman-alpha data sets, including observations from the Sloan Digital Sky Survey. We pay a special attention to systematics, test various possible sources of error, and compare the results of different statistical approaches. Expressed in terms of the mass of a non-resonantly produced sterile neutrino, our bounds read m_NRP > 8 keV (frequentist 99.7% confidence limit) or m_NRP > 12.1 keV (Bayesian 95% credible interval) in the pure Lambda-WDM limit. For the mixed model, we obtain limits on the mass as a function of the warm dark matter fraction F_WDM. Within the mass range studied here (5 keV < m_NRP < infinity), we find that any mass value is allowed when F_WDM < 0.6 (frequentist 99.7% confidence limit); similarly, the Bayesian joint probability on (F_WDM, 1/m_NRP) allows any value of the mass at the 95% confidence level, provided that F_WDM < 0.35.

We discuss the bounds on the mass of Dark Matter (DM) particles, coming from the analysis of DM phase-space distribution in dwarf spheroidal galaxies (dSphs). After reviewing the existing approaches, we choose two methods to derive such a bound. The first one depends on the information about the current phase space distribution of DM particles only, while the second one uses both the initial and final distributions. We discuss the recent data on dSphs as well as astronomical uncertainties in relevant parameters. As an application, we present lower bounds on the mass of DM particles, coming from various dSphs, using both methods. The model-independent bound holds for any type of fermionic DM. Stronger, model-dependent bounds are quoted for several DM models (thermal relics, non-resonantly and resonantly produced sterile neutrinos, etc.). The latter bounds rely on the assumption that baryonic feedback cannot significantly increase the maximum of a distribution function of DM particles. For the scenario in which all the DM is made of sterile neutrinos produced via non-resonant mixing with the active neutrinos (NRP) this gives m_nrp > 1.7 keV. Combining these results in their most conservative form with the X-ray bounds of DM decay lines, we conclude that the NRP scenario remains allowed in a very narrow parameter window only. This conclusion is independent of the results of the Lyman-alpha analysis. The DM model in which sterile neutrinos are resonantly produced in the presence of lepton asymmetry remains viable. Within the minimal neutrino extension of the Standard Model (the nuMSM), both mass and the mixing angle of the DM sterile neutrino are bounded from above and below, which suggests the possibility for its experimental search.

Using the high-resolution spectrometer SPI on board the International Gamma-Ray Astrophysics Laboratory (INTEGRAL), we search for a spectral line produced by a dark matter(DM) particle with a mass in the range 40keV < M_DM < 14MeV, decaying in the DM halo of the Milky Way. To distinguish the DM decay line from numerous instrumental lines found in the SPI background spectrum, we study the dependence of the intensity of the line signal on the offset of the SPI pointing from the direction toward the Galactic Centre. After a critical analysis of the uncertainties of the DM density profile in the inner Galaxy, we find that the intensity of the DM decay line should decrease by at least a factor of 3 when the offset from the Galactic Centre increases from 0 to 180 degrees. We find that such a pronounced variation of the line flux across the sky is not observed for any line, detected with a significance higher than 3 sigma in the SPI background spectrum. Possible DM decay origin is not ruled out only for the unidentified spectral lines, having low (~3 sigma) significance or coinciding in position with the instrumental ones. In the energy interval from 20 keV to 7 MeV, we derive restrictions on the DM decay line flux, implied by the (non-)detection of the DM decay line. For a particular DM candidate, the sterile neutrino of mass MDM, we derive a bound on the mixing angle.

Previous fits of sterile neutrino dark matter models to cosmological data assumed a peculiar production mechanism, which is not representative of the best-motivated particle physics models given current data on neutrino oscillations. These analyses ruled out sterile neutrino masses smaller than 8-10 keV. Here we focus on sterile neutrinos produced resonantly. We show that their cosmological signature can be approximated by that of mixed Cold plus Warm Dark Matter (CWDM). We use recent results on LambdaCWDM models to show that for each mass greater than or equal to 2 keV, there exists at least one model of sterile neutrino accounting for the totality of dark matter, and consistent with Lyman-alpha and other cosmological data. Resonant production occurs in the framework of the nuMSM (the extension of the Standard Model with three right-handed neutrinos). The models we checked to be allowed correspond to parameter values consistent with neutrino oscillation data, baryogenesis and all other dark matter bounds.

In this talk we review existing cosmological and astrophysical bounds on light (with the mass in keV - MeV range) and super-weakly interacting dark matter candidates. A particular attention is paid to the sterile neutrino DM candidate.

We derive constraints on parameters of the radiatively decaying Dark Matter (DM) particles, using XMM-Newton EPIC spectra of the Andromeda galaxy (M31). Using the observations of the outer (5'-13') parts of M31 we improve the existing constraints. For the case of sterile neutrino DM, combining our constraints with the latest computation of abundances of sterile neutrino in the Dodelson-Widrow (DW) scenario, we obtain the lower mass limit m_s < 4 keV, which is stronger than the previous one m_s < 6 kev, obtained recently by Asaka et al. (2007) [hep-ph/0612182]. Comparing this limit with the most recent results on Lyman-alpha forest analysis of Viel et al. (2007) [arXiv:0709.0131] (m_s > 5.6 kev), we argue that the scenario in which all the DM is produced via DW mechanism is ruled out. We discuss however other production mechanisms and note that the sterile neutrino remains a viable candidate of Dark Matter, either warm or cold.

We improve the earlier restrictions on parameters of the dark matter (DM) in the form of a sterile neutrino. The results were obtained from non-observing the DM decay line in the X-ray spectrum of the Milky Way (using the recent XMM-Newton PN blank sky data). We also present a similar constraint coming from the recent XMM-Newton observation of Ursa Minor -- dark, X-ray quiet dwarf spheroidal galaxy. The new Milky way data improve on (by as much as the order of magnitude at masses ~3.5 keV) existing constraints. Although the observation of Ursa Minor has relatively poor statistics, the constraints are comparable to those recently obtained using observations of the Large Magellanic Cloud or M31. This confirms a recent proposal that dwarf satellites of the MW are very interesting candidates for the DM search and dedicated studies should be made to this purpose.

Sterile neutrinos with the mass in the keV range are interesting warm dark matter (WDM) candidates. The restrictions on their parameters (mass and mixing angle) obtained by current X-ray missions (XMM-Newton or Chandra) can only be improved by less than an order of magnitude in the near future. Therefore the new strategy of search is needed. We compare the sensitivities of existing and planned X-ray missions for the detection of WDM particles with the mass ~1-20 keV. We show that existing technology allows an improvement in sensitivity by a factor of 100. Namely, two different designs can achieve such an improvement: [A] a spectrometer with the high spectral resolving power of 0.1%, wide (steradian) field of view, with small effective area of about cm^2 (which can be achieved without focusing optics) or [B] the same type of spectrometer with a smaller (degree) field of view but with a much larger effective area of 10^3 cm^2 (achieved with the help of focusing optics). To illustrate the use of the "type A" design we present the bounds on parameters of the sterile neutrino obtained from analysis of the data taken by an X-ray microcalorimeter. In spite of the very short exposure time (100 sec) the derived bound is comparable to the one found from long XMM-Newton observation.

We derived constraints on parameters of a radiatively decaying warm dark matter particle, e.g., the mass and mixing angle for a sterile neutrino, using Chandra X-ray spectra of a galaxy cluster 1E0657-56 (the ``bullet'' cluster). The constraints are based on nondetection of the sterile neutrino decay emission line. This cluster exhibits spatial separation between the hot intergalactic gas and the dark matter, helping to disentangle their X-ray signals. It also has a very long X-ray observation and a total mass measured via gravitational lensing. This makes the resulting constraints on sterile neutrino complementary to earlier results that used different cluster mass estimates. Our limits are comparable to the best existing constraints.

We propose a strategy of how to look for dark matter (DM) particles possessing a radiative decay channel and derive constraints on their parameters from observations of X-rays from our own Galaxy and its dwarf satellites. When applied to the sterile neutrinos in keV mass range, it allows a significant improvement of restrictions to its parameters, as compared with previous works.

Sterile neutrinos with masses in the keV range are viable candidates for the warm dark matter. We analyze existing data for the extragalactic diffuse X-ray background for signatures of sterile neutrino decay. The absence of detectable signal within current uncertainties of background measurements puts model-independent constraints on allowed values of sterile neutrino mass and mixing angle, which we present in this work.

We find restrictions on the mass and mixing angle of the dark matter sterile neutrinos using X-ray observations of Coma and Virgo galaxy clusters.

In an extention of the Standard Model by three relatively light right-handed neutrinos (the nuMSM model) the role of the dark matter particle is played by the lightest sterile neutrino. We demonstrate that the observations of the extragalactic X-ray background allow to put a strong upper bound on the mass of the lightest active neutrino and predict the absolute values of the mass of the two heavier active neutrinos in the nuMSM, provided that the mass of the dark matter sterile neutrino is larger than 1.8 keV.