Researcher profile

Philipp Mertsch

Philipp Mertsch contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2021arXiv

No longer ballistic, not yet diffusive--the formation of cosmic ray small-scale anisotropies

The arrival directions of TeV-PeV cosmic rays are remarkably uniform due to the isotropization of their directions by scattering on turbulent magnetic fields. Small anisotropies can exist in standard diffusion models, however, only on the largest angular scales. Yet, high-statistics observatories like IceCube and HAWC have found significant deviations from isotropy down to small angular scales. Here, we explain the formation of small-scale anisotropies by considering pairs of cosmic rays that get correlated by their transport through the same realisation of the turbulent magnetic field. We argue that the formation of small-scale anisotropies is the reflection of the particular realisation of the turbulent magnetic field experienced by cosmic rays on time scales intermediate between the early, ballistic regime and the late, diffusive regime. We approach this problem in two different ways: First, we run test particle simulations in synthetic turbulence, covering for the first time the TV rigidities of observations with realistic turbulence parameters. Second, we extend the recently introduced mixing matrix approach and determine the steady-state angular power spectrum. Throughout, we adopt magneto-static, slab-like turbulence. We find excellent agreement between the predicted angular power spectra in both approaches over a large range of rigidities. In the future, measurements of small-scale anisotropies will be valuable in constraining the nature of the turbulent magnetic field in our Galactic neighborhood.

preprint2021arXiv

Self-confinement of low-energy cosmic rays around supernova remnants

Supernova remnants have long been considered as a promising candidate for sources of Galactic cosmic rays. However, modelling cosmic-ray transport around these sources is complicated by the fact that the overdensity of cosmic rays close to their acceleration site can lead to self-confinement, that is the generation of turbulence upon which these particles scatter. Such a highly non-linear problem can be addressed by numerically solving the coupled differential equations describing the evolution in space and time of the escaping particles and of the turbulent plasma waves. In this work, we focus essentially on the propagation of cosmic rays from supernova remnants in the warm ionized and warm neutral phases of the interstellar medium and propose an extended framework to take into account also the effect of energy loss relevant for cosmic rays of energy below 10 GeV. Interestingly, the diffusion coefficient of low-energy cosmic rays could be suppressed by up to 2 orders of magnitude for several tens of kiloyears after the escape from the shock. The cosmic-ray spectrum outside the supernova remnant flattens below 1 GeV at a sufficiently late time reminiscient of the spectral behaviour observed by Voyager. We also find the grammage accumulated around the source to be non-negligible, with important implications for precision fitting of the cosmic-ray spectra.

preprint2020arXiv

Bayesian inference of three-dimensional gas maps: I. Galactic CO

Carbon monoxide (CO) is the best tracer of Galactic molecular hydrogen (H2). Its lowest rotational emission lines are in the radio regime and thanks to Galactic rotation emission at different distances is Doppler shifted. For a given gas flow model the observed spectra can thus be deprojected along the line of sight to infer the gas distribution. We use the CO line survey of Dame et al. (2001) to reconstruct the three-dimensional density of H2. We consider the deprojection as a Bayesian variational inference problem. The posterior distribution of the gas densities allows us to estimate both the mean and uncertainty of the reconstructed density. Unlike most of the previous attempts, we take into account the correlations of gas on a variety of scales which allows curing some of the well-known pathologies, like fingers-of-god effects. Both gas flow models that we adopt incorporate a Galactic bar which induces radial motions in the inner few kiloparsecs and thus offers spectral resolution towards the Galactic centre. We compare our gas maps with those of earlier studies and characterise their statistical properties, e.g. the radial profile of the average surface mass density. We have made our three-dimensional gas maps and their uncertainties available to the community at https://dx.doi.org/10.5281/zenodo.4405437 .

preprint2020arXiv

Test particle simulations of cosmic rays

