
COSMICMAG
Evaluation of the Uncertainties of the Galactic Magnetic Field to Elucidate the Origin of Ultra-High Energy Cosmic Rays
The COSMICMAG project was an International Outgoing Fellowship
funded by European Research Fellowship Programme (Marie Skłodowska-Curie actions) in collaboration with
Prof. Glennys Farrar at NYU.
A brief summary of the project is given below.
The main results of the project were published
here.
Project
COSMICMAG aims at the development of a new analysis approach for charged particle
astronomy to find the astrophysical sources of cosmic particles with
extreme energies.
During the recent years, there has been a rapid progress in our
understanding of Galactic magnetic field (GMF) and the collected
statistics of ultra-high energy cosmic ray events. It should therefore
be possible to approximately correct the measured arrival directions of cosmic rays
on Earth for deflections in the GMF and to significantly improve the
search for their astrophysical sources in this way. The research
objectives of COSMICMAG are therefore: 1) A thorough evaluation of
the uncertainties of the GMF with the aim of producing sky maps of the
deflection uncertainty of cosmic rays given their arrival direction
and rigidity (= energy divided by charge). 2) Development of an optimal
analysis method for correlation studies that takes into account the
deflection uncertainties originating from the uncertainties in both
GMF and the estimated rigidities. 3) Application of the optimal data
analysis method to data of the Pierre Auger Observatory.
Gallery of Results
Below a few pictures illustrating the results and work in progress. Images can be enlarged by clicking on them.

The magnetic field of our Galaxy (GMF) can be estimated using astrophysical
observations, in particular the Farraday rotation of polarized radio emission from astrophysical sources
and synchrotron
radiation of cosmic-ray electrons spiralling in the GMF.
From these observations, parametric models of the GMF can
be constructed as for instance the one of
Jansson and
Farrar (2012) which is visualized in the left figure above.
The main result of the COSMICMAG project is a thorough evaluation of
the uncertainties of estimates of the GMF, resulting from different
functional forms used for the modeling, different data sets used for
the tuning of parameters and different auxillary assumptions used to
interpret the data. The effect of GMF uncertainties on our ability to
estimate the arrival directions of extragalactic ultrahigh-energy
cosmic rays is illustrated in the right panel of the figure above,
where two sky maps are shown
in Galactic
coordinates. A grid of representative arrival directions at Earth
is indicated as colored points. These directions were back-tracked
through different models of the GMF to obtain the arrival direction at
the edge of the Galaxy indicated as lettered squares. The different
letters identify different models of the GMF. Each of them gives a
reasonable description of the available astrophysical data. As can be
seen, for cosmic rays with a large energy-to-charge ratio (60 EV, top
panel) the overall amount of deflection and correspondingly also the
model differences are small, confirming the long-speculated
possibility of charged particle astronomy with protons at ultrahigh
energies. Even for energy-to-charge ratios as low as 20 EV, the
different deflections are mostly confined within well-defined regions
so it seems plausible that a correction for the spatially varying
average deflection based on all models, can still be used to enhance
the capabilities to identify the sources of ultrahigh energy cosmic
rays. More details can be found
in our ICRC
contribution.


The spectral "ankle"-feature in the flux of ultra-high energy cosmic rays as well as the proton-dominated composition observed at the presumed transition from Galactic to extra-galactic cosmic rays can be naturally explained by postulating a suitable source environment around the sources of ultra-high energy cosmic rays. This was shown in collaboration with L. Anchordoqui and G. Farrar. The left picture illustrates our model for extragalactic cosmic rays. Their sources (denoted by yellow stars) are situated in an environment such that the ultrahigh-energy cosmic rays partially photodisnitegrate via interactions with the ambient photon field, before they escape. The escape time depends on the strength of the random magnetic field around the sources. The right picture shows a comparison of the model prediction to the flux measured by the Pierre Auger Observatory. The black solid line is the total model prediction and the colored lines show the contributions for increasingly massive particles starting from protons (red) to silicon (cyan). More details can be found in our paper.

The interstellar medium permeates the space between stars in galaxies. Its properties are of scientific interest on their own, but in the context of the COSMICMAG project we are mainly interested in its influence on the rotation measures that we use to estimate the strength of the Galactic magnetic field. Previous attempts to model the Galactic magnetic field relied on the NE2001 model for the thermal electrons in the interstellar medium. The figure shows the stage of our development of an improved model of thermal electrons. On the left a sky map of the measured Hα emission is shown (composite of VTSS, SHASSA, WHAM data). On the right the emission simulated with our preliminary thermal electron model is shown, including the emission of known HII regions. The gray area masks regions with large obscuration from dust. See this presentation for more information.

The study of the magnetic field of our Galaxy is complicated by the fact that we are observing it
from the inside. To judge whether the modeled magnetic field gives
a fair representation of observations of other galaxies, a simulated view from an imaginary outside observer
is useful. Previous simulation programs such as Hammurabi could only produce simulated sky
maps for an observer sitting at the location of the Earth. This restriction
does not apply to the RUQI
framework developed for the COSMICMAG project. This feature is illustrated in the above figure.
The three panels show simulated observables for the JF12 magnetic field model viewed from 45 degree above the Galactic plane. One the left the
rotation measure is shown and the Stokes parameters Q and U are displayed in the middle and right panel.
In addition to the free choice
of observer position, the RUQI framework supports adaptive line-of-sight integration and dynamic loading of astrophysical models.


