Refine
Has Fulltext
- no (123) (remove)
Year of publication
Document Type
- Article (114)
- Other (5)
- Review (2)
- Monograph/Edited Volume (1)
- Conference Proceeding (1)
Is part of the Bibliography
- yes (123)
Keywords
- acceleration of particles (20)
- ISM: supernova remnants (18)
- cosmic rays (17)
- gamma rays: galaxies (17)
- gamma rays: general (15)
- radiation mechanisms: non-thermal (12)
- shock waves (12)
- galaxies: active (10)
- BL Lacertae objects: general (8)
- galaxies: jets (8)
Institute
- Institut für Physik und Astronomie (123) (remove)
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 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 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.
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.
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.
Discovery of Very-high-energy Emission from RGB J2243+203 and Derivation of Its Redshift Upper Limit
(2017)
Very-high-energy (VHE; > 100 GeV) gamma-ray emission from the blazar RGB J2243+203 was discovered with the VERITAS Cherenkov telescope array, during the period between 2014 December 21 and 24. The VERITAS energy spectrum from this source can be fitted by a power law with a photon index of 4.6 +/- 0.5, and a flux normalization at 0.15 TeV of (6.3 +/- 1.1) x 10(-10) cm(-2) s(-1) TeV-1. The integrated Fermi-LAT flux from 1 to 100 GeV during the VERITAS detection is (4.1 +/- 0.8) x 10(-8) cm(-2) s(-1), which is an order of magnitude larger than the four-year-averaged flux in the same energy range reported in the 3FGL catalog, (4.0 +/- 0.1 x 10(-9) cm(-2) s(-1)). The detection with VERITAS triggered observations in the X-ray band with the Swift-XRT. However, due to scheduling constraints Swift-XRT observations were performed 67 hr after the VERITAS detection, rather than simultaneously with the VERITAS observations. The observed X-ray energy spectrum between 2 and 10 keV can be fitted with a power law with a spectral index of 2.7 +/- 0.2, and the integrated photon flux in the same energy band is (3.6 +/- 0.6) x 10(-13) cm(-2) s(-1). EBL-model-dependent upper limits of the blazar redshift have been derived. Depending on the EBL model used, the upper limit varies in the range from z < 0.9 to z < 1.1.
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.
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.
The F-type star KIC. 8462852 has recently been identified as an exceptional target for search for extraterrestrial intelligence (SETI) observations. We describe an analysis methodology for optical SETI, which we have used to analyze nine hours of serendipitous archival observations of KIC. 8462852 made with the VERITAS gamma-ray observatory between 2009 and 2015. No evidence of pulsed optical beacons, above a pulse intensity at the Earth of approximately 1 photon m(-2), is found. We also discuss the potential use of imaging atmospheric Cherenkov telescope arrays in searching for extremely short duration optical transients in general.
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.
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.
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.
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.
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 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.
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.
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.
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 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%.
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.
Context. RX J1713.7-3946 is the brightest shell-type supernova remnant (SNR) of the TeV gamma-ray sky. Earlier Fermi-LAT results on low energy gamma-ray emission suggested that, despite large uncertainties in the background determination, the spectrum is inconsistent with a hadronic origin.
Aims. We update the GeV-band spectra using improved estimates for the diffuse Galactic gamma-ray emission and more than double the volume of data. We further investigate the viability of hadronic emission models for RX J1713.7-3946.
Methods. We produced a high-resolution map of the diffuse Galactic gamma-ray background corrected for the HI self-absorption and used it in the analysis of more than five years worth of Fermi-LAT data. We used hydrodynamic scaling relations and a kinetic transport equation to calculate the acceleration and propagation of cosmic rays in SNR. We then determined spectra of hadronic gamma-ray emission from RX J1713.7-3946, separately for the SNR interior and the cosmic-ray precursor region of the forward shock, and computed flux variations that would allow us to test the model with observations.
Results. We find that RX J1713.7-3946 is now detected by Fermi-LAT with very high statistical significance, and the source morphology is best described by that seen in the TeV band. The measured spectrum of RX J1713.7-3946 is hard with index gamma = 1.53 +/- 0.07, and the integral flux above 500 MeV is F = (5 : 5 +/- 1 : 1) x 10(-9) photons cm(-2) s(-1). We demonstrate that scenarios based on hadronic emission from the cosmic-ray precursor region are acceptable for RX J1713.7-3946, and we predict a secular flux increase at a few hundred GeV at the level of around 15% over ten years, which may be detectable with the upcoming Cherenkov Telescope Array (CTA) observatory.
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.
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.
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 present deep VERITAS observations of the blazar PKS 1424+240, along with contemporaneous Fermi Large Area Telescope, Swift X-ray Telescope, and Swift UV Optical Telescope data between 2009 February 19 and 2013 June 8. This blazar resides at a redshift of z >= 0.6035, displaying a significantly attenuated gamma-ray flux above 100 GeV due to photon absorption via pair-production with the extragalactic background light. We present more than 100 hr of VERITAS observations over three years, a multiwavelength light curve, and the contemporaneous spectral energy distributions. The source shows a higher flux of (2.1 +/- 0.3) x 10(-7) photons m(-2) s(-1) above 120 GeV in 2009 and 2011 as compared to the flux measured in 2013, corresponding to (1.02 +/- 0.08) x 10-7 photons m(-2) s(-1) above 120 GeV. The measured differential very high energy (VHE; E >= 100 GeV) spectral indices are Gamma = 3.8 +/- 0.3, 4.3 +/- 0.6 and 4.5 +/- 0.2 in 2009, 2011, and 2013, respectively. No significant spectral change across the observation epochs is detected. We find no evidence for variability at gamma-ray opacities of greater than tau = 2, where it is postulated that any variability would be small and occur on timescales longer than a year if hadronic cosmic-ray interactions with extragalactic photon fields provide a secondary VHE photon flux. The data cannot rule out such variability due to low statistics.
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.
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 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.
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.
Prompt emission from the very fluent and nearby (z = 0.34) gamma-ray burst GRB130427A was detected by several orbiting telescopes and by ground-based, wide-field-of-view optical transient monitors. Apart from the intensity and proximity of this GRB, it is exceptional due to the extremely long-lived high-energy (100 MeV to 100 GeV) gamma-ray emission, which was detected by the Large Area Telescope on the Fermi Gamma-Ray Space Telescope for similar to 70 ks after the initial burst. The persistent, hard-spectrum, high-energy emission suggests that the highest-energy gamma rays may have been produced via synchrotron self-Compton processes though there is also evidence that the high-energy emission may instead be an extension of the synchrotron spectrum. VERITAS, a ground-based imaging atmospheric Cherenkov telescope array, began follow-up observations of GRB130427A similar to 71 ks (similar to 20 hr) after the onset of the burst. The GRB was not detected with VERITAS; however, the high elevation of the observations, coupled with the low redshift of the GRB, make VERITAS a very sensitive probe of the emission from GRB130427A for E > 100 GeV. The non-detection and consequent upper limit derived place constraints on the synchrotron self-Compton model of high-energy gamma-ray emission from this burst.
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.
We report the discovery of TeV gamma-ray emission from the Type Ia supernova remnant (SNR) G120.1+1.4, known as Tycho's SNR. Observations performed in the period 2008-2010 with the VERITAS ground-based gamma-ray observatory reveal weak emission coming from the direction of the remnant, compatible with a point source located at 00(h)25(m)27(s).0, +64 degrees 10'50 '' (J2000). The TeV photon spectrum measured by VERITAS can be described with a power law dN/dE = C(E/3.42 TeV)(-Gamma) with Gamma = 1.95 +/- 0.51(stat) +/- 0.30(sys) and C = (1.55 +/- 0.43(stat) +/- 0.47(sys)) x 10(-14) cm(-2) s(-1) TeV-1. The integral flux above 1 TeV corresponds to similar to 0.9% of the steady Crab Nebula emission above the same energy, making it one of the weakest sources yet detected in TeV gamma rays. We present both leptonic and hadronic models that can describe the data. The lowest magnetic field allowed in these models is similar to 80 mu G, which may be interpreted as evidence for magnetic field amplification.
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.
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.
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.
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.
The maximum cosmic-ray energy achievable by acceleration by a relativistic blast wave is derived. It is shown that forward shocks from long gamma-ray bursts (GRBs) in the interstellar medium accelerate protons to large enough energies, and have a sufficient energy budget, to produce the Galactic cosmic-ray component just below the ankle at 4 x 10(18) eV, as per an earlier suggestion. It is further argued that, were extragalactic long GRBs responsible for the component above the ankle as well, the occasional Galactic GRB within the solar circle would contribute more than the observational limits on the outward flux from the solar circle, unless an avoidance scenario, such as intermittency and/or beaming, allows the present-day local flux to be less than 10(-3) of the average. Difficulties with these avoidance scenarios are noted.
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 report the detection of pulsed gamma rays from the Crab pulsar at energies above 100 giga-electron volts (GeV) with the Very Energetic Radiation Imaging Telescope Array System (VERITAS) array of atmospheric Cherenkov telescopes. The detection cannot be explained on the basis of current pulsar models. The photon spectrum of pulsed emission between 100 mega-electron volts and 400 GeV is described by a broken power law that is statistically preferred over a power law with an exponential cutoff. It is unlikely that the observation can be explained by invoking curvature radiation as the origin of the observed gamma rays above 100 GeV. Our findings require that these gamma rays be produced more than 10 stellar radii from the neutron star.
The origin of Galactic cosmic rays is a century-long puzzle. Indirect evidence points to their acceleration by supernova shockwaves, but we know little of their escape from the shock and their evolution through the turbulent medium surrounding massive stars. Gamma rays can probe their spreading through the ambient gas and radiation fields. The Fermi Large Area Telescope (LAT) has observed the star-forming region of Cygnus X. The 1- to 100-gigaelectronvolt images reveal a 50-parsec-wide cocoon of freshly accelerated cosmic rays that flood the cavities carved by the stellar winds and ionization fronts from young stellar clusters. It provides an example to study the youth of cosmic rays in a superbubble environment before they merge into the older Galactic population.
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 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.
Ground-based gamma-ray astronomy has had a major breakthrough with the impressive results obtained using systems of imaging atmospheric Cherenkov telescopes. Ground-based gamma-ray astronomy has a huge potential in astrophysics, particle physics and cosmology. CTA is an international initiative to build the next generation instrument, with a factor of 5-10 improvement in sensitivity in the 100 GeV-10 TeV range and the extension to energies well below 100 GeV and above 100 TeV. CTA will consist of two arrays (one in the north, one in the south) for full sky coverage and will be operated as open observatory. The design of CTA is based on currently available technology. This document reports on the status and presents the major design concepts of CTA.
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.
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.
There is an observational correlation between astrophysical shocks and nonthermal particle distributions extending to high energies. As a first step toward investigating the possible feedback of these particles on the shock at the microscopic level, we perform particle-in-cell (PIC) simulations of a simplified environment consisting of uniform, interpenetrating plasmas, both with and without an additional population of cosmic rays. We vary the relative density of the counterstreaming plasmas, the strength of a homogeneous parallel magnetic field, and the energy density in cosmic rays. We compare the early development of the unstable spectrum for selected configurations without cosmic rays to the growth rates predicted from linear theory, for assurance that the system is well represented by the PIC technique. Within the parameter space explored, we do not detect an unambiguous signature of any cosmic-ray-induced effects on the microscopic instabilities that govern the formation of a shock. We demonstrate that an overly coarse distribution of energetic particles can artificially alter the statistical noise that produces the perturbative seeds of instabilities, and that such effects can be mitigated by increasing the density of computational particles.
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 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.