naima: a Python package for inference of relativistic particle energy distributions from observed nonthermal spectra

The ultimate goal of the observation of nonthermal emission from astrophysical sources is to understand the underlying particle acceleration and evolution processes, and few tools are publicly available to infer the particle distribution properties from the observed photon spectra from X-ray to VHE gamma rays. Here I present naima, an open source Python package that provides models for nonthermal radiative emission from homogeneous distribution of relativistic electrons and protons. Contributions from synchrotron, inverse Compton, nonthermal bremsstrahlung, and neutral-pion decay can be computed for a series of functional shapes of the particle energy distributions, with the possibility of using user-defined particle distribution functions. In addition, naima provides a set of functions that allow to use these models to fit observed nonthermal spectra through an MCMC procedure, obtaining probability distribution functions for the particle distribution parameters. Here I present the models and methods available in naima and an example of their application to the understanding of a galactic nonthermal source. naima's documentation, including how to install the package, is available at this http URL

[1]  F. Aharonian,et al.  SIMPLE ANALYTICAL APPROXIMATIONS FOR TREATMENT OF INVERSE COMPTON SCATTERING OF RELATIVISTIC ELECTRONS IN THE BLACKBODY RADIATION FIELD , 2013, 1310.7971.

[2]  S. Reynolds,et al.  Radio to Gamma-Ray Emission from Shell-Type Supernova Remnants: Predictions from Nonlinear Shock Acceleration Models , 1998, astro-ph/9810158.

[3]  F. Aharonian,et al.  Parametrization of gamma-ray production cross-sections for pp interactions in a broad proton energy range from the kinematic threshold to PeV energies , 2014, 1406.7369.

[4]  Felix Aharonian,et al.  Energy spectra of gamma rays, electrons, and neutrinos produced at proton-proton interactions in the very high energy regime , 2006 .

[5]  G. Blumenthal,et al.  BREMSSTRAHLUNG, SYNCHROTRON RADIATION, AND COMPTON SCATTERING OF HIGH- ENERGY ELECTRONS TRAVERSING DILUTE GASES. , 1970 .

[6]  David J. C. MacKay,et al.  Information Theory, Inference, and Learning Algorithms , 2004, IEEE Transactions on Information Theory.

[7]  F. Aharonian,et al.  Angular, spectral, and time distributions of highest energy protons and associated secondary gamma-rays and neutrinos propagating through extragalactic magnetic and radiation fields , 2010, 1006.1045.

[8]  John D. Hunter,et al.  Matplotlib: A 2D Graphics Environment , 2007, Computing in Science & Engineering.

[9]  T. Koi,et al.  Parameterization of γ, epm, and Neutrino Spectra Produced by p-p Interaction in Astronomical Environments , 2006, astro-ph/0605581.

[10]  F. Aharonian,et al.  Hadronic gamma-rays from RX J1713.7−3946? , 2014, 1406.2322.

[11]  Sebastian Funk,et al.  Primary particle acceleration above 100 TeV in the shell-type supernova remnant RX J1713.7-3946 with deep HESS observations , 2007 .

[12]  Daniel Foreman-Mackey,et al.  emcee: The MCMC Hammer , 2012, 1202.3665.

[13]  Prasanth H. Nair,et al.  Astropy: A community Python package for astronomy , 2013, 1307.6212.