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.
ABSTRACT We present a systematic study of non-thermal electron-proton plasma and its emission processes in starburst galaxies in order to explain the correlation between the luminosity in the radio ...band and the recently observed gamma luminosity. In doing so, a steady state description of the cosmic-ray (CR) electrons and protons within the spatially homogeneous starburst is considered where continuous momentum losses are included as well as catastrophic losses due to diffusion and advection. The primary source of the relativistic CRs, e.g., supernova remnants, provides a quasi-neutral plasma with a power-law spectrum in momentum where we account for rigidity-dependent differences between the electron and proton spectrum. We examine the resulting leptonic and hadronic radiation processes by synchrotron radiation, inverse Compton scattering, Bremsstrahlung, and hadronic pion production. Finally, the observations of NGC 253, M82, NGC 4945, and NGC 1068 in the radio and gamma-ray bands as well as the observed supernova rate are used to constrain a best-fit model. In the case of NGC 253, M82, and NGC 4945 our model is able to accurately describe the data, showing that: (i) supernovae are the dominant particle accelerators for NGC 253, M82, and NGC 4945, but not for NGC 1068; (ii) all considered starburst galaxies are poor proton calorimeters in which for NGC 253 the escape is predominantly driven by the galactic wind, whereas the diffusive escape dominates in NGC 4945 and M82 (at energies >1 TeV); and (iii) secondary electrons from hadronic pion production are important to model the radio flux, but the associated neutrino flux is below the current observation limit.
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.
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.
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.
The recent detection of possible neutrino emission from the blazar TXS 0506+056 was the first high-energy neutrino associated with an astrophysical source, making this special type of active galaxies ...a promising neutrino emitter. The fact that two distinct episodes of neutrino emission were detected with a separation of around 3 yr raises the possibility that emission could be periodic. Periodic emission is expected from supermassive binary black hole systems due to jet precession close to the binary's merger. Here, we show that, if TXS 0506+056 is a binary source, then the next neutrino flare could already have occurred, possibly still hidden in IceCube's not-yet-analyzed data. We derive the binary properties that would lead to the detection of gravitational waves from this system by the Laser Interferometer Space Antenna (LISA) over the next decade. Our results for the first time quantify the timescale of these correlations for the example of TXS 0506+056, connecting the possible neutrino and gravitational-wave signatures of such sources.
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
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.
Context.
Tensions between the diffuse gamma-ray sky observed by the
Fermi
Large Area Telescope (
Fermi
-LAT) and the diffuse, high-energy neutrino sky detected by the IceCube South Pole Neutrino ...Observatory raise questions about our knowledge of high-energy neutrino sources in the gamma-ray regime. While blazars are among the most energetic persistent particle accelerators in the Universe, studies suggest that they could account for up to 10–30% of the neutrino flux measured by IceCube.
Aims.
Our recent results highlight that the associated IceCube neutrinos arrived in a local gamma-ray minimum (dip) of three strong neutrino point-source candidates. Here, we increase the sample of neutrino-source candidates in order to study their gamma-ray light curves.
Methods.
We generated the one-year
Fermi
-LAT light curve for eight neutrino-source candidate blazars (RBS 0958, GB6 J1040+0617, PKS 1313-333, TXS 0506+056, PKS 1454-354, NVSS J042025-374443, PKS 0426-380, and PKS 1502+106), centered on the detection time of the associated IceCube neutrinos. We applied the Bayesian block algorithm to the light curves to characterize their variability.
Results.
Our results indicate that GB6 J1040+0617 was in a phase of high gamma-ray activity, while none of the other seven neutrino-source candidates were statistically bright during the detection of the corresponding neutrinos; indeed, most of the time neutrinos arrived in a faint gamma-ray phase of the light curves. This suggests either that the eight source candidate blazars (associated with seven neutrino events) in our reduced sample are not the sources of the corresponding IceCube neutrinos, or that an in-source effect (e.g., the suppression of gamma rays due to high gamma-gamma opacity) complicates the multimessenger scenario of neutrino emission for these blazars.
Abstract
Motivated by cosmic ray (CR) re-acceleration at a potential Galactic Wind Termination
Shock (GWTS), we present a numerical model for time-dependent Diffusive Shock Acceleration (DSA). We use ...the stochastic differential equation solver (DiffusionSDE) of the cosmic ray
propagation framework CRPropa3.2 with two modifications: An importance sampling module is
introduced to improve statistics at high energies in order to keep the simulation time short. An
adaptive time step is implemented in the DiffusionSDE module. This ensures to efficiently meet
constraints on the time and diffusion step, which is crucial to obtain the correct shock spectra.
The time evolution of the spectrum at a one-dimensional planar shock is verified against the
solution obtained by the grid-based solver VLUGR3 for both energy-independent and energy-dependent
diffusion. We show that the injection of pre-accelerated particles can lead to a broken power law
spectrum in momentum if the incoming spectrum of CRs is harder than the re-accelerated
spectrum. If the injected spectrum is steeper, the shock spectrum dominates at all energies. We
finally apply the developed model to the GWTS by considering a spherically symmetric shock, a
spiral Galactic magnetic field, and anisotropic diffusion. The time-dependent spectrum at the
shock is modeled as a basis for further studies.