Modelling of cosmic ray transport and interpretation of cosmic ray data ultimately rely on a solid understanding of the interactions of charged particles with turbulent magnetic fields. The paradigm over the last 50 years has been the so-called quasi-linear theory, despite some well-known issues. In the absence of a widely accepted extension of quasi-linear theory, wave-particle interactions must also be studied in numerical simulations where the equations of motion are directly solved in a realisation of the turbulent magnetic field. The applications of such test particle simulations of cosmic rays are manifold: testing transport theories, computing parameters like diffusion coefficients or making predictions for phenomena beyond standard diffusion theories, e.g. for cosmic ray small-scale anisotropies. In this review, we seek to give a low-level introduction to test particle simulations of cosmic rays, enabling readers to perform their own test particle simulations. We start with a review of quasi-linear theory, highlighting some of its issues and suggested extensions. Next, we summarise the state-of-the-art in test particle simulations and give concrete recipes for generating synthetic turbulence. We present a couple of examples for applications of such simulations and comment on an important conceptual detail in the backtracking of particles.

preprint2014arXiv

Fingerprints of Galactic Loop I on the Cosmic Microwave Background

We investigate possible imprints of galactic foreground structures such as the "radio loops" in the derived maps of the cosmic microwave background. Surprisingly there is evidence for these not only at radio frequencies through their synchrotron radiation, but also at microwave frequencies where emission by dust dominates. This suggests the mechanism is magnetic dipole radiation from dust grains enriched by metallic iron or ferrimagnetic materials. This new foreground we have identified is present at high galactic latitudes, and potentially dominates over the expected $B$-mode polarization signal due to primordial gravitational waves from inflation.

preprint2013arXiv

Loops and spurs: The angular power spectrum of the Galactic synchrotron background

We present a new model of the diffuse Galactic synchrotron radiation, concentrating on its angular anisotropies. While previous studies have focussed on either the variation of the emissivity on large (kpc) scales, or on fluctuations due to MHD turbulence in the interstellar medium, we unify these approaches to match the angular power spectrum. We note that the usual turbulence cascade calculation ignores spatial correlations at the injection scale due to compression of the interstellar medium by old supernova remnants -- the 'radio loops', only four of which are visible by eye in radio maps. This new component naturally provides the otherwise missing power on intermediate and small scales in the all-sky map at 408 MHz. Our model can enable more reliable subtraction of the synchrotron foreground for studies of CMB anisotropies (both in temperature and polarisation) or searches for dark matter annihilation. We conclude with some remarks on the relevance to modelling of the polarised foreground.

preprint2012arXiv

Loop-induced dark matter direct detection signals from gamma-ray lines

Improved limits as well as tentative claims for dark matter annihilation into gamma-ray lines have been presented recently. We study the direct detection cross section induced from dark matter annihilation into two photons in a model-independent fashion, assuming no additional couplings between dark matter and nuclei. We find a striking non-standard recoil spectrum due to different destructively interfering contributions to the dark matter nucleus scattering cross section. While in the case of s-wave annihilation the current sensitivity of direct detection experiments is insufficient to compete with indirect detection searches, for p-wave annihilation the constraints from direct searches are comparable. This will allow to test dark matter scenarios with p-wave annihilation that predict a large di-photon annihilation cross section in the next generation of experiments.

preprint2011arXiv

2nd-order Fermi acceleration as the origin of the Fermi bubbles

Gamma-ray data from Fermi-LAT show a bi-lobular structure extending up to 50 degrees above and below the Galactic centre, coincident with a possibly related structure in the ROSAT X-ray map which presumably originated in some energy release close to the centre a few million years ago. It has been argued that the gamma-rays arise due to inverse Compton scattering of relativistic electrons accelerated at plasma shocks present in the bubbles. We explore the alternative possibility that the relativistic electrons undergo stochastic 2nd-order Fermi acceleration in the entire volume of the bubbles by plasma wave turbulence. This turbulence is generated behind the outer shock and propagates into the bubble volume, leading to a non-trivial spatial variation of the electron spectral index. Rather than a constant volume emissivity as predicted in other models we find an almost constant surface brightness in gamma-rays and also reproduce the observed sharp edges of the bubbles. We comment on possible cross-checks in other channels.

preprint2011arXiv

A new analytic solution for 2nd-order Fermi acceleration

