Source author record

Martin Pohl

Martin Pohl appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

53works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

53 published item(s)

preprint2022arXiv

Assessing the Impact of Hydrogen Absorption on the Characteristics of the Galactic Center Excess

We present a new reconstruction of the distribution of atomic hydrogen in the inner Galaxy that is based on explicit radiation-transport modelling of line and continuum emission and a gas-flow model in the barred Galaxy that provides distance resolution for lines of sight toward the Galactic Center. The main benefits of the new gas model are, a), the ability to reproduce the negative line signals seen with the H$I$4PI survey and, b), the accounting for gas that primarily manifests itself through absorption. We apply the new model of Galactic atomic hydrogen to an analysis of the diffuse gamma-ray emission from the inner Galaxy, for which an excess at a few GeV was reported that may be related to dark matter. We find with high significance an improved fit to the diffuse gamma-ray emission observed with the Fermi-LAT, if our new H$I$ model is used to estimate the cosmic-ray induced diffuse gamma-ray emission. The fit still requires a nuclear bulge at high significance. Once this is included there is no evidence for a dark-matter signal, be it cuspy or cored. But an additional so-called boxy bulge is still favoured by the data. This finding is robust under the variation of various parameters, for example the excitation temperature of atomic hydrogen, and a number of tests for systematic issues.

preprint2022arXiv

Diffusive shock acceleration at oblique high Mach number shocks

The current paradigm of cosmic ray (CR) origin states that the most part of galactic CRs is produced by supernova remnants. The interaction of supernova ejecta with the interstellar medium after supernova's explosions results in shocks responsible for CR acceleration via diffusive shock acceleration (DSA). We use particle-in-cell (PIC) simulations and a combined PIC-magnetohydrodynamic (PIC-MHD) technique to investigate whether DSA can occur in oblique high Mach number shocks. Using the PIC method, we follow the formation of the shock and determine the fraction of the particles that gets involved in DSA. Then, with this result, we use PIC-MHD simulations to model the large-scale structure of the plasma and the magnetic field surrounding the shock and find out whether or not the reflected particles can generate the upstream turbulence and trigger DSA. We find that the feasibility of this process in oblique shocks depends strongly on the Alfvenic Mach number, and the DSA process is more likely triggered at high Mach number shocks.

preprint2022arXiv

Pre-acceleration in the Electron Foreshock I: Electron Acoustic Waves

To undergo diffusive shock acceleration, electrons need to be pre-accelerated to increase their energies by several orders of magnitude, else their gyro-radii are smaller than the finite width of the shock. In oblique shocks, where the upstream magnetic field orientation is neither parallel or perpendicular to the shock normal, electrons can escape to the shock upstream, modifying the shock foot to a region called the electron foreshock. To determine the pre-acceleration in this region, we undertake PIC simulations of oblique shocks while varying the obliquity and in-plane angles. We show that while the proportion of reflected electrons is negligible for $θ_{\rm Bn} = 74.3^\circ$, it increases to $R \sim 5\%$ for $θ_{\rm Bn} = 30^\circ$, and that, via the electron acoustic instability, these electrons power electrostatic waves upstream with energy density proportional to $R^{0.6}$ and a wavelength $\approx 2 λ_{\rm{se}}$, where $λ_{\rm{se}}$ is the electron skin length. While the initial reflection mechanism is typically a combination of shock surfing acceleration and magnetic mirroring, we show that once the electrostatic waves have been generated upstream they themselves can increase the momenta of upstream electrons parallel to the magnetic field. In $\lesssim 1\%$ of cases, upstream electrons are prematurely turned away from the shock and never injected downstream. In contrast, a similar fraction are re-scattered back towards the shock after reflection, re-interact with the shock with energies much greater than thermal, and cross into the downstream.

preprint2022arXiv

Spectral softening in core-collapse supernova remnant expanding inside wind-blown bubble

Context. Galactic cosmic rays are widely assumed to arise from diffusive shock acceleration, specifically at shocks in supernova remnants (SNRs). These shocks expand in a complex environment, particularly in the core-collapse scenario as these SNRs evolve inside the wind-blown bubbles created by their progenitor stars. The cosmic rays (CRs) at core-collapse SNRs may carry spectral signatures of that complexity. Aims. We study particle acceleration in the core-collapse SNR of a progenitor with initial mass 60 $M_\odot$ and realistic stellar evolution. The SNR shock interacts with discontinuities inside the wind-blown bubble and generates several transmitted and reflected shocks. We analyse their impact on particle spectra and the resulting emission from the remnant. Methods. The hydrodynamic equations for the evolution of SNR inside the pre-supernova circumstellar medium have been solved simultaneously with the transport equation for cosmic rays in test-particle approximation and with the induction equation for the magnetohydrodynamics (MHD) in 1-D spherical symmetry. Results. The evolution of core-collapse SNRs inside complex wind-blown bubbles modifies the spectra of both the particles and their emission. We have found softer particle spectra with spectral indices close to 2.5 during shock propagation inside the shocked wind, and this softness persists at later evolutionary stages. Further, our calculated total production spectrum released into the interstellar medium demonstrates spectral consistency at high energy with the galactic CRs injection spectrum, required in propagation models. The magnetic field structure effectively influences the emission morphology of SNR as it governs the transportation of particles and the synchrotron emissivity.

preprint2022arXiv

Suppression of the TeV pair-beam plasma instability by a tangled weak intergalactic magnetic field

We study the effect of a tangled sub-fG level intergalactic magnetic field (IGMF) on the electrostatic instability of a blazar-induced pair beam. Sufficiently strong IGMF may significantly deflect the TeV pair beams, which would reduce the flux of secondary cascade emission below the observational limits. A similar flux reduction may result from the electrostatic beam-plasma instability, which operates the best in the absence of IGMF. Considering IGMF with correlation lengths smaller than a kpc, we find that weak magnetic fields increase the transverse momentum of the pair beam particles, which dramatically reduces the linear growth rate of the electrostatic instability and hence the energy-loss rate of the pair beam. We show that the beam-plasma instability is eliminated as an effective energy-loss agent at a field strength three orders of magnitude below that needed to suppress the secondary cascade emission by magnetic deflection. For intermediate-strength IGMF, we do not know a viable process to explain the observed absence of GeV-scale cascade emission.

preprint2022arXiv

The electron foreshock at high-Mach-number nonrelativistic oblique shocks

In the Universe matter outside of stars and compact objects is mostly composed of collisionless plasma. The interaction of a supersonic plasma flow with an obstacle results in collisionless shocks that are often associated with intense nonthermal radiation and the production of cosmic ray particles. Motivated by simulations of non-relativistic high-Mach-number shocks in supernova remnants, we investigate the instabilities excited by relativistic electron beams in the extended foreshock of oblique shocks. The phase-space distributions in the inner and outer foreshock regions are derived with a Particle-in-Cell simulation of the shock and used as initial conditions for simulations with periodic boundary conditions to study their relaxation towards equilibrium. We find that the observed electron-beam instabilities agree very well with the predictions of a linear dispersion analysis: the electrostatic electron-acoustic instability dominates in the outer region of the foreshock, while the denser electron beams in the inner foreshock drive the gyroresonant oblique-whistler instability.

