ABSTRACT
Understanding the transport of energetic cosmic rays belongs to the most challenging topics in astrophysics. Diffusion due to scattering by electromagnetic fluctuations is a key process in ...cosmic ray transport. The transition from a ballistic to a diffusive-propagation regime is presented in direct numerical calculations of diffusion coefficients for homogeneous magnetic field lines subject to turbulent perturbations. Simulation results are compared with theoretical derivations of the parallel diffusion coefficient’s dependences on the energy and the fluctuation amplitudes in the limit of weak turbulence. The present study shows that the widely used extrapolation of the energy scaling for the parallel diffusion coefficient to high turbulence levels predicted by quasi-linear theory does not provide a universally accurate description in the resonant-scattering regime. It is highlighted here that the numerically calculated diffusion coefficients can be polluted for low energies due to missing resonant interaction possibilities of the particles with the turbulence. Five reduced-rigidity regimes are established, which are separated by analytical boundaries derived in this work. Consequently, a proper description of cosmic ray propagation can only be achieved by using a turbulence-level-dependent diffusion coefficient and can contribute to solving the Galactic cosmic ray gradient problem.
Cosmic-ray transport in astrophysical environments is often dominated by the diffusion of particles in a magnetic field composed of both a turbulent and a mean component. This process, which is ...two-fold turbulent mixing in that the particle motion is stochastic with respect to the field lines, needs to be understood in order to properly model cosmic-ray signatures. One of the most important aspects in the modeling of cosmic-ray diffusion is that fully resonant scattering, the most effective such process, is only possible if the wave spectrum covers the entire range of propagation angles. By taking the wave spectrum boundaries into account, we quantify cosmic-ray diffusion parallel and perpendicular to the guide field direction at turbulence levels above 5% of the total magnetic field. We apply our results of the parallel and perpendicular diffusion coefficient to the Milky Way. We show that simple purely diffusive transport is in conflict with observations of the inner Galaxy, but that just by taking a Galactic wind into account, data can be matched in the central 5 kpc zone. Further comparison shows that the outer Galaxy at
>
5
kpc, on the other hand, should be dominated by perpendicular diffusion, likely changing to parallel diffusion at the outermost radii of the Milky Way.
Numerical simulations of the propagation of charged particles through magnetic fields solving the equation of motion often leads to the usage of an interpolation in the case of discretely defined ...magnetic fields, typically given on a homogeneous grid structure. However, the interpolation method influences the magnetic field properties and, therefore, also the propagation of particles through these fields. To determine the resulting error, we compare three different interpolation routines-trilinear, tricubic, and nearest neighbor interpolation-in the case of isotropic, turbulent magnetic fields. First, we analyze the impact of the different interpolation methods on the root mean square field strength, the divergence, and the spectrum of the turbulent magnetic field. Here, the nearest neighbor interpolation shows some clear benefits compared with the trilinear method; however, that changes significantly if we consider the particle propagation. In principle, a better interpolation method also yields a better description of the particle transport. In the case of field line random walk, it is shown that none of these methods, especially not the nearest neighbor interpolation, is able to yield an accurate description of the diffusion coefficient, exposing the need for a continuous, grid-less turbulent magnetic field. We optimize the performance of an algorithm that generates such a magnetic field by more than an order of magnitude. Further, the necessary number of wave-modes is determined, so that this continuous method supports realistic simulations over a larger energy range without limitations by the available memory.
ABSTRACT
Understanding the time-scales for diffusive processes and their degree of anisotropy is essential for modelling cosmic ray transport in turbulent magnetic fields. We show that the diffusion ...time-scales are isotropic over a large range of energy and turbulence levels, notwithstanding the high degree of anisotropy exhibited by the components of the diffusion tensor for cases with an ordered magnetic field component. The predictive power of the classical scattering relation as a description for the relation between the parallel and perpendicular diffusion coefficients is discussed and compared to numerical simulations. Very good agreement for a large parameter space is found, transforming classical scattering relation predictions into a computational prescription for the perpendicular component. We discuss and compare these findings, in particular, the time-scales to become diffusive with the time-scales that particles reside in astronomical environments, the so-called escape time-scales. The results show that, especially at high energies, the escape times obtained from diffusion coefficients may exceed the time-scales required for diffusion. In these cases, the escape time cannot be determined by the diffusion coefficients.
Abstract
In astrophysics, the search for sources of the highest-energy cosmic rays continues. For further progress, not only ever better observatories but also ever more realistic numerical ...simulations are needed. We compare different approaches for numerical test simulations of ultra-high-energy cosmic rays in the intergalactic magnetic field and show that all methods provide correct statistical propagation characteristics of the particles in means of their diffusive behaviour. Through convergence tests, we show that the necessary requirements for the methods differ and ultimately reveal significant differences in the required simulation time.
Abstract
Astro-COLIBRI is a novel tool that evaluates alerts of transient observations in real time, filters them by user-specified criteria, and puts them into their multiwavelength and ...multimessenger context. Through fast generation of an overview of persistent sources as well as transient events in the relevant phase space, Astro-COLIBRI contributes to an enhanced discovery potential of both serendipitous and follow-up observations of the transient sky. The software’s architecture comprises a Representational State Transfer Application Programming Interface, both a static and a real-time database, a cloud-based alert system, as well as a website and apps for iOS and Android as clients for users. The latter provide a graphical representation with a summary of the relevant data to allow for the fast identification of interesting phenomena along with an assessment of observing conditions at a large selection of observatories around the world.
Context.
The details of cosmic-ray transport have a strong impact on galaxy evolution. The peak of the cosmic-ray energy distribution is observable in the radio continuum using the electrons as ...proxy.
Aims.
We aim to measure the distance that the cosmic-ray electrons (CREs) are transported during their lifetime in the nearby galaxy M 51 across one order of magnitude in cosmic-ray energy (approximately 1–10 GeV). To this end, we use new ultra-low frequency observations from the LOw Frequency ARay (LOFAR) at 54 MHz and ancillary data between 144 and 8350 MHz.
Methods.
As the CREs originate from supernova remnants, the radio maps are smoothed in comparison to the distribution of the star formation. By convolving the map of the star formation rate (SFR) surface density with a Gaussian kernel, we can linearise the radio–SFR relation. The best-fitting convolution kernel is then our estimate of the CRE transport length.
Results.
We find that the CRE transport length increases at low frequencies, as expected since the CRE have longer lifetimes. The CRE transport length is
l
CRE
= √4
Dt
syn
, where
D
is the isotropic diffusion coefficient and
t
syn
is the CRE lifetime as given by synchrotron and inverse Compton losses. We find that the data can be well fitted by diffusion, where
D
= (2.14 ± 0.13)×10
28
cm
2
s
−1
. With
D
∝
E
0.001 ± 0.185
, the diffusion coefficient is independent of the CRE energy
E
in the range considered.
Conclusions.
Our results suggest that the transport of GeV-cosmic ray electrons in the star-forming discs of galaxies is governed by energy-independent diffusion.
ABSTRACT Understanding the transport of energetic cosmic rays belongs to the most challenging topics in astrophysics. Diffusion due to scattering by electromagnetic fluctuations is a key process in ...cosmic ray transport. The transition from a ballistic to a diffusive-propagation regime is presented in direct numerical calculations of diffusion coefficients for homogeneous magnetic field lines subject to turbulent perturbations. Simulation results are compared with theoretical derivations of the parallel diffusion coefficient’s dependences on the energy and the fluctuation amplitudes in the limit of weak turbulence. The present study shows that the widely used extrapolation of the energy scaling for the parallel diffusion coefficient to high turbulence levels predicted by quasi-linear theory does not provide a universally accurate description in the resonant-scattering regime. It is highlighted here that the numerically calculated diffusion coefficients can be polluted for low energies due to missing resonant interaction possibilities of the particles with the turbulence. Five reduced-rigidity regimes are established, which are separated by analytical boundaries derived in this work. Consequently, a proper description of cosmic ray propagation can only be achieved by using a turbulence-level-dependent diffusion coefficient and can contribute to solving the Galactic cosmic ray gradient problem.
We report a search for nonstandard neutrino interactions (NSI) using eight years of TeV-scale atmospheric muon neutrino data from the IceCube Neutrino Observatory. By reconstructing incident energies ...and zenith angles for atmospheric neutrino events, this analysis presents unified confidence intervals for the NSI parameter epsilon(mu tau). The best-fit value is consistent with no NSI at a p value of 25.2%. With a 90% confidence interval of -0.0041 <= epsilon(mu tau) <= 0.0031 along the real axis and similar strength in the complex plane, this result is the strongest constraint on any NSI parameter from any oscillation channel to date.