Intergalactic Magnetic Fields Research Articles

Prospects for detection of the pair-echo emission from TeV gamma-ray bursts

The intergalactic magnetic field (IGMF) present in the voids of large-scale structures is considered to be the weakest magnetic field in the Universe. Gamma-ray observations of blazars in the GeV-TeV domain have led to lower limits on the IGMF strength based on the search for delayed or extended emission. Nevertheless, these results have been obtained with strong assumptions placed on the unknown source properties. The recent discovery of TeV radiation from gamma-ray bursts (GRBs) has paved the way for IGMF studies with these bright transients. Among the current TeV-detected GRBs, GRB\,190114C, located at a redshift of $z = 0.42$, is the best sampled. Therefore, it can be considered to be representative of the properties of GRBs in the VHE domain. In addition GRB\,221009A ($z = 0.151$) is the brightest event ever detected. We present a phenomenological model based on the intrinsic properties of GRB\,190114C and GRB\,221009A to predict the delayed emission component (pair-echo) in the GeV-TeV band. We investigate the detectability of this component from low-redshift ($z 1$) GRBs for three values of IGMF strength ($10^ $ G, $10^ $ G and $10^ $ G) different observational times ($3$ hrs, $6$ hrs, and $9$ hrs) and source intrinsic properties. We find that for current and future generation gamma -ray instruments extending the observation for at least 3 hours after the GRB detection is a viable strategy for probing the IGMF. We also confirm that GeV-TeV observations of GRBs can probe IGMF strengths on the order of $10^ $ G, representing a competitive alternative to the current studies performed with active galactic nuclei (AGNs)

LOFAR Deep Fields: Probing the sub-mJy regime of polarized extragalactic sources in ELAIS-N1. I. The catalog

Quantifying the number density and physical characteristics of extragalactic polarized sources is important for the successful planning of future studies based on Faraday rotation measure (RM) grids of polarized sources to probe foreground Galactic and intergalactic magnetic fields. However, it is proving very hard to detect polarized signal from the population of very faint (sub-mJy) polarized sources at low radio frequencies, and their properties are mostly unknown. LOFAR can play an important role in such studies thanks to its sensitivity and angular resolution, combined with the precision on the inferred RM values that can be achieved through low-frequency broad-band polarimetry. The aim of this study is to probe the sub-mJy polarized source population with LOFAR. In this first paper, we present the method used to stack LOFAR polarization datasets, the resulting catalog of polarized sources, and the derived polarized source counts. The European Large Area ISO Survey-North 1 (ELAIS-N1) field, one of the deepest of the LOFAR Two-Metre Sky Survey (LoTSS) Deep Fields so far, was selected for a polarimetric study at 114.9–177.4 MHz. A total area of 25 deg$^2$ was imaged at 6"-resolution in the Stokes $Q$ and $U$ parameters. Alignment of polarization angles was done both in frequency and in Faraday space before stacking datasets from 19 eight-hour-long epochs taken in two different LOFAR observing cycles. A search for polarized sources was carried out in the final, stacked dataset, and the properties of the detected sources were examined. The depolarization level of sources known to be polarized at 1.4 GHz was quantified. A one-sigma noise level, $ QU $, of 19 mu Jy/beam was reached in the central part of the field after stacking. Twenty-five polarized sources were detected above 8$ QU $, five of which had not been detected in polarization at any other radio frequencies before. Seven additional polarized components were found by lowering the threshold to 6$ QU $ at positions corresponding to sources known to be polarized at 1.4 GHz. In two radio galaxies, polarization was detected from both radio lobes, so the final number of associated radio continuum sources is 31. The detected sources are weakly polarized, with a median degree of polarization of 1.75<!PCT!> for the sample of sources detected in polarized emission. For the 10 polarized sources previously identified in a pilot LOFAR study of the ELAIS-N1 field at 20"-resolution, the RM values are consistent but the degrees of polarization are higher in the 6"-resolution data. The sources previously detected in polarization at 1.4 GHz are significantly depolarized at 150 MHz. The catalog is used to derive the polarized source counts at 150 MHz. This is the deepest and highest-resolution polarization study at 150 MHz to date. A full characterization of the sources and an analysis of the catalog will be presented in Paper II.

