Galactic Magnetic Field Models Research Articles

Effects of the Galactic Magnetic Field on the UHECR Correlation Studies with Starburst Galaxies

We estimate the biases caused by the coherent deflection of cosmic rays due to the Galactic magnetic field (GMF) in maximum likelihood analyses for searches of ultrahigh-energy cosmic-ray sources in the literature. We simulate mock event data sets with a set of assumptions for the starburst galaxy source model, coherent deflection by a GMF model, and mixed-mass composition. We then conduct a maximum likelihood analysis without accounting for the GMF in the same manner as previous studies. We find that the anisotropic fraction f ani is estimated to be systematically lower than the true value. We estimate the true parameters that are compatible with the best-fit parameters that were reported, and find that except for a narrow region with a large anisotropic fraction and a small separation angular scale, a wide parameter space is still compatible with the experimental results. We also develop a maximum likelihood method that takes the GMF model into account and confirm in Monte Carlo simulations that we can estimate the true parameters within a 1σ contour under the ideal condition that we know the event-by-event mass and the GMF.

Using the “least-deflected” subsamples of ultra-high energy cosmic rays to constrain source population(s)

Context. The source population(s) of ultra-high energy cosmic rays (UHECRs) with energies of ≥57 × 1018 eV (57 EeV) are still unknown and continue to be a subject of debate. The published arrival directions and energies of 303 UHECRs from the Pierre Auger Observatory and Telescope Array enable correlations with astronomical sources. However, deflections of UHECRs in the Galactic Magnetic Field (GMF) can be considerable, especially for heavy UHECR compositions. Aims. We aim to constrain the source(s) of UHECRs by focusing on the subset of UHECRs that experience minimal deflections by the Galactic magnetic field (GMF). Methods. We used Monte Carlo simulations of UHECR trajectories in a total of eight posited GMF models (and six UHECR compositions from H to Fe) to identify two subsamples of about 40 UHECRs. The H compositions of these subsamples suffer deflections of less than a few degrees in one (or most) of the GMFs. Both the total UHECR sample and the “least-deflected” (vs. the remaining “more-deflected”) subsamples of UHECRs were cross-correlated with catalogs of astronomical sources (at D ≲ 100 Mpc) and neutrino arrival directions. Results. For H compositions, the full sample of UHECRs is most closely correlated with Cen A, nearby (< 75 Mpc) radiogalaxies, and optically selected active galactic nuclei (AGNs); the least-deflected sub sample shows a most consistent correlation with nearby (D ≲ 25 Mpc) galaxies, although some GMF models show a preference for the IceCube muon-neutrino arrival detections. For oxygen compositions, the out-of-galaxy arrival directions of least-deflected UHECRs fall on the extragalactic plane (and thus close to nearby galaxies) in many GMFs. For Fe compositions, the out-of-galaxy UHECR arrival directions typically originate in a region ∼30 deg offset from Cen A, except in the case of the JF12 model; this offset can be decreased by re-scaling (by ∼50%) the magnitudes of the dipole and disk fields of the relevant GMFs.

Treasure Maps for Detections of Extreme Energy Cosmic Rays

The origin of ultra-high-energy cosmic rays is a 60 yr old mystery. We show that with more events at the highest energies (above 150 EeV) it may be possible to limit the character of the sources and learn about the intervening magnetic fields. Individual sources become more prominent, relative to the background, as the horizon diminishes. An event-by-event, composition-dependent observatory would allow a “tomography” of the sources as different mass and energy groups probe different GZK horizons. A major goal here is to provide a methodology to distinguish between steady and transient or highly variable sources. Using recent Galactic magnetic field models, we calculate “treasure” sky maps to identify the most promising directions for detecting Extreme Energy Cosmic Rays doublets, events that are close in arrival time and direction. On this basis, we predict the incidence of doublets as a function of the nature of the source host galaxy. Based on the asymmetry in the distribution of time delays, we show that observation of doublets might distinguish source models. In particular, the Telescope Array hotspot could exhibit temporal variability as it is in a “magnetic window” of small time delays. These considerations could improve the use of data with existing facilities and the planning of future ones such as Global Cosmic Ray Observatory (GCOS).