preprint2021arXiv

Leptonic non-thermal emission from supernova remnants evolving in the circumstellar magnetic field

The very-high-energy (VHE; E > 100 GeV) gamma-ray emission observed from a number of Supernova remnants (SNRs) indicates particle acceleration to high energies at the shock of the remnants and a potentially significant contribution to Galactic cosmic rays. It is extremely difficult to determine whether protons (through hadronic interactions and subsequent pion decay) or electrons (through inverse Compton scattering on ambient photon fields) are responsible for this emission. For a successful diagnostic, a good understanding of the spatial and energy distribution of the underlying particle population is crucial. Most SNRs are created in core-collapse explosions and expand into the wind bubble of their progenitor stars. This circumstellar medium features a complex spatial distribution of gas and magnetic field which naturally strongly affects the resulting particle population. In this work, we conduct a detailed study of the spectro-spatial evolution of the electrons accelerated at the forward shock of core-collapse SNRs and their non-thermal radiation, using the RATPaC code that is designed for the time- and spatially dependent treatment of particle acceleration at SNR shocks. We focus on the impact of the spatially inhomogeneous magnetic field through the efficiency of diffusion and synchrotron cooling. It is demonstrated that the structure of the circumstellar magnetic field can leave strong signatures in the spectrum and morphology of the resulting non-thermal emission.

preprint2021arXiv

Mildly relativistic magnetized shocks in electron-ion plasmas -- II. Particle acceleration and heating

Particle acceleration and heating at mildly relativistic magnetized shocks in electron-ion plasma are investigated with unprecedentedly high-resolution two-dimensional particle-in-cell simulations that include ion-scale shock rippling. Electrons are super-adiabatically heated at the shock, and most of the energy transfer from protons to electrons takes place at or downstream of the shock. We are the first to demonstrate that shock rippling is crucial for the energization of electrons at the shock. They remain well below equipartition with the protons. The downstream electron spectra are approximately thermal with a limited supra-thermal power-law component. Our results are discussed in the context of wakefield acceleration and the modelling of electromagnetic radiation from blazar cores.

preprint2020arXiv

Kinetic simulations of nonrelativistic perpendicular shocks of young supernova remnants. III. Magnetic reconnection

Fully kinetic two-dimensional particle-in-cell simulations are used to study electron acceleration at high-Mach-number nonrelativistic perpendicular shocks. SNR shocks are mediated by the Weibel instability which is excited because of an interaction between shock-reflected and upstream ions. Nonlinear evolution of the Weibel instability leads to the formation of current sheets. At the turbulent shock ramp the current sheets decay through magnetic reconnection. The number of reconnection sites strongly depends on the ion-to-electron mass ratio and the Alfvénic Mach number of the simulated shock. Electron acceleration is observed at locations where magnetic reconnection operates. For the highest mass ratios almost all electrons are involved in magnetic reconnection, which makes the magnetic reconnection the dominant acceleration process for electrons at these shocks. We discuss the relevance of our results for 3D systems with realistic ion-to-electron mass ratio.

preprint2020arXiv

Mildly relativistic magnetized shocks in electron-ion plasmas I. Electromagnetic shock structure

Mildly relativistic shocks in magnetized electron-ion plasmas are investigated with 2D kinetic particle-in-cell simulations of unprecedentedly high resolution and large scale for conditions that may be found at internal shocks in blazar cores. Ion-scale effects cause corrugations along the shock surface whose properties somewhat depend on the configuration of the mean perpendicular magnetic field, that is either in or out of the simulation plane. We show that the synchrotron maser instability persists to operate in mildly relativistic shocks in agreement with theoretical predictions and produces coherent emission of upstream-propagating electromagnetic waves. Shock front ripples are excited in both mean-field configurations and they engender effective wave amplification. The interaction of these waves with upstream plasma generates electrostatic wakefields.

preprint2020arXiv

Production of secondary particles in heavy nuclei interactions in supernova remnants

Depending on their type, supernova remnants may have ejecta material with high abundance of heavy elements such as carbon or oxygen. In addition, core-collapse supernovae explode in the wind material of their progenitor star that may also have a high abundance of heavy elements. Hadronic collisions in these enriched media spawn the production of gamma rays, neutrinos, and secondary electrons with spectra that cannot be scaled from those calculated for pp collisions, potentially leading to erroneous results. We used Monte-Carlo event generators to calculate the differential production rate of particles such as gamma rays, neutrinos, and secondary electrons for H, He, C, and O nuclei as projectiles and as target material. The cross sections and the multiplicity spectra are separately computed for each of the 16 combinations of projectile and target. We describe characteristic effects of heavy nuclei in the shape and normalization of the spectra of various particles produced.

preprint2020arXiv

Rapid Particle Acceleration due to Recollimation Shocks and Turbulent Magnetic Fields in Injected Jets with Helical Magnetic Fields

One of the key questions in the study of relativistic jets is how magnetic reconnection occurs and whether it can effectively accelerate electrons in the jet. We performed 3D particle-in-cell (PIC) simulations of a relativistic electron-proton jet of relatively large radius that carries a helical magnetic field. We focussed our investigation on the interaction between the jet and the ambient plasma and explore how the helical magnetic field affects the excitation of kinetic instabilities such as the Weibel instability (WI), the kinetic Kelvin-Helmholtz instability (kKHI), and the mushroom instability (MI). In our simulations these kinetic instabilities are indeed excited, and particles are accelerated. At the linear stage we observe recollimation shocks near the center of the jet. As the electron-proton jet evolves into the deep nonlinear stage, the helical magnetic field becomes untangled due to reconnection-like phenomena, and electrons are repeatedly accelerated as they encounter magnetic-reconnection events in the turbulent magnetic field.

preprint2019arXiv

Cosmic-ray acceleration and escape from post-adiabatic Supernova remnants

Supernova remnants are known to accelerate cosmic rays on account of their non-thermal emission of radio waves, X-rays, and gamma rays. Although there are many models for the acceleration of cosmic rays in Supernova remnants, the escape of cosmic rays from these sources is yet understudied. We use our time-dependent acceleration code RATPaC to study the acceleration of cosmic rays and their escape in post-adiabatic Supernova remnants and calculate the subsequent gamma-ray emission from inverse-Compton scattering and Pion decay. We performed spherically symmetric 1-D simulations in which we simultaneously solve the transport equations for CRs, magnetic turbulence, and the hydrodynamical flow of the thermal plasma in a volume large enough to keep all CRs in the simulation. The transport equations for cosmic-rays and magnetic turbulence are coupled via the cosmic-ray gradient and the spatial diffusion coefficient of the cosmic rays, while the cosmic-ray feedback onto the shock structure can be ignored. Our simulations span 100kyrs, thus covering the remnants evolution until the beginning of the post-adiabatic phase. At later stages of the evolution cosmic rays over a wide range of energy can reside outside of the remnant, creating spectra that are softer than predicted by standard DSA and feature breaks in the 10-100 GeV-range. The total spectrum of cosmic rays released into the interstellar medium has a spectral index of s~2.4 above roughly 10 GeV which is close to that required by Galactic propagation models. We further find the gamma-ray luminosity to peak around an age of 4,000 years for inverse-Compton-dominated high-energy emission. Remnants expanding in low-density media emit generally more inverse-Compton radiation matching the fact that the brightest known supernova remnants - RCW86, Vela Jr, HESSJ1731-347 and RXJ1713.7-3946 - are all expanding in low density environments.

