Institut für Physik und Astronomie
Refine
Year of publication
Document Type
- Article (45)
- Doctoral Thesis (1)
- Other (1)
Language
- English (47)
Is part of the Bibliography
- yes (47) (remove)
Keywords
- gamma rays: general (47) (remove)
Institute
We report the discovery of TeV gamma-ray emission from the Type Ia supernova remnant (SNR) G120.1+1.4, known as Tycho's SNR. Observations performed in the period 2008-2010 with the VERITAS ground-based gamma-ray observatory reveal weak emission coming from the direction of the remnant, compatible with a point source located at 00(h)25(m)27(s).0, +64 degrees 10'50 '' (J2000). The TeV photon spectrum measured by VERITAS can be described with a power law dN/dE = C(E/3.42 TeV)(-Gamma) with Gamma = 1.95 +/- 0.51(stat) +/- 0.30(sys) and C = (1.55 +/- 0.43(stat) +/- 0.47(sys)) x 10(-14) cm(-2) s(-1) TeV-1. The integral flux above 1 TeV corresponds to similar to 0.9% of the steady Crab Nebula emission above the same energy, making it one of the weakest sources yet detected in TeV gamma rays. We present both leptonic and hadronic models that can describe the data. The lowest magnetic field allowed in these models is similar to 80 mu G, which may be interpreted as evidence for magnetic field amplification.
Context. The large jet kinetic power and non-thermal processes occurring in the microquasar SS 433 make this source a good candidate for a very high-energy (VHE) gamma-ray emitter. Gamma-ray fluxes above the sensitivity limits of current Cherenkov telescopes have been predicted for both the central X-ray binary system and the interaction regions of SS 433 jets with the surrounding W50 nebula. Non-thermal emission at lower energies has been previously reported, indicating that efficient particle acceleration is taking place in the system. Aims. We explore the capability of SS 433 to emit VHE gamma rays during periods in which the expected flux attenuation due to periodic eclipses (P-orb similar to 13.1 days) and precession of the circumstellar disk (P-pre similar to 162 days) periodically covering the central binary system is expected to be at its minimum. The eastern and western SS 433/W50 interaction regions are also examined using the whole data set available. We aim to constrain some theoretical models previously developed for this system with our observations. Methods. We made use of dedicated observations from the Major Atmospheric Gamma Imaging Cherenkov telescopes (MAGIC) and High Energy Spectroscopic System (H.E.S.S.) of SS 433 taken from 2006 to 2011. These observation were combined for the first time and accounted for a total effective observation time of 16.5 h, which were scheduled considering the expected phases of minimum absorption of the putative VHE emission. Gamma-ray attenuation does not affect the jet/medium interaction regions. In this case, the analysis of a larger data set amounting to similar to 40-80 h, depending on the region, was employed. Results. No evidence of VHE gamma-ray emission either from the central binary system or from the eastern/western interaction regions was found. Upper limits were computed for the combined data set. Differential fluxes from the central system are found to be less than or similar to 10(-12)-10(-13) TeV-1 cm(-2) s(-1) in an energy interval ranging from similar to few x 100 GeV to similar to few TeV. Integral flux limits down to similar to 10(-12)-10(-13) ph cm(-2) s(-1) and similar to 10(-13)-10(-14) ph cm(-2) s(-1) are obtained at 300 and 800 GeV, respectively. Our results are used to place constraints on the particle acceleration fraction at the inner jet regions and on the physics of the jet/medium interactions. Conclusions. Our findings suggest that the fraction of the jet kinetic power that is transferred to relativistic protons must be relatively small in SS 433, q(p) <= 2.5 x 10(-5), to explain the lack of TeV and neutrino emission from the central system. At the SS 433/W50 interface, the presence of magnetic fields greater than or similar to 10 mu G is derived assuming a synchrotron origin for the observed X-ray emission. This also implies the presence of high-energy electrons with E-e up to 50 TeV, preventing an efficient production of gamma-ray fluxes in these interaction regions.
Relativistic pair beams produced in the cosmic voids by TeV gamma rays from blazars are expected to produce a detectable GeV-scale cascade emission missing in the observations. The suppression of this secondary cascade implies either the deflection of the pair beam by intergalactic magnetic fields (IGMFs) or an energy loss of the beam due to the electrostatic beam-plasma instability. IGMF of femto-Gauss strength is sufficient to significantly deflect the pair beams reducing the flux of secondary cascade below the observational limits. A similar flux reduction may result in the absence of the IGMF from the beam energy loss by the instability before the inverse Compton cooling. This dissertation consists of two studies about the instability role in the evolution of blazar-induced beams.
Firstly, we investigated the effect of sub-fG level IGMF on the beam energy loss by the instability. Considering IGMF with correlation lengths smaller than a few kpc, we found that such fields increase the transverse momentum of the pair beam particles, dramatically reducing the linear growth rate of the electrostatic instability and hence the energy-loss rate of the pair beam. Our results show that the IGMF eliminates beam plasma instability 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 and hence can be excluded.
Secondly, we probed how the beam-plasma instability feeds back on the beam, using a realistic two-dimensional beam distribution. We found that the instability broadens the beam opening angles significantly without any significant energy loss, thus confirming a recent feedback study on a simplified one-dimensional beam distribution. However, narrowing diffusion feedback of the beam particles with Lorentz factors less than 1e6 might become relevant even though initially it is negligible. Finally, when considering the continuous creation of TeV pairs, we found that the beam distribution and the wave spectrum reach a new quasi-steady state, in which the scattering of beam particles persists and the beam opening angle may increase by a factor of hundreds. This new intrinsic scattering of the cascade can result in time delays of around ten years, thus potentially mimicking the IGMF deflection. Understanding the implications on the GeV cascade emission requires accounting for inverse Compton cooling and simulating the beam-plasma system at different points in the IGM.
Context. Blazars are variable sources on various timescales over a broad energy range spanning from radio to very high energy (>100 GeV, hereafter VHE). Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. However, most of the gamma-ray studies performed on Mrk 501 during the past years relate to flaring activity, when the source detection and characterization with the available gamma-ray instrumentation was easier to perform.
Aims. Our goal is to characterize the source gamma-ray emission in detail, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity.
Methods. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. This provided extensive energy and temporal coverage of Mrk 501 throughout the entire campaign.
Results. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%-20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy and a tentative correlation between the X-ray and VHE fluxes. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density.
Conclusions. The one-zone SSC model can adequately describe the broadband spectral energy distribution of the source during the two months covered by the MW campaign. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during such a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.
A SEARCH FOR VERY HIGH ENERGY GAMMA RAYS FROM THE MISSING LINK BINARY PULSAR J1023+0038 WITH VERITAS
(2016)
The binary millisecond radio pulsar PSR J1023+0038 exhibits many characteristics similar to the gamma-ray binary system PSR B1259-63/LS 2883, making it an ideal candidate for the study of high-energy nonthermal emission. It has been the subject of multiwavelength campaigns following the disappearance of the pulsed radio emission in 2013 June, which revealed the appearance of an accretion disk around the neutron star. We present the results of very high energy (VHE) gamma-ray observations carried out by the Very Energetic Radiation Imaging Telescope Array System before and after this change of state. Searches for steady and pulsed emission of both data sets yield no significant gamma-ray signal above 100 GeV, and upper limits are given for both a steady and pulsed gamma-ray flux. These upper limits are used to constrain the magnetic field strength in the shock region of the PSR J1023+0038 system. Assuming that VHE gamma rays are produced via an inverse Compton mechanism in the shock region, we constrain the shock magnetic field to be greater than similar to 2 G before the disappearance of the radio pulsar and greater than similar to 10 G afterward.
We report on very high energy (E > 100 GeV) gamma-ray observations of V407 Cygni, a symbiotic binary that underwent a nova outburst producing 0.1-10 GeV gamma rays during 2010 March 10-26. Observations were made with the Very Energetic Radiation Imaging Telescope Array System during 2010 March 19-26 at relatively large zenith angles due to the position of V407 Cyg. An improved reconstruction technique for large zenith angle observations is presented and used to analyze the data. We do not detect V407 Cygni and place a differential upper limit on the flux at 1.6 TeV of 2.3 x 10(-12) erg cm(-2) s(-1) (at the 95% confidence level). When considered jointly with data from Fermi-LAT, this result places limits on the acceleration of very high energy particles in the nova.
We report on deep observations of the extended TeV gamma-ray source MGRO J1908+06 made with the VERITAS very high energy gamma-ray observatory. Previously, the TeV emission has been attributed to the pulsar wind nebula (PWN) of the Fermi-LAT pulsar PSR J1907+0602. We detect MGRO J1908+06 at a significance level of 14 standard deviations (14 sigma) and measure a photon index of 2.20 +/- 0.10(stat) +/- 0.20(sys). The TeV emission is extended, covering the region near PSR J1907+0602 and also extending toward SNR G40.5-0.5. When fitted with a two-dimensional Gaussian, the intrinsic extension has a standard deviation of sigma(src) = 0 degrees.44 +/- 0 degrees.02. In contrast to other TeV PWNe of similar age in which the TeV spectrum softens with distance from the pulsar, the TeV spectrum measured near the pulsar location is consistent with that measured at a position near the rim of G40.5-0.5, 0 degrees.33 away.
A search for enhanced very high energy GAMMA-RAY emission from the 2013 march crab nebula flare
(2014)
In 2013 March, a flaring episode from the Crab Nebula lasting similar to 2 weeks was detected by Fermi-LAT (Large Area Telescope on board the Fermi Gamma-ray Space Telescope). The Very Energetic Radiation Imaging Telescope Array System (VERITAS) provides simultaneous observations throughout this period. During the flare, Fermi-LAT detected a 20 fold increase in flux above the average synchrotron flux >100 MeV seen from the Crab Nebula. Simultaneous measurements with VERITAS are consistent with the non-variable long-term average Crab Nebula flux at TeV energies. Assuming a linear correlation between the very high energy flux change >1 TeV and the flux change seen in the Fermi-LAT band >100 MeV during the period of simultaneous observations, the linear correlation factor can be constrained to be at most 8.6 x 10(-3) with 95% confidence.
The Galactic Center ridge has been observed extensively in the past by both GeV and TeV gamma-ray instruments revealing a wealth of structure, including a diffuse component and the point sources G0.9+0.1 (a composite supernova remnant) and Sgr A* (believed to be associated with the supermassive black hole located at the center of our Galaxy). Previous very high energy (VHE) gamma-ray observations with the H.E.S.S.. experiment have also detected an extended TeV gamma-ray component along the Galactic plane in the >300 GeV gamma-ray regime. Here we report on observations of the Galactic Center ridge from 2010 to 2014 by the VERITAS telescope array in the >2 TeV energy range. From these observations we (1) provide improved measurements of the differential energy spectrum for Sgr A* in the >2 TeV gamma-ray regime, (2) provide a detection in the >2 TeV gamma-ray emission from the composite SNR G0.9+0.1 and an improved determination of its multi-TeV gamma-ray energy spectrum, and. (3) report on the detection of VER J1746-289, a localized enhancement of >2 TeV gamma-ray emission along the Galactic plane.
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.
We study the frequently used assumption in multi-messenger astrophysics that the gamma-ray and neutrino fluxes are directly connected because they are assumed to be produced by the same photohadronic production chain. An interesting candidate source for this test is the flat-spectrum radio quasar PKS B1424-418, which recently called attention to a potential correlation between an IceCube PeV neutrino event and its burst phase. We simulate both the multi-waveband photon and the neutrino emission from this source using a self-consistent radiation model. We demonstrate that a simple hadronic model cannot adequately describe the spectral energy distribution for this source, but a lepto-hadronic model with a subdominant hadronic component can reproduce the multi-waveband photon spectrum observed during various activity phases of the blazar. As a conclusion, up to about 0.3 neutrino events may coincide with the burst, which implies that the leptonic contribution dominates in the relevant energy band. We also demonstrate that the time-wise correlation between the neutrino event and burst phase is weak.
Context. The High Energy Stereoscopic System Galactic plane survey (HGPS) is to date the most comprehensive census of Galactic gamma -ray sources at very high energies (VHE; 100 GeV <= E <= 100 TeV). As a consequence of the limited sensitivity of this survey, the 78 detected gamma -ray sources comprise only a small and biased subsample of the overall population. The larger part consists of currently unresolved sources, which contribute to large-scale diffuse emission to a still uncertain amount.Aims. We study the VHE gamma -ray source population in the Milky Way. For this purpose population-synthesis models are derived based on the distributions of source positions, extents, and luminosities.Methods. Several azimuth-symmetric and spiral-arm models are compared for spatial source distribution. The luminosity and radius function of the population are derived from the source properties of the HGPS data set and are corrected for the sensitivity bias of the HGPS. Based on these models, VHE source populations are simulated and the subsets of sources detectable according to the HGPS are compared with HGPS sources.Results. The power-law indices of luminosity and radius functions are determined to range between -1.6 and -1.9 for luminosity and -1.1 and -1.6 for radius. A two-arm spiral structure with central bar is discarded as spatial distribution of VHE sources, while azimuth-symmetric distributions and a distribution following a four-arm spiral structure without bar describe the HGPS data reasonably well. The total number of Galactic VHE sources is predicted to be in the range from 800 to 7000 with a total luminosity and flux of (1.6-6.3) x 10(36) ph s(-1) and (3-15) x 10(-10) ph cm(-2) s(-1), respectively.Conclusions. Depending on the model, the HGPS sample accounts for (68-87)% of the emission of the population in the scanned region. This suggests that unresolved sources represent a critical component of the diffuse emission measurable in the HGPS. With the foreseen jump in sensitivity of the Cherenkov Telescope Array, the number of detectable sources is predicted to increase by a factor between 5 and 9.
We report the discovery of an unidentified, extended source of very-high-energy gamma-ray emission, VER J2019+407, within the radio shell of the supernova remnant SNR G78.2+2.1, using 21.4 hr of data taken by the VERITAS gamma-ray observatory in 2009. These data confirm the preliminary indications of gamma-ray emission previously seen in a two-year (2007-2009) blind survey of the Cygnus region by VERITAS. VER J2019+407, which is detected at a post-trials significance of 7.5 standard deviations in the 2009 data, is localized to the northwestern rim of the remnant in a region of enhanced radio and X-ray emission. It has an intrinsic extent of 0 degrees.23 +/- 0 degrees.03(stat-0 degrees.02sys)(+0 degrees.04) and its spectrum is well-characterized by a differential power law (dN/dE = N-0 x (E/TeV)-Gamma) with a photon index of Gamma = 2.37 +/- 0.14(stat) +/- 0.20(sys) and a flux normalization of N-0 = 1.5 +/- 0.2(stat) +/- 0.4(sys) x 10(-12) photon TeV-1 cm(-2) s(-1). This yields an integral flux of 5.2 +/- 0.8(stat) +/- 1.4(sys) x 10(-12) photon cm(-2) s(-1) above 320 GeV, corresponding to 3.7% of the Crab Nebula flux. We consider the relationship of the TeV gamma-ray emission with the GeV gamma-ray emission seen from SNR G78.2+2.1 as well as that seen from a nearby cocoon of freshly accelerated cosmic rays. Multiple scenarios are considered as possible origins for the TeV gamma-ray emission, including hadronic particle acceleration at the SNR shock.
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.
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.