Estimation and reduction of the biases by the galactic magnetic field on the UHECR correlation studies

The correlation studies between ultrahigh-energy cosmic ray (UHECR) anisotropy and source candidates are keys to understand the origin of UHECRs. Especially, studies with consideration of source models, magnetic fields, and mass compositions are necessary. We estimated the biases caused by the coherent deflection due to the galactic magnetic field (GMF) in the previous maximum-likelihood analysis for searching the UHECR sources (Aab et al. 2018, Abbasi et al. 2018). In our work with simulated mock datasets, we find that the anisotropic fraction fani is estimated systematically lower than the true value when we ignore the effect caused by the GMF. We also develop the maximum-likelihood method which includes the GMF model and confirm that the estimated parameters would be improved. We apply the method to the observational datasets obtained from the Telescope Array and Auger experiments.

Galactic halo bubble magnetic fields and UHECR deflections

ABSTRACT We consider the synchrotron emission from electrons out in the Galactic halo bubble region where the Fermi bubble structures reside. Utilizing a simple analytical expression for the non-thermal electron distribution and a toy magnetic field model, we simulate polarized synchrotron emission maps at a frequency of 30 GHz. Comparing these maps with the observational data, we obtain constraints on the parameters of our toy Galactic halo bubble magnetic field model. Utilizing this parameter value range for the toy magnetic field model, we determine the corresponding range of arrival directions and suppression factors of ultra high energy cosmic rays (UHECRs) from potential local source locations. We find that high levels of flux suppression (down to 2 per cent) and large deflection angles (≥80°) are possible for source locations whose line of sight pass through the Galactic halo bubble region. We conclude that the magnetic field out in the Galactic halo bubble region can strongly dominate the level of deflection UHECRs experience whilst propagating from local sources to Earth.

Orion A’s complete 3D magnetic field morphology

Magnetic fields permeate the interstellar medium and are important in the star formation process. Determining the three-dimensional (3D) magnetic fields of molecular clouds will allow us to better understand their role in the evolution of these clouds and the formation of stars. We fully reconstruct the approximate 3D magnetic field morphology of the Orion A molecular cloud (on scales of a few to ∼100 pc) using Galactic magnetic field models, as well as available line-of-sight and plane-of-sky magnetic field observations. While previous studies identified the 3D magnetic field morphology of the Orion A cloud as an arc shape, in this study we provide the orientation of this arc-shaped field and its plane-of-sky direction for the first time. We find that this 3D field is a tilted, semi-convex (from our point of view) structure and mostly points in the direction of decreasing latitude and longitude on the plane of the sky from our vantage point. The previously identified bubbles and events in this region were key in shaping this arc-shaped magnetic field morphology.

3D magnetic-field morphology of the Perseus molecular cloud

Context. Despite recent observational and theoretical advances in mapping the magnetic fields associated with molecular clouds, their three-dimensional (3D) morphology remains unresolved. Multi-wavelength and multi-scale observations will allow us to paint a comprehensive picture of the magnetic fields of these star-forming regions. Aims. We reconstructed the 3D magnetic field morphology associated with the Perseus molecular cloud and compared it with predictions of cloud-formation models. These cloud-formation models predict a bending of magnetic fields associated with filamentary molecular clouds. We compared the orientation and direction of this field bending with our 3D magnetic-field view of the Perseus cloud. Methods. We used previous line-of-sight and plane-of-sky magnetic field observations as well as Galactic magnetic field models to reconstruct the complete 3D magnetic field vectors and morphology associated with the Perseus cloud. Results. We approximated the 3D magnetic field morphology of the cloud as a concave arc that points in the decreasing longitude direction in the plane of the sky (from our point of view). This field morphology preserves a memory of the Galactic magnetic field. In order to compare this morphology to cloud-formation model predictions, we assume that the cloud retains a memory of its most recent interaction. After incorporating velocity observations, we find that the line-of-sight magnetic field observations are consistent with predictions of shock-cloud-interaction models. Conclusions. To our knowledge, this is the first time that the 3D magnetic fields of a molecular cloud have been reconstructed. We find the 3D magnetic field morphology of the Perseus cloud to be consistent with the predictions of the shock-cloud-interaction model that describes the formation mechanism of filamentary molecular clouds.