preprint2019arXiv

PIC Simulation Methods for Cosmic Radiation and Plasma Instabilities

Particle acceleration in collisionless plasma systems is a central question in astroplasma and astroparticle physics. The structure of the acceleration regions, electron-ion energy equilibration, preacceleration of particles at shocks to permit further energization by diffusive shock acceleration, require knowledge of the distribution function of particles besides the structure and dynamic of electromagnetic fields, and hence a kinetic description is desirable. Particle-in-cell simulations offer an appropriate, if computationally expensive method of essentially conducting numerical experiments that explore kinetic phenomena in collisionless plasma. We review recent results of PIC simulations of astrophysical plasma systems, particle acceleration, and the instabilities that shape them.

preprint2016arXiv

Microscopic Processes in Global Relativistic Jets Containing Helical Magnetic Fields

In the study of relativistic jets one of the key open questions is their interaction with the environment on the microscopic level. Here, we study the initial evolution of both electron$-$proton ($e^{-}-p^{+}$) and electron$-$positron ($e^{\pm}$) relativistic jets containing helical magnetic fields, focusing on their interaction with an ambient plasma. We have performed simulations of "global" jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). In our initial simulation study these kinetic instabilities are suppressed and new types of instabilities can grow. In the $e^{-}-p^{+}$ jet simulation a recollimation-like instability occurs and jet electrons are strongly perturbed. In the $e^{\pm}$ jet simulation a recollimation-like instability occurs at early times followed by a kinetic instability and the general structure is similar to a simulation without helical magnetic field. Simulations using much larger systems are required in order to thoroughly follow the evolution of global jets containing helical magnetic fields.

preprint2016arXiv

Non-relativistic perpendicular shocks modeling young supernova remnants: nonstationary dynamics and particle acceleration at forward and reverse shocks

For parameters that are applicable to the conditions at young supernova remnants, we present results of 2D3V particle-in-cell simulations of a non-relativistic plasma shock with a large-scale perpendicular magnetic field inclined at 45-deg angle to the simulation plane to approximate 3D physics. We developed an improved clean setup that uses the collision of two plasma slabs with different density and velocity, leading to the development of two distinctive shocks and a contact discontinuity. The shock formation is mediated by Weibel-type filamentation instabilities that generate magnetic turbulence. Cyclic reformation is observed in both shocks with similar period, for which we note global variations on account of shock rippling and local variations arising from turbulent current filaments. The shock rippling occurs on spatial and temporal scales given by gyro-motions of shock-reflected ions. The drift motion of electrons and ions is not a gradient drift, but commensurates with E x B drift. We observe a stable suprathermal tail in the ion spectra, but no electron acceleration because the amplitude of Buneman modes in the shock foot is insufficient for trapping relativistic electrons. We see no evidence of turbulent reconnection. A comparison with other 2D simulation results suggests that the plasma beta and the ion-to-electron mass ratio are not decisive for efficient electron acceleration, but pre-acceleration efficacy might be reduced with respect to the 2D results once three-dimensional effects are fully accounted for. Other microphysical factors may also be at play to limit the amplitude of Buneman waves or prevent return of electrons to the foot region.

preprint2016arXiv

Particle diffusion and localized acceleration in inhomogeneous AGN jets - Part II: stochastic variation

We study the stochastic variation of blazar emission under a 2-D spatially resolved leptonic jet model we previously developed. Random events of particle acceleration and injection in small zones within the emission region are assumed to be responsible for flux variations. In addition to producing spectral energy distributions that describe the observed flux of Mrk 421, we further analyze the timing properties of the simulated light curves, such as the power spectral density (PSD) at different bands, flux-flux correlations, as well as the cross-correlation function between X-rays and TeV γ-rays. We find spectral breaks in the PSD at a timescale comparable to the dominant characteristic time scale in the system, which is usually the pre-defined decay time scale of an acceleration event. Cooling imposes a delay, and so PSDs taken at lower energy bands in each emission component (synchrotron or inverse Compton) generally break at longer timescales. The flux-flux correlation between X-rays and TeV γ-rays can be either quadratic or linear, depending on whether or not there are large variation of the injection into the particle acceleration process. When the relationship is quadratic, the TeV flares lag the X-ray flares, and the optical & GeV flares are large enough to be comparable to the ones in X-ray. When the relationship is linear, the lags are insignificant, and the optical & GeV flares are small.

preprint2016arXiv

The Quiescent Intracluster Medium in the Core of the Perseus Cluster

Clusters of galaxies are the most massive gravitationally-bound objects in the Universe and are still forming. They are thus important probes of cosmological parameters and a host of astrophysical processes. Knowledge of the dynamics of the pervasive hot gas, which dominates in mass over stars in a cluster, is a crucial missing ingredient. It can enable new insights into mechanical energy injection by the central supermassive black hole and the use of hydrostatic equilibrium for the determination of cluster masses. X-rays from the core of the Perseus cluster are emitted by the 50 million K diffuse hot plasma filling its gravitational potential well. The Active Galactic Nucleus of the central galaxy NGC1275 is pumping jetted energy into the surrounding intracluster medium, creating buoyant bubbles filled with relativistic plasma. These likely induce motions in the intracluster medium and heat the inner gas preventing runaway radiative cooling; a process known as Active Galactic Nucleus Feedback. Here we report on Hitomi X-ray observations of the Perseus cluster core, which reveal a remarkably quiescent atmosphere where the gas has a line-of-sight velocity dispersion of 164+/-10 km/s in a region 30-60 kpc from the central nucleus. A gradient in the line-of-sight velocity of 150+/-70 km/s is found across the 60 kpc image of the cluster core. Turbulent pressure support in the gas is 4% or less of the thermodynamic pressure, with large scale shear at most doubling that estimate. We infer that total cluster masses determined from hydrostatic equilibrium in the central regions need little correction for turbulent pressure.

preprint2016arXiv

Transport of magnetic turbulence in Supernova remnants

