Transport of positrons in the interstellar medium

This work investigates some aspects of the transport of low-energy positrons in the interstellar medium (ISM). We consider resonance interactions with magnetohydrodynamic waves above the resonance threshold. Below the threshold, collisions take over and deflect positrons in their motion parallel to magnetic-field lines. Using Monte-Carlo simulations, we model the propagation and energy losses of positrons in the different phases of the ISM until they annihilate. We suggest that positrons produced in the disk by an old population of stars, with initial kinetic energies below 1 MeV, and propagating in the spiral magnetic field of the disk, can probably not penetrate the Galactic bulge.Comment: 4 pages, 3 figures, accepted for publication in the proceeding of the 6th INTEGRAL worksho

Lipschitz optimization methods for fitting a sum of damped sinusoids to a series of observations

A general nonlinear regression model is considered in the form of fitting a sum of damped sinusoids to a series of non-uniform observations. The problem of parameter estimation in this model is important in many applications like signal processing. The corresponding continuous optimization problem is typically difficult due to the high multiextremal character of the objective function. It is shown how Lipschitz-based deterministic methods can be well-suited for studying these challenging global optimization problems, when a limited computational budget is given and some guarantee of the found solution is required

Jackknife resampling technique on mocks: an alternative method for covariance matrix estimation

We present a fast and robust alternative method to compute covariance matrix in case of cosmology studies. Our method is based on the jackknife resampling applied on simulation mock catalogues. Using a set of 600 BOSS DR11 mock catalogues as a reference, we find that the jackknife technique gives a similar galaxy clustering covariance matrix estimate by requiring a smaller number of mocks. A comparison of convergence rates show that $\sim$7 times fewer simulations are needed to get a similar accuracy on variance. We expect this technique to be applied in any analysis where the number of available N-body simulations is low.Comment: 11 pages, 11 figures, 2 table

The all-sky distribution of 511 keV electron-positron annihilation emission

We present a map of 511 keV electron-positron annihilation emission, based on data accumulated with the SPI spectrometer aboard ESA's INTEGRAL gamma-ray observatory, that covers approximately 95% of the celestial sphere. 511 keV line emission is significantly detected towards the galactic bulge region and, at a very low level, from the galactic disk. The bulge emission is highly symmetric and is centred on the galactic centre with an extension of 8 deg. The emission is equally well described by models that represent the stellar bulge or halo populations. The disk morphology is only weakly constrained by the present data, being compatible with both the distribution of young and old stellar populations. The 511 keV line flux from the bulge and disk components is 1.05e-3 ph cm-2 s-1 and 0.7e-3 ph cm-2 s-1, respectively, corresponding to a bulge-to-disk flux ratio in the range 1-3. Assuming a positronium fraction of 0.93 this translates into annihilation rates of 1.5e43 s-1 and 3e42 s-1, respectively. The ratio of the bulge luminosity to that of the disk is in the range 3-9. We find no evidence for a point-like source in addition to the diffuse emission, down to a typical flux limit of 1e-4 ph cm-2 s-1. We also find no evidence for the positive latitude enhancement that has been reported from OSSE measurements; our 3 sigma upper flux limit for this feature is 1.5e-4 ph cm-2 s-1. The disk emission can be attributed to the beta+ decay of the radioactive species 26Al and 44Ti. The bulge emission arises from a different source which has only a weak or no disk component. We suggest that Type Ia supernovae and/or low-mass X-ray binaries are the prime candidates for the source of the galactic bulge positrons. Light dark matter annihilation could also explain the observed 511 keV bulge emission characteristics.Comment: accepted for publication in A&

HKT Geometry and Fake Five Dimensional Supergravity

Recent results on the relation between hyper-Kahler geometry with torsion and solutions admitting Killing spinors in minimal de sitter supergravity are extended to more general supergravity models with vector multiplets.Comment: 14 pages, latex. Minor typos corrected, references adde

Stochastic algorithms for solving structured low-rank approximation problems

In this paper, we investigate the complexity of the numerical construction of the Hankel structured low-rank approximation (HSLRA) problem, and develop a family of algorithms to solve this problem. Briefly, HSLRA is the problem of finding the closest (in some pre-defined norm) rank r approximation of a given Hankel matrix, which is also of Hankel structure. We demonstrate that finding optimal solutions of this problem is very hard. For example, we argue that if HSLRA is considered as a problem of estimating parameters of damped sinusoids, then the associated optimization problem is basically unsolvable. We discuss what is known as the orthogonality condition, which solutions to the HSLRA problem should satisfy, and describe how any approximation may be corrected to achieve this orthogonality. Unlike many other methods described in the literature the family of algorithms we propose has the property of guaranteed convergence

Cosmological inference including massive neutrinos from the matter power spectrum: biases induced by uncertainties in the covariance matrix

Data analysis from upcoming large galaxy redshift surveys, such as Euclid and DESI will significantly improve constraints on cosmological parameters. To optimally extract the information from these galaxy surveys, it is important to control with a high level of confidence the uncertainty and bias arising from the estimation of the covariance that affects the inference of cosmological parameters. In this work, we are addressing two different but closely related issues: (i) the sampling noise present in a covariance matrix estimated from a finite set of simulations and (ii) the impact on cosmological constraints of the non-Gaussian contribution to the covariance matrix of the power spectrum. We focus on the parameter estimation obtained from fitting the matter power spectrum in real space, using the DEMNUni N-body simulations. Regarding the first issue, we adopt two different approaches to reduce the sampling noise in the precision matrix that propagates in the parameter space: on the one hand using an alternative estimator of the covariance matrix based on a non-linear shrinkage, NERCOME; and on the other hand employing a method of fast generation of approximate mock catalogs, COVMOS. We find that NERCOME can significantly reduce the noise induced on the posterior distribution of parameters, but at the cost of a systematic overestimation of the error bars on the cosmological parameters. We show that using a COVMOS covariance matrix estimated from a large number of realisations (10~000) results in unbiased cosmological constraints. Regarding the second issue, we quantify the impact on cosmological constraints of the non-Gaussian part of the power spectrum covariance purely coming from non-linear clustering. We find that when this term is neglected, both the errors and central values of the estimated parameters are affected for a scale cut \kmax > 0.2\ \invMpc.Comment: 21 pages, 2 appendices, 20 figures. Submitted to A&