Institut für Physik und Astronomie
Refine
Has Fulltext
- no (112)
Year of publication
Document Type
- Article (112) (remove)
Language
- English (112)
Is part of the Bibliography
- yes (112)
Keywords
- acceleration of particles (20)
- ISM: supernova remnants (17)
- cosmic rays (17)
- gamma rays: galaxies (17)
- gamma rays: general (15)
- radiation mechanisms: non-thermal (12)
- shock waves (11)
- galaxies: active (10)
- BL Lacertae objects: general (8)
- galaxies: jets (8)
Institute
- Institut für Physik und Astronomie (112) (remove)
Thermal electrons have gyroradii many orders of magnitude smaller than the finite width of a shock, thus need to be pre-accelerated before they can cross it and be accelerated by diffusive shock acceleration. One region where pre-acceleration may occur is the inner foreshock, which upstream electrons must pass through before any potential downstream crossing. In this paper, we perform a large-scale particle-in-cell simulation that generates a single shock with parameters motivated from supernova remnants. Within the foreshock, reflected electrons excite the oblique whistler instability and produce electromagnetic whistler waves, which comove with the upstream flow and as nonlinear structures eventually reach radii of up to 5 ion-gyroradii. We show that the inner electromagnetic configuration of the whistlers evolves into complex nonlinear structures bound by a strong magnetic field around four times the upstream value. Although these nonlinear structures do not in general interact with cospatial upstream electrons, they resonate with electrons that have been reflected at the shock. We show that they can scatter, or even trap, reflected electrons, confining around 0.8% of the total upstream electron population to the region close to the shock where they can undergo substantial pre-acceleration. This acceleration process is similar to, yet approximately three times more efficient than, stochastic shock drift acceleration.
Assessing the impact of hydrogen absorption on the characteristics of the Galactic center excess
(2022)
We present a new reconstruction of the distribution of atomic hydrogen in the inner Galaxy that is based on explicit radiation transport modeling 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 HI4PI 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 of a dark-matter signal, be it cuspy or cored. But an additional so-called boxy bulge is still favored 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.
The current paradigm of cosmic-ray (CR) origin states that the greater part of galactic CRs is produced by supernova remnants. The interaction of supernova ejecta with the interstellar medium after a 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. 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 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 to be triggered at high Mach number shocks.
Core-collapse supernova remnants are the gaseous nebulae of galactic interstellar media (ISM) formed after the explosive death of massive stars. Their morphology and emission properties depend both on the surrounding circumstellar structure shaped by the stellar wind-ISM interaction of the progenitor star and on the local conditions of the ambient medium. In the warm phase of the Galactic plane (n approximate to 1 cm(-3), T approximate to 8000 K), an organized magnetic field of strength 7 mu G has profound consequences on the morphology of the wind bubble of massive stars at rest. In this paper, we show through 2.5D magnetohydrodynamical simulations, in the context of a Wolf-Rayet-evolving 35 M 0 star, that it affects the development of its supernova remnant. When the supernova remnant reaches its middle age (15-20 kyr), it adopts a tubular shape that results from the interaction between the isotropic supernova ejecta and the anisotropic, magnetized, shocked stellar progenitor bubble into which the supernova blast wave expands. Our calculations for non-thermal emission, i.e. radio synchrotron and inverse-Compton radiation, reveal that such supernova remnants can, due to projection effects, appear as rectangular objects in certain cases. This mechanism for shaping a supernova remnant is similar to the bipolar and elliptical planetary nebula production by wind-wind interaction in the low-mass regime of stellar evolution. If such a rectangular core-collapse supernova remnant is created, the progenitor star must not have been a runaway star. We propose that such a mechanism is at work in the shaping of the asymmetric core-collapse supernova remnant Puppis A.
To undergo diffusive shock acceleration, electrons need to be preaccelerated to increase their energies by several orders of magnitude, else their gyroradii will be smaller than the finite width of the shock. In oblique shocks, where the upstream magnetic field orientation is neither parallel nor 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 preacceleration in this region, we undertake particle-in-cell simulations of oblique shocks while varying the obliquity and in-plane angles. We show that while the proportion of reflected electrons is negligible for theta (Bn) = 74.degrees 3, it increases to R similar to 5% for theta (Bn) = 30 degrees, and that, via the electron acoustic instability, these electrons power electrostatic waves upstream with energy density proportional to R (0.6) and a wavelength approximate to 2 lambda (se), where lambda (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 less than or similar to 1% of cases, upstream electrons are prematurely turned away from the shock and never injected downstream. In contrast, a similar fraction is rescattered back toward the shock after reflection, reinteracts with the shock with energies much greater than thermal, and crosses into the downstream.
Cosmic-ray acceleration at non-relativistic shocks relies on scattering by turbulence that the cosmic rays drive upstream of the shock. We explore the rate of energy transfer from cosmic rays to non-resonant Bell modes and the spectral softening it implies. Accounting for the finite time available for turbulence driving at supernova-remnant shocks yields a smaller spectral impact than found earlier with steady-state considerations. Generally, for diffusion scaling with the Bohm rate by a factor eta, the change in spectral index is at most eta divided by the Alfvenic Mach number of the thermal sub-shock. For M (A) less than or similar to 50 it is well below this limit. Only for very fast shocks and very efficient cosmic-ray acceleration can the change in spectral index reach 0.1. For standard SNR parameters, it is negligible. Independent confirmation is derived by considering the synchrotron energy losses of electrons: if intense nonthermal multi-keV emission is produced, the energy loss, and hence the spectral steepening, is very small for hadronic cosmic rays that produce TeV-band gamma-ray emission.
Suppression of the TeV Pair-beam-Plasma Instability by a Tangled Weak Intergalactic Magnetic Field
(2022)
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 kiloparsec, 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.
A very small fraction of (runaway) massive stars have masses exceeding 60-70 M-circle dot and are predicted to evolve as luminous blue variable and Wolf-Rayet stars before ending their lives as core-collapse supernovae. Our 2D axisymmetric hydrodynamical simulations explore how a fast wind (2000 km s(-1)) and high mass-loss rate (10(-5)M(circle dot) yr(-1)) can impact the morphology of the circumstellar medium. It is shaped as 100 pc-scale wind nebula that can be pierced by the driving star when it supersonically moves with velocity 20-40 km s(-1) through the interstellar medium (ISM) in the Galactic plane. The motion of such runaway stars displaces the position of the supernova explosion out of their bow shock nebula, imposing asymmetries to the eventual shock wave expansion and engendering Cygnus-loop-like supernova remnants. We conclude that the size (up to more than 200 pc) of the filamentary wind cavity in which the chemically enriched supernova ejecta expand, mixing efficiently the wind and ISM materials by at least 10 per cent in number density, can be used as a tracer of the runaway nature of the very massive progenitors of such 0.1Myr old remnants. Our results motivate further observational campaigns devoted to the bow shock of the very massive stars BD+43 degrees 3654 and to the close surroundings of the synchrotron-emitting Wolf-Rayet shell G2.4+1.4.
Context.
Tycho's supernova remnant (SNR) is associated with the historical supernova (SN) event SN 1572 of Type Ia. The explosion occurred in a relatively clean environment, and was visually observed, providing an age estimate. This SNR therefore represents an ideal astrophysical test-bed for the study of cosmic-ray acceleration and related phenomena. A number of studies suggest that shock acceleration with particle feedback and very efficient magnetic-field amplification combined with Alfvenic drift are needed to explain the rather soft radio spectrum and the narrow rims observed in X-rays.
Aims.
We show that the broadband spectrum of Tycho's SNR can alternatively be well explained when accounting for stochastic acceleration as a secondary process. The re-acceleration of particles in the turbulent region immediately downstream of the shock should be efficient enough to impact particle spectra over several decades in energy. The so-called Alfvenic drift and particle feedback on the shock structure are not required in this scenario. Additionally, we investigate whether synchrotron losses or magnetic-field damping play a more profound role in the formation of the non-thermal filaments.
Methods.
We solved the full particle transport equation in test-particle mode using hydrodynamic simulations of the SNR plasma flow. The background magnetic field was either computed from the induction equation or follows analytic profiles, depending on the model considered. Fast-mode waves in the downstream region provide the diffusion of particles in momentum space.
Results.
We show that the broadband spectrum of Tycho can be well explained if magnetic-field damping and stochastic re-acceleration of particles are taken into account. Although not as efficient as standard diffusive shock acceleration, stochastic acceleration leaves its imprint on the particle spectra, which is especially notable in the emission at radio wavelengths. We find a lower limit for the post-shock magnetic-field strength similar to 330 mu G, implying efficient amplification even for the magnetic-field damping scenario. Magnetic-field damping is necessary for the formation of the filaments in the radio range, while the X-ray filaments are shaped by both the synchrotron losses and magnetic-field damping.
B2 1215+30 is a BL-Lac-type blazar that was first detected at TeV energies by the MAGIC atmospheric Cherenkov telescopes and subsequently confirmed by the Very Energetic Radiation Imaging Telescope Array System (VERITAS) observatory with data collected between 2009 and 2012. In 2014 February 08, VERITAS detected a large-amplitude flare from B2. 1215+30 during routine monitoring observations of the blazar 1ES. 1218+304, located in the same field of view. The TeV flux reached 2.4 times the Crab Nebula flux with a variability timescale of <3.6 hr. Multiwavelength observations with Fermi-LAT, Swift, and the Tuorla Observatory revealed a correlated high GeV flux state and no significant optical counterpart to the flare, with a spectral energy distribution where the gamma-ray luminosity exceeds the synchrotron luminosity. When interpreted in the framework of a onezone leptonic model, the observed emission implies a high degree of beaming, with Doppler factor delta > 10, and an electron population with spectral index p < 2.3.
Imaging atmospheric Cherenkov telescopes (IACTs) are equipped with sensitive photomultiplier tube (PMT) cameras. Exposure to high levels of background illumination degrades the efficiency of and potentially destroys these photo-detectors over time, so IACTs cannot be operated in the same configuration in the presence of bright moonlight as under dark skies. Since September 2012, observations have been carried out with the VERITAS IACTs under bright moonlight (defined as about three times the night-sky-background (NSB) of a dark extragalactic field, typically occurring when Moon illumination > 35%) in two observing modes, firstly by reducing the voltage applied to the PMTs and, secondly, with the addition of ultra-violet (UV) bandpass filters to the cameras. This has allowed observations at up to about 30 times previous NSB levels (around 80% Moon illumination), resulting in 30% more observing time between the two modes over the course of a year. These additional observations have already allowed for the detection of a flare from the 1ES 1727 + 502 and for an observing program targeting a measurement of the cosmic-ray positron fraction. We provide details of these new observing modes and their performance relative to the standard VERITAS observations. (C) 2017 Elsevier B.V. All rights reserved.
Prospects for Cherenkov Telescope Array Observations of the Young Supernova Remnant RX J1713.7-3946
(2017)
We perform simulations for future Cherenkov Telescope Array (CTA) observations of RX J1713.7-3946, a young supernova remnant (SNR) and one of the brightest sources ever discovered in very high energy (VHE) gamma rays. Special attention is paid to exploring possible spatial (anti) correlations of gamma rays with emission at other wavelengths, in particular X-rays and CO/H I emission. We present a series of simulated images of RX J1713.7-3946 for CTA based on a set of observationally motivated models for the gamma-ray emission. In these models, VHE gamma rays produced by high-energy electrons are assumed to trace the nonthermal X-ray emission observed by XMM-Newton, whereas those originating from relativistic protons delineate the local gas distributions. The local atomic and molecular gas distributions are deduced by the NANTEN team from CO and H I observations. Our primary goal is to show how one can distinguish the emission mechanism(s) of the gamma rays (i.e., hadronic versus leptonic, or a mixture of the two) through information provided by their spatial distribution, spectra, and time variation. This work is the first attempt to quantitatively evaluate the capabilities of CTA to achieve various proposed scientific goals by observing this important cosmic particle accelerator.
Combined with measurements made by very-long-baseline interferometry, the observations of fast TeV gamma-ray flares probe the structure and emission mechanism of blazar jets. However, only a handful of such flares have been detected to date, and only within the last few years have these flares been observed from lower-frequency-peaked BL. Lac objects and flat-spectrum radio quasars. We report on a fast TeV gamma-ray flare from the blazar BL. Lacertae observed by the Very Energetic Radiation Imaging Telescope Array System (VERITAS). with a rise time of similar to 2.3 hr and a decay time of similar to 36 min. The peak flux above 200 GeV is (4.2 +/- 0.6) x 10(-6) photon m(-2) s(-1) measured with a 4-minute-binned light curve, corresponding to similar to 180% of the flux that is observed from the Crab Nebula above the same energy threshold. Variability contemporaneous with the TeV gamma-ray flare was observed in GeV gamma-ray, X-ray, and optical flux, as well as in optical and radio polarization. Additionally, a possible moving emission feature with superluminal apparent velocity was identified in Very Long Baseline Array observations at 43 GHz, potentially passing the radio core of the jet around the time of the gamma-ray flare. We discuss the constraints on the size, Lorentz factor, and location of the emitting region of the flare, and the interpretations with several theoretical models that invoke relativistic plasma passing stationary shocks.
On 2015 March 23, the Very Energetic Radiation Imaging Telescope Array System (VERITAS) responded to a Swift-Burst Alert Telescope (BAT) detection of a gamma-ray burst, with observations beginning 270 s after the onset of BAT emission, and only 135 s after the main BAT emission peak. No statistically significant signal is detected above 140 GeV. The VERITAS upper limit on the fluence in a 40-minute integration corresponds to about 1% of the prompt fluence. Our limit is particularly significant because the very-high-energy (VHE) observation started only similar to 2 minutes after the prompt emission peaked, and Fermi-Large Area Telescope observations of numerous other bursts have revealed that the high-energy emission is typically delayed relative to the prompt radiation and lasts significantly longer. Also, the proximity of GRB 150323A (z = 0.593) limits the attenuation by the extragalactic background light to similar to 50% at 100-200 GeV. We conclude that GRB 150323A had an intrinsically very weak high-energy afterglow, or that the GeV spectrum had a turnover below similar to 100 GeV. If the GRB exploded into the stellar wind of a massive progenitor, the VHE non-detection constrains the wind density parameter to be A greater than or similar to 3 x 10(11) g . cm(-1), consistent with a standard Wolf-Rayet progenitor. Alternatively, the VHE emission from the blast wave would be weak in a very tenuous medium such as the interstellar medium, which therefore cannot be ruled out as the environment of GRB 150323A.
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.
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.
On 2017 September 22, the IceCube Neutrino Observatory reported the detection of the high-energy neutrino event IC 170922A, of potential astrophysical origin. It was soon determined that the neutrino direction was consistent with the location of the gamma-ray blazar TXS 0506+056. (3FGL J0509.4+ 0541), which was in an elevated gamma-ray emission state as measured by the Fermi satellite. Very Energetic Radiation Imaging Telescope Array System (VERITAS) observations of the neutrino/blazar region started on 2017 September 23 in response to the neutrino alert and continued through 2018 February 6. While no significant very-high-energy (VHE; E > 100 GeV) emission was observed from the blazar by VERITAS in the two-week period immediately following the IceCube alert, TXS 0506+ 056 was detected by VERITAS with a significance of 5.8 standard deviations (sigma) in the full 35 hr data set. The average photon flux of the source during this period was (8.9 +/- 1.6). x. 10(-12) cm(-2) s(-1), or 1.6% of the Crab Nebula flux, above an energy threshold of 110 GeV, with a soft spectral index of 4.8. +/-. 1.3.
We present results from deep observations toward the Cygnus region using 300 hr of very high energy (VHE)gamma-ray data taken with the VERITAS Cerenkov telescope array and over 7 yr of high-energy.-ray data taken with the Fermi satellite at an energy above 1 GeV. As the brightest region of diffuse gamma-ray emission in the northern sky, the Cygnus region provides a promising area to probe the origins of cosmic rays. We report the identification of a potential Fermi-LAT counterpart to VER J2031+415 (TeV J2032+4130) and resolve the extended VHE source VER J2019+368 into two source candidates (VER J2018+367* and VER J2020+368*) and characterize their energy spectra. The Fermi-LAT morphology of 3FGL J2021.0+4031e (the Gamma Cygni supernova remnant) was examined, and a region of enhanced emission coincident with VER J2019+407 was identified and jointly fit with the VERITAS data. By modeling 3FGL J2015.6+3709 as two sources, one located at the location of the pulsar wind nebula CTB 87 and one at the quasar QSO J2015+371, a continuous spectrum from 1 GeV to 10 TeV was extracted for VER J2016+371 (CTB 87). An additional 71 locations coincident with Fermi-LAT sources and other potential objects of interest were tested for VHE gamma-ray emission, with no emission detected and upper limits on the differential flux placed at an average of 2.3% of the Crab Nebula flux. We interpret these observations in a multiwavelength context and present the most detailed gamma-ray view of the region to date.
HESS J1943+213
(2018)
HESS J1943+213 is a very high energy (VHE; > 100 GeV) gamma-ray source in the direction of the Galactic plane. Studies exploring the classification of the source are converging toward its identification as an extreme synchrotron BL Lac object. Here we present 38 hr of VERITAS observations of HESS J1943+213 taken over 2 yr. The source is detected with a significance of similar to 20 standard deviations, showing a remarkably stable flux and spectrum in VHE gamma-rays. Multifrequency Very Long Baseline Array (VLBA) observations of the source confirm the extended, jet-like structure previously found in the 1.6 GHz band with the European VLBI Network and detect this component in the 4.6 and 7.3 GHz bands. The radio spectral indices of the core and the jet and the level of polarization derived from the VLBA observations are in a range typical for blazars. Data from VERITAS, Fermi-LAT, Swift-XRT, the FLWO 48 ' telescope, and archival infrared and hard X-ray observations are used to construct and model the spectral energy distribution (SED) of the source with a synchrotron self-Compton model. The well-measured gamma-ray peak of the SED with VERITAS and Fermi-LAT provides constraining upper limits on the source redshift. Possible contribution of secondary gamma-rays from ultra-high-energy cosmic-ray-initiated electromagnetic cascades to the gamma-ray emission is explored, finding that only a segment of the VHE spectrum can be accommodated with this process. A variability search is performed across X-ray and gamma-ray bands. No statistically significant flux or spectral variability is detected.