Context. Supernova remnants are known as sources of galactic cosmic rays for their non-thermal emission of radio waves, X-rays, and gamma-rays. However, the observed soft broken power-law spectra are hard to reproduce within standard acceleration theory based on the assumption of Bohm diffusion and steady-state calculations. Aims. We point out that a time-dependent treatment of the acceleration process together with a self-consistent treatment of the scattering turbulence amplification is necessary. Methods. We numerically solve the coupled system of transport equations for cosmic rays and isotropic Alfvenic turbulence. The equations are coupled through the growth rate of turbulence determined by the cosmic-ray gradient and the spatial diffusion coefficient of cosmic rays determined by the energy density of the turbulence. The system is solved on a co-moving expanding grid extending upstream for dozens of shock radii, allowing for the self-consistent study of cosmic-ray diffusion in the vicinity of their acceleration site. The transport equation for cosmic rays is solved in a test-particle approach. Results. We demonstrate that the system is typically not in a steady state. In fact, even after several thousand years of evolution, no equilibrium situation is reached. The resulting time-dependent particle spectra strongly differ from those derived assuming a steady state and Bohm diffusion. Our results indicate that proper accounting for the evolution of the scattering turbulence and hence the particle diffusion coefficient is crucial for the formation of the observed soft spectra. In any case, the need to continuously develop magnetic turbulence upstream of the shock introduces non-linearity in addition to that imposed by cosmic-ray feedback.

preprint2015arXiv

AMS tracking in-orbit performance

AMS-02 is a high precision magnetic spectrometer for cosmic rays in the GeV to TeV energy range. Its tracker consists of nine layers of double-sided silicon microstrip sensors. They are used to locate the trajectories of cosmic rays in the 0.14 T field of a cylindrical magnet, thus measuring their rigidity $p/Z$ and charge sign. In addition, they deliver a high resolution measurement of the absolute charge $|Z|$. The detector has been designed to operate in space with a position resolution of about 10 $μ$m for each hit and charge identification capabilities up to $Z=26$. In this talk I describe the performance in orbit of this detector component and its impact on the overall performance of the spectrometer.

preprint2015arXiv

Indirect and direct search for dark matter

The majority of the matter in the universe is still unidentified and under investigation by both direct and indirect means. Many experiments searching for the recoil of dark-matter particles off target nuclei in underground laboratories have established increasingly strong constraints on the mass and scattering cross sections of weakly interacting particles, and some have even seen hints at a possible signal. Other experiments search for a possible mixing of photons with light scalar or pseudo-scalar particles that could also constitute dark matter. Furthermore, annihilation or decay of dark matter can contribute to charged cosmic rays, photons at all energies, and neutrinos. Many existing and future ground-based and satellite experiments are sensitive to such signals. Finally, data from the Large Hadron Collider at CERN are scrutinized for missing energy as a signature of new weakly interacting particles that may be related to dark matter. In this review article we summarize the status of the field with an emphasis on the complementarity between direct detection in dedicated laboratory experiments, indirect detection in the cosmic radiation, and searches at particle accelerators.

preprint2014arXiv

Latest Results from the Alpha Magnetic Spectrometer on the International Space Station

I review latest results from AMS on electrons and positrons in primary cosmic rays in the GeV to TeV energy range. Separate fluxes for electrons and positrons as well as their sum are presented. Neither of the fluxes is compatible with a simple power law. New data on the fraction of positrons in the joint electron and positron flux are also presented, which extend the energy range of our previous observation and increase its precision. The new results show, for the first time, that above about 200 GeV the positron fraction no longer exhibits an increase with energy. The results confirm that a common new source of electrons and positrons exists.

preprint2014arXiv

Magnetic Field Amplification and Flat Spectrum Radio Quasars

We perform time-dependent, spatially-resolved simulations of blazar emission to evaluate several flaring scenarios related to magnetic-field amplification and enhanced particle acceleration. The code explicitly accounts for light-travel-time effects and is applied to flares observed in the flat spectrum radio quasar (FSRQ) PKS 0208-512, which show optical/γ-ray correlation at some times, but orphan optical flares at other times. Changes in both the magnetic field and the particle acceleration efficiency are explored as causes of flares. Generally, external Compton emission appears to describe the available data better than a synchrotron self-Compton scenario, and in particular orphan optical flares are difficult to produce in the SSC framework. X-ray soft-excesses, γ-ray spectral hardening, and the detections at very high energies of certain FSRQs during flares find natural explanations in the EC scenario with particle acceleration change. Likewise, optical flares with/without γ-ray counterparts can be explained by different allocations of energy between the magnetization and particle acceleration, which may be related to the orientation of the magnetic field relative to the jet flow. We also calculate the degree of linear polarization and polarization angle as a function of time for a jet with helical magnetic field. Tightening of the magnetic helix immediately downstream of the jet perturbations, where flares occur, can be sufficient to explain the increases in the degree of polarization and a rotation by >= 180 degree of the observed polarization angle, if light-travel-time effects are properly considered.

preprint2014arXiv

Magnetic Field Amplification and Saturation in Turbulence Behind a Relativistic Shock

We have investigated via two-dimensional relativistic MHD simulations the long-term evolution of turbulence created by a relativistic shock propagating through an inhomogeneous medium. In the postshock region, magnetic field is strongly amplified by turbulent motions triggered by preshock density inhomogeneities. Using a long-simulation box we have followed the magnetic-field amplification until it is fully developed and saturated. The turbulent velocity is sub-relativistic even for a strong shock. Magnetic-field amplification is controled by the turbulent motion and saturation occurs when the magnetic energy is comparable to the turbulent kinetic energy. Magnetic-field amplification and saturation depend on the initial strength and direction of the magnetic field in the preshock medium, and on the shock strength. If the initial magnetic field is perpendicular to the shock normal, the magnetic field is first compressed at the shock and then can be amplified by turbulent motion in the postshock region. Saturation occurs when the magnetic energy becomes comparable to the turbulent kinetic energy in the postshock region. If the initial magnetic field in the preshock medium is strong, the postshock region becomes turbulent but significant field amplification does not occur. If the magnetic energy after shock compression is larger than the turbulent kinetic energy in the postshock region, significant field amplification does not occur. We discuss possible applications of our results to gamma-ray bursts and active galactic nuclei.

preprint2014arXiv

PANGU: A High Resolution Gamma-ray Space Telescope

We describe the instrument concept of a high angular resolution telescope dedicated to the sub-GeV (from $\gtrsim$10 MeV to $\gtrsim$1 GeV) gamma-ray photon detection. This mission, named PANGU (PAir-productioN Gamma-ray Unit), has been suggested as a candidate for the joint small mission between the European Space Agency (ESA) and the Chinese Academy of Science (CAS). A wide range of topics of both astronomy and fundamental physics can be attacked with PANGU, covering Galactic and extragalactic cosmic-ray physics, extreme physics of a variety of extended (e.g. supernova remnants, galaxies, galaxy clusters) and compact (e.g. black holes, pulsars, gamma-ray bursts) objects, solar and terrestrial gamma-ray phenomena, and searching for dark matter decay and/or annihilation signature etc. The unprecedented point spread function can be achieved with a pair-production telescope with a large number of thin active tracking layers to precisely reconstruct the pair-produced electron and positron tracks. Scintillating fibers or thin silicon micro-strip detectors are suitable technology for such a tracker. The energy measurement is achieved by measuring the momentum of the electrons and positrons through a magnetic field. The innovated spectrometer approach provides superior photon pointing resolution, and is particular suitable in the sub-GeV range. The level of tracking precision makes it possible to measure the polarization of gamma rays, which would open up a new frontier in gamma-ray astronomy. The frequent full-sky survey at sub-GeV with PANGU's large field of view and significantly improved point spread function would provide crucial information to GeV-TeV astrophysics for current/future missions including Fermi, DAMPE, HERD, and CTA, and other multi-wavelength telescopes.

