Refine
Year of publication
Document Type
- Article (110)
- Other (5)
- Postprint (3)
- Review (2)
- Monograph/Edited Volume (1)
- Conference Proceeding (1)
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 (9)
- galaxies: jets (8)
Institute
X-RAY OBSERVATIONS OF BOW SHOCKS AROUND RUNAWAY O STARS. THE CASE OF zeta OPH AND BD+43 degrees 3654
(2016)
Non-thermal radiation has been predicted within bow shocks around runaway stars by recent theoretical works. We present X-ray observations toward the runaway stars zeta Oph by Chandra and Suzaku and of BD+43 degrees 3654 by XMM-Newton to search for the presence of non-thermal X-ray emission. We found no evidence of non-thermal emission spatially coincident with the bow shocks; nonetheless, diffuse emission was detected in the vicinity of zeta Oph. After a careful analysis of its spectral characteristics, we conclude that this emission has a thermal nature with a plasma temperature of T approximate to 2 x 10(6) K. The cometary shape of this emission seems to be in line with recent predictions of radiation-hydrodynamic models of runaway stars. The case of BD+43 degrees 3654 is puzzling, as non-thermal emission has been reported in a previous work for this source.
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.
We present very-high-energy gamma-ray observations of the BL Lac object 1ES 2344+514 taken by the Very Energetic Radiation Imaging Telescope Array System between 2007 and 2015. 1ES 2344+514 is detected with a statistical significance above the background of 20.8 sigma in 47.2 h (livetime) of observations, making this the most comprehensive very-high-energy study of 1ES 2344+514 to date. Using these observations, the temporal properties of 1ES 2344+514 are studied on short and long times-scales. We fit a constant-flux model to nightly and seasonally binned light curves and apply a fractional variability test to determine the stability of the source on different time-scales. We reject the constant-flux model for the 2007-2008 and 2014-2015 nightly binned light curves and for the long-term seasonally binned light curve at the > 3 sigma level. The spectra of the time-averaged emission before and after correction for attenuation by the extragalactic background light are obtained. The observed time-averaged spectrum above 200 GeV is satisfactorily fitted (x(2)/NDF = 7.89/6) by a power-law function with an index Gamma = 2.46 +/- 0.06(stat) +/- 0.20(sys) and extends to at least 8 TeV. The extragalactic-backgroundlight-deabsorbed spectrum is adequately fit (x(2)/NDF = 6.73/6) by a power-law function with an index Gamma = 2.15 +/- 0.06(stat) +/- 0.20(sys) while an F-test indicates that the power law with an exponential cut-off function provides a marginally better fit (x(2)/NDF = 2.56/5) at the 2.1 sigma level. The source location is found to be consistent with the published radio location and its spatial extent is consistent with a point source.
The very high energy (VHE; E > 100 GeV) blazar Markarian 501 was observed between April 17 and May 5 (MJD 54 938-54 956), 2009, as part of an extensive multiwavelength campaign from radio to VHE. Strong VHE yray activity was detected on May 1st with Whipple and VERITAS, when the flux (E > 400 GeV) increased to 10 times the preflare baseline flux (3.9 x 10(-11) ph cm(-2) s(-1)), reaching five times the flux of the Crab Nebula. This coincided with a decrease in the optical polarization and a rotation of the polarization angle by 15. This VHE flare showed a fast flux variation with an increase of a factor similar to 4 in 25 min, and a falling time of similar to 50 min. We present the observations of the quiescent state previous to the flare and of the high state after the flare, focusing on the flux and spectral variability from Whipple, VERITAS, Fermi-LAT, RXTE, and Swift combined with optical and radio data.
We report results from TeV gamma-ray observations of the microquasar Cygnus X-3. The observations were made with the Very Energetic Radiation Imaging Telescope Array System (VERITAS) over a time period from 2007 June 11 to 2011 November 28. VERITAS is most sensitive to gamma rays at energies between 85 GeV and 30 TeV. The effective exposure time amounts to a total of about 44 hr, with the observations covering six distinct radio/X-ray states of the object. No significant TeV gamma-ray emission was detected in any of the states, nor with all observations combined. The lack of a positive signal, especially in the states where GeV gamma rays were detected, places constraints on TeV gamma-ray production in Cygnus X-3. We discuss the implications of the results.
We report on very high energy (>100 GeV) gamma-ray observations of Swift J164449.3+573451, an unusual transient object first detected by the Swift Observatory and later detected by multiple radio, optical, and X-ray observatories. A total exposure of 28 hr was obtained on Swift J164449.3+573451 with the Very Energetic Radiation Imaging Telescope Array System ( VERITAS) during 2011 March 28-April 15. We do not detect the source and place a differential upper limit on the emission at 500 GeV during these observations of 1.4 x 10(-12) erg cm(-2) s(-1) (99% confidence level). We also present time-resolved upper limits and use a flux limit averaged over the X-ray flaring period to constrain various emission scenarios that can accommodate both the radio-through-X-ray emission detected from the source and the lack of detection by VERITAS.
We present the results of observations of the TeV binary LS I + 61 degrees 303 with the VERITAS telescope array between 2008 and 2010, at energies above 300 GeV. In the past, both ground-based gamma-ray telescopes VERITAS and MAGIC have reported detections of TeV emission near the apastron phases of the binary orbit. The observations presented here show no strong evidence for TeV emission during these orbital phases; however, during observations taken in late 2010, significant emission was detected from the source close to the phase of superior conjunction (much closer to periastron passage) at a 5.6 standard deviation (5.6 sigma) post-trials significance. In total, between 2008 October and 2010 December a total exposure of 64.5 hr was accumulated with VERITAS on LS I + 61 degrees 303, resulting in an excess at the 3.3 sigma significance level for constant emission over the entire integrated data set. The flux upper limits derived for emission during the previously reliably active TeV phases (i.e., close to apastron) are less than 5% of the Crab Nebula flux in the same energy range. This result stands in apparent contrast to previous observations by both MAGIC and VERITAS which detected the source during these phases at 10% of the Crab Nebula flux. During the two year span of observations, a large amount of X-ray data were also accrued on LS I + 61 degrees 303 by the Swift X-ray Telescope and the Rossi X-ray Timing Explorer Proportional Counter Array. We find no evidence for a correlation between emission in the X-ray and TeV regimes during 20 directly overlapping observations. We also comment on data obtained contemporaneously by the Fermi Large Area Telescope.
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 VERITAS observations of the BL Lac object PG 1553+113 spanning the years 2010, 2011, and 2012. The time-averaged spectrum, measured between 160 and 560 GeV, is well described by a power law with a spectral index of 4.33 +/- 0.09. The time-averaged integral flux above 200 GeV measured for this period was (1.69 +/- 0.06) x 10(-11) photons cm(-2) s(-1), corresponding to 6.9% of the Crab Nebula flux. We also present the combined gamma-ray spectrum from the Fermi Large Area Telescope and VERITAS covering an energy range from 100 MeV to 560 GeV. The data are well fit by a power law with an exponential cutoff at 101.9 +/- 3.2 GeV. The origin of the cutoff could be intrinsic to PG 1553+113 or be due to the gamma-ray opacity of our universe through pair production off the extragalactic background light (EBL). Given lower limits to the redshift of z > 0.395 based on optical/UV observations of PG 1553+113, the cutoff would be dominated by EBL absorption. Conversely, the small statistical uncertainties of the VERITAS energy spectrum have allowed us to provide a robust upper limit on the redshift of PG 1553+113 of z <= 0.62. A strongly elevated mean flux of (2.50 +/- 0.14) x10(-11) photons cm(-2) s(-1) (10.3% of the Crab Nebula flux) was observed during 2012, with the daily flux reaching as high as (4.44 +/- 0.71) x10(-11) photons cm(-2) s(-1) (18.3% of the Crab Nebula flux) on MJD 56048. The light curve measured during the 2012 observing season is marginally inconsistent with a steady flux, giving a chi(2) probability for a steady flux of 0.03%.
We report on VERITAS very high energy (VHE; E >= 100 GeV) observations of six blazars selected from the Fermi Large Area Telescope First Source Catalog (1FGL). The gamma-ray emission from 1FGL sources was extrapolated up to the VHE band, taking gamma-ray absorption by the extragalactic background light into account. This allowed the selection of six bright, hard-spectrum blazars that were good candidate TeV emitters. Spectroscopic redshift measurements were attempted with the Keck Telescope for the targets without Sloan Digital Sky Survey spectroscopic data. No VHE emission is detected during the observations of the six sources described here. Corresponding TeV upper limits are presented, along with contemporaneous Fermi observations and non-concurrent Swift UVOT and X-Ray Telescope data. The blazar broadband spectral energy distributions (SEDs) are assembled and modeled with a single-zone synchrotron self-Compton model. The SED built for each of the six blazars shows a synchrotron peak bordering between the intermediate-and high-spectrum-peak classifications, with four of the six resulting in particle-dominated emission regions.
We present the results of 16 Swift-triggered Gamma-ray burst (GRB) follow-up observations taken with the Very Energetic Radiation Imaging Telescope Array System (VERITAS) telescope array from 2007 January to 2009 June. The median energy threshold and response time of these observations were 260 GeV and 320 s, respectively. Observations had an average duration of 90 minutes. Each burst is analyzed independently in two modes: over the whole duration of the observations and again over a shorter timescale determined by the maximum VERITAS sensitivity to a burst with a t(-1.5) time profile. This temporal model is characteristic of GRB afterglows with high-energy, long-lived emission that have been detected by the Large Area Telescope on board the Fermi satellite. No significant very high energy (VHE) gamma-ray emission was detected and upper limits above the VERITAS threshold energy are calculated. The VERITAS upper limits are corrected for gamma-ray extinction by the extragalactic background light and interpreted in the context of the keV emission detected by Swift. For some bursts the VHE emission must have less power than the keV emission, placing constraints on inverse Compton models of VHE emission.
VERITAS has been monitoring the very-high-energy (VHE; > 100 GeV) gamma-ray activity of the radio galaxy M87 since 2007. During 2008, flaring activity on a timescale of a few days was observed with a peak flux of (0.70 +/- 0.16) x 10(-11) cm(-2) s(-1) at energies above 350 GeV. In 2010 April, VERITAS detected a flare from M 87 with peak flux of (2.71 +/- 0.68) x 10(-11) cm(-2) s(-1) for E > 350 GeV. The source was observed for six consecutive nights during the flare, resulting in a total of 21 hr of good-quality data. The most rapid flux variation occurred on the trailing edge of the flare with an exponential flux decay time of 0.90(-0.15)(+0.22) days. The shortest detected exponential rise time is three times as long, at 2.87(+1.65)(-0.99) days. The quality of the data sample is such that spectral analysis can be performed for three periods: rising flux, peak flux, and falling flux. The spectra obtained are consistent with power-law forms. The spectral index at the peak of the flare is equal to 2.19 +/- 0.07. There is some indication that the spectrum is softer in the falling phase of the flare than the peak phase, with a confidence level corresponding to 3.6 standard deviations. We discuss the implications of these results for the acceleration and cooling rates of VHE electrons in M 87 and the constraints they provide on the physical size of the emitting region.
The VERITAS array of Cherenkov telescopes has carried out a deep observational program on the nearby dwarf spheroidal galaxy Segue 1. We report on the results of nearly 48 hours of good quality selected data, taken between January 2010 and May 2011. No significant gamma-ray emission is detected at the nominal position of Segue 1, and upper limits on the integrated flux are derived. According to recent studies, Segue 1 is the most dark matter-dominated dwarf spheroidal galaxy currently known. We derive stringent bounds on various annihilating and decaying dark matter particle models. The upper limits on the velocity-weighted annihilation cross-section are <sigma upsilon >(95%) (CL) less than or similar to 10(-23) cm(3) s(-1), improving our limits from previous observations of dwarf spheroidal galaxies by at least a factor of 2 for dark matter particle masses m(chi) greater than or similar to 300 GeV. The lower limits on the decay lifetime are at the level of tau(95%) (CL) greater than or similar to 10(24) s. Finally, we address the interpretation of the cosmic ray lepton anomalies measured by ATIC and PAMELA in terms of dark matter annihilation, and show that the VERITAS observations of Segue 1 disfavor such a scenario.
We present results from multiwavelength observations of the BL Lacertae object 1ES 1741 + 196, including results in the very high energy gamma-ray regime using the Very Energetic Radiation Imaging Telescope Array System (VERITAS). The VERITAS time-averaged spectrum, measured above 180 GeV, is well modelled by a power law with a spectral index of 2.7 +/- 0.7(stat) +/- 0.2(syst). The integral flux above 180 GeV is (3.9 +/- 0.8(stat) +/- 1.0(syst)) x 10(-8) m(-2) s(-1), corresponding to 1.6 per cent of the Crab nebula flux on average. The multiwavelength spectral energy distribution of the source suggests that 1ES 1741+196 is an extreme-high-frequency-peaked BL Lacertae object. The observations analysed in this paper extend over a period of six years, during which time no strong flares were observed in any band. This analysis is therefore one of the few characterizations of a blazar in a non-flaring state.
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.
Between the beginning of its full-scale scientific operations in 2007 and 2012, the VERITAS Cherenkov telescope array observed more than 130 blazars; of these, 26 were detected as very-high-energy (VHE; E > 100 GeV) gamma-ray sources. In this work, we present the analysis results of a sample of 114 undetected objects. The observations constitute a total live-time of similar to 570 hr. The sample includes several unidentified Fermi-Large Area Telescope (LAT) sources (located at high Galactic latitude) as well as all the sources from the second Fermi-LAT catalog that are contained within the field of view of the VERITAS observations. We have also performed optical spectroscopy measurements in order to estimate the redshift of some of these blazars that do not have spectroscopic distance estimates. We present new optical spectra from the Kast instrument on the Shane telescope at the Lick observatory for 18 blazars included in this work, which allowed for the successful measurement or constraint on the redshift of four of them. For each of the blazars included in our sample, we provide the flux upper limit in the VERITAS energy band. We also study the properties of the significance distributions and we present the result of a stacked analysis of the data set, which shows a 4s excess.
We investigate the temporal and spectral correlations between flux and anisotropy fluctuations of TeV-band cosmic rays in light of recent data taken with IceCube. We find that for a conventional distribution of cosmic-ray sources, the dipole anisotropy is higher than observed, even if source discreteness is taken into account. Moreover, even for a shallow distribution of galactic cosmic-ray sources and a reacceleration model, fluctuations arising from source discreteness provide a probability only of the order of 10% that the cosmic-ray anisotropy limits of the recent IceCube analysis are met. This probability estimate is nearly independent of the exact choice of source rate, but generous for a large halo size. The location of the intensity maximum far from the Galactic Center is naturally reproduced.
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.
Context. Supernova remnants are known as sources of Galactic cosmic rays for their nonthermal emission of radio waves, X-rays, and gamma rays. However, the observed soft broken power-law spectra are hard to reproduce within standard acceleration theory based on the assumption of Bohm diffusion and steady-state calculations. Aims. We point out that a time-dependent treatment of the acceleration process together with a self-consistent treatment of the scattering turbulence amplification is necessary. Methods. We numerically solve the coupled system of transport equations for cosmic rays and isotropic Alfvenic turbulence. The equations are coupled through the growth rate of turbulence determined by the cosmic-ray gradient and the spatial diffusion coefficient of cosmic rays determined by the energy density of the turbulence. The system is solved on a comoving expanding grid extending upstream for dozens of shock radii, allowing for the self-consistent study of cosmic-ray diffusion in the vicinity of their acceleration site. The transport equation for cosmic rays is solved in a test-particle approach. Results. We demonstrate that the system is typically not in a steady state. In fact, even after several thousand years of evolution, no equilibrium situation is reached. The resulting time-dependent particle spectra strongly differ from those derived assuming a steady state and Bohm diffusion. Our results indicate that proper accounting for the evolution of the scattering turbulence and hence the particle diffusion coefficient is crucial for the formation of the observed soft spectra. In any case, the need to continuously develop magnetic turbulence upstream of the shock introduces nonlinearity in addition to that imposed by cosmic-ray feedback.
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.
Time-dependent escape of cosmic rays from supernova remnants, and their interaction with dense media
(2012)
Context. Supernova remnants (SNRs) are thought to be the main source of Galactic cosmic rays (CRs) up to the "knee" in CR spectrum. During the evolution of a SNR, the bulk of the CRs are confined inside the SNR shell. The highest-energy particles leave the system continuously, while the remaining adiabatically cooled particles are released when the SNR has expanded sufficiently and decelerated so that the magnetic field at the shock is no longer able to confine them. Particles escaping from the parent system may interact with nearby molecular clouds, producing.-rays in the process via pion decay. The soft gamma-ray spectra observed for a number of SNRs interacting with molecular clouds, however, challenge current theories of non-linear particle acceleration that predict harder spectra.
Aims. We study how the spectrum of escaped particles depends on the time-dependent acceleration history in both Type Ia and core-collapse SNRs, as well as on different assumptions about the diffusion coefficient in the vicinity of the SNR.
Methods. We solve the CR transport equation in a test-particle approach combined with numerical simulations of SNR evolution.
Results. We extend our method for calculating the CR acceleration in SNRs to trace the escaped particles in a large volume around SNRs. We calculate the evolution of the spectra of CRs that have escaped from a SNR into a molecular cloud or dense shell for two diffusion models. We find a strong confinement of CRs in a close region around the SNR, and a strong dilution effect for CRs that were able to propagate out as far as a few SNR radii.
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.
The statistical distribution of energies among particles responsible for long gamma-ray burst (GRB) emission is analyzed in light of recent results of the Fermi Observatory. The all-sky flux, F., recorded by the Gamma-Ray Burst Monitor (GBM) is shown, despite its larger energy range, to be not significantly larger than that reported by the Burst and Transient Explorer, suggesting a relatively small flux in the 3-30MeV energy range. The present-day energy input rate in gamma-rays recorded by the GBM from long GRBs is found, assuming star formation rates in the literature, to be W(0) = 0.5F gamma H/c = 5x10(42) erg Mpc(-3) yr(-1). The Large Area Telescope fluence, when observed, is about 5%-10% per decade of the total, in good agreement with the predictions of saturated, nonlinear shock acceleration. The high- energy component of long GRBs, as measured by Fermi, is found to contain only similar to 10-2.5 of the energy needed to produce ultrahigh-energy cosmic rays (UHECRs) above 4 EeV, assuming the latter to be extragalactic, when various numerical factors are carefully included, if the cosmic-ray source spectrum has a spectral index of -2. The observed. - ray fraction of the required UHECR energy is even smaller if the source spectrum is softer than E-2. The AMANDA II limits rule out such a GRB origin for UHECRs if much more than 10(-2) of the cosmic-ray energy goes into neutrinos that are within, and simultaneous with, the gamma-ray beam. It is suggested that "orphan" neutrinos out of the gamma-ray beam might be identifiable via orphan afterglow or other wide angle signatures of GRBs in lieu of coincidence with prompt gamma-rays, and it is recommended that feasible single neutrino trigger criteria be established to search for such coincidences.
Context. The Galactic gamma-ray diffuse emission is currently observed in the GeV-TeV energy range with unprecedented accuracy by the Fermi satellite. Understanding this component is crucial because it provides a background to many different signals, such as extragalactic sources or annihilating dark matter. It is timely to reinvestigate how it is calculated and to assess the various uncertainties that are likely to affect the accuracy of the predictions.
Aims. The Galactic gamma-ray diffuse emission is mostly produced above a few GeV by the interactions of cosmic ray primaries impinging on the interstellar material. The theoretical error on that component is derived by exploring various potential sources of uncertainty. Particular attention is paid to cosmic ray propagation. Nuclear cross sections, the proton and helium fluxes at the Earth's position, the Galactic radial profile of supernova remnants, and the hydrogen distribution can also severely affect the signal.
Methods. The propagation of cosmic ray species throughout the Galaxy is described in the framework of a semi-analytic two-zone diffusion/convection model. The gamma-ray flux is reliably and quickly determined. This allows conversion of the constraints set by the boron-to-carbon data into a theoretical uncertainty on the diffuse emission. New deconvolutions of the HI and CO sky maps are also used to get the hydrogen distribution within the Galaxy.
Results. The thickness of the cosmic ray diffusive halo is found to have a significant effect on the Galactic gamma-ray diffuse emission, while the interplay between diffusion and convection has little influence on the signal. The uncertainties related to nuclear cross sections and to the primary cosmic ray fluxes at the Earth are significant. The radial distribution of supernova remnants along the Galactic plane turns out to be a key ingredient. As expected, the predictions are extremely sensitive to the spatial distribution of hydrogen within the Milky Way.
Conclusions. Most of the sources of uncertainty are likely to be reduced in the near future. The stress should be put (i) on better determination of the thickness of the cosmic ray diffusive halo; and (ii) on refined observations of the radial profile of supernova remnants.
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.
e-ASTROGAM (‘enhanced ASTROGAM’) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV – the lower energy limit can be pushed to energies as low as 150 keV, albeit with rapidly degrading angular resolution, for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with 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 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 provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and the promise of eLISA.
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.
Aims. We analyze the emission in the 0.3-30 GeV energy range of gamma-ray bursts detected with the Fermi Gamma-ray Space Telescope. We concentrate on bursts that were previously only detected with the Gamma-Ray Burst Monitor in the keV energy range. These bursts will then be compared to the bursts that were individually detected with the Large Area Telescope at higher energies.
Methods. To estimate the emission of faint GRBs we used nonstandard analysis methods and sum over many GRBs to find an average signal that is significantly above background level. We used a subsample of 99 GRBs listed in the Burst Catalog from the first two years of observation.
Results. Although most are not individually detectable, the bursts not detected by the Large Area Telescope on average emit a significant flux in the energy range from 0.3 GeV to 30 GeV, but their cumulative energy fluence is only 8% of that of all GRBs. Likewise, the GeV-to-MeV flux ratio is less and the GeV-band spectra are softer. We confirm that the GeV-band emission lasts much longer than the emission found in the keV energy range. The average allsky energy flux from GRBs in the GeV band is 6.4 x 10(-4) erg cm(-2) yr(-1) or only similar to 4% of the energy flux of cosmic rays above the ankle at 10(18.6) eV.
The 2010 very high energy gamma-ray flare and 10 years ofmulti-wavelength oservations of M 87
(2012)
The giant radio galaxy M 87 with its proximity (16 Mpc), famous jet, and very massive black hole ((3-6) x 10(9) M-circle dot) provides a unique opportunity to investigate the origin of very high energy (VHE; E > 100 GeV) gamma-ray emission generated in relativistic outflows and the surroundings of supermassive black holes. M 87 has been established as a VHE gamma-ray emitter since 2006. The VHE gamma-ray emission displays strong variability on timescales as short as a day. In this paper, results from a joint VHE monitoring campaign on M 87 by the MAGIC and VERITAS instruments in 2010 are reported. During the campaign, a flare at VHE was detected triggering further observations at VHE (H.E.S.S.), X-rays (Chandra), and radio (43 GHz Very Long Baseline Array, VLBA). The excellent sampling of the VHE gamma-ray light curve enables one to derive a precise temporal characterization of the flare: the single, isolated flare is well described by a two-sided exponential function with significantly different flux rise and decay times of tau(rise)(d) = (1.69 +/- 0.30) days and tau(decay)(d) = (0.611 +/- 0.080) days, respectively. While the overall variability pattern of the 2010 flare appears somewhat different from that of previous VHE flares in 2005 and 2008, they share very similar timescales (similar to day), peak fluxes (Phi(>0.35 TeV) similar or equal to (1-3) x 10(-11) photons cm(-2) s(-1)), and VHE spectra. VLBA radio observations of 43 GHz of the inner jet regions indicate no enhanced flux in 2010 in contrast to observations in 2008, where an increase of the radio flux of the innermost core regions coincided with a VHE flare. On the other hand, Chandra X-ray observations taken similar to 3 days after the peak of the VHE gamma-ray emission reveal an enhanced flux from the core (flux increased by factor similar to 2; variability timescale <2 days). The long-term (2001-2010) multi-wavelength (MWL) light curve of M 87, spanning from radio to VHE and including data from Hubble Space Telescope, Liverpool Telescope, Very Large Array, and European VLBI Network, is used to further investigate the origin of the VHE gamma-ray emission. No unique, common MWL signature of the three VHE flares has been identified. In the outer kiloparsec jet region, in particular in HST-1, no enhanced MWL activity was detected in 2008 and 2010, disfavoring it as the origin of the VHE flares during these years. Shortly after two of the three flares (2008 and 2010), the X-ray core was observed to be at a higher flux level than its characteristic range (determined from more than 60 monitoring observations: 2002-2009). In 2005, the strong flux dominance of HST-1 could have suppressed the detection of such a feature. Published models for VHE gamma-ray emission from M 87 are reviewed in the light of the new data.
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.
We present the results of a multi-wavelength campaign targeting the blazar 1ES 1218+30.4 with observations with the 1.3 m McGraw-Hill optical telescope, the Rossi X-ray Timing Explorer (RXTE), the Fermi Gamma-Ray Space Telescope, and the Very Energetic Radiation Imaging Telescope Array System (VERITAS). The RXTE and VERITAS observations were spread over a 13 day period and revealed clear evidence for flux variability, and a strong X-ray and gamma-ray flare on 2009 February 26 (MJD 54888). The campaign delivered a well-sampled broadband energy spectrum with simultaneous RXTE and VERITAS very high energy (VHE, > 100 GeV) observations, as well as contemporaneous optical and Fermi observations. The 1ES 1218+30.4 broadband energy spectrum-the first with simultaneous X-ray and VHE gamma-ray energy spectra-is of particular interest as the source is located at a high cosmological redshift for a VHE source (z = 0.182), leading to strong absorption of VHE gamma rays by photons from the optical/infrared extragalactic background light (EBL) via gamma VHE +gamma EBL -> e(+) e(-)pair-creation processes. We model the data with a one-zone synchrotron self-Compton (SSC) emission model and with the extragalactic absorption predicted by several recent EBL models. We find that the observations are consistent with the SSC scenario and all the EBL models considered in this work. We discuss observational and theoretical avenues to improve on the EBL constraints.
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.
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.
A non-resonant cosmic ray (CR) current-driven instability may operate in the shock precursors of young supernova remnants and be responsible for magnetic-field amplification, plasma heating and turbulence. Earlier simulations demonstrated magnetic-field amplification, and in kinetic studies a reduction of the relative drift between CRs and thermal plasma was observed as backreaction. However, all published simulations used periodic boundary conditions, which do not account for mass conservation in decelerating flows and only allow the temporal development to be studied. Here we report results of fully kinetic particle-in-cell simulations with open boundaries that permit inflow of plasma on one side of the simulation box and outflow at the other end, hence allowing an investigation of both the temporal and the spatial development of the instability. Magnetic-field amplification proceeds as in studies with periodic boundaries and, observed here for the first time, the reduction of relative drifts causes the formation of a shock-like compression structure at which a fraction of the plasma ions are reflected. Turbulent electric field generated by the non-resonant instability inelastically scatters CRs, modifying and anisotropizing their energy distribution. Spatial CR scattering is compatible with Bohm diffusion. Electromagnetic turbulence leads to significant non-adiabatic heating of the background plasma maintaining bulk equipartition between ions and electrons. The highest temperatures are reached at sites of large-amplitude electrostatic fields. Ion spectra show supra-thermal tails resulting from stochastic scattering in the turbulent electric field. Together, these modifications in the plasma flow will affect the properties of the shock and particle acceleration there.
We present very high energy (VHE) imaging of MGRO J2019+37 obtained with the VERITAS observatory. The bright extended (similar to 2 degrees) unidentified Milagro source is located toward the rich star formation region Cygnus-X. MGRO J2019+37 is resolved into two VERITAS sources. The faint, point-like source VER J2016+371 overlaps CTB 87, a filled-center remnant (SNR) with no evidence of a supernova remnant shell at the present time. Its spectrum is well fit in the 0.65-10 TeV energy range by a power-law model with photon index 2.3 +/- 0.4. VER J2019+378 is a bright extended (similar to 1 degrees) source that likely accounts for the bulk of the Milagro emission and is notably coincident with PSR J2021+3651 and the star formation region Sh 2-104. Its spectrum in the range 1-30 TeV is well fit with a power-law model of photon index 1.75 +/- 0.3, among the hardest values measured in the VHE band, comparable to that observed near Vela-X. We explore the unusual spectrum and morphology in the radio and X-ray bands to constrain possible emission mechanisms for this source.
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.
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.
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.
We report on the detection of a very rapid TeV gamma-ray flare from BL Lacertae on 2011 June 28 with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). The flaring activity was observed during a 34.6 minute exposure, when the integral flux above 200 GeV reached (3.4 +/- 0.6) x 10(-6) photons m(-2) s(-1), roughly 125% of the Crab Nebula flux measured by VERITAS. The light curve indicates that the observations missed the rising phase of the flare but covered a significant portion of the decaying phase. The exponential decay time was determined to be 13 +/- 4 minutes, making it one of the most rapid gamma-ray flares seen from a TeV blazar. The gamma-ray spectrum of BL Lacertae during the flare was soft, with a photon index of 3.6 +/- 0.4, which is in agreement with the measurement made previously by MAGIC in a lower flaring state. Contemporaneous radio observations of the source with the Very Long Baseline Array revealed the emergence of a new, superluminal component from the core around the time of the TeV gamma-ray flare, accompanied by changes in the optical polarization angle. Changes in flux also appear to have occurred at optical, UV, and GeV gamma-ray wavelengths at the time of the flare, although they are difficult to quantify precisely due to sparse coverage. A strong flare was seen at radio wavelengths roughly four months later, which might be related to the gamma-ray flaring activities. We discuss the implications of these multiwavelength results.
Recent PIC simulations of relativistic electron-positron (electron-ion) jets injected into a stationary medium show that particle acceleration occurs in the shocked regions. Simulations show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields and for particle acceleration. These magnetic fields contribute to the electron’s transverse eflection behind the shock. The “jitter” radiation from deflected electrons in turbulent magnetic fields has properties different from synchrotron radiation calculated in a uniform magnetic field. This jitter radiation may be important for understanding the complex time evolution and/or spectral structure of gamma-ray bursts, relativistic jets in general, and supernova remnants. In order to calculate radiation from first principles and go beyond the standard synchrotron model, we have used PIC simulations. We present synthetic spectra to compare with the spectra obtained from Fermi observations.
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.
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.
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.
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.
We present very-high-resolution 1D MHD simulations of the late-stage supernova remnants (SNRs). In the post-adiabatic stage, the magnetic field has an important and significant dynamical effect on the shock dynamics, the flow structure, and hence the acceleration and emission of cosmic rays. We find that the tangential component of the magnetic field provides pressure support that to a fair degree prevents the collapse of the radiative shell and thus limits the total compression ratio of the partially or fully radiative forward shock. A consequence is that the spectra of cosmic rays would not be as hard as in hydrodynamic simulations. We also investigated the effect on the flow profiles of the magnetic-field inclination and a large-scale gradient in the gas density and/or the magnetic field. A positive density gradient shortens the evolutionary stages, whereas a shock obliquity lowers the shock compression. The compression of the tangential component of the magnetic field leads to its dominance in the downstream region of post-adiabatic shocks for a wide range of orientation of the upstream field, which may explain why one preferentially observes tangential radio polarization in old SNRs. As most cosmic rays are produced at late stages of SNR evolution, the post-adiabatic phase and the influence of the magnetic field during it are most important for modeling the cosmic-ray acceleration at old SNRs and the gamma-ray emission from late-stage SNRs interacting with clouds.
Physik für alle
(2014)
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.
Particle spectra from acceleration at forward and reverse shocks of young Type Ia Supernova Remnants
(2012)
We study cosmic-ray acceleration in young Type Ia Supernova Remnants (SNRs) by means of test-particle diffusive shock acceleration theory and 1-D hydrodynamical simulations of their evolution. In addition to acceleration at the forward shock, we explore the particle acceleration at the reverse shock in the presence of a possible substantial magnetic field, and consequently the impact of this acceleration on the particle spectra in the remnant. We investigate the time evolution of the spectra for various time-dependent profiles of the magnetic field in the shocked region of the remnant. We test a possible influence on particle spectra of the Alfvenic drift of scattering centers in the precursor regions of the shocks. In addition, we study the radiation spectra and morphology in a broad band from radio to gamma-rays. It is demonstrated that the reverse shock contribution to the cosmic-ray particle population of young Type la SNRs may be significant, modifying the spatial distribution of particles and noticeably affecting the volume-integrated particle spectra in young SNRs. In particular spectral structures may arise in test-particle calculations that are often discussed as signatures of non-linear cosmic-ray modification of shocks. Therefore, the spectrum and morphology of emission, and their time evolution, differ from pure forward-shock solutions.
Particle diffusion and localized acceleration in inhomogeneous AGN jets - II. Stochastic variation
(2016)
We study the stochastic variation of blazar emission under a 2D spatially resolved leptonic jet model we previously developed. Random events of particle acceleration and injection in small zones within the emission region are assumed to be responsible for flux variations. In addition to producing spectral energy distributions that describe the observed flux of Mrk 421, we further analyse the timing properties of the simulated light curves, such as the power spectral density (PSD) at different bands, flux-flux correlations, aswell as the cross-correlation function between X-rays and TeV gamma-rays. We find spectral breaks in the PSD at a time-scale comparable to the dominant characteristic time-scale in the system, which is usually the predefined decay time-scale of an acceleration event. Cooling imposes a delay, and so PSDs taken at lower energy bands in each emission component (synchrotron or inverse Compton) generally break at longer time-scales. The flux-flux correlation between X-rays and TeV gamma-rays can be either quadratic or linear, depending on whether or not there are large variation of the injection into the particle acceleration process. When the relationship is quadratic, the TeV flares lag the X-ray flares, and the optical and GeV flares are large enough to be comparable to the ones in X-ray. When the relationship is linear, the lags are insignificant, and the optical and GeV flares are small.
Particle diffusion and localized acceleration in inhomogeneous AGN jets - I. Steady-state spectra
(2015)
We study the acceleration, transport, and emission of particles in relativistic jets. Localized stochastic particle acceleration, spatial diffusion, and synchrotron as well as synchrotron self-Compton (SSC) emission are considered in a leptonic model. To account for inhomogeneity, we use a 2D axisymmetric cylindrical geometry for both relativistic electrons and magnetic field. In this first phase of our work, we focus on steady-state spectra that develop from a time-dependent model. We demonstrate that small isolated acceleration region in a much larger emission volume are sufficient to accelerate particles to high energy. Diffusive escape from these small regions provides a natural explanation for the spectral form of the jet emission. The location of the acceleration regions within the jet is found to affect the cooling break of the spectrum in this diffusive model. Diffusion-caused energy-dependent inhomogeneity in the jets predicts that the SSC spectrum is harder than the synchrotron spectrum. There can also be a spectral hardening towards the high-energy section of the synchrotron spectrum, if particle escape is relatively slow. These two spectral hardening effects indicate that the jet inhomogeneity might be a natural explanation for the unexpected hard. gamma-ray spectra observed in some blazars.
We develop a simple one-zone model of the steady-state Crab nebula spectrum encompassing both the radio/soft X-ray and the GeV/multi-TeV observations. By solving the transport equation for GeV-TeV electrons injected at the wind termination shock as a log-parabola momentum distribution and evolved via energy losses, we determine analytically the resulting differential energy spectrum of photons. We find an impressive agreement with the observed spectrum of synchrotron emission, and the synchrotron self-Compton component reproduces the previously unexplained broad 200-GeV peak that matches the Fermi/Large Area Telescope (LAT) data beyond 1 GeV with the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) data. We determine the parameters of the single log-parabola electron injection distribution, in contrast with multiple broken power-law electron spectra proposed in the literature. The resulting photon differential spectrum provides a natural interpretation of the deviation from power law customarily fitted with empirical multiple broken power laws. Our model can be applied to the radio-to-multi-TeV spectrum of a variety of astrophysical outflows, including pulsar wind nebulae and supernova remnants, as well as to interplanetary shocks.
In this review we confront the current theoretical understanding of particle acceleration at relativistic outflows with recent observational results on various source classes thought to involve such outflows, e.g. gamma-ray bursts, active galactic nuclei, and pulsar wind nebulae. We highlight the possible contributions of these sources to ultra-high-energy cosmic rays.
Nonrelativistic electrostatic unmagnetized shocks are frequently observed in laboratory plasmas and they are likely to exist in astrophysical plasmas. Their maximum speed, expressed in units of the ion acoustic speed far upstream of the shock, depends only on the electron-to-ion temperature ratio if binary collisions are absent. The formation and evolution of such shocks is examined here for a wide range of shock speeds with particle-in-cell simulations. The initial temperatures of the electrons and the 400 times heavier ions are equal. Shocks form on electron time scales at Mach numbers between 1.7 and 2.2. Shocks with Mach numbers up to 2.5 form after tens of inverse ion plasma frequencies. The density of the shock-reflected ion beam increases and the number of ions crossing the shock thus decreases with an increasing Mach number, causing a slower expansion of the downstream region in its rest frame. The interval occupied by this ion beam is on a positive potential relative to the far upstream. This potential pre-heats the electrons ahead of the shock even in the absence of beam instabilities and decouples the electron temperature in the foreshock ahead of the shock from the one in the far upstream plasma. The effective Mach number of the shock is reduced by this electron heating. This effect can potentially stabilize nonrelativistic electrostatic shocks moving as fast as supernova remnant shocks.
We study the propagation of ultra-high-energy cosmic rays (UHECRs) in the Galaxy, concentrating on the energy range below the ankle in the spectrum at 4 EeV. A Monte Carlo method, based on analytical solutions to the time-dependent diffusion problem, is used to account for intermittency by placing sources at random locations. Assuming a source population that scales with baryon mass density or star formation (e.g., long GRB), we derive constraints arising from intermittency and the observational limits on the composition and anisotropy. It is shown that the composition and anisotropy at 10(18) eV are difficult to reproduce and require that either (1) the particle mean free path is much smaller than a gyroradius, implying the escape time is very long, (2) the composition is heavier than suggested by recent Auger data, (3) the ultra-high-energy sub-ankle component is mostly extragalactic, or (4) we are living in a rare lull in the UHECR production, and the current UHECR intensity is far below the Galactic time average. We therefore recommend a strong observational focus on determining the UHECR composition around 10(18) eV.
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.
Observations of the young supernova remnant RX J1713.7-3946 with the fermi large area telescope
(2011)
We present observations of the young supernova remnant (SNR) RX J1713.7-3946 with the Fermi Large Area Telescope (LAT). We clearly detect a source positionally coincident with the SNR. The source is extended with a best-fit extension of 0 degrees.55 +/- 0 degrees.04 matching the size of the non-thermal X-ray and TeV gamma-ray emission from the remnant. The positional coincidence and the matching extended emission allow us to identify the LAT source with SNR RX J1713.7-3946. The spectrum of the source can be described by a very hard power law with a photon index of Gamma = 1.5 +/- 0.1 that coincides in normalization with the steeper H. E. S. S.-detected gamma-ray spectrum at higher energies. The broadband gamma-ray emission is consistent with a leptonic origin as the dominant mechanism for the gamma-ray emission.
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.
Since the detection of nonthermal radio emission from the bow shock of the massive runaway star BD +43 degrees 3654, simple models have predicted high-energy emission, at X-rays and gamma-rays, from these Galactic sources. Observational searches for this emission so far give no conclusive evidence but a few candidates at gamma-rays. In this work we aim at developing a more sophisticated model for the nonthermal emission from massive runaway star bow shocks. The main goal is to establish whether these systems are efficient nonthermal emitters, even if they are not strong enough yet to be detected. For modeling the collision between the stellar wind and the interstellar medium we use 2D hydrodynamic simulations. We then adopt the flow profile of the wind and the ambient medium obtained with the simulation as the plasma state for solving the transport of energetic particles injected in the system, as well as the nonthermal emission they produce. For this purpose we solve a 3D (two spatial vertical bar energy) advection-diffusion equation in the test-particle approximation. We find that a massive runaway star with a powerful wind converts 0.16%-0.4% of the power injected in electrons into nonthermal emission, mostly produced by inverse Compton scattering of dust-emitted photons by relativistic electrons, and second by synchrotron radiation. This represents a fraction of similar to 10(-5) to 10(-4) of the wind kinetic power. Given the better sensibility of current instruments at radio wavelengths, these systems are more prone to be detected at radio through the synchrotron emission they produce rather than at gamma energies.
We present results of 2D3V particle-in-cell simulations of nonrelativistic plasma collisions with absent or parallel large-scale magnetic field for parameters applicable to the conditions at young supernova remnants. We study the collision of plasma slabs of different density, leading to two different shocks and a contact discontinuity. Electron dynamics play an important role in the development of the system. While nonrelativistic shocks in both unmagnetized and magnetized plasmas can be mediated by Weibel-type instabilities, the efficiency of shock-formation processes is higher when a large-scale magnetic field is present. The electron distributions downstream of the forward and reverse shocks are generally isotropic, whereas that is not always the case for the ions. We do not see any significant evidence of pre-acceleration, neither in the electron population nor in the ion distribution.
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.
Multiwavelength observations of the AGN 1ES 0414+009 with veritas, FERMI-LAT, SWIFT-XRT, AND MDM
(2012)
We present observations of the BL Lac object 1ES 0414+009 in the >200 GeV gamma-ray band by the VERITAS array of Cherenkov telescopes. 1ES 0414+009 was observed by VERITAS between 2008 January and 2011 February, resulting in 56.2 hr of good quality pointed observations. These observations resulted in a detection of 822 events from the source corresponding to a statistical significance of 6.4 standard deviations (6.4 sigma) above the background. The source flux, showing no evidence for variability, is measured as (5.2 +/- 1.1(stat) +/- 2.6(sys)) x 10(-12) photons cm(-2) s(-1) above 200 GeV, equivalent to approximately 2% of the Crab Nebula flux above this energy. The differential photon spectrum from 230 GeV to 850 GeV is well fit by a power law with a photon index of Gamma = 3.4 +/- 0.5(stat) +/- 0.3(sys) and a flux normalization of (1.6 +/- 0.3(stat) +/- 0.8(sys)) x 10(-11) photons cm(-2) s(-1) at 300 GeV. We also present multiwavelength results taken in the optical (MDM), x-ray (Swift-XRT), and GeV (Fermi-LAT) bands and use these results to construct a broadband spectral energy distribution (SED). Modeling of this SED indicates that homogenous one-zone leptonic scenarios are not adequate to describe emission from the system, with a lepto-hadronic model providing a better fit to the data.
We report on the VERITAS discovery of very high energy (VHE) gamma-ray emission above 200 GeV from the high-frequency-peaked BL Lac (HBL) object RX J0648.7+1516 (GB J0648+1516), associated with 1FGL J0648.8+1516. The photon spectrum above 200 GeV is fitted by a power law dN/dE = F-0(E/E-0)(-Gamma) with a photon index Gamma of 4.4 +/- 0.8(stat) +/- 0.3(syst) and a flux normalization F-0 of (2.3 +/- 0.5(stat) +/- 1.2(sys)) x 10(-11) TeV-1 cm(-2) s(-1) with E-0 = 300 GeV. No VHE variability is detected during VERITAS observations of RX J0648.7+1516 between 2010 March 4 and April 15. Following the VHE discovery, the optical identification and spectroscopic redshift were obtained using the Shane 3 m Telescope at the Lick Observatory, showing the unidentified object to be a BL Lac type with a redshift of z = 0.179. Broadband multiwavelength observations contemporaneous with the VERITAS exposure period can be used to subclassify the blazar as an HBL object, including data from the MDM observatory, Swift-UVOT, and X-Ray Telescope, and continuous monitoring at photon energies above 1 GeV from the Fermi Large Area Telescope (LAT). We find that in the absence of undetected, high-energy rapid variability, the one-zone synchrotron self-Compton (SSC) model overproduces the high-energy gamma-ray emission measured by the Fermi-LAT over 2.3 years. The spectral energy distribution can be parameterized satisfactorily with an external-Compton or lepto-hadronic model, which have two and six additional free parameters, respectively, compared to the one-zone SSC model.
We report on the VERITAS observations of the high-frequency peaked BL Lac object 1ES 1959+650 in the period 2007-2011. This source is detected at TeV energies by VERITAS at 16.4 standard deviation (sigma) significance in 7.6 hr of observation in a low flux state. A multiwavelength spectral energy distribution (SED) is constructed from contemporaneous data from VERITAS, Fermi-LAT, RXTE PCA, and Swift UVOT. Swift XRT data is not included in the SED due to a lack of simultaneous observations with VERITAS. In contrast to the orphan gamma-ray flare exhibited by this source in 2002, the X-ray flux of the source is found to vary by an order of magnitude, while other energy regimes exhibit less variable emission. A quasi-equilibrium synchrotron self-Compton model with an additional external radiation field is used to describe three SEDs corresponding to the lowest, highest, and average X-ray states. The variation in the X-ray spectrum is modeled by changing the electron injection spectral index, with minor adjustments of the kinetic luminosity in electrons. This scenario produces small-scale flux variability of the order of less than or similar to 2 in the high energy (E > 1MeV) and very high energy (E > 100 GeV) gamma-ray regimes, which is corroborated by the Fermi-LAT, VERITAS, and Whipple 10 m telescope light curves.
We present the results of a multiwavelength observational campaign on the TeV binary system LS I +61 degrees 303 with the VERITAS telescope array (>200 GeV), Fermi-LAT (0.3-300 GeV), and Swift/XRT (2-10 keV). The data were taken from 2011 December through 2012 January and show a strong detection in all three wavebands. During this period VERITAS obtained 24.9 hr of quality selected livetime data in which LS I +61 degrees 303 was detected at a statistical significance of 11.9 sigma. These TeV observations show evidence for nightly variability in the TeV regime at a post-trial significance of 3.6 sigma. The combination of the simultaneously obtained TeV and X-ray fluxes do not demonstrate any evidence for a correlation between emission in the two bands. For the first time since the launch of the Fermi satellite in 2008, this TeV detection allows the construction of a detailed MeV-TeV spectral energy distribution from LS I +61 degrees 303. This spectrum shows a distinct cutoff in emission near 4 GeV, with emission seen by the VERITAS observations following a simple power-law above 200 GeV. This feature in the spectrum of LS I +61 degrees 303, obtained from overlapping observations with Fermi-LAT and VERITAS, may indicate that there are two distinct populations of accelerated particles producing the GeV and TeV emission.
The formation of unmagnetized electrostatic shock-like structures with a high Mach number is examined with one-and two-dimensional particle-in-cell (PIC) simulations. The structures are generated through the collision of two identical plasma clouds, which consist of equally hot electrons and ions with a mass ratio of 250. The Mach number of the collision speed with respect to the initial ion acoustic speed of the plasma is set to 4.6. This high Mach number delays the formation of such structures by tens of inverse ion plasma frequencies. A pair of stable shock-like structures is observed after this time in the 1D simulation, which gradually evolves into electrostatic shocks. The ion acoustic instability, which can develop in the 2D simulation but not in the 1D one, competes with the nonlinear process that gives rise to these structures. The oblique ion acoustic waves fragment their electric field. The transition layer, across which the bulk of the ions change their speed, widens and their speed change is reduced. Double layer-shock hybrid structures develop.
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.
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.
In this study, we investigate the interaction of jets with their environment at a microscopic level, which is a key open question in the study of relativistic jets. Using small simulation systems during past research, we initially studied the evolution of both electron-proton and electron-positron relativistic jets containing helical magnetic fields, by focusing on their interactions with an ambient plasma. Here, using larger jet radii, we have performed simulations of global jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities, such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the mushroom instability (MI). We found that the evolution of global jets strongly depends on the size of the jet radius. For example, phase bunching of jet electrons, in particular in the electron-proton jet, is mixed with a larger jet radius as a result of the more complicated structures of magnetic fields with excited kinetic instabilities. In our simulation, these kinetic instabilities led to new types of instabilities in global jets. In the electron-proton jet simulation, a modified recollimation occurred, and jet electrons were strongly perturbed. In the electron-positron jet simulation, mixed kinetic instabilities occurred early, followed by a turbulence-like structure. Simulations using much larger (and longer) systems are required in order to further thoroughly investigate the evolution of global jets containing helical magnetic fields.
In the study of relativistic jets one of the key open questions is their interaction with the environment on the microscopic level. Here, we study the initial evolution of both electron-proton (e(-)-p(+)) and electron-positron (e(+/-)) relativistic jets containing helical magnetic fields, focusing on their interaction with an ambient plasma. We have performed simulations of "global" jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). In our initial simulation study these kinetic instabilities are suppressed and new types of instabilities can grow. In the e(-)-p(+) jet simulation a recollimation-like instability occurs and jet electrons are strongly perturbed. In the e(+/-) jet simulation a recollimation-like instability occurs at early times followed by a kinetic instability and the general structure is similar to a simulation without helical magnetic field. Simulations using much larger systems are required in order to thoroughly follow the evolution of global jets containing helical magnetic fields.
In the study of relativistic jets one of the key open questions is their interaction with the environment on the microscopic level. Here, we study the initial evolution of both electron-proton (e(-)-p(+)) and electron-positron (e(+/-)) relativistic jets containing helical magnetic fields, focusing on their interaction with an ambient plasma. We have performed simulations of "global" jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). In our initial simulation study these kinetic instabilities are suppressed and new types of instabilities can grow. In the e(-)-p(+) jet simulation a recollimation-like instability occurs and jet electrons are strongly perturbed. In the e(+/-) jet simulation a recollimation-like instability occurs at early times followed by a kinetic instability and the general structure is similar to a simulation without helical magnetic field. Simulations using much larger systems are required in order to thoroughly follow the evolution of global jets containing helical magnetic fields.
We present a new measurement of the energy spectrum of iron nuclei in cosmic rays from 20 TeV to 500 TeV; The measurement makes use of a template-based analysis method, which, for the first time, is applied to the energy reconstruction of iron-induced air showers recorded by the VERITAS array of imaging atmospheric Cherenkov telescopes. The event selection makes use of the direct Cherenkov light which is emitted by charged particles before the first interaction, as well as other parameters related to the shape of the recorded air shower images. The measured spectrum is well described by a power law dF/dE = f(0) center dot (E/E-0)(-gamma) over the full energy range, with gamma = 2.82 +/- 0.30(stat)(-0.27)(+0.24)(syst) and f(0) = (4.82 +/- 0.98(stat)(-2.70)(+2.12)(syst)) x 10(-7) m(-2) s(-1) TeV-1 sr(-1) at E-0 = 50 TeV, with no indication of a cutoff or spectral break. The measured differential flux is compatible with previous results, with improved statistical uncertainty at the highest energies.
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.
We perform two-dimensional relativistic magnetohydrodynamic simulations of a mildly relativistic shock propagating through an inhomogeneous medium. We show that the postshock region becomes turbulent owing to preshock density inhomogeneity, and the magnetic field is strongly amplified due to the stretching and folding of field lines in the turbulent velocity field. The amplified magnetic field evolves into a filamentary structure in two-dimensional simulations. The magnetic energy spectrum is flatter than the Kolmogorov spectrum and indicates that a so-called small-scale dynamo is occurring in the postshock region. We also find that the amount of magnetic-field amplification depends on the direction of the mean preshock magnetic field, and the timescale of magnetic-field growth depends on the shock strength.
Magnetic field generation in a jet-sheath plasma via the kinetic Kelvin-Helmholtz instability
(2013)
We have investigated the generation of magnetic fields associated with velocity shear between an unmagnetized relativistic jet and an unmagnetized sheath plasma. We have examined the strong magnetic fields generated by kinetic shear (Kelvin-Helmholtz) instabilities. Compared to the previous studies using counter-streaming performed by Alves et al. (2012), the structure of the kinetic Kelvin-Helmholtz instability (KKHI) of our jet-sheath configuration is slightly different, even for the global evolution of the strong transverse magnetic field. In our simulations the major components of growing modes are the electric field E-z, perpendicular to the flow boundary, and the magnetic field B-y, transverse to the flow direction. After the B-y component is excited, an induced electric field E-x, parallel to the flow direction, becomes significant. However, other field components remain small. We find that the structure and growth rate of KKHI with mass ratios m(i)/m(e) = 1836 and m(i)/m(e) = 20 are similar. In our simulations in the nonlinear stage is not as clear as in counter-streaming cases. The growth rate for a mildly-relativistic jet case (gamma(j) = 1.5) is larger than for a relativistic jet case (gamma(j) = 15).
We have investigated via 2D relativistic magnetohydrodynamic simulations the long-term evolution of turbulence created by a relativistic shock propagating through an inhomogeneous medium. In the post-shock region, magnetic field is strongly amplified by turbulent motions triggered by pre-shock density inhomogeneities. Using a long-simulation box we have followed the magnetic field amplification until it is fully developed and saturated. The turbulent velocity is subrelativistic even for a strong shock. Magnetic field amplification is controlled by the turbulent motion and saturation occurs when the magnetic energy is comparable to the turbulent kinetic energy. Magnetic field amplification and saturation depend on the initial strength and direction of the magnetic field in the pre-shock medium, and on the shock strength. If the initial magnetic field is perpendicular to the shock normal, the magnetic field is first compressed at the shock and then can be amplified by turbulent motion in the post-shock region. Saturation occurs when the magnetic energy becomes comparable to the turbulent kinetic energy in the post-shock region. If the initial magnetic field in the pre-shock medium is strong, the post-shock region becomes turbulent but significant field amplification does not occur. If the magnetic energy after shock compression is larger than the turbulent kinetic energy in the post-shock region, significant field amplification does not occur. We discuss possible applications of our results to gamma-ray bursts and active galactic nuclei.
We perform time-dependent, spatially resolved simulations of blazar emission to evaluate several flaring scenarios related to magnetic-field amplification and enhanced particle acceleration. The code explicitly accounts for light-travel-time effects and is applied to flares observed in the flat spectrum radio quasar (FSRQ) PKS 0208-512, which show optical/gamma-ray correlation at some times, but orphan optical flares at other times. Changes in both the magnetic field and the particle acceleration efficiency are explored as causes of flares. Generally, external Compton (EC) emission appears to describe the available data better than a synchrotron self-Compton (SSC) scenario, and in particular orphan optical flares are difficult to produce in the SSC framework. X-ray soft-excesses, gamma-ray spectral hardening, and the detections at very high energies of certain FSRQs during flares find natural explanations in the EC scenario with particle acceleration change. Likewise, optical flares with/without gamma-ray counterparts can be explained by different allocations of energy between the magnetization and particle acceleration, which may be related to the orientation of the magnetic field relative to the jet flow. We also calculate the degree of linear polarization and polarization angle as a function of time for a jet with helical magnetic field. Tightening of the magnetic helix immediately downstream of the jet perturbations, where flares occur, can be sufficient to explain the increases in the degree of polarization and a rotation by a parts per thousand yen180A degrees of the observed polarization angle, if light-travel-time effects are properly considered.
Magnetic agnetic field generation in core-sheath jets via the kinetic Kelvin-Helmholtz instability
(2014)
We have investigated magnetic field generation in velocity shears via the kinetic Kelvin-Helmholtz instability (kKHI) using a relativistic plasma jet core and stationary plasma sheath. Our three-dimensional particle-in-cell simulations consider plasma jet cores with Lorentz factors of 1.5, 5, and 15 for both electron-proton and electron-positron plasmas. For electron-proton plasmas, we find generation of strong large-scale DC currents and magnetic fields that extend over the entire shear surface and reach thicknesses of a few tens of electron skin depths. For electron-positron plasmas, we find generation of alternating currents and magnetic fields. Jet and sheath plasmas are accelerated across the shear surface in the strong magnetic fields generated by the kKHI. The mixing of jet and sheath plasmas generates a transverse structure similar to that produced by the Weibel instability.
We report on VERITAS observations of the BL Lac object B2 1215+30 between 2008 and 2012. During this period, the source was detected at very high energies (VHEs; E > 100 GeV) by VERITAS with a significance of 8.9s and showed clear variability on timescales larger than months. In 2011, the source was found to be in a relatively bright state and a power-law fit to the differential photon spectrum yields a spectral index of 3.6 +/- 0.4(stat) +/- 0.3(syst) with an integral flux above 200 GeV of (8.0 +/- 0.9(stat) +/- 3.2(syst)) x 10(-12) cm(-2) s(-1). No short term variability could be detected during the bright state in 2011. Multi-wavelength data were obtained contemporaneously with the VERITAS observations in 2011 and cover optical (Super-LOTIS, MDM, Swift/UVOT), X-ray (Swift/XRT), and gamma-ray (Fermi-LAT) frequencies. These were used to construct the spectral energy distribution (SED) of B2 1215+30. A one-zone leptonic model is used to model the blazar emission and the results are compared to those of MAGIC from early 2011 and other VERITAS-detected blazars. The SED can be reproduced well with model parameters typical for VHE-detected BL Lac 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.
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.
The hypothesis is considered that the present, local Galactic cosmic-ray spectrum is, due to source intermittency, softer than average over time and over the Galaxy. Measurements of muogenic nuclides underground could provide an independent measurement of the time-averaged spectrum. Source intermittency could also account for the surprising low anisotropy reported by the IceCube Collaboration. Predictions for Galactic emission of ultrahigh-energy (UHE) quanta, such as UHE gamma rays and neutrinos, might be higher or lower than previously estimated.
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.
The majority of the matter in the universe is still unidentified and under investigation by both direct and indirect means. Many experiments searching for the recoil of dark-matter particles off target nuclei in underground laboratories have established increasingly strong constraints on the mass and scattering cross sections of weakly interacting particles, and some have even seen hints at a possible signal. Other experiments search for a possible mixing of photons with light scalar or pseudo-scalar particles that could also constitute dark matter. Furthermore, annihilation or decay of dark matter can contribute to charged cosmic rays, photons at all energies, and neutrinos. Many existing and future ground-based and satellite experiments are sensitive to such signals. Finally, data from the Large Hadron Collider at CERN are scrutinized for missing energy as a signature of new weakly interacting particles that may be related to dark matter. In this review article we summarize the status of the field with an emphasis on the complementarity between direct detection in dedicated laboratory experiments, indirect detection in the cosmic radiation, and searches at particle accelerators. (C) 2015 Elsevier B.V. All rights reserved.
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.
We report on an approach to produce quasi continuous-wave (cw) electron beams with an average beam current of milliamperes and a mean beam energy of a few MeV in a pulsed RF injector. Potential applications are in the planned laboratory astrophysics programs at DESY. The beam generation is based on field emission from a specially designed metallic field emitter. A quasi cw beam profile is formed over subsequent RF cycles at the resonance frequency of the gun cavity. This is realized by debunching in a cut disk structure accelerating cavity (booster) downstream of the gun. The peak and average beam currents can be tuned in beam dynamics simulations by adjusting operation conditions of the booster cavity. Optimization of the transverse beam size at specific positions (e.g., entrance of the plasma experiment) is performed by applying magnetic focusing fields provided by solenoids along the beam line. In this paper, the design of a microtip field emitter is introduced and characterized in electromagnetic field simulations in the gun cavity. A series of particle tracking simulations are conducted for multi-parametric optimization of the parameters of the produced quasi cw electron beams. The obtained results will be presented and discussed. In addition, measurements of the parasitic field emission (PFE) current (dark current) in the PITZ gun will be exemplarily shown to distinguish its order of magnitude from the produced beam current by the designed field emitter.
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.
Giant X-ray outbursts, with luminosities of about 10(37) erg s(-1), are observed roughly every five years from the nearby Be/pulsar binary 1A 0535+262. In this article, we present observations of the source with VERITAS at very high energies (VHEs; E > 100 GeV) triggered by the X-ray outburst in 2009 December. The observations started shortly after the onset of the outburst and provided comprehensive coverage of the episode, as well as the 111 day binary orbit. No VHE emission is evident at any time. We also examined data from the contemporaneous observations of 1A 0535+262 with the Fermi/Large Area Telescope at high-energy photons (E > 0.1 GeV) and failed to detect the source at GeV energies. The X-ray continua measured with the Swift/X-Ray Telescope and the RXTE/PCA can be well described by the combination of blackbody and Comptonized emission from thermal electrons. Therefore, the gamma-ray and X-ray observations suggest the absence of a significant population of non-thermal particles in the system. This distinguishes 1A 0535+262 from those Be X-ray binaries (such as PSR B1259-63 and LS I +61 degrees 303) that have been detected at GeV-TeV energies. We discuss the implications of the results on theoretical models.
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.
We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 2013 April 1 and August 10, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope, Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Metsahovi, and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a Light Detection and Ranging (LIDAR) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution (SED) between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) show evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton (SSC) model to five simultaneous broadband SEDs. We find that the SSC model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission.
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.
The TeV binary system LS I +61 degrees 303 is known for its regular, non-thermal emission pattern that traces the orbital period of the compact object in its 26.5 day orbit around its B0 Ve star companion. The system typically presents elevated TeV emission around apastron passage with flux levels between 5% and 15% of the steady flux from the Crab Nebula (> 300 GeV). In this article, VERITAS observations of LS I + 61 degrees. 303 taken in late 2014 are presented, during which bright TeV flares around apastron at flux levels peaking above 30% of the Crab Nebula flux were detected. This is the brightest such activity from this source ever seen in the TeV regime. The strong outbursts have rise and fall times of less than a day. The short timescale of the flares, in conjunction with the observation of 10 TeV photons from LS I + 61 degrees 303 during the flares, provides constraints on the properties of the accelerator in the source.
In the study of relativistic jets one of the key open questions is their interaction with the environment. Here. we study the initial evolution of both electron-proton (e(-) - p(+)) and electron-positron (e(+/-)) relativistic jets, focusing on their lateral interaction with ambient plasma. We follow the evolution of toroidal magnetic fields generated by both the kinetic Kelvin-Helmholtz and Mushroom instabilities. For an e(-) - p(+) jet, the induced magnetic field collimates the jet and electrons are perpendicularly accelerated. As the instabilities saturate and subsequently weaken, the magnetic polarity switches from clockwise to counterclockwise in the middle of the jet. For an e(+/-) jet, we find strong mixing of electrons and positrons with the ambient plasma, resulting in the creation of a bow shock. The merging of current filaments generates density inhomogeneities that. initiate a forward shock. Strong jet-ambient plasma mixing prevents a full development of the jet (on the scale studied), revealing evidence for both jet collimation and particle acceleration in the forming bow shock. Differences in the magnetic field structure generated by e(-) - p(+) and e(+/-) jets may contribute to the polarization properties of the observed emission in AGN jets and gamma-ray bursts.
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.
We perform particle-in-cell simulations of perpendicular nonrelativistic collisionless shocks to study electron heating and pre-acceleration for parameters that permit the extrapolation to the conditions at young supernova remnants. Our high-resolution large-scale numerical experiments sample a representative portion of the shock surface and demonstrate that the efficiency of electron injection is strongly modulated with the phase of the shock reformation. For plasmas with low and moderate temperature (plasma beta beta p =5.10(-4) and 0.5 beta p =), we explore the nonlinear shock structure and electron pre-acceleration for various orientations of the large-scale magnetic field with respect to the simulation plane, while keeping it at 90 degrees to the shock normal. Ion reflection off of the shock leads to the formation of magnetic filaments in the shock ramp, resulting from Weibel-type instabilities, and electrostatic Buneman modes in the shock foot. In all of the cases under study, the latter provides first-stage electron energization through the shock-surfing acceleration mechanism. The subsequent energization strongly depends on the field orientation and proceeds through adiabatic or second-order Fermi acceleration processes for configurations with the out-of-plane and in-plane field components, respectively. For strictly out-of-plane field, the fraction of suprathermal electrons is much higher than for other configurations, because only in this case are the Buneman modes fully captured by the 2D simulation grid. Shocks in plasma with moderate bp provide more efficient pre-acceleration. The relevance of our results to the physics of fully 3D systems is discussed.
We report on the discovery of high-energy (HE; E > 0.1 GeV) and very high energy (VHE; E > 100 GeV) gamma-ray emission from the high-frequency-peaked BL Lac object RBS 0413. VERITAS, a ground-based gamma-ray observatory, detected VHE. rays from RBS 0413 with a statistical significance of 5.5 standard deviations (sigma) and a gamma-ray flux of (1.5 +/- 0.6(stat) +/- 0.7(syst)) x 10(-8) photons m(-2) s(-1) (similar to 1% of the Crab Nebula flux) above 250 GeV. The observed spectrum can be described by a power law with a photon index of 3.18 +/- 0.68(stat) +/- 0.30(syst). Contemporaneous observations with the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope detected HE gamma rays from RBS 0413 with a (stat)istical significance of more than 9 sigma, a power-law photon index of 1.57 +/- 0.12(stat-0.12sys')(+0.11) and a gamma-ray flux between 300 MeV and 300 GeV of (1.64 +/- 0.43(stat-0.22sys)(+ 0.31)) x 10(-5) photons m(-2) s(-1). We present the results from Fermi-LAT and VERITAS, including a spectral energy distribution modeling of the gamma-ray, quasi-simultaneous X-ray (Swift-XRT), ultraviolet (Swift-UVOT), and R-band optical (MDM) data. We find that, if conditions close to equipartition are required, both the combined synchrotron self-Compton/external-Compton and the lepto-hadronic models are preferred over a pure synchrotron self-Compton model.