Refine
Has Fulltext
- no (122) (remove)
Year of publication
Document Type
- Article (114)
- Other (5)
- Review (2)
- Conference Proceeding (1)
Language
- English (122) (remove)
Is part of the Bibliography
- yes (122)
Keywords
- acceleration of particles (20)
- ISM: supernova remnants (18)
- cosmic rays (17)
- gamma rays: galaxies (17)
- gamma rays: general (15)
- radiation mechanisms: non-thermal (12)
- shock waves (12)
- galaxies: active (10)
- BL Lacertae objects: general (8)
- galaxies: jets (8)
Institute
- Institut für Physik und Astronomie (122) (remove)
Core-collapse supernova remnants are structures of the interstellar medium (ISM) left behind the explosive death of most massive stars ( ?40 M-?). Since they result in the expansion of the supernova shock wave into the gaseous environment shaped by the star's wind history, their morphology constitutes an insight into the past evolution of their progenitor star. Particularly, fast-mo ving massiv e stars can produce asymmetric core-collapse superno va remnants. We inv estigate the mixing of materials in core-collapse supernova remnants generated by a moving massive 35 M-? star, in a magnetized ISM. Stellar rotation and the wind magnetic field are time-dependently included into the models which follow the entire evolution of the stellar surroundings from the zero-age main-sequence to 80 kyr after the supernova explosion. It is found that very little main-sequence material is present in remnants from moving stars, that the Wolf-Rayet wind mixes very efficiently within the 10 kyr after the explosion, while the red supergiant material is still unmixed by 30 per cent within 50 kyr after the supernova. Our results indicate that the faster the stellar motion, the more complex the internal organization of the supernova remnant and the more ef fecti ve the mixing of ejecta therein. In contrast, the mixing of stellar wind material is only weakly affected by progenitor motion, if at all.
Leptonic nonthermal emission from supernova remnants evolving in the circumstellar magnetic field
(2022)
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 nonthermal 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 nonthermal emission.
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.
We develop a simple two-zone interpretation of the broadband baseline Crab nebula spectrum between 10(-5) eV and similar to 100 TeV by using two distinct log-parabola energetic electrons distributions. We determine analytically the very-high energy photon spectrum as originated by inverse-Compton scattering of the far-infrared soft ambient photons within the nebula off a first population of electrons energized at the nebula termination shock. The broad and flat 200 GeV peak jointly observed by Fermi/LAT and MAGIC is naturally reproduced. The synchrotron radiation from a second energetic electron population explains the spectrum from the radio range up to similar to 10 keV. We infer from observations the energy dependence of the microscopic probability of remaining in proximity of the shock of the accelerating electrons.
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.
e-ASTROGAM is a concept for a breakthrough observatory space mission carrying a gamma-ray telescope dedicated to the study of the non-thermal Universe in the photon energy range from 0.15 MeV to 3 GeV. The lower energy limit can be pushed down to energies as low as 30 keV for gamma-ray burst detection with the calorimeter. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with remarkable polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous and current generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will be a major player of the multiwavelength, multimessenger time-domain astronomy of the 2030s, and provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LISA, LIGO, Virgo, KAGRA, the Einstein Telescope and the Cosmic Explorer, IceCube-Gen2 and KM3NeT, SKA, ALMA, JWST, E-ELT, LSST, Athena, and the Cherenkov Telescope Array.
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.