preprint2014arXiv

Particle detection technology for space-borne astroparticle experiments

I review the transfer of technology from accelerator-based equipment to space-borne astroparticle detectors. Requirements for detection, identification and measurement of ions, electrons and photons in space are recalled. The additional requirements and restrictions imposed by the launch process in manned and unmanned space flight, as well as by the hostile environment in orbit, are analyzed. Technology readiness criteria and risk mitigation strategies are reviewed. Recent examples are given of missions and instruments in orbit, under construction or in the planning phase.

preprint2014arXiv

Particle diffusion and localized acceleration in inhomogeneous AGN jets - Part I: Steady-state spectra

We study the acceleration, transport, and emission of particles in relativistic jets. Localized stochastic particle acceleration, spatial diffusion, and synchrotron as well as synchrotron self-Compton emission are considered in a leptonic model. To account for inhomogeneity, we use a 2D axi-symmetric cylindrical geometry for both relativistic electrons and magnetic field. In this first phase of our work, we focus on steady-state spectra that develop from a time-dependent model. We demonstrate that small isolated acceleration region in a much larger emission volume are sufficient to accelerate particles to high energy. Diffusive escape from these small regions provides a natural explanation for the spectral form of the jet emission. The location of the acceleration regions within the jet is found to affect the cooling break of the spectrum in this diffusive model. Diffusion-caused energy-dependent inhomogeneity in the jets predicts that the SSC spectrum is harder than the synchrotron spectrum. There can also be a spectral hardening towards the high-energy section of the synchrotron spectrum, if particle escape is relatively slow. These two spectral hardening effects indicate that the jet inhomogeneity might be a natural explanation for the unexpected hard γ-ray spectra observed in some blazars.

preprint2014arXiv

The ASTRO-H X-ray Astronomy Satellite

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions developed by the Institute of Space and Astronautical Science (ISAS), with a planned launch in 2015. The ASTRO-H mission is equipped with a suite of sensitive instruments with the highest energy resolution ever achieved at E > 3 keV and a wide energy range spanning four decades in energy from soft X-rays to gamma-rays. The simultaneous broad band pass, coupled with the high spectral resolution of Delta E < 7 eV of the micro-calorimeter, will enable a wide variety of important science themes to be pursued. ASTRO-H is expected to provide breakthrough results in scientific areas as diverse as the large-scale structure of the Universe and its evolution, the behavior of matter in the gravitational strong field regime, the physical conditions in sites of cosmic-ray acceleration, and the distribution of dark matter in galaxy clusters at different redshifts.

preprint2013arXiv

Is the Galactic Cosmic Ray Spectrum Constant in Time?

The hypothesis is considered that the present Galactic cosmic ray spectrum is at present softer than its time average due to source intermittency. Measurements of muogenic nuclides underground could provide an independent measurement of the time averaged spectrum. Source intermittency could also account for the surprising low anisotropy reported by the IceCube collaboration. Predictions for Galactic emission of ultrahigh-energy quanta, such as UHE gamma rays and neutrinos, might be higher or lower than previously estimated.

preprint2013arXiv

Magnetic Field Amplification and Blazar Flares

Recent multiwavelength observations of PKS 0208-512 by SMARTS, Fermi, and Swift revealed that gamma-ray and optical light curves of this flat spectrum radio quasars are highly correlated, but with an exception of one large optical flare having no corresponding gamma-ray activity or even detection. On the other hand, recent advances in SNRs observations and plasma simulations both reveal that magnetic field downstream of astrophysical shocks can be largely amplified beyond simple shock compression. These amplifications, along with their associated particle acceleration, might contribute to blazar flares, including the peculiar flare of PKS 0208-512. Using our time dependent multizone blazar emission code, we evaluate several scenarios that may represent such phenomena. This code combines Monte Carlo method that tracks the radiative processes including inverse Compton scattering, and Fokker-Planck equation that follows the cooling and acceleration of particles. It is a comprehensive time dependent code that fully takes into account the light travel time effects. In this study, both the changes of the magnetic field and acceleration efficiency are explored as the cause of blazar flares. Under these assumption, synchrotron self-Compton and external Compton scenarios produce distinct features that favor the external Compton scenario. The optical flares with/without gamma-ray counterparts can be explained by different allocations of energy between the magnetization and particle acceleration, which in turn can be affected by the relative orientation between the magnetic field and the shock flow. We compare the details of the observations and simulation, and highlight what implications this study has on our understanding of relativistic jets.

preprint2013arXiv

Modification of the formation of high-Mach number electrostatic shock-like structures by the ion acoustic instability

The formation of unmagnetized electrostatic shock-like structures with a high Mach number is examined with one- and two-dimensional particle-in-cell (PIC) simulations. The structures are generated through the collision of two identical plasma clouds, which consist of equally hot electrons and ions with a mass ratio of 250. The Mach number of the collision speed with respect to the initial ion acoustic speed of the plasma is set to 4.6. This high Mach number delays the formation of such structures by tens of inverse ion plasma frequencies. A pair of stable shock-like structures is observed after this time in the 1D simulation, which gradually evolve into electrostatic shocks. The ion acoustic instability, which can develop in the 2D simulation but not in the 1D one, competes with the nonlinear process that gives rise to these structures. The oblique ion acoustic waves fragment their electric field. The transition layer, across which the bulk of the ions change their speed, widens and their speed change is reduced. Double layer-shock hybrid structures develop.

preprint2013arXiv

The average GeV-band Emission from Gamma-Ray Bursts

We analyze the emission in the 0.3-30 GeV energy range of Gamma-Ray Bursts detected with the Fermi Gamma-ray Space Telescope. We concentrate on bursts that were previously only detected with the Gamma-Ray Burst Monitor in the keV energy range. These bursts will then be compared to the bursts that were individually detected with the Large Area Telescope at higher energies. To estimate the emission of faint GRBs we use non-standard analysis methods and sum over many GRBs to find an average signal which is significantly above background level. We use a subsample of 99 GRBs listed in the Burst Catalog from the first two years of observation. Although mostly not individually detectable, the bursts not detected by the Large Area Telescope on average emit a significant flux in the energy range from 0.3 GeV to 30 GeV, but their cumulative energy fluence is only 8% of that of all GRBs. Likewise, the GeV-to-MeV flux ratio is less and the GeV-band spectra are softer. We confirm that the GeV-band emission lasts much longer than the emission found in the keV energy range. The average allsky energy flux from GRBs in the GeV band is 6.4*10^-4 erg cm^-2 yr^-1 or only 4% of the energy flux of cosmic rays above the ankle at 10^18.6 eV.

preprint2012arXiv

Calibration strategies for the LAD instrument on-board LOFT

