GeV Band Research Articles

MAGNETIC PAIR CREATION TRANSPARENCY IN GAMMA-RAY PULSARS

Magnetic pair creation $\gamma \to e^+e^-$ has been at the core of radio pulsar paradigms and central to polar cap models of gamma-ray pulsars for over three decades. The Fermi gamma-ray pulsar population now exceeds 140 sources and has defined an important part of Fermi's science legacy. Among the population characteristics well established is the common occurrence of exponential turnovers in their spectra in the 1--10 GeV range. These turnovers are too gradual to arise from magnetic pair creation in the strong magnetic fields of pulsar inner magnetospheres. By demanding insignificant photon attenuation precipitated by such single-photon pair creation, the energies of these turnovers for Fermi pulsars can be used to compute lower bounds for the typical altitude of GeV band emission. This paper explores such pair transparency constraints below the turnover energy, and updates earlier altitude bound determinations of that have been deployed in various Fermi pulsar papers. For low altitude emission locales, general relativistic influences are found to be important. Rotational aberration influences are also explored, and are found to be small at low altitudes, except near the magnetic pole. The altitude bounds are typically in the range of 2-7 stellar radii for the young Fermi pulsar population, and provide key information on the emission altitude in radio quiet pulsars that do not possess double-peaked pulse profiles. The bound for the Crab pulsar is at a much higher altitude, with the putative detection by MAGIC out to 350-400 GeV implying a lower bound of 310km to the emission region, i.e., approximately 20% of the light cylinder radius. These results are also extended to the super-critical field domain, where it is found that emission in magnetars originating below around 10 stellar radii will not appear in the Fermi-LAT band.

HIGH-ENERGY EMISSIONS FROM THE GAMMA-RAY BINARY LS 5039

We study mechanisms of multi-wavelength emissions (X-ray, GeV and TeV gamma-rays) from the gamma-ray binary LS~5039. This paper is composed of two parts. In the first part, we report on results of observational analysis using four year data of \fermi\ Large Area Telescope. Due to the improvement of instrumental response function and increase of the statistics, the observational uncertainties of the spectrum in $\sim$100-300 MeV bands and $>10$GeV bands are significantly improved. The present data analysis suggests that the 0.1-100GeV emissions from LS~5039 contain three different components; (i) the first component contributes to $<$1GeV emissions around superior conjunction, (ii) the second component dominates in 1-10GeV energy bands and (iii) the third component is compatible to lower energy tail of the TeV emissions. In the second part, we develop an emission model to explain the properties of the phase-resolved emissions in multi-wavelength observations. Assuming that LS~5039 includes a pulsar, we argue that both emissions from magnetospheric outer gap and inverse-Compton scattering process of cold-relativistic pulsar wind contribute to the observed GeV emissions. We assume that the pulsar is wrapped by two kinds of termination shock; Shock-I due to the interaction between the pulsar wind and the stellar wind and Shock-II due to the effect of the orbital motion. We propose that the X-rays are produced by the synchrotron radiation at Shock-I region and the TeV gamma-rays are produced by the inverse-Compton scattering process at Shock-II region.

ORIGIN OF MULTI-BAND EMISSION FROM THE MICROQUASAR CYGNUS X-1

We study the origin of non-thermal emissions from the Galactic black hole X-ray binary Cygnus X--1 which is a confirmed high mass microquasar. By analogy with those methods used in studies of active galactic nuclei, we propose a two-dimensional, time-dependent radiation model from the microquasar Cygnus X--1. In this model, the evolution equation for relativistic electrons in a conical jet are numerically solved by including escape, adiabatic and various radiative losses. The radiative processes involved are synchrotron emission, its self-Compton scattering, and inverse Compton scatterings of an accretion disk and its surrounding stellar companion. This model also includes an electromagnetic cascade process of an anisotropic $\gamma$-$\gamma$ interaction. We study the spectral properties of electron evolution and its emission spectral characteristic at different heights of the emission region located in the jet. We find that radio data from Cygnus X--1 are reproduced by the synchrotron emission, the \emph{Fermi} LAT measurements by the synchrotron emission and Comptonization of photons of the stellar companion, the TeV band emission fluxes by the Comptonization of the stellar photons. Our results show that: (1) Radio emission region extends from the binary system scales to the termination of the jet. (2) The GeV band emissions should originate from the distance close to the binary system scales. (3) The TeV band emissions could be inside the binary system, and these emissions could be probed by the upcoming CTA telescope. (4) The MeV tail emissions, which produces a strongly linearly polarized signal, are emitted inside the binary system. The location of the emissions is very close to the inner region of the jet.

