Refine
Year of publication
Document Type
- Article (114)
- Other (5)
- Postprint (2)
- Review (2)
- Monograph/Edited Volume (1)
- Conference Proceeding (1)
Is part of the Bibliography
- yes (125)
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 (125) (remove)
VERITAS and Fermi-LAT Observations of TeV Gamma-Ray Sources Discovered by HAWC in the 2HWC Catalog
(2018)
The High Altitude Water Cherenkov (HAWC) collaboration recently published their 2HWC catalog, listing 39 very high energy (VHE; >100 GeV) gamma-ray sources based on 507 days of observation. Among these, 19 sources are not associated with previously known teraelectronvolt (TeV) gamma-ray sources. We have studied 14 of these sources without known counterparts with VERITAS and Fermi-LAT. VERITAS detected weak gamma-ray emission in the 1 TeV-30 TeV band in the region of DA 495, a pulsar wind nebula coinciding with 2HWC J1953+294, confirming the discovery of the source by HAWC. We did not find any counterpart for the selected 14 new HAWC sources from our analysis of Fermi-LAT data for energies higher than 10 GeV. During the search, we detected gigaelectronvolt (GeV) gamma-ray emission coincident with a known TeV pulsar wind nebula, SNR G54.1+0.3 (VER J1930+188), and a 2HWC source, 2HWC J1930+188. The fluxes for isolated, steady sources in the 2HWC catalog are generally in good agreement with those measured by imaging atmospheric Cherenkov telescopes. However, the VERITAS fluxes for SNR G54.1+0.3, DA 495, and TeV J2032+4130 are lower than those measured by HAWC, and several new HAWC sources are not detected by VERITAS. This is likely due to a change in spectral shape, source extension, or the influence of diffuse emission in the source region.
We report on observations of the pulsar/Be star binary system PSR J2032+4127/MT91 213 in the energy range between 100 GeV and 20 TeV with the Very Energetic Radiation Imaging Telescope Array and Major Atmospheric Gamma Imaging Cherenkov telescope arrays. The binary orbit has a period of approximately 50 years, with the most recent periastron occurring on 2017 November 13. Our observations span from 18 months prior to periastron to one month after. A new point-like gamma-ray source is detected, coincident with the location of PSR J2032+4127/MT91 213. The gamma-ray light curve and spectrum are well characterized over the periastron passage. The flux is variable over at least an order of magnitude, peaking at periastron, thus providing a firm association of the TeV source with the pulsar/Be star system. Observations prior to periastron show a cutoff in the spectrum at an energy around 0.5 TeV. This result adds a new member to the small population of known TeV binaries, and it identifies only the second source of this class in which the nature and properties of the compact object are firmly established. We compare the gamma-ray results with the light curve measured with the X-ray Telescope on board the Neil Gehrels Swift Observatory and with the predictions of recent theoretical models of the system. We conclude that significant revision of the models is required to explain the details of the emission that we have observed, and we discuss the relationship between the binary system and the overlapping steady extended source, TeV J2032+4130.
Vela Jr. (RX J0852.0-4622) is one of just a few known supernova remnants (SNRs) with a resolved shell across the whole electromagnetic spectrum from radio to very-high-energy (>100 GeV; VHE) gamma-rays. Its proximity and large size allow for detailed spatially resolved observations of the source, making Vela Jr. one of the primary sources used for the study of particle acceleration and emission mechanisms in SNRs. High-resolution X-ray observations reveal a steepening of the spectrum toward the interior of the remnant. In this study we aim for a self-consistent radiation model of Vela Jr. which at the same time would explain the broadband emission from the source and its intensity distribution. We solve the full particle transport equation combined with the high-resolution one-dimensional (1D) hydrodynamic simulations (using Pluto code) and subsequently calculate the radiation from the remnant. The equations are solved in the test particle regime. We test two models for the magnetic field profile downstream of the shock: damped magnetic field, which accounts for the damping of strong magnetic turbulence downstream, and transported magnetic field. Neither of these scenarios can fully explain the observed radial dependence of the X-ray spectrum under spherical symmetry. We show, however, that the softening of the spectrum and the X-ray intensity profile can be explained under the assumption that the emission is enhanced within a cone.
This is a correction notice for ‘Post-adiabatic supernova remnants in an interstellar magnetic field: oblique shocks and non-uniform environment’ (DOI: https://doi.org/10.1093/mnras/sty1750), which was published in MNRAS 479, 4253–4270 (2018). The publisher regrets to inform that the colour was missing from the colour scales in Figs 8(a)–(d) and Figs 9(a) and (b). This has now been corrected online. The publisher apologizes for this error.
In September 2017, the IceCube Neutrino Observatory recorded a very-high-energy neutrino in directional coincidence with a blazar in an unusually bright gamma-ray state, TXS0506 + 056 (refs(1,2)). Blazars are prominent photon sources in the Universe because they harbour a relativistic jet whose radiation is strongly collimated and amplified. High-energy atomic nuclei known as cosmic rays can produce neutrinos; thus, the recent detection may help in identifying the sources of the diffuse neutrino flux(3) and the energetic cosmic rays. Here we report a self-consistent analysis of the physical relation between the observed neutrino and the blazar, in particular the time evolution and spectral behaviour of neutrino and photon emission. We demonstrate that a moderate enhancement in the number of cosmic rays during the flare can yield a very strong increase in the neutrino flux, which is limited by co-produced hard X-rays and teraelectronvolt gamma rays. We also test typical radiation models(4,5) for compatibility and identify several model classes(6,7) as incompatible with the observations. We investigate to what degree the findings can be generalized to the entire population of blazars, determine the relation between their output in photons, neutrinos and cosmic rays, and suggest how to optimize the strategy of future observations.
Many solar wind observations at 1 au indicate that the proton (as well as electron) temperature anisotropy is limited. The data distribution in the (A(a), beta(a),(parallel to))-plane have a rhombic-shaped form around beta(a),(parallel to) similar to 1. The boundaries of the temperature anisotropy at beta(a),(parallel to) > 1 can be well explained by the threshold conditions of the mirror (whistler) and oblique proton (electron) firehose instabilities in a bi-Maxwellian plasma, whereas the physical mechanism of the similar restriction at beta(a),(parallel to) < 1 is still under debate. One possible option is Coulomb collisions, which we revisit in the current work. We derive the relaxation rate nu(A)(aa) of the temperature anisotropy in a bi-Maxwellian plasma that we then study analytically and by observed proton data from WIND. We found that nu(A)(pp) increases toward small beta(p),(parallel to) < 1. We matched the data distribution in the (A(p), beta(p),(parallel to))-plane with the constant contour nu(A)(pp) = 2.8 . 10(-6) s(-1), corresponding to the minimum value for collisions to play a role. This contour fits rather well the left boundary of the rhombic-shaped data distribution in the (A(p), beta(p),(parallel to))-plane. Thus, Coulomb collisions are an interesting candidate for explaining the limitations of the temperature anisotropy in the solar wind with small beta(a),(parallel to) < 1 at 1 au.
The particle-in-cell (PIC) method was developed to investigate microscopic phenomena, and with the advances in computing power, newly developed codes have been used for several fields, such as astrophysical, magnetospheric, and solar plasmas. PIC applications have grown extensively, with large computing powers available on supercomputers such as Pleiades and Blue Waters in the US. For astrophysical plasma research, PIC methods have been utilized for several topics, such as reconnection, pulsar dynamics, non-relativistic shocks, relativistic shocks, and relativistic jets. PIC simulations of relativistic jets have been reviewed with emphasis placed on the physics involved in the simulations. This review summarizes PIC simulations, starting with the Weibel instability in slab models of jets, and then focuses on global jet evolution in helical magnetic field geometry. In particular, we address kinetic Kelvin-Helmholtz instabilities and mushroom instabilities.
The variation in the intensity of cosmic rays at small angular scales is attributed to the interstellar turbulence in the vicinity of the Solar system. We show that a turbulent origin of the small-scale structures implies that the morphology of the observed cosmic ray intensity skymap varies with our location in the interstellar turbulence. The gyroradius of cosmic rays is shown to be the length scale associated with an observable change in the skymap over a radian angular scale. The extent to which the intensity at a certain angular scale varies is proportional to the change in our location with a maximum change of about the amplitude of intensity variation at that scale in the existing skymap. We suggest that for TeV cosmic rays a measurable variation could occur over a time-scale of a decade due to the Earth’s motion through the interstellar medium, if interstellar turbulence persists down to the gyroradius, about 300 μpc for TeV-ish cosmic rays. Observational evidence of the variability, or an absence of it, could provide a useful insight into the physical origin of the small-scale anisotropy.
We revisit the effect of nonlinear Landau (NL) damping on the electrostatic instability of blazar-induced pair beams, using a realistic pair-beam distribution. We employ a simplified 2D model in k-space to study the evolution of the electric-field spectrum and to calculate the relaxation time of the beam. We demonstrate that the 2D model is an adequate representation of the 3D physics. We find that nonlinear Landau damping, once it operates efficiently, transports essentially the entire wave energy to small wave numbers where wave driving is weak or absent. The relaxation time also strongly depends on the intergalactic medium temperature, T-IGM, and for T-IGM << 10 eV, and in the absence of any other damping mechanism, the relaxation time of the pair beam is longer than the inverse Compton (IC) scattering time. The weak late-time beam energy losses arise from the accumulation of wave energy at small k, that nonlinearly drains the wave energy at the resonant k of the pair-beam instability. Any other dissipation process operating at small k would reduce that wave-energy drain and hence lead to stronger pair-beam energy losses. As an example, collisions reduce the relaxation time by an order of magnitude, although their rate is very small. Other nonlinear processes, such as the modulation instability, could provide additional damping of the nonresonant waves and dramatically reduce the relaxation time of the pair beam. An accurate description of the spectral evolution of the electrostatic waves is crucial for calculating the relaxation time of the pair beam.
A Search for Pulsed Very High-energy Gamma-Rays from 13 Young Pulsars in Archival VERITAS Data
(2019)
We conduct a search for periodic emission in the very high-energy (VHE) gamma-ray band (E > 100 GeV) from a total of 13 pulsars in an archival VERITAS data set with a total exposure of over 450 hr. The set of pulsars includes many of the brightest young gamma-ray pulsars visible in the Northern Hemisphere. The data analysis resulted in nondetections of pulsed VHE gamma-rays from each pulsar. Upper limits on a potential VHE gamma-ray flux are derived at the 95% confidence level above three energy thresholds using two methods. These are the first such searches for pulsed VHE emission from each of the pulsars, and the obtained limits constrain a possible flux component manifesting at VHEs as is seen for the Crab pulsar.
The angular size of a star is a critical factor in determining its basic properties1. Direct measurement of stellar angular diameters is difficult: at interstellar distances stars are generally too small to resolve by any individual imaging telescope. This fundamental limitation can be overcome by studying the diffraction pattern in the shadow cast when an asteroid occults a star2, but only when the photometric uncertainty is smaller than the noise added by atmospheric scintillation3. Atmospheric Cherenkov telescopes used for particle astrophysics observations have not generally been exploited for optical astronomy due to the modest optical quality of the mirror surface. However, their large mirror area makes them well suited for such high-time-resolution precision photometry measurements4. Here we report two occultations of stars observed by the Very Energetic Radiation Imaging Telescope Array System (VERITAS)5 Cherenkov telescopes with millisecond sampling, from which we are able to provide a direct measurement of the occulted stars’ angular diameter at the ≤0.1 mas scale. This is a resolution never achieved before with optical measurements and represents an order of magnitude improvement over the equivalent lunar occultation method6. We compare the resulting stellar radius with empirically derived estimates from temperature and brightness measurements, confirming the latter can be biased for stars with ambiguous stellar classifications.
Electron injection at high Mach number nonrelativistic perpendicular shocks is studied here for parameters that are applicable to young SNR shocks. Using high-resolution large-scale two-dimensional fully kinetic particle-in-cell simulations and tracing individual particles, we in detail analyze the shock-surfing acceleration (SSA) of electrons at the leading edge of the shock foot. The central question is to what degree the process can be captured in 2D3V simulations. We find that the energy gain in SSA always arises from the electrostatic field of a Buneman wave. Electron energization is more efficient in the out-of-plane orientation of the large-scale magnetic field because both the phase speed and the amplitude of the waves are higher than for the in-plane scenario. Also, a larger number of electrons is trapped by the waves compared to the in-plane configuration. We conclude that significant modifications of the simulation parameters are needed to reach the same level of SSA efficiency as in simulations with out-of-plane magnetic field or 3D simulations.
Preliminary study for the laboratory experiment of cosmic-rays driven magnetic field amplification
(2019)
To understand astrophysical magnetic-field amplification, we conducted a feasibility study for a laboratory experiment of a non-resonant streaming instability at the Photo Injector Test Facility at DESY, Zeuthen site (PITZ). This non-resonant streaming instability, also known as Bell’s instability, is generally regarded as a candidate for the amplification of interstellar magnetic field in the upstream region of supernova-remnant shocks, which is crucial for the efficiency of diffusive shock acceleration. In the beam-plasma system composed of a radio-frequency electron gun and a gas-discharge plasma cell, the goal of our experiment is to demonstrate the development of the non-resonant streaming instability and to find its saturation level in the laboratory environment. Since we find that the electron beam will be significantly decelerated on account of an electrostatic streaming instability, which will decrease the growth rate of desired non-resonant streaming instability, we discuss possible ways to suppress the electrostatic streaming instability by considering the characteristics of a field-emission-based quasi continuous-wave electron beam.
Nonthermal emission from the reverse shock of the youngest galactic supernova remnant G1.9+0.3
(2019)
Context. The youngest Galactic supernova remnant G1.9+0.3 is an interesting target for next-generation gamma-ray observatories. So far, the remnant is only detected in the radio and the X-ray bands, but its young age of approximate to 100 yr and inferred shock speed of approximate to 14 000 km s(-1) could make it an efficient particle accelerator. Aims. We aim to model the observed radio and X-ray spectra together with the morphology of the remnant. At the same time, we aim to estimate the gamma-ray flux from the source and evaluate the prospects of its detection with future gamma-ray experiments. Methods. We performed spherical symmetric 1D simulations with the RATPaC code, in which we simultaneously solved the transport equation for cosmic rays, the transport equation for magnetic turbulence, and the hydro-dynamical equations for the gas flow. Separately computed distributions of the particles accelerated at the forward and the reverse shock were then used to calculate the spectra of synchrotron, inverse Compton, and pion-decay radiation from the source. Results. The emission from G1.9+0.3 can be self-consistently explained within the test-particle limit. We find that the X-ray flux is dominated by emission from the forward shock while most of the radio emission originates near the reverse shock, which makes G1.9+0.3 the first remnant with nonthermal radiation detected from the reverse shock. The flux of very-high-energy gamma-ray emission from G1.9+0.3 is expected to be close to the sensitivity threshold of the Cherenkov Telescope Array. The limited time available to grow large-scale turbulence limits the maximum energy of particles to values below 100 TeV, hence G1.9+0.3 is not a PeVatron.
The extragalactic background light (EBL), a diffuse photon field in the optical and infrared range, is a record of radiative processes over the universe?s history. Spectral measurements of blazars at very high energies (>100 GeV) enable the reconstruction of the spectral energy distribution (SED) of the EBL, as the blazar spectra are modified by redshift- and energy-dependent interactions of the gamma-ray photons with the EBL. The spectra of 14 VERITAS-detected blazars are included in a new measurement of the EBL SED that is independent of EBL SED models. The resulting SED covers an EBL wavelength range of 0.56?56 ?m, and is in good agreement with lower limits obtained by assuming that the EBL is entirely due to radiation from cataloged galaxies.
Precursor Wave Amplification by Ion-Electron Coupling through Wakefield in Relativistic Shocks
(2019)
We investigated electromagnetic precursor wave emission in relativistic shocks by using two-dimensional particle-in-cell simulations. We found that the wave amplitude is significantly enhanced by a positive feedback process associated with ion-electron coupling through the wakefields for high magnetization. The wakefields collapse during the nonlinear process of the parametric decay instability in the near-upstream region, where nonthermal electrons and ions are generated. The intense coherent emission and the particle acceleration may operate in high-energy astrophysical objects.
We explore electron preacceleration at high-Mach-number nonrelativistic perpendicular shocks at, e.g., young supernova remnants, which are a prerequisite of further acceleration to very high energies via diffusive shock acceleration. Using fully kinetic particle-in-cell simulations of shocks and electron dynamics in them, we investigate the influence of shock-surfing acceleration (SSA) at the shock foot on the nonthermal population of electrons downstream of the shock. The SSA is followed by further energization at the shock ramp where the Weibel instability spawns a type of second-order Fermi acceleration. The combination of these two processes leads to the formation of a nonthermal electron population, but the importance of SSA becomes smaller for larger ion-to-electron mass ratios in the simulation. We discuss the resulting electron spectra and the relevance of our results to the physics of systems with real ion-to-electron mass ratios and fully three-dimensional behavior.
Aims. We aim to characterize the multiwavelength emission from Markarian 501 (Mrk 501), quantify the energy-dependent variability, study the potential multiband correlations, and describe the temporal evolution of the broadband emission within leptonic theoretical scenarios. Methods. We organized a multiwavelength campaign to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Results. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of similar to 0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was similar to 3 CU, and the peak of the high-energy spectral component was found to be at similar to 2 TeV. Both the X-ray and VHE gamma-ray spectral slopes were measured to be extremely hard, with spectral indices <2 during most of the observing campaign, regardless of the X-ray and VHE flux. This study reports the hardest Mrk 501 VHE spectra measured to date. The fractional variability was found to increase with energy, with the highest variability occurring at VHE. Using the complete data set, we found correlation between the X-ray and VHE bands; however, if the June 9 flare is excluded, the correlation disappears (significance <3 sigma) despite the existence of substantial variability in the X-ray and VHE bands throughout the campaign. Conclusions. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency-peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The data set acquired shows that the broadband spectral energy distribution (SED) of Mrk 501, and its transient evolution, is very complex, requiring, within the framework of synchrotron self-Compton (SSC) models, various emission regions for a satisfactory description. Nevertheless the one-zone SSC scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behavior seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.