The Scientific objectives of the LOFT mission, e.g., the study of the Neutron Star equation of state and of the Strong Gravity, require accurate energy, time and flux calibration for the 500k channels of the SDD detectors, as well as the knowledge of the detector dead-time and of the detector response with respect to the incident angle of the photons. We report here the evaluations made to asses the calibration issues for the LAD instrument. The strategies for both ground and on-board calibrations, including astrophysical observations, show that the goals are achievable within the current technologies.

preprint2012arXiv

Could Cosmic Rays Affect Instabilities in the Transition Layer of Nonrelativistic Collisionless Shocks?

There is an observational correlation between astrophysical shocks and non-thermal particle distributions extending to high energies. As a first step toward investigating the possible feedback of these particles on the shock at the microscopic level, we perform particle-in-cell (PIC) simulations of a simplified environment consisting of uniform, interpenetrating plasmas, both with and without an additional population of cosmic rays. We vary the relative density of the counterstreaming plasmas, the strength of a homogeneous parallel magnetic field, and the energy density in cosmic rays. We compare the early development of the unstable spectrum for selected configurations without cosmic rays to the growth rates predicted from linear theory, for assurance that the system is well represented by the PIC technique. Within the parameter space explored, we do not detect an unambiguous signature of any cosmic-ray-induced effects on the microscopic instabilities that govern the formation of a shock. We demonstrate that an overly coarse distribution of energetic particles can artificially alter the statistical noise that produces the perturbative seeds of instabilities, and that such effects can be mitigated by increasing the density of computational particles.

preprint2012arXiv

Nonrelativistic parallel shocks in unmagnetized and weakly magnetized plasmas

We present results of 2D3V particle-in-cell simulations of non-relativistic plasma collisions with absent or parallel large-scale magnetic field for parameters applicable to the conditions at young supernova remnants. We study the collision of plasma slabs of different density, leading to two different shocks and a contact discontinuity. Electron dynamics play an important role in the development of the system. While non-relativistic shocks in both unmagnetized and magnetized plasmas can be mediated by Weibel-type instabilities, the efficiency of shock-formation processes is higher when a large-scale magnetic field is present. The electron distributions downstream of the forward and reverse shocks are generally isotropic, whereas that is not always the case for the ions. We do not see any significant evidence of pre-acceleration, neither in the electron population nor in the ion distribution.

preprint2012arXiv

The ASTRO-H X-ray Observatory

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions initiated by the Institute of Space and Astronautical Science (ISAS). ASTRO-H will investigate the physics of the high-energy universe via a suite of four instruments, covering a very wide energy range, from 0.3 keV to 600 keV. These instruments include a high-resolution, high-throughput spectrometer sensitive over 0.3-2 keV with high spectral resolution of Delta E < 7 eV, enabled by a micro-calorimeter array located in the focal plane of thin-foil X-ray optics; hard X-ray imaging spectrometers covering 5-80 keV, located in the focal plane of multilayer-coated, focusing hard X-ray mirrors; a wide-field imaging spectrometer sensitive over 0.4-12 keV, with an X-ray CCD camera in the focal plane of a soft X-ray telescope; and a non-focusing Compton-camera type soft gamma-ray detector, sensitive in the 40-600 keV band. The simultaneous broad bandpass, coupled with high spectral resolution, will enable the pursuit of a wide variety of important science themes.

preprint2012arXiv

The properties of non-thermal X-ray filaments in young supernova remnants

Context. Young supernova remnants (SNRs) exhibit narrow filaments of non-thermal X-ray emission whose widths can be limited either by electron energy losses or damping of the magnetic field. Aims. We want to investigate whether or not different models of these filaments can be observationally tested. Methods. Using observational parameters of four historical remnants, we calculate the filament profiles and compare the spectra of the filaments with those of the total non-thermal emission. For that purpose, we solve an one-dimensional stationary transport equation for the isotropic differential number density of the electrons. Results. We find that the difference between the spectra of filament and total non-thermal emission above 1 keV is more pronounced in the damping model than in the energy-loss model. Conclusions. A considerable damping of the magnetic field can result in an observable difference between the spectra of filament and total non-thermal emission, thus potentially permitting an observational discrimination between the energy-loss model and the damping model of the X-ray filaments.

preprint2012arXiv

Understanding TeV-band cosmic-ray anisotropy

We investigate the temporal and spectral correlations between flux and anisotropy fluctuations of TeV-band cosmic rays in the light of recent data taken with IceCube. We find that for a conventional distribution of cosmic-ray sources the dipole anisotropy is higher than observed, even if source discreteness is taken into account. Moreover, even for a shallow distribution of galactic cosmic-ray sources and a reacceleration model, fluctuations arising from source discreteness provide a probability only of the order of 10% that the cosmic-ray anisotropy limits of the recent IceCube analysis are met. This probability estimate is nearly independent of the exact choice of source rate, but generous for a large halo size. The location of the intensity maximum far from the Galactic Center is naturally reproduced.

preprint2011arXiv

Can Ultrahigh Energy Cosmic Rays Come from Gamma-Ray Bursts? II: Cosmic Rays Below the Ankle and Galactic GRB

The maximum cosmic ray energy achievable by acceleration by a relativistic blast wave is derived. It is shown that forward shocks from long GRB in the interstellar medium are powerful enough to produce the Galactic cosmic-ray component up to the ankle at $4\times 10^{18}$eV, as per an earlier suggestion (Levinson and Eichler, 1993). It is further argued that, were extragalactic long GRB responsible for the component {\it above} the ankle as well, the contribution from an occasional Galactic GRB within the solar circle would yield more than the observational limits on the outward flux from the solar circle, unless intermittency and/or beaming causes the present-day contribution to be less than $10^{-3}$ times the time average, and difficulties with these avoidance scenarios are also noted.

preprint2011arXiv

Origin of Ultrahigh Energy Galactic Cosmic Rays: The Isotropy Problem

We study the propagation of ultra-high-energy cosmic rays (UHECR) in the Galaxy, concentrating on the energy range below the ankle in the spectrum at 4 EeV. A Monte-Carlo method, based on analytical solutions to the time-dependent diffusion problem, is used to account for intermittency by placing sources at random locations. Assuming a source population that scales with baryon mass density or star formation (e.g. long GRB), we derive constraints arising from intermittency and the observational limits on the composition and anisotropy. It is shown that the composition and anisotropy at 1e18 eV are difficult to reproduce and require that either a) the particle mean free path is much smaller than a gyroradius, implying the escape time is very long, b) the composition is heavier than suggested by recent Auger data, c) the ultrahigh-energy sub-ankle component is mostly extragalactic, or d) we are living in a rare lull in the ultrahigh-energy cosmic-ray production, and the current ultrahigh-energy cosmic-ray intensity is far below the Galactic time average. We therefore recommend a strong observational focus on determining the UHECR composition around 1e18 eV.

preprint2011arXiv

Particle Spectra from Acceleration at Forward and Reverse Shocks of Young Type Ia Supernova Remnants