Extragalactic magnetic field constraints from ultrahigh-energy cosmic rays from local galaxies

ABSTRACT We interpret the correlation between local star-forming galaxy positions and ultrahigh-energy cosmic-ray (UHECR) directions, recently detected by the Pierre Auger Observatory (PAO), in terms of physical parameters: the local density of sources and the magnetic fields governing the UHECR propagation. We include a Galactic magnetic field model on top of a random extragalactic magnetic field description to determine the level of UHECR deflections expected from an ensemble of source positions. Besides deflections in magnetic fields, we also take into account energy losses with background photon fields as well as spectrum and composition measurements by the PAO. We find consistency between the PAO anisotropy measurement and the local star-forming galaxy density for large extragalactic magnetic field strengths with $B \gt 0.2 \ \rm nG$ (for a coherence length of $1 \ \rm Mpc$) at the 5σ confidence level. Larger source densities lead to more isotropic background and consequently allow for weaker extragalactic magnetic fields. However, the acceleration of UHECR by such abundant sources is more challenging to motivate. Too large source densities and extragalactic magnetic field strengths, on the other hand, are also disfavoured as that decreases the expected level of anisotropy. This leads to upper limits of $B \lt 22 \ \rm nG$ and $\rho _0 \lt 8.4 \times 10^{-2} \ \rm Mpc^{-3}$ at the 90 per cent confidence level.

The Imprint of Large-scale Structure on the Ultrahigh-energy Cosmic-Ray Sky

Abstract Ultrahigh-energy cosmic rays (UHECRs) are atomic nuclei from space with vastly higher energies than any other particles ever observed. Their origin and chemical composition remain a mystery. As we show here, the large and intermediate angular scale anisotropies observed by the Pierre Auger Observatory are a powerful tool for understanding the origin of UHECRs. Without specifying any particular production mechanism but only postulating that the source distribution follows the matter distribution of the local universe, a good accounting of the magnitude, direction, and energy dependence of the dipole anisotropy at energies above 8 × 1018 eV is obtained after taking into account the impact of energy losses during propagation (the “GZK horizon”), diffusion in the extragalactic magnetic field, and deflections in the Galactic magnetic field (GMF). This is a major step toward the long-standing hope of using UHECR anisotropies to constrain UHECR composition and magnetic fields. The observed dipole anisotropy is incompatible with a pure proton composition in this scenario. With a more accurate treatment of energy losses, it should be possible to further constrain the cosmic-ray composition and properties of the extragalactic magnetic field, self-consistently improve the GMF model, and potentially expose individual UHECR sources.

X-Ray Searches for Axions from Super Star Clusters

Axions may be produced in abundance inside stellar cores and then convert into observable x rays in the Galactic magnetic fields. We focus on the Quintuplet and Westerlund 1 super star clusters, which host large numbers of hot, young stars including Wolf-Rayet stars; these stars produce axions efficiently through the axion-photon coupling. We use Galactic magnetic field models to calculate the expected x-ray flux locally from axions emitted from these clusters. We then combine the axion model predictions with archival Nuclear Spectroscopic Telescope Array (NuSTAR) data from 10-80keV to search for evidence of axions. We find no significant evidence for axions and constrain the axion-photon coupling g_{aγγ}≲3.6×10^{-12} GeV^{-1} for masses m_{a}≲5×10^{-11} eV at 95% confidence.

