Relativistic jets from blazars are known to be sources of very-high-energy gamma rays (VHEGRs). During their propagation in the intergalactic space, VHEGRs interact with pervasive cosmological photon fields such as the extragalactic background light (EBL) and the cosmic microwave background (CMB), producing electron-positron pairs. These pairs can upscatter CMB/EBL photons to high energies via inverse Compton scattering, thereby continuing the cascade process. This is often used to set limits on intergalactic magnetic fields (IGMFs). However, the picture may change if plasma instabilities, arising due to the interaction of the pairs with the intergalactic medium (IGM), cool down the electrons/positrons faster than inverse Compton scattering. As a consequence, the kinetic energy lost by the pairs to the IGM could cause a hardening in the observed gamma-ray spectrum at energies below $\sim$100 GeV. Here we study several types and models of instabilities and assess their impact for interpreting observations of distant blazars. Our results suggest that plasma instabilities can describe the spectra of some blazars and mimic the effects of IGMFs, depending on parameters such as intrinsic spectrum of the object, the density and temperature of the IGM, and the luminosity of the beam. On the other hand, we find that for our fiducial set of parameters plasma instabilities do not have a major impact on the spectra of some of the blazars studied. Therefore, they may be used for constraining IGMFs.

Primordial magnetic fields (PMFs), being present before the epoch of cosmic recombination, induce small-scale baryonic density fluctuations. These inhomogeneities lead to an inhomogeneous recombination process that alters the peaks and heights of the large-scale anisotropies of the cosmic microwave backround (CMB) radiation. Utilizing numerical compressible MHD calculations and a Monte Carlo Markov chain analysis, which compares calculated CMB anisotropies with those observed by the \textit{WMAP} and \textit{Planck} satellites, we derive limits on the magnitude of putative PMFs. We find that the {\it total} remaining present day field, integrated over all scales, cannot exceed 47 pG for scale-invariant PMFs and 8.9 pG for PMFs with a violet Batchelor spectrum at 95\% confidence level. These limits are more than one order of magnitude more stringent than any prior stated limits on PMFs from the CMB which have not accounted for this effect.

The origin of the first magnetic fields in the Universe is a standing problem in cosmology. Intergalactic magnetic fields (IGMFs) may be an untapped window to the primeval Universe, providing further constrains on magnetogenesis. We demonstrate the feasibility of using ultra-high-energy cosmic rays (UHECRs) to constrain the helicity of IGMFs by performing simulations of cosmic-ray propagation in simple magnetic field configurations. We show that the first harmonic moments of the arrival distribution of UHECRs may be used to measure the absolute value of the helicity and its sign.

At TeV energies and above gamma rays can induce electromagnetic cascades, whose charged component is sensitive to intervening intergalactic magnetic fields (IGMFs). When interpreting gamma-ray measurements in the energy range between a few GeV and hundreds of TeV, one has to carefully account for effects due to IGMFs, which depend on their strength and power spectrum. Therefore, gamma-ray-induced electromagnetic cascades can be used as probes of cosmic magnetism, since their arrival distribution as well as spectral and temporal properties can provide unique information about IGMFs, whose origin and properties are currently poorly understood. In this contribution we present an efficient three-dimensional Monte Carlo code for simulations of gamma-ray propagation. We focus on the effects of different configurations of IGMFs, in particular magnetic helicity and the power spectrum of stochastic fields, on the morphology of the arrival directions of gamma rays, and discuss the prospects for detecting pair haloes around distant blazars.

We determine the effect of intergalactic magnetic fields on the distribution of high energy gamma rays by performing three-dimensional Monte Carlo simulations of the development of gamma-ray-induced electromagnetic cascades in the magnetized intergalactic medium. We employ the so-called 'Large Sphere Observer' method to efficiently simulate blazar gamma ray halos. We study magnetic fields with a Batchelor spectrum and with maximal left- and right-handed helicities. We also consider the case of sources whose jets are tilted with respect to the line of sight. We verify the formation of extended gamma ray halos around the source direction, and observe spiral-like patterns if the magnetic field is helical. We apply the $Q$-statistics to the simulated halos to extract their spiral nature and also propose an alternative method, the $S$-statistics. Both methods provide a quantative way to infer the helicity of the intervening magnetic fields from the morphology of individual blazar halos for magnetic field strengths $B \gtrsim 10^{-15}\,{\rm G}$ and magnetic coherence lengths $L_{\rm c} \gtrsim 100\,{\rm Mpc}$. We show that the $S$-statistics has a better performance than the $Q$-statistics when assessing magnetic helicity from the simulated halos.