We study cosmic-ray acceleration in young Type Ia Supernova Remnants (SNRs) by means of test-particle diffusive shock acceleration theory and 1-D hydrodynamical simulations of their evolution. In addition to acceleration at the forward shock, we explore the particle acceleration at the reverse shock in the presence of a possible substantial magnetic field, and consequently the impact of this acceleration on the particle spectra in the remnant. We investigate the time evolution of the spectra for various time-dependent profiles of the magnetic field in the shocked region of the remnant. We test a possible influence on particle spectra of the Alfvénic drift of scattering centers in the precursor regions of the shocks. In addition, we study the radiation spectra and morphology in a broad band from radio to gamma-rays. It is demonstrated that the reverse shock contribution to the cosmic-ray particle population of young Type Ia SNRs may be significant, modifying the spatial distribution of particles and noticeably affecting the volume-integrated particle spectra in young SNRs. In particular spectral structures may arise in test-particle calculations that are often discussed as signatures of non-linear cosmic-ray modification of shocks. Therefore, the spectrum and morphology of emission, and their time evolution, differ from pure forward-shock solutions.

preprint2011arXiv

Systematic effects in the estimate of the local gamma-ray emissivity

We show in this letter that estimates of the local emissivity of γ-rays in the GeV-TeV range suffer uncertainties which are of the same order of magnitude as the current Fermi results. Primary cosmic-ray fluxes, cosmic-ray propagation, interstellar helium abundance and γ-ray production crosssections all affect the estimate of this quantity. We also show that the so-called nuclear enhancement factor -- though widely used so far to model the γ-ray emissivity -- is no longer a relevant quantity given the latest measurements of the primary cosmic ray proton and helium spectra.

preprint2011arXiv

The GeV-TeV Galactic gamma-ray diffuse emission I. Uncertainties in the predictions of the hadronic component

The Galactic gamma-ray diffuse emission is currently observed in the GeV-TeV energy range with unprecedented accuracy by the Fermi satellite. Understanding this component is crucial as it provides a background to many different signals such as extragalactic sources or annihilating dark matter. It is timely to reinvestigate how it is calculated and to assess the various uncertainties which are likely to affect the accuracy of the predictions. The Galactic gamma-ray diffuse emission is mostly produced above a few GeV by the interactions of cosmic ray primaries impinging on the interstellar material. The theoretical error on that component is derived by exploring various potential sources of uncertainty. Particular attention is paid to cosmic ray propagation. Nuclear cross sections, the proton and helium fluxes at the Earth, the Galactic radial profile of supernova remnants and the hydrogen distribution can also severely affect the signal. The propagation of cosmic ray species throughout the Galaxy is described in the framework of a semi-analytic two-zone diffusion/convection model. This allows to convert the constraints set by the boron-to-carbon data into a theoretical uncertainty on the diffuse emission. New deconvolutions of the HI and CO sky maps are also used to get the hydrogen distribution within the Galaxy. The thickness of the cosmic ray diffusive halo is found to have a significant effect on the Galactic gamma-ray diffuse emission while the interplay between diffusion and convection has little influence on the signal. The uncertainties related to nuclear cross sections and to the primary cosmic ray fluxes at the Earth are significant. The radial distribution of supernova remnants along the Galactic plane turns out to be a key ingredient. As expected, the predictions are extremely sensitive to the spatial distribution of hydrogen within the Milky Way.

preprint2010arXiv

Fermi Constrains Dark Matter Origin of High Energy Positron Anomaly

Fermi measurements of the high-latitude gamma-ray background strongly constrain a decaying-dark-matter origin for the 1--100 GeV Galactic positron anomaly measured with PAMELA. Inverse-Compton scattering of the microwave background by the emergent positrons produces a bump in the diffuse 100-200 MeV gamma-ray background that would protrude from the observed background at these energies. The positrons are thus constrained to emerge from the decay process at a typical energy between ~100 GeV and ~250 GeV. By considering only gamma-ray emission of the excess positrons and electrons, we derive a minimum diffuse gamma-ray flux that, apart from the positron spectrum assumed, is independent of the actual decay modes. Any gamma-rays produced directly by the dark-matter decay leads to an additional signal that make the observational limits more severe. A similar constraint on the energy of emergent positrons from annihilation in dark-matter substructures is argued to exist, according to recent estimates of enhancement in low-mass dark-matter substructures, and improved simulations of such substructure will further sharpen this constraint.

preprint2010arXiv

Magnetic Field Amplification by Turbulence in A Relativistic Shock Propagating through An Inhomogeneous Medium

We perform two-dimensional relativistic magnetohydrodynamic simulations of a mildly relativistic shock propagating through an inhomogeneous medium. We show that the postshock region becomes turbulent owing to preshock density inhomogeneity, and the magnetic field is strongly amplified due to the stretching and folding of field lines in the turbulent velocity field. The amplified magnetic field evolves into a filamentary structure in two-dimensional simulations. The magnetic energy spectrum is flatter than the Kolmogorov spectrum and indicates that a so-called small-scale dynamo is occurring in the postshock region. We also find that the amount of magnetic-field amplification depends on the direction of the mean preshock magnetic field, and the time scale of magnetic-field growth depends on the shock strength.

preprint2010arXiv

The ASTRO-H Mission

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions initiated by the Institute of Space and Astronautical Science (ISAS). ASTRO-H will investigate the physics of the high-energy universe by performing high-resolution, high-throughput spectroscopy with moderate angular resolution. ASTRO-H covers very wide energy range from 0.3 keV to 600 keV. ASTRO-H allows a combination of wide band X-ray spectroscopy (5-80 keV) provided by multilayer coating, focusing hard X-ray mirrors and hard X-ray imaging detectors, and high energy-resolution soft X-ray spectroscopy (0.3-12 keV) provided by thin-foil X-ray optics and a micro-calorimeter array. The mission will also carry an X-ray CCD camera as a focal plane detector for a soft X-ray telescope (0.4-12 keV) and a non-focusing soft gamma-ray detector (40-600 keV) . The micro-calorimeter system is developed by an international collaboration led by ISAS/JAXA and NASA. The simultaneous broad bandpass, coupled with high spectral resolution of Delta E ~7 eV provided by the micro-calorimeter will enable a wide variety of important science themes to be pursued.

preprint2010arXiv

The High Energy Budget Allocations in Shocks and GRB

The statistical distribution of energies among particles responsible for long Gamma Ray Burst (GRB) emission is analyzed in light of recent results of the Fermi Observatory. The allsky flux, $F_γ$, recorded by the Gamma Ray Burst Monitor (GBM) is shown, despite its larger energy range, to be not significantly larger than that reported by the Burst and Transient Explorer (BATSE), suggesting a relatively small flux in the 3 - 30 MeV energy range. The present-day energy input rate in $γ$-rays recorded by the GBM from long GRB is found, assuming star-formation rates in the literature, to be $\dot W(0)=0.5 F_γ H/c = 5 \times 10^{42}\ \rm{erg/Mpc^3 yr}$. The Large Area Telescope (LAT) fluence, when observed, is about 5-10\% per decade of the total, in good agreement with the predictions of saturated, non-linear shock acceleration. The high-energy component of long GRBs, as measured by Fermi, is found to contain only $\sim 10^{-2.5}$ of the energy needed to produce ultrahigh-energy cosmic rays (UHECR) above 4 Eev, assuming the latter to be extragalactic, when various numerical factors are carefully included, if the cosmic ray source spectrum has a spectral index of -2. The observed $γ$-ray fraction of the required UHECR energy is even smaller if the source spectrum is softer than $E^{-2}$. The AMANDA II limits rule out such a GRB origin for UHECR if much more than $10^{-2}$ of the cosmic ray energy goes into neutrinos that are within, and simultaneous with, the $γ$-ray beam. It is suggested that "orphan" neutrinos out of the $γ$-ray beam might be identifiable via orphan afterglow { or other wide angle signatures of GRB in lieu of coincidence with prompt $γ$-rays}, and it is recommended that feasible single neutrino trigger criteria be established to search for such coincidences.