The idea to detect cosmic rays with a cheap citizen-science approach was put forward by the DECO and CRAYFIS groups. On the left the population density is shown using data from the Global Rural-Urban Mapping Project . On the right the amount of user participation needed to obtain a similar area as "professional" cosmic ray observatories can be seen. From these figures it can be concluded that it may indeed be possible to detect a few events with a network of smartphones, but a large enough area to collect enough events to perform high-quality science seems difficult to reach, even under favorable assumptions. See our preprint for more information.
Publications
Publications related to the COSMICMAG project:-
Uncertainties in the Magnetic Field of the Milky Way
M. Unger and G.R. Farrar, PoS(ICRC2017) 558 (preprint) -
Origin of the ankle in the ultrahigh energy cosmic ray spectrum, and
of the extragalactic protons below it
M. Unger, G.R. Farrar and L.A. Anchordoqui, Phys.Rev. D92 (2015) no.12, 123001 (preprint) -
Detailed simulations of Fermi-LAT constraints on UHECR production scenarios
M. Muzio, G.R. Farrar and M. Unger, PoS(ICRC2017) 557 - (In)Feasibility of
Studying Ultra-High-Energy Cosmic Rays with Smartphones
M. Unger and G.R. Farrar, arXiv:1505.04777 [astro-ph.IM] (preprint)
Talks and Posters
Galactic Magnetic Field
- Uncertainties in the Magnetic Field of the Milky Way
poster presented at the 35th International Cosmic Ray Conference, Busan, Korea, 2017 - Ultrahigh Energy Cosmic Rays and
the Galactic Magnetic Field
talk given at the ULB Seminar, University of Brussels, Belgium, 2017 - Ultrahigh Energy Cosmic Rays and
the Galactic Magnetic Field
talk given at the HEP seminar, Penn State University, USA, 2017 - Uncertainties of Models of the Galactic Magnetic Field
talk given at the Pierre Auger Collaboration meeting, Malargüe, Argentina, 2017 - Global Modeling of the Galactic Magnetic Field and Thermal Electrons
talk given at the Crutcher and Heiles Conference (Star formation, magnetic fields, and diffuse matter in the Galaxy), Madison, USA, 2016 - Air Showers and Magnetic Field Deflections
talk given at the 1st ISSI International Team Meeting on Bayesian modeling of the Galactic magnetic field, Bern, Switzerland, 2014
Ultra-High Energy Cosmic Rays
- Highlights from the Pierre Auger Observatory
invited plenary talk at the 35th International Cosmic Ray Conference, Busan, Korea, 2017 - Origin of the Ankle and the Light Composition at EeV Energies
talk given at the NEUCOS workshop, DESY Zeuthen, Germany, 2017. - Mass Composition of Ultra-high Energy Cosmic Rays
talk given at HEP seminar, Institute of Physics of the Czech Academy of Sciences, Prague, Czech Republic, 2016. - Mass Composition of Cosmic Rays: Introduction and Overview
invited talk at the HAP Workshop on the Non-Thermal Universe, Erlangen, Germany, 2016 - Ultra-High Energy Cosmic Rays
invited talk at the workshop on Multi-Messenger Approaches to Cosmic Rays (MACROS), Penn State, USA, 2016 - Mass Composition and Hadronic Interactions
invited talk at the Very High Energy Particle Astronomy Workshop (VHEPA), Honolulu, USA, 2016 - Mass Composition Studies
with the Pierre Auger Observatory
invited talk at the workshop on Sources of Galactic cosmic rays (SOCOR), Paris, France, 2016 - Report of the Working Group on the Composition of
Ultra-High Energy Cosmic Rays
contributed talk at the 34th International Cosmic Ray Conference (ICRC), The Hague, Netherlands, 2015 - Ultra-High Energy Cosmic Rays
invited talk at the IceCube Particle Astrophysics Symposium (IPA), Madison, USA, 2015
Teaching
As part of the workplan of the project, a graduate course on High Energy Astrophysics was prepared for the fall semester 2015 at NYU.Software
The software developed for COSMICMAG is hosted on the IAP subversion server. A first public release will be made available after the end of the project (second half of 2017).- RUQI: a C++ framework for the calculation rotation measures, dispersion measures, emission measures, scattering measures and synchrotron emission in galaxies. Supports adaptive line-of-sight integration, dynamic loading of astrophysical models, arbitrary observer positions etc.
- AstroTools: A collection of C++ utilities for astrophysical calculations. Includes streaming and visualization of HEALPix data using the ROOT data analysis framework.
Useful Links
- ATNF pulsar catalogue
- NVSS rotation measures catalogue
- WHAM survey of Hα emission
- Planck legacy archive
- WISE catalogue of HII regions
- Interstellar dust reddening with Pan-STARRS 1
- Kinematic Distance Calculator
- Lutz-Kelker Bias Tool
- NE2001 thermal electron model
- YMW16 thermal electron model
- HEALPix software for equal-area pixelization on a sphere
- CRT software for tracking of cosmic rays through magnetic fields
- CRPropa software for propagation of ultra-high energy cosmic rays
- Hammurabi software for the simulation of magnetic-field related observables
- DRAGON software for transport of cosmic rays in the galaxy
- GALPROP software for transport of cosmic rays in the galaxy
- ISSI team on Bayesian modeling of the galactic magnetic field