The observation in the GeV band of distant blazars has been recently used to put constraints on the Extragalactic Background Light (EBL) and Extragalactic Magnetic Fields (EGMF). To support such claims one has to assume that the leptonic component of the electromagnetic cascade initiated by blazar gamma-rays is deflected away by strong enough EGMF, suppressing the signal in the Fermi window. Apart from magnetic fields, the development of such a cascade might be affected by plasma instabilities due to interactions with the ionized component of the Intergalactic Medium (IGM). In this paper we model the electromagnetic cascade through a Monte Carlo simulation in which both effects are taken into account separately, and we derive constraints on these scenarios from the combined Fermi-HESS data set. In the specific case of 1ES 0229+200 observations, we show that both explanations of the GeV flux suppression are compatible with the available data, specifically by assuming a magnetic field of $B \gtrsim 10^{-16}\,\rm{G}$ or an IGM temperature of $T \lesssim 5 \times 10^{4}\,\rm{K}$ along the line of sight. Future observations of the spectra of high redshift ($z\lesssim 1$) TeV objects will help to distinguish magnetic field and plasma effects on electromagnetic cascades in the IGM.

We extend our recent derivation of the time evolution equations for the energy content of magnetic fields and turbulent motions for incompressible, homogeneous, and isotropic turbulence to include the case of nonvanishing helicity. These equations are subsequently numerically integrated in order to predict the present day primordial magnetic field strength and correlation length, depending on its initial helicity and magnetic energy density. We find that all prior analytic predictions for helical magnetic fields, such as the epoch when they become maximally helical and their subsequent growth of correlation length L ~ a^{1/3} and decrease of magnetic field strength B ~ a^{-1/3} with scale factor a, are well confirmed by the simulations. An initially fully helical primordial magnetic field is a factor 4 10^4 stronger at the present epoch then its nonhelical counterpart when generated during the electroweak epoch.

We present a derivation of the time evolution equations for the energy content of nonhelical magnetic fields and the accompanying turbulent flows from first principles of incompressible magnetohydrodynamics in the general framework of homogeneous and isotropic turbulence. This is then applied to the early Universe, i.e., the evolution of primordial magnetic fields. Numerically integrating the equations, we find that most of the energy is concentrated at an integral wavenumber scale k_I where the turbulence turn over time equals the Hubble time. At larger length scales L, i.e., smaller wavenumbers q = 2 \pi / L << k_I, independent of the assumed turbulent flow power spectrum, mode-mode coupling tends to develop a small q magnetic field tail with a Batchelor spectrum proportional to the fourth inverse power of L and therefore a scaling for the magnetic field of B ~ L^(-5/2).

Motivated by experimental indications of a significant presence of heavy nuclei in the cosmic ray flux at ultra high energies ($\gtrsim 10^{19} \eV$), we consider the effects of Planck scale suppressed Lorentz Invariance Violation (LIV) on the propagation of cosmic ray nuclei. In particular we focus on LIV effects on the photodisintegration of nuclei onto the background radiation fields. After a general discussion of the behavior of the relevant quantities, we apply our formalism to a simplified model where the LIV parameters of the various nuclei are assumed to kinematically result from a single LIV parameter for the constituent nucleons, $\eta$, and we derive constraints on $\eta$. Assuming a nucleus of a particular species to be actually present at $10^{20}$ eV the following constraints can be placed: $-3\times10^{-2} \lesssim \eta \lesssim 4$ for $^{56}$Fe, $-2\times10^{-3} \lesssim \eta \lesssim 3\times10^{-2}$ for $^{16}$O and $-7\times10^{-5} \lesssim \eta \lesssim 1\times10^{-4}$ for $^{4}$He, respectively.