A delayed 400 GeV photon from GRB 221009A and implication on the intergalactic magnetic field

Large High Altitude Air Shower Observatory has detected 0.2 − 13 TeV emission of GRB 221009A within 2000 s since the trigger. Here we report the detection of a 400 GeV photon, without accompanying prominent low-energy emission, by Fermi Large Area Telescope in this direction with a 0.4 days’ delay. Given an intergalactic magnetic field strength of about 4 × 10−17 G, which is comparable to limits from TeV blazars, the delayed 400 GeV photon can be explained as the cascade emission of about 10 TeV gamma rays. We estimate the probabilities of the cascade emission that can result in one detectable photon beyond 100 GeV by Fermi Large Area Telescope within 0.3 − 1 days is about 2% whereas it is about 20.5% within 0.3 − 250 days. Our results show that Synchrotron Self-Compton explanation is less favored with probabilities lower by a factor of about 3 − 30 than the cascade scenario.

Constraining the Astrophysical Origin of Intergalactic Magnetic Fields

High-energy photons can produce electron–positron pairs upon interacting with the extragalactic background light. These pairs will in turn be deflected by the intergalactic magnetic field (IGMF), before possibly up-scattering photons of the cosmic microwave background, thereby initiating an electromagnetic cascade. The nonobservation of an excess of GeV photons and an extended halo around individual blazars due to this electromagnetic cascade can be used to constrain the properties of the IGMF. In this work, we use publicly available data of 1ES 0229+200 obtained with the Fermi Large Area Telescope and the High Energy Stereoscopic System to constrain cosmological MHD simulations of various magnetogenesis scenarios, and find that all models without a strong space-filling primordial component or overoptimistic dynamo amplifications can be excluded at the 95% confidence level. In fact, we find that the fraction of space filled by a strong IGMF has to be at least f ≳ 0.67, thus excluding most astrophysical production scenarios. Moreover, we set lower limits of B 0 > 5.1 × 10−15 G (B 0 > 1.0 × 10−14 G) for a space-filling primordial IGMF for a blazar activity time of Δt = 104 yr (Δt = 107 yr).

Nonlinear Feedback of the Electrostatic Instability on the Blazar-induced Pair Beam and GeV Cascade

Relativistic pair beams produced in the cosmic voids by TeV gamma-rays from blazars are expected to produce a detectable GeV-scale cascade that is missing in the observations. The suppression of this secondary cascade implies either the deflection of the pair beam by intergalactic magnetic fields or, alternatively, an energy loss of the beam due to the beam-plasma instability. Here, we study how the beam-plasma instability feeds back on the beam, using a realistic two-dimensional beam distribution. We find that the instability broadens the beam opening angles significantly without any significant energy loss, thus confirming a recent feedback study on a simplified one-dimensional beam distribution. However, narrowing diffusion feedback of the beam particles with Lorentz factors less than 106 might become relevant, even though initially it is negligible. Finally, when considering the continuous creation of TeV pairs, we find that the beam distribution and the wave spectrum reach a new quasi-steady state, in which the scattering of beam particles persists and the beam opening angle may increase by a factor of hundreds. Understanding the implications on the GeV cascade emission requires accounting for inverse-Compton cooling.

A Detectable Ultra-high-energy Cosmic-Ray Outburst from GRB 221009A

Gamma-ray bursts (GRBs) have been proposed as one of the promising sources of ultra-high-energy cosmic rays (UHECRs), but observational evidence is still lacking. The nearby brightest of all time GRB 221009A, a once in 1000 yr event, is able to accelerate protons to ∼103 EeV, and then generate extremly energetic neutrons via the photomeson production interaction. Protons arriving at the Milky Way are dominated by neutron-decay-induced protons. The intergalactic magnetic fields would not yield a sizable delay of the ≥10 EeV cosmic rays if its strength is ≲10−13 G, while Galactic magnetic fields would cause a significant time delay. We predict that a UHECR burst from GRB 221009A would be detectable by the Pierre Auger Observatory and the TA×4, within ∼10 yr. The detection of such a UHECR outburst will provide the direct evidence for UHECR acceleration in GRBs.