preprint2009arXiv

Aperiodic magnetic turbulence produced by relativistic ion beams

Magnetic-field generation by a relativistic ion beam propagating through an electron-ion plasma along a homogeneous magnetic field is investigated with 2.5D high-resolution particle-in-cell (PIC) simulations. The studies test predictions of a strong amplification of short-wavelength modes of magnetic turbulence upstream of nonrelativistic and relativistic parallel shocks associated with supernova remnants, jets of active galactic nuclei, and gamma-ray bursts. We find good agreement in the properties of the turbulence observed in our simulations compared with the dispersion relation calculated for linear waves with arbitrary orientation of ${\vec k}$. Depending on the parameters, the backreaction on the ion beam leads to filamentation of the ambient plasma and the beam, which in turn influences the properties of the magnetic turbulence. For mildly- and ultra-relativistic beams, the instability saturates at field amplitudes a few times larger than the homogeneous magnetic field strength. This result matches our recent studies of nonrelativistically drifting, hot cosmic-ray particles upstream of supernova-remnant shocks which indicated only a moderate magnetic-field amplification by nonresonant instabilities. We also demonstrate that the aperiodic turbulence generated by the beam can provide efficient particle scattering with a rate compatible with Bohm diffusion. Representing the ion beam as a constant external current, i.e. excluding a backreaction of the magnetic turbulence on the beam, we observe non-resonant parallel modes with wavelength and growth rate as predicted by analytic calculations. In this unrealistic setup the magnetic field is amplified to amplitudes far exceeding the homogeneous field, as observed in recent MHD and PIC simulations.

preprint2009arXiv

Kinetic simulations of turbulent magnetic-field growth by streaming cosmic rays

Efficient acceleration of cosmic rays (via the mechanism of diffusive shock acceleration) requires turbulent, amplified magnetic fields in the shock's upstream region. We present results of multidimensional particle-in-cell simulations aimed at observing the magnetic field amplification that is expected to arise from the cosmic-ray current ahead of the shock, and the impact on the properties of the upstream interstellar medium. We find that the initial structure and peak strength of the amplified field is somewhat sensitive to the choice of parameters, but that the field growth saturates in a similar manner in all cases: the back-reaction on the cosmic rays leads to modification of their rest-frame distribution and also a net transfer of momentum to the interstellar medium, substantially weakening their relative drift while also implying the development of a modified shock. The upstream medium becomes turbulent, with significant spatial fluctuations in density and velocity, the latter in particular leading to moderate upstream heating; such fluctuations will also have a strong influence on the shock structure.

preprint2009arXiv

Radio polarimetry signatures of strong magnetic turbulence in Supernova Remnants

We discuss the emission and transport of polarized radio-band synchrotron radiation near the forward shocks of young shell-type supernova remnants, for which X-ray data indicate a strong amplification of turbulent magnetic field. Modeling the magnetic turbulence through the superposition of waves, we calculate the degree of polarization and the magnetic polarization direction which is at $90^\circ$ to the conventional electric polarization direction. We find that isotropic strong turbulence will produce weakly polarized radio emission even in the absence of internal Faraday rotation. If anisotropy is imposed on the magnetic-field structure, the degree of polarization can be significantly increased, provided internal Faraday rotation is inefficient. Both for shock compression and a mixture with a homogeneous field, the increase in polarization degree goes along with a fairly precise alignment of the magnetic-polarization angle with the direction of the dominant magnetic-field component, implying tangential magnetic polarization at the rims in the case of shock compression. We compare our model with high-resolution radio polarimetry data of Tycho's remnant. Using the absence of internal Faraday rotation we find a soft limit for the amplitude of magnetic turbulence, $δB \lesssim 200\ {\rm μG}$. The data are compatible with a turbulent magnetic field superimposed on a radial large-scale field of similar amplitude, $δB\simeq B_0$. An alternative viable scenario involves anisotropic turbulence with stronger amplitudes in the radial direction, as was observed in recent MHD simulations of shocks propagating through a medium with significant density fluctuations.

preprint2008arXiv

Cosmic-ray electron signatures of dark matter

There is evidence for an excess in cosmic-ray electrons at about 500 GeV energy, that may be related to dark-matter annihilation. I have calculated the expected electron contributions from a pulsar and from Kaluza-Klein dark matter, based on a realistic treatment of the electron propagation in the Galaxy. Both pulsars and dark-matter clumps are quasi-pointlike and few, and therefore their electron contributions at Earth generally have spectra that deviate from the average spectrum one would calculate for a smooth source distribution. I find that pulsars younger than about 10^5 years naturally cause a narrow peak at a few hundred GeV in the locally observed electron spectrum, similar to that observed. On the other hand, for a density n_c = 10 /kpc^3 of dark-matter clumps the sharp cut-off in the contribution from Kaluza-Klein particles is sometimes more pronounced, but often smoothed out and indistinguishable from a pulsar source, and therefore the spectral shape of the electron excess is insufficient to discriminate a dark-matter origin from more conventional astrophysical explanations. The amplitude of variations in the spectral feature caused by dark matter predominantly depends on the density of dark-matter clumps, which is not well known.

preprint1996arXiv

The galactic center arc as source of high-energy gamma-rays

We discuss the radio {\it arc} at the galactic center to be the counterpart of the high-energy \gr source 2EG J1746-2852. The observed \gr flux is constant within the statistical limits and the spectrum is very hard. The lack of variability makes it highly unlikely that any of the compact sources in the vicinity of the Galactic Center is the counterpart of 2EG J1746-2852. This includes the peculiar source Sgr A$^\ast$ at the very center of the Galaxy, which is often discussed to harbour a black hole of $10^6\, M_\odot$. Existing radio data on the {\it arc} support the view that its synchrotron emission originates from cooling, initially monoenergetic electrons which diffuse and convect from their sources to the outer extensions of the {\it arc}. If the source of high-energy electrons coincides with the sickle region (G0.18-0.04), as indicated by the radio data, then the ambient far- infrared photons can be up-scattered to \grs by inverse-Compton inter- action with the young high-energy electrons. We solve the continuity equation for the electrons including terms for diffusion, convection, monoenergetic injection, and the full energy loss. We show that both the flux and the spectrum of 2EG J1746-2852 can be well explained.