Normalized to: Pumpe, D.

**4** article(s) in total. **14** co-authors, from **1** to **4** common article(s). Median position in authors list is **1,0**.

[1] oai:arXiv.org:1802.02013 [pdf] - 1789403

Denoising, deconvolving and decomposing multi-domain photon
observations- The D4PO algorithm

Comments:

Submitted: **2018-02-06**, last modified: **2018-08-27**

Astronomical imaging based on photon count data is a non-trivial task. In
this context we show how to denoise, deconvolve, and decompose multi-domain
photon observations. The primary objective is to incorporate accurate and well
motivated likelihood and prior models in order to give reliable estimates about
morphologically different but superimposed photon flux components present in
the data set. Thereby we denoise and deconvolve photon counts, while
simultaneously decomposing them into diffuse, point-like and uninteresting
background radiation fluxes. The decomposition is based on a probabilistic
hierarchical Bayesian parameter model within the framework of information field
theory (IFT). In contrast to its predecessor D$^3$PO, D$^4$PO reconstructs
multi-domain components. Thereby each component is defined over its own direct
product of multiple independent domains, for example location and energy.
D$^4$PO has the capability to reconstruct correlation structures over each of
the sub-domains of a component separately. Thereby the inferred correlations
implicitly define the morphologically different source components, except for
the spatial correlations of the point-like flux. Point-like source fluxes are
spatially uncorrelated by definition. The capabilities of the algorithm are
demonstrated by means of a synthetic, but realistic, mock data set, providing
spectral and spatial information about each detected photon. D$^4$PO
successfully denoised, deconvolved, and decomposed a photon count image into
diffuse, point-like and background flux, each being functions of location as
well as energy. Moreover, uncertainty estimates of the reconstructed fields as
well as of their correlation structure are provided employing their posterior
density function and accounting for the manifolds the domains reside on.

[2] oai:arXiv.org:1708.05702 [pdf] - 1644122

Search for quasi-periodic signals in magnetar giant flares

Comments:

Submitted: **2017-08-18**, last modified: **2018-03-05**

Quasi-periodic oscillations (QPOs) discovered in the decaying tails of giant
flares of magnetars are believed to be torsional oscillations of neutron stars.
These QPOs have a high potential to constrain properties of high-density
matter. In search for quasi-periodic signals, we study the light curves of the
giant flares of SGR 1806-20 and SGR 1900+14, with a non-parametric Bayesian
signal inference method called D$^3$PO. The D$^3$PO algorithm models the raw
photon counts as a continuous flux and takes the Poissonian shot noise as well
as all instrument effects into account. It reconstructs the logarithmic flux
and its power spectrum from the data. Using this fully noise-aware method, we
do not confirm previously reported frequency lines at $\nu\gtrsim17\,$Hz
because they fall into the noise-dominated regime. However, we find two new
potential candidates for oscillations at $9.2\,$Hz (SGR 1806-20) and $7.7\,$Hz
(SGR 1900+14). If these are real and the fundamental magneto-elastic
oscillations of the magnetars, current theoretical models would favour
relatively weak magnetic fields $\bar B\sim 6\times10^{13} - 3\times10^{14}\,$G
(SGR 1806-20) and a relatively low shear velocity inside the crust compared to
previous findings.

[3] oai:arXiv.org:1708.01073 [pdf] - 1586716

NIFTy 3 - Numerical Information Field Theory - A Python framework for
multicomponent signal inference on HPC clusters

Comments: 18 pages, 7 figures, 1 table, available at
https://gitlab.mpcdf.mpg.de/ift/NIFTy/

Submitted: **2017-08-03**

NIFTy, "Numerical Information Field Theory", is a software framework designed
to ease the development and implementation of field inference algorithms. Field
equations are formulated independently of the underlying spatial geometry
allowing the user to focus on the algorithmic design. Under the hood, NIFTy
ensures that the discretization of the implemented equations is consistent.
This enables the user to prototype an algorithm rapidly in 1D and then apply it
to high-dimensional real-world problems. This paper introduces NIFTy 3, a major
upgrade to the original NIFTy framework. NIFTy 3 allows the user to run
inference algorithms on massively parallel high performance computing clusters
without changing the implementation of the field equations. It supports
n-dimensional Cartesian spaces, spherical spaces, power spaces, and product
spaces as well as transforms to their harmonic counterparts. Furthermore, NIFTy
3 is able to treat non-scalar fields. The functionality and performance of the
software package is demonstrated with example code, which implements a real
inference algorithm from the realm of information field theory. NIFTy 3 is
open-source software available under the GNU General Public License v3 (GPL-3)
at https://gitlab.mpcdf.mpg.de/ift/NIFTy/

[4] oai:arXiv.org:1601.07901 [pdf] - 1443808

Dynamic system classifier

Comments: 11 pages, 8 figures

Submitted: **2016-01-28**, last modified: **2016-07-26**

Stochastic differential equations describe well many physical, biological and
sociological systems, despite the simplification often made in their
derivation. Here the usage of simple stochastic differential equations to
characterize and classify complex dynamical systems is proposed within a
Bayesian framework. To this end, we develop a dynamic system classifier (DSC).
The DSC first abstracts training data of a system in terms of time dependent
coefficients of the descriptive stochastic differential equation. Thereby the
DSC identifies unique correlation structures within the training data. For
definiteness we restrict the presentation of DSC to oscillation processes with
a time dependent frequency {\omega}(t) and damping factor {\gamma}(t). Although
real systems might be more complex, this simple oscillator captures many
characteristic features. The {\omega} and {\gamma} timelines represent the
abstract system characterization and permit the construction of efficient
signal classifiers. Numerical experiments show that such classifiers perform
well even in the low signal-to-noise regime.