Magnetic moments of astrophysical neutrinos

We study the impact of neutrino magnetic moments on astrophysical neutrinos, in particular supernova neutrinos and ultra-high energy neutrinos from extragalactic sources. We show that magnetic moment-induced conversion of Dirac neutrinos from left-handed states into unobservable right-handed singlet states can substantially change the flux and flavour composition of these neutrinos at Earth. Notably, neutrinos from a supernova's neutronisation burst, whose flux can be predicted with 𝒪(10%) accuracy, offer a discovery reach to neutrino magnetic moments ∼ few × 10-13 μB , up to one order of magnitude below current limits. For high-energy neutrinos from distant sources, for which no robust flux prediction exists, we show how the flavour composition at Earth can be used as a handle to establish the presence of non-negligible magnetic moments, potentially down to few× 10-17 μB if the measurement can be performed on neutrinos from a single source. In both cases, the sensitivity strongly depends on the galactic (intergalactic) magnetic field profiles along the line of sight. Therefore, while a discovery is possible down to very small values of the magnetic moment, the absence of a discovery does not imply an equally strong limit. We also comment on the dependence of our results on the right-handed neutrino mass, paying special attention to the transition from coherent deflection by a classical magnetic field to incoherent scattering on individual scattering targets. Finally, we show that a measurement of Standard Model Dirac neutrino magnetic moments, of order 10-19 μB , could be possible under rather optimistic, but not completely outrageous, assumptions using flavour ratios of high-energy astrophysical neutrinos.

Constraints on the intergalactic magnetic field from Fermi/LAT observations of the ‘pair echo’ of GRB 221009A

Delayed ‘pair-echo’ signal from interactions of very-high-energy γ rays in the intergalactic medium can be used for the detection of the intergalactic magnetic field (IGMF). We used the data of the Fermi/LAT telescope coupled with LHAASO observatory measurements to confirm the presence of IGMF along the line of sight to the γ-ray burst GRB 221009A. Comparing the Fermi/LAT measurements with the expected level of the pair-echo flux, set by the multi-TeV LHAASO detection, we derived a lower bound 10−19 G on the IGMF with correlation length l larger than 1 Mpc, improving as l−1/2 for shorter correlation lengths. This provides an independent verification of the existence of a lower bound on IGMF in the voids of the large-scale structure, previously derived from the observations of active galactic nuclei.

Graviton to Photon Conversion in Curved Space-Time and External Magnetic Field

The suppression of relic gravitational waves due to their conversion into electromagnetic radiation in a cosmological magnetic field is studied. The coupled system of equations describing gravitational and electromagnetic wave propagation in an arbitrary curved space-time and in external magnetic field is derived. The subsequent elimination of photons from the beam due to their interaction with the primary plasma is taken into account. The resulting system of equations is solved numerically in the Friedman–LeMaitre–Robertson–Walker metric for the upper limit of the intergalactic magnetic field strength of 1 nGs. We conclude that the gravitational wave conversion into photons in the intergalactic magnetic field cannot significantly change the amplitude of the relic gravitational wave and their frequency spectrum.

Limit on the intergalactic magnetic field from the ultrahigh-energy cosmic ray hotspot in the Perseus-Pisces region

Limit on the intergalactic magnetic field from the ultrahigh-energy cosmic ray hotspot in the Perseus-Pisces region

Strong gravitational lensing of blazar gamma-radiation and intergalactic magnetic fields