EXTREME BLAZARS STUDIED WITHFERMI-LAT ANDSUZAKU: 1ES 0347–121 AND BLAZAR CANDIDATE HESS J1943+213

We report on our study of high-energy properties of two peculiar TeV emitters: the "extreme blazar" 1ES 0347-121 and the "extreme blazar candidate" HESS J1943+213 located near the Galactic Plane. Both objects are characterized by quiescent synchrotron emission with flat spectra extending up to the hard X-ray range, and both were reported to be missing GeV counterparts in the Fermi-LAT 2-year Source Catalog. We analyze a 4.5 year accumulation of the Fermi-LAT data, resulting in the detection of 1ES 0347-121 in the GeV band, as well as in improved upper limits for HESS J1943+213. We also present the analysis results of newly acquired Suzaku data for HESS J1943+213. The X-ray spectrum is well represented by a single power law extending up to 25 keV with photon index 2.00+/-0.02 and a moderate absorption in excess of the Galactic value, in agreement with previous X-ray observations. No short-term X-ray variability was found over the 80 ks duration of the Suzaku exposure. Under the blazar hypothesis, we modeled the spectral energy distributions of 1ES 0347-121 and HESS J1943+213, and derived constraints on the intergalactic magnetic field strength and source energetics. We conclude that although the classification of HESS J1943+213 has not yet been determined, the blazar hypothesis remains the most plausible option, since in particular the broad-band spectra of the two analyzed sources along with the source model parameters closely resemble each other, and the newly available WISE and UKIDSS data for HESS J1943+213 are consistent with the presence of an elliptical host at the distance of approximately ~600 Mpc.

Gamma-ray luminosity function of BL Lac objects

The gamma-ray luminosity function (GLF) of BL Lac objects is constructed by using a sample of BL Lac objects with redshifts selected from the Second LAT AGN catalog. The GLFs of BL Lacs in the frame of the pure density evolution (PDE), the pure luminosity evolution (PLE), and the luminosity-dependent density (LDDE) models are determined by using the Markov Chain Monte Carlo (MCMC) technique, respectively. Our results suggest that the PDE model can give best description for BL Lac GLF based on the combination of constraints of model parameters and good fits to the observed data of Fermi-Large Area Telescope (LAT) BL Lacs, but other two models (PLE and LDDE) cannot be excluded. Based on our constructed GLFs, the contribution to the extragalactic diffuse gamma-ray background (EGRB) from BL Lacs is estimated, and ∼1–5 per cent of the EGRB in the 0.1–100 GeV band is found to come from unresolved BL Lacs (including the cascade emission). In addition, it is found that the BL Lac GLF is very different from flat spectrum radio quasar GLF and then the contribution of blazars to the EGRB should be estimated separately.

DETAILED INVESTIGATION OF THE GAMMA-RAY EMISSION IN THE VICINITY OF SNR W28 WITHFERMI-LAT