A novel analytical model of the magnetic field configuration in the Galactic center

Context. Cosmic-ray propagation is strongly dependent on the large-scale configuration of the Galactic magnetic field. In particular, the Galactic center region provides highly interesting cosmic-ray data from gamma-ray maps and it is clear that a large fraction of the cosmic rays detected at Earth originate in this region of the Galaxy. Yet because of confusion from line-of-sight integration, the magnetic field structure in the Galactic center is not well known and no large-scale magnetic field model exists at present. Aims. In this paper, we develop a magnetic field model, derived from observational data on the diffuse gas, nonthermal radio filaments, and molecular clouds. Methods. We derive an analytical description of the magnetic field structure in the central molecular zone by combining observational data with the theoretical modeling of the basic properties of magnetic fields. Results. We provide a first description of the large-scale magnetic field in the Galactic center region. We present first test simulations of cosmic-ray propagation and the impact of the magnetic field structure on the cosmic-ray distribution in the three dimensions. Conclusions. Our magnetic field model is able to describe the main features of polarization maps; it is particularly important to note that they are significantly better than standard global Galactic magnetic field models. It can also be used to model cosmic-ray propagation in the Galactic center region more accurately.

Galactic cosmic ray propagation models using Picard

We will give an overview of recent developments in numerical modelling of CR transport in our Galaxy with the PICARD code. Corresponding numerical models aim at reproducing CR spectra and also diffuse gamma-ray emission from the Galaxy by solving the cosmic-ray transport problem in a three-dimensional model of the Galaxy. We address the impact of different transport physics processes on the flux and distribution of diffuse Galactic gamma rays: we investigate CR transport using different cosmic-ray source distribution models ranging from smooth source distributions to localised sources related to high-energy gamma-ray observations. Additionally, we consider anisotropic cosmic-ray diffusion in an improved Galactic magnetic field model. The choice of changing the different transport parameters is most readily visible in the inverse-Compton channel, which shows features not present in commonly-used axisymmetric transport models.

Dust polarization modelling at large scale over the northern Galactic cap using EBHIS and Planck data

The primary source of systematic uncertainty in the quest for the B-mode polarization of the Cosmic Microwave Background (CMB) introduced by primordial gravitational waves is polarized thermal emission from Galactic dust. Therefore, accurate characterization and separation of the polarized thermal dust emission is an essential step in distinguishing such a faint CMB B-mode signal. We provide a modelling framework to simulate polarized thermal dust emission based on the model described in Ghosh et al. (2017, A&A, 601, A71), making use of both the Planck dust and Effelsberg-Bonn HI surveys over the northern Galactic cap. Our seven-parameter dust model, incorporating both HI gas in three different column density templates as a proxy for spatially variable dust intensity and a phenomenological model of Galactic magnetic field, is able to reproduce both one- and two-point statistics of the observed dust polarization maps seen by Planck at 353 GHz over a selected low-column density region in the northern Galactic cap. This work has important applications in assessing the accuracy of component separation methods and in quantifying the confidence level of separating polarized Galactic emission and the CMB B-mode signal, as is needed for ongoing and future CMB missions.

The complex large-scale magnetic fields in the first Galactic quadrant as revealed by the Faraday depth profile disparity

ABSTRACT The Milky Way is one of the very few spiral galaxies known to host large-scale magnetic field reversals. The existence of the field reversal in the first Galactic quadrant near the Sagittarius spiral arm has been well established, yet poorly characterized due to the insufficient number of reliable Faraday depths (FDs) from extragalactic radio sources (EGSs) through this reversal region. We have therefore performed broad-band (1–$2\, {\rm GHz}$) spectropolarimetric observations with the Karl G. Jansky Very Large Array (VLA) to determine the FD values of 194 EGSs in the Galactic longitude range of 20°–52° within ±5° from the Galactic mid-plane, covering the Sagittarius arm tangent. This factor of five increase in the EGS FD density has led to the discovery of a disparity in FD values across the Galactic mid-plane in the Galactic longitude range of 40°–52°. Combined with existing pulsar FD measurements, we suggest that the Sagittarius arm can host an odd-parity disc field. We further compared our newly derived EGS FDs with the predictions of three major Galactic magnetic field models, and concluded that none of them can adequately reproduce our observational results. This has led to our development of new, improved models of the Milky Way disc magnetic field that will serve as an important step towards major future improvements in Galactic magnetic field models.