The influence of intergalactic magnetic fields on the strong gravitational lensing of blazar secondary gamma radiation is discussed. Currently, two cases of strong gravitational lensing of blazar gamma-radiation are known, where radiation is deflected by galaxies on the line of sight between the blazars and Earth. The magnetic field can affect the movements of electron-positron pairs generated by primary radiation, thereby changing the directions of secondary gamma radiation. It modifies the equation of the gravitational lens and leads to the dependence of the observed signal in the secondary gamma radiation on the energy of photons and magnetic field. Accordingly, it is possible, in principle, to estimate the intergalactic magnetic fields from the time delay of signals, from the angular position of images (for future high-resolution gamma-ray telescopes) or from the shape of the observed energy spectrum. This method is demonstrated by the example of the blazar B0218+357. In this case, however, it is not possible to obtain useful constraints due to the large distances to the blazar and lens galaxy. The result is only a lower limit on the magnetic field G, which is weaker than other existing constraints. However, future discoveries of lensed blazars may provide more favorable opportunities for measuring the magnetic fields, especially with the help of a new generation of gamma-ray telescopes such as e-ASTROGAM, GECAM, and SVOM as well as future gamma-ray telescopes with a high angular resolution, .

Sampling the Faraday rotation sky of TNG50: imprint of the magnetized circumgalactic medium around Milky Way-like galaxies

ABSTRACT Faraday rotation measure (RM) is arguably the most practical observational tracer of magnetic fields in the diffuse circumgalactic medium (CGM). We sample synthetic Faraday rotation skies of Milky Way-like galaxies in TNG50 of the IllustrisTNG project by placing an observer inside the galaxies at a solar circle-like position. Our synthetic RM grids emulate specifications of current and upcoming surveys; the NRAO VLA Sky Survey (NVSS), the Polarisation Sky Survey of the Universe’s Magnetism (POSSUM), and a future Square Kilometre Array (SKA1-mid) polarization survey. It has been suggested that magnetic fields regulate the survival of high-velocity clouds. However, there is only a small number of observational detections of magnetized clouds thus far. In the first part of the paper, we test conditions for the detection of magnetized circumgalactic clouds. Based on the synthetic RM samplings of clouds in the simulations, we predict upcoming polarimetric surveys will open opportunities for the detection of even low-mass and distant clouds. In the second part of the paper, we investigate the imprint of the CGM in the all-sky RM distribution. We test whether the RM variation produced by the CGM is correlated with global galaxy properties, such as distance to a satellite, specific star formation rate, neutral hydrogen covering fraction, and accretion rate to the supermassive black hole. We argue that the observed fluctuation in the RM measurements on scales less than 1○, which has been considered an indication of intergalactic magnetic fields, might in fact incorporate a significant contribution of the Milky Way CGM.

Constraints on the Intergalactic Magnetic Field Strength from γ-Ray Observations of GRB 221009A

Characteristics of the cascade gamma-ray signal resulting from very-high-energy gamma-ray sources, such as gamma-ray bursts, can be used to constrain the strength and structure of intergalactic magnetic fields (IGMFs). There has been a debate on whether GRB 190114C, the first gamma-ray burst with observed TeV photons, can constrain the IGMF. Recently, LHAASO detected the brightest-of-all-time GRB 221009A, which has much larger energy in the TeV band, and the spectrum extends to energy above 10 TeV, providing an unprecedented opportunity to study IGMF. We perform a Monte Carlo simulation of the cascade process with the public ELMAG code, considering the TeV data of GRB 221009A observed by LHAASO. By comparing the resulting cascade emission with the flux limit obtained from Fermi-LAT observations, we infer a limit of B ≥ 10−18.5 G for IGMF. Though this limit may not be as strong as the limit from blazars, it serves as an independent constraint on IGMF from a new class of TeV sources.

Modeling the propagation of very-high-energy γ-rays with the CRbeam code: Comparison with CRPropa and ELMAG codes