A new analytic solution for 2nd-order Fermi acceleration is presented. In particular, we consider time-dependent rates for stochastic acceleration, diffusive and convective escape as well as adiabatic losses. The power law index q of the turbulence spectrum is unconstrained and can therefore account for Kolmogorov (q = 5/3) and Kraichnan (q = 3/2) turbulence, Bohm diffusion (q = 1) as well as the hard-sphere approximation (q = 2). This considerably improves beyond solutions known to date and will prove a useful tool for more realistic modelling of 2nd-order Fermi acceleration in a variety of astrophysical environments.

preprint2011arXiv

Fermi gamma-ray `bubbles' from stochastic acceleration of electrons

Gamma-ray data from Fermi-LAT reveal a bi-lobular structure extending up to 50 degrees above and below the galactic centre, which presumably originated in some form of energy release there less than a few million years ago. It has been argued that the gamma-rays arise from hadronic interactions of high energy cosmic rays which are advected out by a strong wind, or from inverse-Compton scattering of relativistic electrons accelerated at plasma shocks present in the bubbles. We explore the alternative possibility that the relativistic electrons are undergoing stochastic 2nd-order Fermi acceleration by plasma wave turbulence through the entire volume of the bubbles. The observed gamma-ray spectral shape is then explained naturally by the resulting hard electron spectrum and inverse Compton losses. Rather than a constant volume emissivity as in other models, we predict a nearly constant surface brightness, and reproduce the observed sharp edges of the bubbles.

preprint2011arXiv

The high energy neutrino cross-section in the Standard Model and its uncertainty

Updated predictions are presented for high energy neutrino and antineutrino charged and neutral current cross-sections within the conventional DGLAP formalism of NLO QCD using modern PDF fits. PDF uncertainties from model assumptions and parametrization bias are considered in addition to the experimental uncertainties. Particular attention is paid to assumptions and biases which could signal the need for extension of the conventional formalism to include effects such as ln(1/x) resummation or non-linear effects of high gluon density.

preprint2010arXiv

Cosmic ray backgrounds for dark matter indirect detection

Recently, dark matter indirect searches have gained a lot of attention, mostly due to the possibility of recent anomalies in cosmic rays and microwave sky maps being due to the annihilation or decay of dark matter. In this thesis, we argue however that these signals are plagued by irreducible astrophysical backgrounds and show how plausible conventional physics can mimic the alleged dark matter signals. In particular, we consider the possibility that the rise in the positron fraction observed by the PAMELA satellite is due to the production through (hadronic) cosmic ray spallation and subsequent acceleration of positrons, in the same sources as the primary cosmic rays. We present a new analytical estimate of the range of possible fluctuations in the high energy electron flux due to the discreteness of plausible cosmic ray sources. Fitting our result for the total electron-positron flux measured by the Fermi satellite allows us to fix the only free parameter of the model and make an independent prediction for the positron fraction. This model can be tested by considering similar effects expected for nuclear secondary-to-primary ratios such as B/C. A rise predicted above O(100) GeV/n would be an unique confirmation of our explanation for a rising positron fraction and rule out the dark matter explanation. Furthermore, we review the assumptions made in the extraction of the `WMAP haze' which has also been claimed to be due to electrons and positrons from dark matter annihilation in the Galactic centre region. We argue that the energy-dependence of their diffusion makes the extraction of the haze through template fitting unreliable. The systematic effects introduced by this can, under specific circumstances, reproduce the residual, suggesting that the `haze' may be just an artefact of the template subtraction.

preprint2010arXiv

Systematic effects in the extraction of the 'WMAP haze'

The extraction of a 'haze' from the WMAP microwave skymaps is based on subtraction of known foregrounds, viz. free-free (bremsstrahlung), thermal dust and synchrotron, each traced by other skymaps. While the 408 MHz all-sky survey is used for the synchrotron template, the WMAP bands are at tens of GHz where the spatial distribution of the radiating cosmic ray electrons ought to be quite different because of the energy-dependence of their diffusion in the Galaxy. The systematic uncertainty this introduces in the residual skymap is comparable to the claimed haze and can, for certain source distributions, have a very similar spectrum and latitudinal profile and even a somewhat similar morphology. Hence caution must be exercised in interpreting the 'haze' as a physical signature of, e.g., dark matter annihilation in the Galactic centre.