We present a detailed investigation of the $\gamma$-ray emission in the vicinity of the supernova remnant (SNR) W28 (G6.4$-$0.1) observed by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope. We detected significant $\gamma$-ray emission spatially coincident with TeV sources HESS J1800$-$240A, B, and C, located outside the radio boundary of the SNR. Their spectra in the 2-100 GeV band are consistent with the extrapolation of the power-law spectra of the TeV sources. We also identified a new source of GeV emission, dubbed Source W, which lies outside the boundary of TeV sources and coincides with radio emission from the western part of W28. All of the GeV $\gamma$-ray sources overlap with molecular clouds in the velocity range from 0 to 20 km s$^{-1}$. Under the assumption that the $\gamma$-ray emission towards HESS J1800-240A, B, and C comes from $\pi^0$ decay due to the interaction between the molecular clouds and cosmic rays (CRs) escaping from W28, they can be naturally explained by a single model in which the CR diffusion coefficient is smaller than the theoretical expectation in the interstellar space. The total energy of the CRs escaping from W28 is constrained through the same modeling to be larger than $\sim$ 2 $\times$ 10$^{49}$ erg. The emission from Source W can also be explained with the same CR escape scenario.

Probing the cosmic-ray content of galaxy clusters by stackingFermi-LAT count maps

Observations in radio have shown that galaxy clusters are giant reservoirs of cosmic rays (CR). Although a gamma- ray signal from the cluster volume is expected to arise through interactions of CR protons with the ambient plasma, a confirming observation is still missing. We search for a cumulative gamma-ray emission in direction of galaxy clusters by analysing a collection of stacked Fermi-LAT count maps. Additionally, we investigate possible systematic differences in the emission between cool-core and non-cool core cluster populations. Making use of a sample of 53 clusters selected from the HIFLUGCS catalog, we do not detect a significant signal from the stacked sample. The upper limit on the average flux per cluster derived for the total stacked sample is at the level of a few 1e-11 ph cm-2 s-1 at 95% confidence level in the 1-300 GeV band, assuming power-law spectra with photon indices 2.0, 2.4, 2.8 and 3.2. Separate stacking of the cool core and non-cool core clusters in the sample lead to similar values of around 5e-11 ph cm-2 s-1 and 2e-11 ph cm-2 s-1, respectively. Under the assumption that decaying pi0, produced in collisions between CRs and the ambient thermal gas, are responsible for the gamma-ray emission, we set upper limits on the average CR content in galaxy clusters. For the entire cluster population, our upper limit on the gamma-ray flux translates into an upper limit on the average CR-to-thermal energy ratio of 4.6% for a photon index of 2.4, although it is possible for individual systems to exceed this limit. Our 95% upper limits are at the level expected from numerical simulations, which likely suggests that the injection of CR at cosmological shocks is less efficient than previously assumed.

Constraining the parameter space of the one-zone synchrotron-self-Compton model for GeV-TeV detected BL Lacertae objects

The one-zone synchrotron-self-Compton (SSC) model aims to describe the spectral energy distribution (SED) of BL Lac objects via synchrotron emission by a non-thermal population of electrons and positrons in a single homogeneous emission region, partially upscattered to gamma-rays by the particles themselves. The model is usually considered as degenerate, given that the number of free parameters is higher than the number of observables. It is thus common to model the SED by choosing a single set of values for the SSC-model parameters that provide a good description of the data, without studying the entire parameter space. We present here a new numerical algorithm which permits us to find the complete set of solutions, using the information coming from the detection in the GeV and TeV energy bands. The algorithm is composed of three separate steps: we first prepare a grid of simulated SEDs and extract from each SED the values of the observables; we then parametrize each observable as a function of the SSC parameters; we finally solve the system for a given set of observables. We iteratively solve the system to take into account uncertainties in the values of the observables, producing a family of solutions. We present a first application of our algorithm to the typical high-frequency-peaked BL Lac object 1RXS J101015.9-311909, provide constraints on the SSC parameters, and discuss the result in terms of our understanding of the blazar emitting region.

THE FERMI BUBBLES AS A SCALED-UP VERSION OF SUPERNOVA REMNANTS