Context. Very-high-energy γ-rays produce electron positron pairs in interactions with low-energy photons of extragalactic background light during propagation through the intergalactic medium. The electron-positron pairs generate secondary γ-rays detectable by γ-ray telescopes. This secondary emission can be used to detect intergalactic magnetic fields (IGMF) in the voids of large-scale structure. Aims. A new γ-ray observatory, namely, Cherenkov Telescope Array (CTA), will provide an increase in sensitivity for detections of these secondary γ-ray emission and enable the measurement of its properties for sources at cosmological distances. The interpretation of the CTA data, including the detection of IGMF and study of its properties and origins, will require precision modeling of the primary and secondary γ-ray fluxes. Methods. We assess the precision of the modeling of the secondary γ-ray emission using model calculations with publicly available Monte Carlo codes CRPropa and ELMAG and compare their predictions with theoretical expectations and with model calculations of a newly developed CRbeam code. Results. We find that model predictions of different codes differ by up to 50% for low-redshift sources, with discrepancies increasing up to order-of-magnitude level with the increasing source redshifts. We identify the origin of these discrepancies and demonstrate that after eliminating the inaccuracies found, the discrepancies between the three codes are reduced to 10% when modeling nearby sources with z ~ 0.1. We argue that the new CRbeam code provides reliable predictions for the spectral, timing, and imaging properties of the secondary γ-ray signal for both nearby and distant sources with z ~ 1. Thus, it can be used to study gamma-ray sources and IGMF with a level of precision that is appropriate for the prospective CTA study of the effects of γ-ray propagation through the intergalactic medium.

Constraints on the Intergalactic Magnetic Field Using Fermi-LAT and H.E.S.S. Blazar Observations

Magnetic fields in galaxies and galaxy clusters are believed to be the result of the amplification of intergalactic seed fields during the formation of large-scale structures in the universe. However, the origin, strength, and morphology of this intergalactic magnetic field (IGMF) remain unknown. Lower limits on (or indirect detection of) the IGMF can be obtained from observations of high-energy gamma rays from distant blazars. Gamma rays interact with the extragalactic background light to produce electron−positron pairs, which can subsequently initiate electromagnetic cascades. The gamma-ray signature of the cascade depends on the IGMF since it deflects the pairs. Here we report on a new search for this cascade emission using a combined data set from the Fermi Large Area Telescope and the High Energy Stereoscopic System. Using state-of-the-art Monte Carlo predictions for the cascade signal, our results place a lower limit on the IGMF of B > 7.1 × 10−16 G for a coherence length of 1 Mpc even when blazar duty cycles as short as 10 yr are assumed. This improves on previous lower limits by a factor of 2. For longer duty cycles of 104 (107) yr, IGMF strengths below 1.8 × 10−14 G (3.9 × 10−14 G) are excluded, which rules out specific models for IGMF generation in the early universe.

Effects of large-scale magnetic fields on the observed composition of ultrahigh-energy cosmic rays

Ultra high-energy (UHE) cosmic rays (CRs) from distant sources interact with intergalactic radiation fields, leading to their spallation and attenuation. They are also deflected in intergalactic magnetic fields (IGMFs), particularly those associated with Mpc-scale structures. These deflections extend the propagation times of CR particles, forming a magnetic horizon for each CR species. The cumulative cooling and interactions of a CR ensemble also modifies their spectral shape and composition observed on Earth. We construct a transport formulation to calculate the observed UHE CR spectral composition for 4 classes of source population. The effects on CR propagation brought about by IGMFs are modeled as scattering processes during transport, by centers associated with cosmic filaments. Our calculations demonstrate that IGMFs can have a marked effect on observed UHE CRs, and that source population models are degenerate with IGMF properties. Interpretation of observations, including the endorsement or rejection of any particular source classes, thus needs careful consideration of the structural properties and evolution of IGMFs. Future observations providing tighter constraints on IGMF properties will significantly improve confidence in assessing UHE CR sources and their intrinsic CR production properties.

Intergalactic magnetic field studies by means of γ -ray emission from GRB 190114C

