Green Open-Access
Refine
Language
- English (25)
Is part of the Bibliography
- yes (25)
Keywords
- gamma rays: general (6)
- ISM: supernova remnants (4)
- instabilities (4)
- acceleration of particles (3)
- Laboratory astrophysics (2)
- astroparticle physics (2)
- cosmic rays (2)
- galaxies: active (2)
- gamma rays: galaxies (2)
- magnetic fields (2)
Institute
Many solar wind observations at 1 au indicate that the proton (as well as electron) temperature anisotropy is limited. The data distribution in the (A(a), beta(a),(parallel to))-plane have a rhombic-shaped form around beta(a),(parallel to) similar to 1. The boundaries of the temperature anisotropy at beta(a),(parallel to) > 1 can be well explained by the threshold conditions of the mirror (whistler) and oblique proton (electron) firehose instabilities in a bi-Maxwellian plasma, whereas the physical mechanism of the similar restriction at beta(a),(parallel to) < 1 is still under debate. One possible option is Coulomb collisions, which we revisit in the current work. We derive the relaxation rate nu(A)(aa) of the temperature anisotropy in a bi-Maxwellian plasma that we then study analytically and by observed proton data from WIND. We found that nu(A)(pp) increases toward small beta(p),(parallel to) < 1. We matched the data distribution in the (A(p), beta(p),(parallel to))-plane with the constant contour nu(A)(pp) = 2.8 . 10(-6) s(-1), corresponding to the minimum value for collisions to play a role. This contour fits rather well the left boundary of the rhombic-shaped data distribution in the (A(p), beta(p),(parallel to))-plane. Thus, Coulomb collisions are an interesting candidate for explaining the limitations of the temperature anisotropy in the solar wind with small beta(a),(parallel to) < 1 at 1 au.
This work revisits the electrostatic instability for blazar-induced pair beams propagating through the intergalactic medium (IGM) using linear analysis and PIC simulations. We study the impact of the realistic distribution function of pairs resulting from the interaction of high-energy gamma-rays with the extragalactic background light. We present analytical and numerical calculations of the linear growth rate of the instability for the arbitrary orientation of wave vectors. Our results explicitly demonstrate that the finite angular spread of the beam dramatically affects the growth rate of the waves, leading to the fastest growth for wave vectors quasi-parallel to the beam direction and a growth rate at oblique directions that is only a factor of 2-4 smaller compared to the maximum. To study the nonlinear beam relaxation, we performed PIC simulations that take into account a realistic wide-energy distribution of beam particles. The parameters of the simulated beam-plasma system provide an adequate physical picture that can be extrapolated to realistic blazar-induced pairs. In our simulations, the beam looses only 1% of its energy, and we analytically estimate that the beam would lose its total energy over about 100 simulation times. An analytical scaling is then used to extrapolate the parameters of realistic blazar-induced pair beams. We find that they can dissipate their energy slightly faster by the electrostatic instability than through inverse-Compton scattering. The uncertainties arising from, e.g., details of the primary gamma-ray spectrum are too large to make firm statements for individual blazars, and an analysis based on their specific properties is required.
We revisit the effect of nonlinear Landau (NL) damping on the electrostatic instability of blazar-induced pair beams, using a realistic pair-beam distribution. We employ a simplified 2D model in k-space to study the evolution of the electric-field spectrum and to calculate the relaxation time of the beam. We demonstrate that the 2D model is an adequate representation of the 3D physics. We find that nonlinear Landau damping, once it operates efficiently, transports essentially the entire wave energy to small wave numbers where wave driving is weak or absent. The relaxation time also strongly depends on the intergalactic medium temperature, T-IGM, and for T-IGM << 10 eV, and in the absence of any other damping mechanism, the relaxation time of the pair beam is longer than the inverse Compton (IC) scattering time. The weak late-time beam energy losses arise from the accumulation of wave energy at small k, that nonlinearly drains the wave energy at the resonant k of the pair-beam instability. Any other dissipation process operating at small k would reduce that wave-energy drain and hence lead to stronger pair-beam energy losses. As an example, collisions reduce the relaxation time by an order of magnitude, although their rate is very small. Other nonlinear processes, such as the modulation instability, could provide additional damping of the nonresonant waves and dramatically reduce the relaxation time of the pair beam. An accurate description of the spectral evolution of the electrostatic waves is crucial for calculating the relaxation time of the pair beam.
We present very-high-resolution 1D MHD simulations of the late-stage supernova remnants (SNRs). In the post-adiabatic stage, the magnetic field has an important and significant dynamical effect on the shock dynamics, the flow structure, and hence the acceleration and emission of cosmic rays. We find that the tangential component of the magnetic field provides pressure support that to a fair degree prevents the collapse of the radiative shell and thus limits the total compression ratio of the partially or fully radiative forward shock. A consequence is that the spectra of cosmic rays would not be as hard as in hydrodynamic simulations. We also investigated the effect on the flow profiles of the magnetic-field inclination and a large-scale gradient in the gas density and/or the magnetic field. A positive density gradient shortens the evolutionary stages, whereas a shock obliquity lowers the shock compression. The compression of the tangential component of the magnetic field leads to its dominance in the downstream region of post-adiabatic shocks for a wide range of orientation of the upstream field, which may explain why one preferentially observes tangential radio polarization in old SNRs. As most cosmic rays are produced at late stages of SNR evolution, the post-adiabatic phase and the influence of the magnetic field during it are most important for modeling the cosmic-ray acceleration at old SNRs and the gamma-ray emission from late-stage SNRs interacting with clouds.
An anomalous gamma-ray excess emission has been found in the Fermi Large Area Telescope data1 covering the centre of the Galaxy2,3. Several theories have been proposed for this ‘Galactic centre excess’. They include self-annihilation of dark-matter particles4, an unresolved population of millisecond pulsars5, an unresolved population of young pulsars6, or a series of burst events7. Here, we report on an analysis that exploits hydrodynamical modelling to register the position of interstellar gas associated with diffuse Galactic gamma-ray emission. We find evidence that the Galactic centre excess gamma rays are statistically better described by the stellar over-density in the Galactic bulge and the nuclear stellar bulge, rather than a spherical excess. Given its non-spherical nature, we argue that the Galactic centre excess is not a dark-matter phenomenon but rather associated with the stellar population of the Galactic bulge and the nuclear bulge.
The variation in the intensity of cosmic rays at small angular scales is attributed to the interstellar turbulence in the vicinity of the Solar system. We show that a turbulent origin of the small-scale structures implies that the morphology of the observed cosmic ray intensity skymap varies with our location in the interstellar turbulence. The gyroradius of cosmic rays is shown to be the length scale associated with an observable change in the skymap over a radian angular scale. The extent to which the intensity at a certain angular scale varies is proportional to the change in our location with a maximum change of about the amplitude of intensity variation at that scale in the existing skymap. We suggest that for TeV cosmic rays a measurable variation could occur over a time-scale of a decade due to the Earth’s motion through the interstellar medium, if interstellar turbulence persists down to the gyroradius, about 300 μpc for TeV-ish cosmic rays. Observational evidence of the variability, or an absence of it, could provide a useful insight into the physical origin of the small-scale anisotropy.
Preliminary study for the laboratory experiment of cosmic-rays driven magnetic field amplification
(2019)
To understand astrophysical magnetic-field amplification, we conducted a feasibility study for a laboratory experiment of a non-resonant streaming instability at the Photo Injector Test Facility at DESY, Zeuthen site (PITZ). This non-resonant streaming instability, also known as Bell’s instability, is generally regarded as a candidate for the amplification of interstellar magnetic field in the upstream region of supernova-remnant shocks, which is crucial for the efficiency of diffusive shock acceleration. In the beam-plasma system composed of a radio-frequency electron gun and a gas-discharge plasma cell, the goal of our experiment is to demonstrate the development of the non-resonant streaming instability and to find its saturation level in the laboratory environment. Since we find that the electron beam will be significantly decelerated on account of an electrostatic streaming instability, which will decrease the growth rate of desired non-resonant streaming instability, we discuss possible ways to suppress the electrostatic streaming instability by considering the characteristics of a field-emission-based quasi continuous-wave electron beam.
Precursor Wave Amplification by Ion-Electron Coupling through Wakefield in Relativistic Shocks
(2019)
We investigated electromagnetic precursor wave emission in relativistic shocks by using two-dimensional particle-in-cell simulations. We found that the wave amplitude is significantly enhanced by a positive feedback process associated with ion-electron coupling through the wakefields for high magnetization. The wakefields collapse during the nonlinear process of the parametric decay instability in the near-upstream region, where nonthermal electrons and ions are generated. The intense coherent emission and the particle acceleration may operate in high-energy astrophysical objects.
In September 2017, the IceCube Neutrino Observatory recorded a very-high-energy neutrino in directional coincidence with a blazar in an unusually bright gamma-ray state, TXS0506 + 056 (refs(1,2)). Blazars are prominent photon sources in the Universe because they harbour a relativistic jet whose radiation is strongly collimated and amplified. High-energy atomic nuclei known as cosmic rays can produce neutrinos; thus, the recent detection may help in identifying the sources of the diffuse neutrino flux(3) and the energetic cosmic rays. Here we report a self-consistent analysis of the physical relation between the observed neutrino and the blazar, in particular the time evolution and spectral behaviour of neutrino and photon emission. We demonstrate that a moderate enhancement in the number of cosmic rays during the flare can yield a very strong increase in the neutrino flux, which is limited by co-produced hard X-rays and teraelectronvolt gamma rays. We also test typical radiation models(4,5) for compatibility and identify several model classes(6,7) as incompatible with the observations. We investigate to what degree the findings can be generalized to the entire population of blazars, determine the relation between their output in photons, neutrinos and cosmic rays, and suggest how to optimize the strategy of future observations.
Ground-based gamma-ray astronomy has had a major breakthrough with the impressive results obtained using systems of imaging atmospheric Cherenkov telescopes. Ground-based gamma-ray astronomy has a huge potential in astrophysics, particle physics and cosmology. CTA is an international initiative to build the next generation instrument, with a factor of 5-10 improvement in sensitivity in the 100 GeV-10 TeV range and the extension to energies well below 100 GeV and above 100 TeV. CTA will consist of two arrays (one in the north, one in the south) for full sky coverage and will be operated as open observatory. The design of CTA is based on currently available technology. This document reports on the status and presents the major design concepts of CTA.