In this study, we treat the Fermi bubbles as a scaled-up version of supernova remnants (SNRs). The bubbles are created through activities of the super-massive black hole (SMBH) or starbursts at the Galactic center (GC). Cosmic-rays (CRs) are accelerated at the forward shocks of the bubbles like SNRs, which means that we cannot decide whether the bubbles were created by the SMBH or starbursts from the radiation from the CRs. We follow the evolution of CR distribution by solving a diffusion-advection equation, considering the reduction of the diffusion coefficient by CR streaming. In this model, gamma-rays are created through hadronic interaction between CR protons and the gas in the Galactic halo. In the GeV band, we can well reproduce the observed flat distribution of gamma-ray surface brightness, because some amount of gas is left behind the shock. The edge of the bubbles is fairly sharp owing to the high gas density behind the shock and the reduction of the diffusion coefficient there. The latter also contributes the hard gamma-ray spectrum of the bubbles. We find that the CR acceleration at the shock has started when the bubbles were small, and the time-scale of the energy injection at the GC was much smaller than the age of the bubbles. We predict that if CRs are accelerated to the TeV regime, the apparent bubble size should be larger in the TeV band, which could be used to discriminate our hadronic model from other leptonic models. We also present neutrino fluxes.

What caused the GeV flare of PSR B1259-63?

PSR B1259-63 is a gamma-ray binary system composed of a high spindown pulsar and a massive star. Non-thermal emission up to TeV energies is observed near periastron passage, attributed to emission from high energy e+e- pairs accelerated at the shock with the circumstellar material from the companion star, resulting in a small-scale pulsar wind nebula. Weak gamma-ray emission was detected by the Fermi/LAT at the last periastron passage, unexpectedly followed 30 days later by a strong flare, limited to the GeV band, during which the luminosity nearly reached the spindown power of the pulsar. The origin of this GeV flare remains mysterious. We investigate whether the flare could have been caused by pairs, located in the vicinity of the pulsar, up-scattering X-ray photons from the surrounding pulsar wind nebula rather than UV stellar photons, as usually assumed. Such a model is suggested by the geometry of the interaction region at the time of the flare. We compute the gamma-ray lightcurve for this scenario, based on a simplified description of the interaction region, and compare it to the observations. The GeV lightcurve peaks well after periastron with this geometry. The pairs are inferred to have a Lorentz factor ~500. They also produce an MeV flare with a luminosity ~1e34 erg/s prior to periastron passage. A significant drawback is the very high energy density of target photons required for efficient GeV emission. We propose to associate the GeV-emitting pairs with the Maxwellian expected at shock locations corresponding to high pulsar latitudes, while the rest of the non-thermal emission arises from pairs accelerated in the equatorial region of the pulsar wind termination shock.

H.E.S.S. discovery of VHEγ-rays from the quasar PKS 1510−089