The presence of delayed GeV emission after a strong transient, such as a GRB (Gamma-Ray Burst), in the VHE (Very-High Energy, $E>100$ GeV) band can be the signature of a non-zero magnetic field in the intergalactic medium. We used a synchrotron self-Compton multiwavelength model to infer an analytical description of the intrinsic VHE spectrum (corrected for absorption by the Extragalactic Background Light, EBL) of GRB$\,$190114C to predict the lightcurves and SEDs of the delayed emission with Monte Carlo simulations for different IGMF (Intergalactic Magnetic Field) configurations (strengths $B=8\times10^{-21}$ G, $10^{-20}$ G, $3\times 10^{-20}$G and correlation length $\lambda>1$ Mpc), and compared them with the Fermi-LAT (Fermi Large Area Telescope) limits computed for several exposure times. We found that Fermi LAT is not sensitive enough to constrain any IGMF strengths using GRB$\,$190114C.

Cosmic neutrino flux and spin flavor oscillations in intergalactic medium

The ultra high energy (UHE) cosmic neutrinos are expected to play a pivotal role in the disquisition of physics beyond the standard model of particle physics as well as serve as an ideal cosmic messengers. This epitomizes the selling point of several currently running or planned neutrino telescopes. The UHE cosmic neutrinos usually perambulate gargantuan scales in the extragalactic universe having a magnetic field. If neutrinos have a finite magnetic moment (μν) owing to quantum loop corrections, this may result in spin-flavor oscillations, which can affect the cosmic neutrino flux. Using the current limit and assuming neutrinos to be Dirac particles, we show that the flux of cosmic neutrinos will reduce by half if they traverse few Mpcs through the intergalactic magnetic field, in the range of μG to nG. Moreover, one can safely neglect the effect of μν if the current upper bound is improved by a few orders of magnitude even if the neutrinos travel through the size of the visible universe.

Search of the pair echo signatures in the high-energy light curve of GRB190114C

A model of the time delayed electromagnetic cascade "echo" is applied to the bright gamma-ray burst GRB190114C - the first gamma-ray burst to be contemporaneously detected in high and very high energy gamma-ray bands. It is shown that the internal spread of the cascade in the absence of the intervening magnetic fields dilutes the "echo" emission over $10^3-10^5$ seconds depending on the energy. Accounting for the measured source flux in the $0.3-1$ TeV gamma-ray band, the prediction of the "echo" model is shown to match the detected lower-energy gamma-ray emission $10^4$ seconds after the burst. However, the "echo" emission remains indistinguishable from the intrinsic GRB190114C flux within the measurement uncertainties. Implications of this in the context of the intergalactic magnetic field measurement are discussed.

A lower bound on intergalactic magnetic fields from time variability of 1ES 0229+200 from MAGIC and Fermi/LAT observations

Context. Extended and delayed emission around distant TeV sources induced by the effects of propagation of γ ray s through the intergalactic medium can be used for the measurement of the intergalactic magnetic field (IGMF). Aims. We search for delayed GeV emission from the hard-spectrum TeV γ-ray emitting blazar 1ES 0229+200, with the goal of detecting or constraining the IGMF-dependent secondary flux generated during the propagation of TeV γ rays through the intergalactic medium. Methods. We analysed the most recent MAGIC observations over a 5 year time span, and complemented them with historic data of the H.E.S.S. and VERITAS telescopes, along with a 12-year-long exposure of the Fermi/LAT telescope. We used them to trace source evolution in the GeV–TeV band over a decade and a half. We used Monte Carlo simulations to predict the delayed secondary γ-ray flux, modulated by the source variability, as revealed by TeV-band observations. We then compared these predictions for various assumed IGMF strengths to all available measurements of the γ-ray flux evolution. Results. We find that the source flux in the energy range above 200 GeV experiences variations around its average on the 14-year time span of observations. No evidence for the flux variability is found in the 1 − 100 GeV energy range accessible to Fermi/LAT. The non-detection of variability due to delayed emission from electromagnetic cascade developing in the intergalactic medium imposes a lower bound of B > 1.8 × 10−17 G for the long-correlation-length IGMF and B > 10−14 G for an IGMF of cosmological origin. Though weaker than the one previously derived from the analysis of Fermi/LAT data, this bound is more robust, being based on a conservative intrinsic source spectrum estimate and accounting for the details of source variability in the TeV energy band. We discuss implications of this bound for cosmological magnetic fields that might explain the baryon asymmetry of the Universe.