The mass and anisotropy profiles of galaxy clusters from the projected phase space density: testing the method on simulated data

We present a new method of constraining the mass and velocity anisotropy profiles of galaxy clusters from kinematic data. The method is based on a model of the phase space density which allows the anisotropy to vary with radius between two asymptotic values. The characteristic scale of transition between these asymptotes is fixed and tuned to a typical anisotropy profile resulting from cosmological simulations. The model is parametrized by two values of anisotropy, at the centre of the cluster and at infinity, and two parameters of the NFW density profile, the scale radius and the scale mass. In order to test the performance of the method in reconstructing the true cluster parameters we analyze mock kinematic data for 20 relaxed galaxy clusters generated from a cosmological simulation of the standard LCDM model. We use Bayesian methods of inference and the analysis is carried out following the Markov Chain Monte Carlo approach. The parameters of the mass profile are reproduced quite well, but we note that the mass is typically underestimated by 15 percent, probably due to the presence of small velocity substructures. The constraints on the anisotropy profile for a single cluster are in general barely conclusive. Although the central asymptotic value is determined accurately, the outer one is subject to significant systematic errors caused by substructures at large clustercentric distance. The anisotropy profile is much better constrained if one performs joint analysis of at least a few clusters. In this case it is possible to reproduce the radial variation of the anisotropy over two decades in radius inside the virial sphere.Comment: 11 pages, 10 figures, accepted for publication in MNRA

A critical reassessment of particle Dark Matter limits from dwarf satellites

Dwarf satellite galaxies are ideal laboratories for identifying particle Dark Matter signals. When setting limits on particle Dark Matter properties from null searches, it becomes however crucial the level at which the Dark Matter density profile within these systems is constrained by observations. In the limit in which the spherical Jeans equation is assumed to be valid for a given tracer stellar population, we study the solution of this equation having the Dark Matter mass profile as an output rather than as a trial parametric input. Within our new formulation, we address to what level dwarf spheroidal galaxies feature a reliable mass estimator. We assess then possible extrapolation of the density profiles in the inner regions and -- keeping explicit the dependence on the orbital anisotropy profile of the tracer population -- we derive general trends on the line-of-sight integral of the density profile squared, a quantity commonly dubbed J-factor and crucial to estimate fluxes from prompt Dark Matter pair annihilations. Taking Ursa Minor as a study case among Milky Way satellites, we perform Bayesian inference using the available kinematical data for this galaxy. Contrary to all previous studies, we avoid marginalization over quantities poorly constrained by observations or by theoretical arguments. We find minimal J-factors to be about 2 to 4 times smaller than commonly quoted estimates, approximately relaxing by the same amount the limit on Dark Matter pair annihilation cross section from gamma-ray surveys of Ursa Minor. At the same time, if one goes back to a fixed trial parametric form for the density, e.g. using a NFW or Burkert profile, we show that the minimal J can hardly be reduced by more than a factor of 1.5. \ua9 2016 IOP Publishing Ltd and Sissa Medialab srl

A lower bound on the mass of Dark Matter particles

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.Comment: 20 pages, published in JCA

Spike pattern recognition by supervised classification in low dimensional embedding space

© The Author(s) 2016. This article is published with open access at Springerlink.com under the terms of the Creative Commons Attribution License 4.0, (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.Epileptiform discharges in interictal electroencephalography (EEG) form the mainstay of epilepsy diagnosis and localization of seizure onset. Visual analysis is rater-dependent and time consuming, especially for long-term recordings, while computerized methods can provide efficiency in reviewing long EEG recordings. This paper presents a machine learning approach for automated detection of epileptiform discharges (spikes). The proposed method first detects spike patterns by calculating similarity to a coarse shape model of a spike waveform and then refines the results by identifying subtle differences between actual spikes and false detections. Pattern classification is performed using support vector machines in a low dimensional space on which the original waveforms are embedded by locality preserving projections. The automatic detection results are compared to experts’ manual annotations (101 spikes) on a whole-night sleep EEG recording. The high sensitivity (97 %) and the low false positive rate (0.1 min−1), calculated by intra-patient cross-validation, highlight the potential of the method for automated interictal EEG assessment.Peer reviewedFinal Published versio

MEASUREMENTS OF ELECTRIC DIPOLE MOMENTS OF A1Σ+A^{1}\Sigma^{+} BAO BY STARK-HANLE-EFFECT AND QUANTUM BEATS

Author Institution:Electric dipole moments of a series of vibrational levels of the excited electronic state $A^{1} \Sigma^{+} (v^{\prime} = 0, 1,2, 3)$ of BaO are measured by two different experimental techniques Using a single mode cw dye laser, rotationally resolved excitation spectra of the bandheads of the transitions $A^{1} \Sigma^{+} (v^{\prime} = 0, 1, 2, 3) - X^{1} \Sigma^{+} (v^{\prime\prime} = 0, 1)$ are obatained. By applying electric fields (0-500 V $cm^{-1}$) Stark-Hanle effect signals on the R(O) lines of these bands are observed. A pulsed $N_{2}$-laser pumped dye laser is used to measure the lifetimes of the same levels by the method of delayed coincidences. With these lifetime values the dipole moments are extracted from the halfeidts of the Hanle signals. In addition to this, dipole moment values are obtained from quantum beat modulations which are superimposed on the exponential decay when an electric field is applied. The similarities and differences of the two methods are discussed. It turns out that the quantum beat technique is more reliable since saturation influences are negligible in this case. The measured dipole moments are $\mu_{0} = 2.970, \mu_{1} = 2.660, \mu_{2} = 3.150 \mu_{3} = 3.180$