The quasar PKS 1510-089 (z=0.361) was observed with the H.E.S.S. array of imaging atmospheric Cherenkov telescopes during high states in the optical and GeV bands, to search for very high energy (VHE, defined as E >= 0.1 TeV) emission. VHE \gamma-rays were detected with a statistical significance of 9.2 standard deviations in 15.8 hours of H.E.S.S. data taken during March and April 2009. A VHE integral flux of I(0.15 TeV < E < 1.0 TeV) = (1.0 +- 0.2 (stat) +- 0.2 (sys) x 10^{-11} cm^{-2}s^{-1} is measured. The best-fit power law to the VHE data has a photon index of \Gamma=5.4 +- 0.7 (stat) +- 0.3 (sys). The GeV and optical light curves show pronounced variability during the period of H.E.S.S. observations. However, there is insufficient evidence to claim statistically significant variability in the VHE data. Because of its relatively high redshift, the VHE flux from PKS 1510-089 should suffer considerable attenuation in the intergalactic space due to the extragalactic background light (EBL). Hence, the measured \gamma-ray spectrum is used to derive upper limits on the opacity due to EBL, which are found to be comparable with the previously derived limits from relatively-nearby BL Lac objects. Unlike typical VHE-detected blazars where the broadband spectrum is dominated by non-thermal radiation at all wavelengths, the quasar PKS 1510-089 has a bright thermal component in the optical to UV frequency band. Among all VHE detected blazars, PKS 1510-089 has the most luminous broad line region (BLR). The detection of VHE emission from this quasar indicates a low level of \gamma-\gamma absorption on the internal optical to UV photon field.

A NEW RESULT ON THE ORIGIN OF THE EXTRAGALACTIC GAMMA-RAY BACKGROUND

In the paper, we continually use the method of image stacking to study the origin of the extragalactic gamma-ray background (EGB) at GeV bands, and find that the Faint Images of the Radio Sky at Twenty centimeters (FIRST) sources undetected by the Large Area Telescope on the Fermi Gamma-ray Space Telescope can contribute about (56{plus minus}6)% of the EGB. Because the FIRST is a flux limited sample of radio sources with incompleteness at the faint limit, we consider that the point-sources, including blazars, non-blazar AGNs,starburst galaxies, could produce a much larger fraction of the EGB.

VARIABILITY OF GAMMA-RAY EMISSION FROM BLAZARS ON BLACK HOLE TIMESCALES

We investigate the variability properties of blazars in the GeV band using the data of the Fermi/LAT telescope. We find that blazars exhibit variability on the scales down to the minimal timescale resolvable by Fermi, which is a function of the peak photon count rate in the LAT. This implies that the real minimal variability timescales for the majority of blazars are typically shorter than those resolvable by the LAT. We find that for several blazars these minimal variability timescales reach those associated to the blazar central engine, the supermassive black hole. At the same time, none of the blazars exhibits variability on the timescale shorter than the black hole horizon light crossing time and/or the period of rotation around the last stable circular orbit. Based on this fact, we argue that the timing properties of the gamma-ray signal could be determined by the processes in the direct vicinity of the supermassive black hole.

The average GeV-band emission from gamma-ray bursts

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. To estimate the emission of faint GRBs we use non-standard analysis methods and sum over many GRBs to find an average signal which is significantly above background level. We use a subsample of 99 GRBs listed in the Burst Catalog from the first two years of observation. Although mostly 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*10^-4 erg cm^-2 yr^-1 or only 4% of the energy flux of cosmic rays above the ankle at 10^18.6 eV.

INVERSE COMPTON SCATTERING MODEL FOR X-RAY EMISSION OF THE GAMMA-RAY BINARY LS 5039

We propose a model for the gamma-ray binary LS 5039 in which the X-ray emission is due to the inverse Compton (IC) process instead of the synchrotron radiation. Although the synchrotron model has been discussed in previous studies, it requires a strong magnetic field which leads to a severe suppression of the TeV gamma-ray flux in conflict with H.E.S.S. observations. In this paper, we calculate the IC emission by low energy electrons (\gamma_e \lesssim 10^3) in the Thomson regime. We find that IC emission of the low energy electrons can explain the X-ray flux and spectrum observed with Suzaku if the minimum Lorentz factor of injected electrons \gamma_min is around 10^3. In addition, we show that the Suzaku light curve is well reproduced if \gamma_min varies in proportion to the Fermi flux when the distribution function of injected electrons at higher energies is fixed. We conclude that the emission from LS 5039 is well explained by the model with the IC emission from electrons whose injection properties are dependent on the orbital phase. Since the X-ray flux is primarily determined by the total number of cooling electrons, this conclusion is rather robust, although some mismatches between the model and observations at the GeV band remain in the present formulation.

Exploring the extreme universe with the Fermi Gamma-Ray Space Telescope

The Fermi orbiter, mapping the entire gamma-ray sky every three hours, monitors the cosmos for high-energy phenomena both fleeting and enduring.

MULTIWAVELENGTH OBSERVATIONS OF THE AGN 1ES 0414+009 WITH VERITAS,FERMI-LAT,SWIFT-XRT, AND MDM

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 January 2008 and February 2011, resulting in 56.2 hours 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 * 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 an 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 * 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.

THE LUMINOSITY FUNCTION OFFERMI-DETECTED FLAT-SPECTRUM RADIO QUASARS

Fermi has provided the largest sample of {\gamma}-ray selected blazars to date. In this work we use a complete sample of FSRQs detected during the first year of operation to determine the luminosity function (LF) and its evolution with cosmic time. The number density of FSRQs grows dramatically up to redshift \sim0.5-2.0 and declines thereafter. The redshift of the peak in the density is luminosity dependent, with more luminous sources peaking at earlier times; thus the LF of {\gamma}-ray FSRQs follows a luminosity-dependent density evolution similarly to that of radio-quiet AGN. Also using data from the Swift Burst Alert Telescope we derive the average spectral energy distribution of FSRQs in the 10 keV-100 GeV band and show that there is no correlation of the peak {\gamma}-ray luminosity with {\gamma}-ray peak frequency. The coupling of the SED and LF allows us to predict that the contribution of FSRQs to the Fermi isotropic {\gamma}-ray background is 9.3(+1.6/-1.0) (\pm3% systematic uncertainty) in the 0.1-100GeV band. Finally we determine the LF of unbeamed FSRQs, finding that FSRQs have an average Lorentz factor of {\gamma} = 11.7(+3.3/-2.2), that most are seen within 5\circ of the jet axis, and that they represent only ~0.1 % of the parent population.

γ-RAY LUMINOSITY FUNCTION AND THE CONTRIBUTION TO EXTRAGALACTIC γ-RAY BACKGROUND FORFERMI-DETECTED BLAZARS

We study γ-ray luminosity functions (GLFs) and the contribution to the extragalactic diffuse γ-rays for blazars. After using a clean blazar sample in which 485 blazars are observed at both the 8.4 GHz and ∼ GeV bands, we investigate the radio/γ-ray luminosity correlation and find out that a strong correlation exists between the radio and γ-ray luminosities using the partial correlation analysis method. Based on such a correlation, we obtain the GLF, while also considering the GLF based on studies of the X-ray luminosity function of active galactic nuclei and the GLF of the EGRET blazars. K-S tests indicate that these three kinds of GLFs are consistent with the observed data of our sample in a reasonable parameter range. Using these GLFs, we further investigate the blazar contribution to the extragalactic diffuse γ-ray background (EGRB). Our results show that (1) the contribution of blazars are ∼21.3% ±2.7%, ∼12.9% ±1.7%, and ∼37.9% ±4.8% of unresolved Fermi EGRB flux for the three GLFs and (2) the contribution of flat-spectrum radio quasars are ∼11.2% ±1.6%, ∼6.9% ±1.0%, and ∼9.0% ±1.2% of unresolved Fermi EGRB flux for the three GLFs.

JITTER SELF-COMPTON PROCESS: GeV EMISSION OF GRB 100728A

Jitter radiation, the emission of relativistic electrons in a random and small-scale magnetic field, has been applied to explain the gamma-ray burst (GRB) prompt emission. The seed photons produced from jitter radiation can be scattered by thermal/nonthermal electrons to the high-energy bands. This mechanism is called jitter self-Compton (JSC) radiation. GRB 100728A, which was simultaneously observed by the Swift and Fermi, is a great example to constrain the physical processes of jitter and JSC. In our work, we utilize jitter/JSC radiation to reproduce the multiwavelength spectrum of GRB 100728A. In particular, due to JSC radiation, the powerful emission above the GeV band is the result of those jitter photons in X-ray band scattered by the relativistic electrons with a mixed thermal-nonthermal energy distribution. We also combine the geometric effect of microemitters to the radiation mechanism, such that the "jet-in-jet" scenario is considered. The observed GRB duration is the result of summing up all of the contributions from those microemitters in the bulk jet.