A more detailed look at Galactic magnetic field models: using free–free absorption in HII regions

Context. Cosmic rays (CRs) and the Galactic magnetic field (GMF) are fundamental actors in many processes in the Milky Way. The observed interaction product of these actors is Galactic synchrotron emission integrated over the line of sight (LOS). A comparison to simulations can be made with this tracer using existing GMF models and CR density models. This probes the GMF strength and morphology and the CR density. Aims. Our aim is to provide insight into the Galactic CR density and the distribution and morphology of the GMF strength by exploring and explaining the differences between the simulations and observations of synchrotron intensity. Methods. At low radio frequencies HII regions become opaque due to free–free absorption. Using these HII regions we can measure the synchrotron intensity over a part of the LOS through the Galaxy. The measured intensity per unit path length, that is, the emissivity, for HII regions at different distances, allows us to probe the variation in synchrotron emission not only across the sky but also in the third dimension of distance. Performing these measurements on a large scale is one of the new applications of the window opened by current low-frequency arrays. Using a number of existing GMF models in conjunction with the Galactic CR modeling code GALPROP, we can simulate these synchrotron emissivities. Results. We present an updated catalog, compiled from the literature, of low-frequency absorption measurements of HII regions, their distances, and electron temperatures. We report a simulated emissivity that shows a compatible trend for HII regions that are near the observer. However, we observe a systematically increasing synchrotron emissivity for HII regions that are far from the observer, which is not compatible with the values simulated by the GMF models and GALPROP. Conclusions. Current GMF models plus a GALPROP generated CR density model cannot explain low-frequency absorption measurements. One possibility is that distances to all HII regions catalogued at the kinematic “far” distance are erroneously determined, although this is unlikely since it ignores all evidence for far distances in the literature. However, a detection bias due to the nature of this tracer requires us to keep in mind that certain sources may be missed in an observation. The other possibilities are an enhanced emissivity in the outer Galaxy or a diminished emissivity in the inner Galaxy.

Modulations in Spectra of Galactic Gamma-ray sources as a result of photon-ALPs mixing.

Axion like particles (ALPs) are fundamental pseudo scalar particles with properties similar to Axions which are a well-known extensions of the standard model to solve the strong CP problem in Quantum Chromodynamics. ALPs can oscillate into photons and vice versa in the presence of an external tranversal magnetic field. This oscillation of photon and ALPs could have important implications for astronomical observations, i.e. a characteristic energy dependent attenuation in Gamma ray spectra for astrophysical sources. Here we have revisited the opportunity to search Photon-ALPs coupling in the disappearance channel. We use eight years of Fermi Pass 8 data of a selection of promising galactic Gamma-ray source candidates and study the modulation in the spectra in accordance with Photon-ALPs mixing and estimate best fit values of the parameters i.e. Photon-ALPs coupling constant (gʱγγ) and ALPs mass(mα). For the magnetic field we use large scale galactic magnetic field models based on Faraday rotation measurements and we have also studied the survival probability of photons in the Galactic plane.

Constraints on axion-like particle properties with TeV gamma-ray observations of Galactic sources

Axion-like particles (ALPs) can oscillate into photons and vice versa in the magnetic fields, which provides a feasible explanation of the observed transparency of the Universe for TeV gamma rays. The coupling gaγ for the ALP mass ma⩽ 10−7 eV has been tightly constrained by the Fermi-LAT and H.E.S.S. observations of some extragalactic sources. In this work, we search for the ALP oscillation effect with the H.E.S.S. observation data of some TeV sources in the Galactic plane for the first time. No significant ALP signal is found. The constraints on the ALP parameters, though sensitively depend on the Galactic magnetic field (GMF) model, are more stringent than the limits set by the CAST experiment at ma ∼ 0.1μ eV. The upcoming CTA observations of these Galactic TeV sources are able to improve the constraints substantially.

Solenoidal Improvements for the JF12 Galactic Magnetic Field Model

Abstract The popular JF12 analytic model by Jansson & Farrar provides a quantitative description of the Galaxy’s large-scale magnetic field, which is widely used in various astrophysical applications. However, both the poloidal X-type component and the spiral disk component of JF12 exhibit regions in which the magnetic divergence constraint is violated. We first propose a cure for this problem, resulting in a truly solenoidal large-scale spiral field. Second, the otherwise straight field lines of the X-type component exhibit kinks in the Galactic plane that, in addition to implying the presence of a singular current sheet, may pose difficulties for e.g., numerical tracing of cosmic-ray particles. We propose and discuss two possible strategies to mitigate this problem. Although all corrections are kept as minimal as possible, the extended set of model parameters will have to be carefully readjusted in order to fully restore the agreement to observational data that the unmodified JF12 field is based on. Furthermore, the performance of our improved version of the field model is quantitatively assessed by test simulations using the CRPropa Galactic cosmic-ray propagation code.

Arrival Directions of Ultrahigh-Energy Cosmic Rays

We discuss some important properties of the arrival directions of ultrahigh-energy cosmic rays detected by the Telescope Array and the Pierre Auger experiments using data released by both collaborations. We also discuss the potential sources and analyse the effect of correcting the arrival directions to take into account the deflections of the cosmic rays by the magnetic field of our galaxy for one specific model of the galactic magnetic field under several assumptions about the composition of the primary cosmic rays. The results show that for a proton-dominated composition the maxima deflections, due to the galactic magnetic field model considered in this paper, are about 20° for energies higher the 60 EeV, allowing the approximation of doing astronomy with charged particles, while for medium and heavy compositions the galactic magnetic deflections are much higher and cannot be excluded from any search of the sources of ultrahigh-energy cosmic rays. We also study the validity of the backpropagation method by means of forward propagation of the particles from potential astrophysical sources to the vicinity of the Earth.

Faraday Rotation of Extended Emission as a Probe of the Large-Scale Galactic Magnetic Field

The Galactic magnetic field is an integral constituent of the interstellar medium (ISM), and knowledge of its structure is crucial to understanding Galactic dynamics. The Rotation Measures (RM) of extragalactic (EG) sources have been the basis of comprehensive Galactic magnetic field models. Polarised extended emission (XE) is also seen along lines of sight through the Galactic disk, and also displays the effects of Faraday rotation. Our aim is to investigate and understand the relationship between EG and XE RMs near the Galactic plane, and to determine how the XE RMs, a hitherto unused resource, can be used as a probe of the large-scale Galactic magnetic field. We used polarisation data from the Canadian Galactic Plane Survey (CGPS), observed near 1420 MHz with the Dominion Radio Astrophysical Observatory (DRAO) Synthesis Telescope. We calculated RMs from a linear fit to the polarisation angles as a function of wavelength squared in four frequency channels, for both the EG sources and the XE. Across the CGPS area, 55 ∘ < ℓ < 193 ∘ , − 3 ∘ < b < 5 ∘ , the RMs of the XE closely track the RMs of the EG sources, with XE RMs about half the value of EG-source RMs. The exceptions are places where large local HII complexes heavily depolarise more distant emission. We conclude that there is valuable information in the XE RM dataset. The factor of 2 between the two types of RM values is close to that expected from a Burn slab model of the ISM. This result indicates that, at least in the outer Galaxy, the EG and XE sources are likely probing similar depths, and that the Faraday rotating medium and the synchrotron emitting medium have similar variation with galactocentric distance.