Fermi LAT Observations of the Supernova Remnant W28 (G6.4-0.1)

To be re- submitted to ApJ: v5.7, May 26, 2010

Fermi LAT Observations of the Supernova Remnant W28 (G6.4-0.1)

^ : A. A. Abdo^ '^ M. Ackermann\ M. Ajello'^, A. AUafort^, L. Baldini^ J. Ballet^ G. Barbiellini^ ^

^ ; D. Bastieri''"', K. Bechtol\ R. Bellazzini^ B. Berenji"^, R. D. Blandford'^, E. D. Bloom\

>^: E. Bonamente"'!^ A. W. Borgland'^, A. Bouvier"^, T. J. Brandt^^-^\ J. Bregeon^ M. Brigida^^'^^

^ : R Bruel^^ R. Buehler^, S. Buson'', G. A. Caliandrol^ R. A. Cameron^, R A. Caraveoi'',

2 ; S. Carrigan^°, J. M. Casandjian^ C. Cecchi'^''^, O. gelik20'2i.22^ A. Chekhtman^^^ j. Chiang\

^ : S. Ciprini^^^ R. Claus^, J. Cohen-Tanugi2\ J. Conrad^^^e^?^ q D Dermer^, R de Palma'^.ie^

: E. do Couto e Silva\ R S. Drell^, R. Dubois'*, D. Dumora^^^'', C. Farnier^^, C. Favuzzi^^■^^

[ij; S. J. Fegan", Y. Fukazawa^°, Y. Fukui^\ S. Funk"*, R Fusco*^■*^ F. Gargano*^ N. Gehrels^",

X : S. Germani*''*^ N. Giglietto^5'^^ E Giordano'5''^ T. Glanzman\ G. Godfrey\ I. A. Grenier^

^ : J. E. Grovel S. Guiriec^^ D. Hadasch^^ Y Hanabata^o, A. K. Harding^", E. Hays^", D. Horan",

7^1 R. E. Hughes^"*, G. Johannesson'*, A. S. Johnson"*, W. N. Johnson^, T. Kamae"*, H. Katagiri^"*,

2 ! J- Kataoka^"*, J. Knodlseder", M. Kuss^, J. Lande"*, L. Latronico^, S.-H. Lee"*,

c/3 \ M. Lemoine-Goumard^^'^^, M. Llena Garde^^'^^, F. Longo^*^, F. Loparco*^*'', M. N. Lovellette^,

I — \\ P. Lubrano*^'*^, A. Makeev^'^\ M. N. Mazziotta^*", R F. Michelson"*, W. Mitthumsiri"*,

^^ ; T. Mizuno^°, A. A. Moiseev^*'^^ C. Monte*^'^^ M. E. Monzani'*, A. Morselli^^

> ; I. V. Moskalenko^, S. Murgia^ T. Nakamori^^ R L. Nolan\ J. R Norris^^ E. Nuss^^, M. Ohno^^

^ ; T. Ohsugi^^ N. Omodei^, E. Orlando^i, J. F Ormes3^ M. Ozaki^^ J. H. Panetta\

^ ; D. Parent2'23.28,29^ Y peiassa^"*, M. Pepe*''*^ M. Pesce-Rollins^ F. Piron^^, T. A. Porter"*,

^ ; S. Rain6'^''^ R. Rando''"', M. Razzano^ A. Reimer"*^^ O. Reimer^^'^ T. Reposeur^^'^^,

O ; A. Y Rodriguez'^ M. Roth^^, H. R-W. Sadrozinski^^, A. Sander^^^ R M. Saz Parkinson4^

O

X

C. Sgr6^ E. J. Siskind4^ D. A. Smith^^^^ R D. Smith'^^ G. Spandre^ R Spinelli^^^'^

M. S. Strickman^, D. J. Suson^^, H. Tajima^^ H. Takahashi^o, T. Takahashi^^ T. Tanaka"*'\

J. B. Thayer^, J. G. Thayer^, D. J. Thompson^", L. Tibaldo'^'^°''^'"*\ O. Tibolla4^ D. E Torres33'l^

H : G. Tosti"'i2, Y Uchiyama^'i, T. Uehara^^, T. L. Usher^ V. Vasileiou^i 22^ N. Vilchezi\

V. Vitale^^'^^ A. R Waite^, R Wang^, B. L. Winer'^, K. S. Wood^, H. Yamamoto^*,

R. Yamazaki^^", Z. Yang^^^e^ t. Ylinen^^^^i-^e^ y^ Ziegler'

.44

'Corresponding authors: H. Katagiri, katagiri@hep01.hepl.hiroshima-u.ac.jp; H. Tajima, hta-

jima@slac.stanford.edu; T. Tanaka, ttanaka@slac.stanford.edu; Y. Uchiyama, uchiyama@slac.stanford.edu.

^Space Science Division, Naval Research Laboratory, Washington, DC 20375, USA

^National Research Council Research Associate, National Academy of Sciences, Washington, DC 20001, USA

"*W. W. Hansen Experimental Physics Laboratory, Kavli Institute for Particle Astrophysics and Cosmology, De-
partment of Physics and SLAC National Accelerator Laboratory, Stanford University, Stanford, CA 94305, USA

^Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, 1-56127 Pisa, Italy

^Laboratoire AIM, CEA-IRFU/CNRS/Universite Paris Diderot, Service d' Astrophysique, CEA Saclay, 91 191 Gif
sur Yvette, France

^Istituto Nazionale di Fisica Nucleare, Sezione di Trieste, 1-34127 Trieste, Italy

^Dipartimento di Fisica, Universita di Trieste, 1-34127 Trieste, Italy

'istituto Nazionale di Fisica Nucleare, Sezione di Padova, 1-35131 Padova, Italy

'"Dipartimento di Fisica "G. Galilei", Universita di Padova, 1-35131 Padova, Italy

^'istituto Nazionale di Fisica Nucleare, Sezione di Perugia, 1-06123 Perugia, Italy

'^Dipartimento di Fisica, Universita degli Studi di Perugia, 1-06123 Perugia, Italy

'^Centre d'Etude Spatiale des Rayonnements, CNRS/UPS, BP 44346, F-30128 Toulouse Cedex 4, France

'"^Department of Physics, Center for Cosmology and Astro-Particle Physics, The Ohio State University, Columbus,
OH 43210, USA

'^Dipartimento di Fisica "M. Merlin" dell' Universita e del Politecnico di Bari, 1-70126 Bari, Italy

'^Istituto Nazionale di Fisica Nucleare, Sezione di Bari, 70126 Bari, Italy

'^Laboratoire Leprince-Ringuet, Ecole polytechnique, CNRS/IN2P3, Palaiseau, France

'^Institut de Ciencies de I'Espai (lEEC-CSIC), Campus UAB, 08193 Barcelona, Spain

''INAF-Istituto di Astrofisica Spaziale e Fisica Cosmica, 1-20133 Milano, Italy

^^NASA Goddard Space FUght Center, Greenbelt, MD 20771, USA

^'Center for Research and Exploration in Space Science and Technology (CRESST) and NASA Goddard Space
Fhght Center, Greenbelt, MD 20771, USA

^^Department of Physics and Center for Space Sciences and Technology, University of Maryland Baltimore County,
Baltimore, MD 21250, USA

^^George Mason University, Fairfax, VA 22030, USA

^"^Laboratoire de Physique Theorique et Astroparticules, Universite Montpellier 2, CNRS/IN2P3, Montpellier,
France

^^Department of Physics, Stockholm University, AlbaNova, SE-106 91 Stockholm, Sweden

ABSTRACT

^^The Oskar Klein Centre for Cosmoparticle Physics, AlbaNova, SE-106 91 Stockholm, Sweden

^^Royal Swedish Academy of Sciences Research Fellow, funded by a grant from the K. A. Wallenberg Foundation

2**CNRS/IN2P3, Centre d'Etudes Nucleaires Bordeaux Gradignan, UMR 5797, Gradignan, 33175, France

^^Universite de Bordeaux, Centre d'Etudes Nucleaires Bordeaux Gradignan, UMR 5797, Gradignan, 33 175, France

^"Department of Physical Sciences, Hiroshima University, Higashi-Hiroshima, Hiroshima 739-8526, Japan

^'Department of Physics and Astrophysics, Nagoya University, Chikusa-ku Nagoya 464-8602, Japan

^^Center for Space Plasma and Aeronomic Research (CSPAR), University of Alabama in Huntsville, Huntsville,
AL 35899, USA

^^Institucio Catalana de Recerca i Estudis Avan9ats (ICREA), Barcelona, Spain

^"^Research Institute for Science and Engineering, Waseda University, 3-4-1, Okubo, Shinjuku, Tokyo, 169-8555
Japan

^^Department of Physics and Department of Astronomy, University of Maryland, College Park, MD 20742, USA

^^Istituto Nazionale di Fisica Nucleare, Sezione di Roma "Tor Vergata", 1-00133 Roma, Italy

^^Department of Physics, Tokyo Institute of Technology, Meguro City, Tokyo 152-8551, Japan

^^Department of Physics and Astronomy, University of Denver, Denver, CO 80208, USA

^^Institute of Space and Astronautical Science, JAXA, 3-1-1 Yoshinodai, Sagamihara, Kanagawa 229-8510, Japan

^•^Hiroshima Astrophysical Science Center, Hiroshima University, Higashi-Hiroshima, Hiroshima 739-8526, Japan

"^'Max-Planck Institut fiir extraterrestrische Physik, 85748 Garching, Germany

"'^Institut fiir Astro- und Teilchenphysik and Institut fiir Theoretische Physik, Leopold-Franzens-Universitat Inns-
bruck, A-6020 Innsbruck, Austria

"^^Department of Physics, University of Washington, Seatde, WA 98195-1560, USA

"'"^Santa Cruz Institute for Particle Physics, Department of Physics and Department of Astronomy and Astrophysics,
University of CaHfornia at Santa Cruz, Santa Cruz, CA 95064, USA

"♦^NYCB Real-Time Computing Inc., Lattingtown, NY 1 1560-1025, USA

"'''Department of Chemistry and Physics, Purdue University Calumet, Hammond, IN 46323-2094, USA

''^Partially supported by the International Doctorate on Astroparticle Physics (IDAPP) program

"'^Institut fiir Theoretische Physik and Astrophysik, Universitat Wiirzburg, D-97074 Wiirzburg, Germany

"'^Dipartimento di Fisica, Universita di Roma "Tor Vergata", 1-00133 Roma, Italy

^"^Department of Physics, Royal Institute of Technology (KTH), AlbaNova, SE-106 91 Stockholm, Sweden

^'School of Pure and Applied Natural Sciences, University of Kalmar, SE-391 82 Kalmar, Sweden

We present detailed analysis of the two gamma-ray sources, IFGL J1801.3-2322c
and IFGL J1800.5-2359c, that have been found toward the supernova remnant (SNR)
W28 with the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space
Telescope. IFGL J1801.3-2322c is found to be an extended source within the
boundary of SNR W28, and to extensively overlap with the TeV gamma-ray source
HESS J1801-233, which is associated with a dense molecular cloud interacting
with the supernova remnant. The gamma-ray spectrum measured with LAT from
0.2-100 GeV can be described by a broken power-law function with a break of
~ 1 GeV, and photon indices of 2.09 ± 0.08 (stat) ± 0.28 (sys) below the break and
2.74 ± 0.06 (stat) ± 0.09 (sys) above the break. Given the clear association between
HESS Jl 801-233 and the shocked molecular cloud and a smoothly connected spec-
trum in the GeV-TeV band, we consider the origin of the gamma-ray emission in
both GeV and TeV ranges to be the interaction between particles accelerated in the
SNR and the molecular cloud. The decay of neutral pions produced in interactions
between accelerated hadrons and dense molecular gas provide a reasonable explana-
tion for the broadband gamma-ray spectrum. IFGL J1800.5-2359c located outside
the southern boundary of SNR W28, cannot be resolved. An upper limit on the size
of the gamma-ray emission was estimated to be ~ 16' using events above ~ 2 GeV
under the assumption of a circular shape with uniform surface brightness. It appears
to coincide with the TeV source HESS J1800-240B, which is considered to be associ-
ated with a dense molecular cloud that contains the ultra compact H II region W28A2
(G5. 89-0. 39). We found no significant gamma-ray emission in the LAT energy band
at the positions of TeV sources HESS J1800-230A and HESS J1800-230C. The LAT
data for HESS J1800-230A combined with the TeV data points indicate a spectral
break between 10 GeV and 100 GeV.

Subject headings: cosmic rays — acceleration of particles — ISM: individual objects
(W28, G6. 4-0.1) — ISM: supernova remnants — gamma rays: ISM

1. Introduction

Diffusive shock acceleration operating at supernova shock waves can distribute particles to

very h igh energies with a power-law form having number index about 2 (e.g. jBlandford & Eichlei

19871). It is generally expected that if a dense molecular cloud is overtaken by a supernova blast

wave, the sho cked molecular cloud c an be illuminated by relativistic particles accelerated at super-

nova shocks (|Aharonian et al.lll994|) . If the accelerated particles are comprised mostly of protons.

say > 100 times more abundant than electrons like the observed Galactic cosmic rays, decays of

neutral pions produced in inelastic collisions of the accelerated protons with dense gas are expected
to be a dominant radiation component in the gamma-ray spectrum of the cosmic-ray-illuminated
molecular cloud. Although an earlier attempt to detect TeV gamma-ray emission from supernova
remna nts (SNRs) that have evidence for molecular cloud interactions with the Whipple telescope
failed (jBu cklev et al.lll998b. two arch etvpical S NRs interacting with m olecular clouds, IC 443
jAlbert et al...2007l : lAcciari et alJboog I) and W28 jAharonian et allboosl) . have been detected with
the current generation of imaging atmospheric Cherenkov telescopes. However, the identities of
the particles responsible for the TeV sources remains elusive.

The advent of the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Tele-
scope has brought a new opportunity to study the gamma-ray emission from SNRs at G eV energies.

LAT observa tions of middle-ag ed SNRs inte racting with molec ular clouds, W51C (|Abdo et al

20091), W44 (lAbdo et al.ll2010ah . and IC 443 (lAbdoet all 1201 Obh . have revealed bright extended
gamma-ray sources coincident with the SNRs. The gamma-ray luminosity reaches ~ 10^^ erg s~\
which immediately rules out an inverse-Compton origin of the GeV gamma rays since it requires
total electron energy comparable to or larger than the typical kinetic energy released by a supernova
explosion, ~ 10^^ erg. The gamma-ray spectra of the three remnants exhibit remarkable spectral
breaks at an energy of several GeV, making these SNRs much less luminous at TeV energies. This
characteristic demonstrates the importance of observations at GeV energies.

W28 is a mixed-morphology SNR, characterized by center-filled thermal X-ray emission and
shell-like radio morpholo gy. In addition X-ray ob servations show limb-brightened shells in the
northeast and southwest (IRho & Borkowskil |2002|) . The shell-like radio emission is prominent

in the northeastern region with slightly fainter emission a t the northern b oundaries (|Dubner et al.

2000|) . Interactions of the SNR with molecular clouds dWoottenlll98lh along its no rthern anc
northeastern boundaries are traced by the high concentration of 1720 MHz OH mas ers (IFrail et al.

1994: Claussen et al.

1997l . ll999h . and high density (hh > 10^ cm ^) shocked gas (lArikawa et al.

19991 : lReachetal.ll2005b . The overall shape of W28 is elliptical with a size of 50' x 45'. W28
is located within a complex, star-forming region along the Gal actic plane to ward the large H II
regions (M8 and M20) and young clusters (e.g., NGC 6530) (iGoudisI 19761). The o bservations

of molecular lines place SNR W28 at a distanc e of ~ 2 kpc (I Velazquez et al. Il2002() . Estimates
for its age vary between 3.5 and 15 x 10"* yrs (IKaspi et al.lll993h. W28 is considered to be an

evolved remnant in th e radiative stage o f SNR evolution (|Lozinskay all 1 992() . which is supported by
optical observations dLozinskay al 1 1 974|) . Measurements with Energetic Gamma-Ray Experiment
Telescope (EGRE T) onboard the Comp ton Gamma-ray Observatory found a gamma-ray source,
3EG J 1800-233 8 faartman et al.lll999l) in the W28 field. However, its association to SNR W28
remained unclear mainly due to large source location uncertainties from EGRET. A gamma-ray
source is listed i n the W28 field in the AGILE (Astro-rivelatore Gamma a Immagini LEggero)
one year catalog (|Pittori et al.ll2009|) . However, detailed analysis of this field is not published by

-6

AGILE yet.

H.E.S.S. observations of the W2 8 field have revealed fo ur TeV gamma-ray sources position-
ally coincident with molecular clouds (lAharonian et alJlOOSh : HESS Jl 801-233, located along the

northeastern boundary of W28, and a complex of sources, HESS J1800-240A, B and C, located
~ 30' south of SNR W28. HESS J1801-233 coincides with a molecular cloud interacting with
SNR W28, providing one of the best examples of a cosmic -ray-illuminated cloud. Understanding
the origins of TeV emission in HESS J 1800-240 ABC is of particular interest; they may be due to
cosmic rays that have diffused from W28.

In this paper, we report Fermi LAT observations of gamma-ray sources in the W28 field in
the GeV domain. First, we give a brief description of the observation and gamma-ray selection
in Section [2l The analysis procedures and results are explained in Section [3l where the spatial
extension and spectra of the LAT sources in the W28 field are described. Discussion is given in
Section in followed by conclusions in Section [51

2. OBSERVATION AND DATA REDUCTION

The LAT is the main instrument on Fermi sensitive to gamma rays. The energy range spans
from ~ 20 MeV to > 300 GeV, although as noted below in the present analysis we use only events
with energies > 200 MeV. It is an electron-positron pair production telescope, built with tungsten
foils and silicon microstrip detectors to measure the arrival directions of incoming gamma rays, and
a hodoscopic cesium iodide calorimeter to determine the gamma-ray energies. They are surrounded
by 89 segmented plastic scintillators that serve as an anticoincidence detector to reject charged
particle eve nts. Details of the LAT instrument and pre-launch expectations of the performance can

be found in lAtwood et al.l (120091) . Relative to earlier gamma-ray missions, the LAT has a large
~ 2.4 sr field of view, a large effective area (~ 8000 cm^ for >1 GeV if on-axis) and improved
angular resolution or point-spread function (PSF, better than 1° for 68% containment at 1 GeV).

Routine science operations of the LAT began on August 4, 2008, after the conclusion of a
commissioning period. We have analyzed events in the W28 field, collected from August 4, 2008,
to July 4, 2009, with a total exposure of ~ 2.8 x 10^° cm^ s (at 1 GeV). During this time interval,
the LAT was operated in sky survey mode nearly all of the time. In this observing mode the LAT
scans the sky, obtaining complete sky coverage every 2 orbits (~ 3 hr) and relatively uniform
exposures over time.

We used the standard LAT analysis software, ScienceTools version v9rl5, which is available

-7

from the Fermi Science Support Center (FSSCU and applied the following event selection criteria:
a) events have the highest probability of being gamma rays, i.e., they should be classified as so-

called Pass 6 diffuse class (JAtwood et al.ll2009|) . b) the reconstructed zenith angles of the arrival

direction of gamma rays should be less than 105°, to minimize contamination from Earth albedo
gamma rays, c) the center of the LAT field of view should be within 39° from the zenith in order
to exclude data from the short time intervals when the field of view can be partly occulted by
the earth. There are no gamma-ray bursts detected by the LAT within 15° of the W28 field, thus
we did not need to apply any additional time cut. The energy range analyzed here is restricted
to > 200 MeV to avoid possible large systematic uncertainties due to the strong Galactic diffuse
emission near the Galactic center, smaller effective area, and much broader PSF at lower energies.

3. ANALYSIS AND RESULTS

3.1. Source position and Spatial Extension

Figure [H shows a smoothed count map in the 2-10 GeV energy band in a 10° x 10° region
around W28. Figure [2]gives comparisons with images of other wavebands in close-up view. Black
contours indicate the H.E.S.S. significance map in (a), CO (J=l-0) line intensity taken by NAN-
TEN for the velocity range from km s"^ to 20 km s "' (corresponding to kinematic distances of
approximately to 4 kpc) in (b), and a VLA image in (c). Correlations between GeV gamma rays
observed by LAT and some of the H.E.S.S. sources are evident. Gamma rays are also bright in the
brightest spots in the CO contours, and bright regions in the gamma-ray image extensively overlap
with bright reg ions in the VLA co ntours. There are two LAT sources in the vicinity of W28 in the

IFGL catalog (lAbdo et al.ll2010dl) : IFGL J1801.3-2322c, and IFGL J1800.5-2359c. No obvious
pulsations of gamma rays are found in these sources. Hereafter we refer to IFGL J1801.3-2322c
as Source N, and to IFGL J1800.5-2359c as Source S.

In order to quantitatively evaluate the extension and location of these two sources, we apply
the maximum likelihood tool, gt like, which is publicly available as part of Fermi ScienceTools.
The likelihood is the product of the probability of observing the gamma-ray counts of each spatial
and energy bin given the emission model, and the best parameter values are estimated by maximiz-
ing the likelihood of the data given the model (IMattox et al.ll996b . The probability density function

for the likelihood analysis included a) individual sources detected in the preliminary LAT 1-year
catalog, b) the Galactic diffuse emission resulting from cosmic-ray interactions with the interstellar

'Software and documentation of the Fermi ScienceTools are distributed by Fermi Science Support Center at
[http://fermi.gsfc.nasa.gov/ssc

-8

medium and radiation based on the LAT standard diffuse background model gll_iem_v02 available
from FSSCq and c) the isotropic component to represent extragalactic and residual cosmic -ray
backgrounds using the isotropic spectral template isotropic_iem_v02 from FSSC. Note that we
make energy dependent corrections of the the Galactic diffuse model by multiplying a power law
function with the spectral index free to vary in the fit. This correction gives better spectral fits by
taking account of local systematic discrepancies between the data and the Galactic diffuse model.
The region of interest for the binned maximum likelihood analysis was a square region of 20° x 20°
centered on W28 with a pixel size of 0.°1. The instrument response functions (IRFs) used in
our work were the "Pass 6 v3" (P6_V3) IRFs, which were developed following launch to address
gamma-ray detection inefficiencies that are correlated with background rates.

Using the tools described above, we investigated spatial extensions of the two LAT W28
sources. Here we used only events above 2.15 GeV to take advantage of narrower PSF in the higher
energy band. Before investigating the detailed extension, we first determined the flux and spectral
index for all components except for Source N and Source S. In this process, the normalization of
the Galactic diffuse emission and the flux and spectral index of power-law model for the sources
within 5 degrees of Source N were set free to account for the effects of sources around W28
on the background flux in the fit. The flux and spectral index except for those of Source N and
Source S are fixed hereafter. We used a radially symmetric uniform disk to evaluate the extension.
We varied the radius and location of the disk while holding the position of Source S fixed at the
catalog position, and evaluated the resulting maximum likelihood value (Lex) with respect to the
maximum likelihood for no source hypothesis (Lq) and the point source hypothesis (Lps). A two-
source hypothesis (L2s) is also evaluated by scanning positions of two point sources.

The best likelihood ratio -21n(Lps/Lex) ~ 76 is obtained for Source N with a disk radius of
a = 0.°39, which rejects a point source hypothesis at more than 8 a. Note that we v alidated this

proced ure by applying this method on a nearby gamma-ray pulsar, PSR J 1809-2332 (lAbdo et al.

2010cl) (~ 2° away from W28), where we find the extension to be consistent with a point source.
In addition, the best likelihood ratio -21n(Lo/Lex) ~ 336 for the disk shape is substantially better
than that for two point sources, -21n(Lo/L2s) ~ 276, where the positions of the two sources were
free in the optimization. Therefore we use the best-fit disk model for further analysis of Source N.
The best location of the disk model is found at (R.A., Dec.) = (18''0r"2P, -23°26'26") with an
error radius of 0.°03 at 68 % confidence level.

The extension of Source S is also investigated using the same procedure as above. We did not
find significant extension. An upper limit on the size of the gamma-ray emission was obtained by
investigating the decrease of the likelihood with increasing source size. Under the assumption of a

^The model can be downloaded from|http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.litniil

disk shape, the upper limit amounts to 16' at the 68% confidence level. We treat Source S as a point
source hereafter. The best-fit location of Source S was estimated to be (R.A., Dec.) = (18''00'"59\
-24°11'31") using gamma rays above 2.15 GeV. This location is 13' away from the location of
IFGL J1800.5-2359c using gamma rays above 0.1 GeV and assuming IFGL J1801.3-2322c is
a point source. In order to investigate the origin of the position difference, we performed local-
ization with the same energy cut as the catalog analysis assuming a point source hypothesis and
the extended source hypothesis. We found that the localization with a point source hypothesis is
consistent with the catalog position while the extended source hypothesis yields the location con-
sistent with our localization with £ > 2.15 GeV cut. From this comparison, we conclude that the
localization difference is mostly due to the extension of Source N. We use the average of two lo-
calizations, (R.A., Dec.) = (18''00'"47'*, -24°05'45"), for the following analyses, with a systematic
uncertainty of half of the positional difference, ~ 6'.

We find no strong gamma-ray emission in the LAT energy band at the locations of HESS J 1 800-240 A
and C. We evaluated the likelihood of point sources at these locations and found likelihood ratios
-21n(Lo/Lps) ^ 10 (corresponding to a ~ 3 a) for HESS J1800-240A and -21n(Lo/Lps) ^ 9 (cor-
responding to a ~ 3 cr) for HESS J1800-240C, which confirms no significant emission at these
locations.

3.2. Energy Spectrum

We used the maximum likelihood fit tool, gtlike, for the spectral analysis of LAT sources
in this region. Flux points were obtained by performing the maximum likelihood analysis in each
energy bin. We used eight logarithmically spaced energy bins from 215 MeV to 100 GeV for
Source N. We did not divide the energy range from 0.2 GeV-100 GeV into bins for Source S due to
poor statistics. Fig.[3]shows resulting SEDs for (a) Source N and (b) Source S. The 68% confidence
region is illustrated for Source S up to ~ 20 GeV. We do not have sufficient number of events above
~ 20 GeV to constrain the spectrum. Upper limits at 90% confidence level are calculated assuming
a photon index of 2 if the detection is not significant in an energy bin, i.e., the likelihood ratio with
respect to no source is less than 9. Note that the value of the spectral index has negligible effect
on the upper limits. We take into account systematic errors due to uncertainties of the extension,
the Galactic diffuse model, the LAT effective area, and the effect of the nearby gamma-ray pulsar
PSR J 1809-2332. Systematic errors associated with the extension are estimated by varying the
source size by =b 1 a. We also evaluated the effect of the shape by comparing the disk shape
used in this paper and a circular two-dimensional Gaussian. Systematic errors due to the Galactic
diffuse model are estimated by using the residual gamma-ray data with respect to the best fit model
in the region where no LAT source is present, specifically, / = 2.°3-4.°7 andZ7 = -l.°25-0.°75. The

-10-

observed residual is energy dependent and can be modeled as ~ (116 {El\ GeV)^"*^^'" - 100) %
of the total Galactic diffuse flux. The normalization of the Galactic diffuse model is adjusted
according to the above equation to estimate the systematic error on the source flux. We evaluated
the position dependence of this residual and found a dispersion of ~5%, which is not a large effect.
Uncertainties of the Galactic diffuse model make a dominant contribution to systematic errors in
all the energy bins. Uncertainties of the LAT effective area is 10% at 100 MeV , decreasing to 5%

at 500 MeV, and increasing to 20% at 10 GeV and above (IRando et al. Il2009|) . The effect of the

nearby gamma-ray pulsar was estimated by varying the flux by ± 1 a.

We evaluated a possibility of a spectral break in the LAT energy band by comparing the
likelihood of the spectral fit for the LAT data between a simple power law and a broken power law
as a spectral model of Source N and Source S. The fits yield the likelihood ratio -21n(LpL/LBPL) ~
62 for Source N, where LpL and Lbpl are the likelihoods for the simple power-law model and
the broken power-law model, respectively. The likelihood ratio slightly decreases to at least ~
48 (corresponding to 6.6 a with two degrees of freedom) in the worst case accounting for 1 a
systematic uncertainties. Thus we conclude that Source N has a spectral break at 1.0 ± 0.2 GeV,
where the error is dominated by statistics. Photon indices are 2.09 ± 0.08 (stat) ± 0.28 (sys) below
the break and 2.74 ± 0.06 (stat) ± 0.09 (sys) above the break. Note t hat the photon index a bove

the break is consistent with the H.E.S.S. measurement of 2.66 ±0.27 (|Aharonian et al.ll2008[) . On
the other hand, we do not find any evidence of a spectral break for Source S. The photon index is
found to be 2.19 ± 0.14 (stat) ± 0.41 (sys).

We placed upper limits on the gamma-ray flux in the LAT band at the positions of HESS J 1 800-240A
and HESS J1800-240C. The upper limits are compared with the H.E.S.S. spectra in Fig.|4](a) and
(b). The gamma-ray upper limits in the GeV band at the location of HESS J1800-240A appear
to be inconsistent if the source spectrum is a simple power law. We evaluated the possibility of
a break between the LAT and H.E.S.S. energy band by comparing the likelihood of the spectral
fit for the LAT data between a broken power law and a simple power law as a spectral model of
this source. We fixed the spectral index and the flux values using the H.E.S.S. measurements for
the simple power-law model while we fixed the spectral index above the break and the flux at the
same value as the simple power law and varied the break energy and the spectral index below the
break for the broken power-law model. We used the parameter values at 1 a away from the best
fit values from H.E.S.S. in the direction where we expect less flux in the LAT energy band, i.e.
harder spectral index and lower flux than the best fit values. The fits yield the likelihood ratio
-21n(LpL/LBPL) ~ 35 (corresponding to a 5.5 a significance with two degrees of freedom) for a
break energy between 10 GeV and 100 GeV. The likelihood ratio decreases to ~ 8 correspond-
ing to a ~ 2.4 cr significance if we use parameter values 2 a away from the best-fit values from
H.E.S.S. We conclude that HESS J1800-240A is likely to have a spectral break between 10 GeV
and 100 GeV.

-11-

4. DISCUSSION

4.1. IFGL J1801.3-2322C (Source N)

IFGL J1801.3-2322c (Source N) is found to be extended and positionally coincident with
HESS J1801-233. Given the clear spatial match between the TeV source and the molecular cloud
interacting with SNR W28 and a smoothly connected spectrum in the GeV-TeV band, we assume
that the bulk of GeV and TeV gamma-ray emission comes from the molecular cloud illuminated by
particles accelerated in supernova shocks. Below, we adopt the simplest assumption, that GeV and
TeV gamma rays are emitted by a population of accelerated protons and electrons distributed in
the same region characterized by constant density and magnetic field strength. It should be noted
that the imaging resolution of the LAT is not good enough to differentiate the GeV emission region
from the TeV region. We assume the injected electrons have the same momentum distribution
as the protons. This assumption implies a break in the particle momentum spectrum because
the spectral index of the radio data, corresponding to lower particle momenta, is much harder
than that for the gamma-rays, which correspond to higher particle momenta. Therefore, we use a
broken power-law to describe the particle momentum distribution in the region. Electrons suffer
energy losses due to ionization (or Coulomb scattering), bremsstrahlung, synchrotron processes,
and inverse Compton scattering. The modification of t he electron spectral distribution due to such

losses was calculated according to lAtoyan et al.l (| 19951) . where electrons are assumed to be injected

at t = from an impulsive source. Since diffusive shock acceleration theory generally predicts
particle accelerations in the Sedov phase with a typical duration 10^-10"* yrs, the assumption of an
impulsive source would be a good approximation for SNRs with the age of lO'^-lO^ yrs. We adopt
4 X 10^ yrs for the age of W28 in this modeling. Note that here we consider the minimum momenta
of protons and electrons to be 100 MeV c"^ since the details of the proton/electron injection process
are poorly known. The gamma-ray spectrum from tt" decay produced by the interaction of protons
with ambient hydrogen is scaled by a facto r of 1.84 to account for helium and heavy nuclei in

target material and cosmic-ray composition (|Morill2009|) .

First, we consider a rr^-decay model to account for the broadband gamma-ray spectrum. In-
deed, the hadronic scenario gives the most satisfactory explanation f or the GeV gamrn a rays ob-

served in other SN Rs interacting with molecular gas such as W51C (lAbdo et al.ll2009|) and W44

(|Abdo et al.ll2010a|) . The number index of protons in the high-energy regime is constrained to be
Oh ~ 2.7 from the gamma-ray spectral slope. The observed gamma-ray luminosity requires the
gas density to be much larger than hh ~ 1 cm"^ averaged over the entire SNR shell in order not
to exceed the typical kinetic energy of a supernova explosion (~ 10^' erg). The resulting total
proton energy, Wp ~ 1.3 x 10"^^ • (10^ cnT^ /nu) ■ {d/2 kpc)^ erg, is less than 10% of the typical
kinetic energy of supernova explosions and quite reasonable. Note that Wp is not the total energy
of accelerated protons but that of the in-situ protons in the molecular clouds. Using the parame-

12

ters summarized in Table [H we calculated radiation model curves as shown in Figure [5] (a). As
described above, the spectral index of the proton momentum below the break is determined to be
Ol ~ 1.7 by modeling the radio spectrum as synchrotron radiation by relativistic electrons because
we assume that protons and electrons have identical injection spectra. It is difficult to derive the
break point of the proton momentum spectrum from the break of the gamma-ray spectrum since it
lies in the region where we expect a gamma-ray spectral curvature due to kinematics of tt" produc-
tion and decays. Because of this, the gamma-ray spectrum gives an upper bound for the momentum
break at ~ 5 GeV c"^ . The momentum break cannot be lower than ~ 2 GeV c"^ to avoid conflict
with the radio data. Here we adopt 2 GeV c~^ . The magnetic field strength is constrained to be
B ~ 160 fiG, for an electron-to-proton ratio of K^p = 0.01, which is the ratio found in the local
cosmic-ray abundance. Here the ratio is defined at a particle momentum of 1 GeV c"^ . This choice
of Kep results in minor contribution from the electron bremsstrahlung emission in the gamma-ray
band.

On the other hand, it is difficult to model the GeV-TeV spectrum by the electron bremsstrahlung
component alone because the break in the electron spectrum corresponding to the gamma-ray spec-
trum will appear in the radio data as sh own in Fig. [5](b) although it might be contaminated by a
thermal component (IDubner et al.l 120001) . Moreover, the magnetic field strength is constrained to
be less than 5^4 /iG for the age of 4 x lO'* yrs since no apparent synchrotron cooling is observed
at TeV energies. The low magnetic field requires the gas density to be smaller than hh ~ 5 cm"^
from the flux ratio of the synchrotron component in the radio band and the bremsstrahlung in the
gamma-ray band. These low values for the magnetic field and the gas density are in disagreement
with our assumption that the gamma-ray emission comes from the dense molecular cloud.

The gamma-ray spectrum is formally reproduced by inverse Compton (IC) emission from
accelerated electrons in Figl5](c). The interstellar radiation field for the IC scattering (see Table [B
is comprised of infrared, optical and the cosmic microwave background (CMB). The infrared and
optical componen ts are taken from the interstellar radiation field at the location of W28 in the
GALPROP code ([Porter et al.l 120081) . Since the spectral shape of the non-CMB component is
very complex, it is approximated by two infrared and two optical blackbody components. The
flux ratio between the IC and the synchrotron components constrains the magnetic field to be less
than 0.6/iG. The total energy in electrons is calculated to be W^ ~ 9.0 x 10^° ■ (d/2 kpc)^ erg for
an energy density of ~ 1.8 eV cm"^ for the interstellar radiation field. This model requires a
rather low gas density of hh ~ 2 x 10"^ cm"^ to suppress the electron bremsstrahlung, which is in
contradiction with our assumption.

Observations with Fermi have demonstrated that bright, extended gamma-ray sources co

incident with middle-aged SNRs interacting w ith molecular clou ds, such as W44 (|Abdo et al

2010ah . W51C (lAbdoetal.ll2009h . and IC 443 (lAbdo et all 1201 Obi) , exhibit spectral breaks from

13

9 GeV c"^ to a few tens of GeV c"' in the proton momentum spectrum. Our observations of
IFGL J1801.3-2322c in the vicinity of W28 combined with the radio data constrain the proton
momentum break to be in the range, 2-5 GeV c"^ . The observed energy distribution of relativistic
particles (most likely protons) could be influenced greatly by diffusive transport of particles. If
so, the relationship between the observed particle spectrum and the acceleration spectrum should
be rather complicated. The steep particle spectrum, a ~ 2.7, deduced for IFGL J 1801.3-2322c
does not necessarily represent the acceleration index. iGabici & AharonianI (|2007h discussed the
time evolution of non-thermal emission from molecular clouds illuminated by cosmic rays from a
nearby SNR and predicted a steep gamma-ray spectrum for an old SNR due to energy-dependent
diffusion of cosmic rays. More detailed studies of the properties of interactions between SNRs
and molecular clouds are required for a comprehensive description of the differences in the break
momenta and photon indices above the break among the above SNRs.

4.2. IFGL J1800.5-2359C (Source S)

IFGL J1800.5-2359c (Source S) was found to have no significant extension and spatially
coincides with the TeV source HESS J1800-240B, which is associated with molecular clouds
that contain the ultra-compact H II region W28A2 (G5. 89-0. 39). The observed GeV-TeV spec-
trum for Source S can be formally described by 7r° dominated, bremsstrahlung dominated, and
IC dominated mode ls in a similar manners to Source N, as shown in Fig. [6l Free-free emission
(|G6mez et al.lll99l|) is responsible for the radio spectrum of W28A2, which provides only upper
limits in the SED. Thus v alues of ^ and hh are le ss constrained by pion and bremsstrahlung models.

We adopt hh = 10^ cm ^ (lAharonian et al.

2008|) for both models. From the positional coincidence.

W28A2 can be invoked as a possible source of high-energy particles resp onsible for the gamma-ray

emission. W28A2 exhibits very ener getic bipolar molecular outflows (IHarvey & Forceillelll988

Acord et al. 1997: SoUins et al.

2004|) which would arise from the accretion of matter by a stellar

progenitor. iKlaassen et al.l (|2006|) estimated the total kinetic energy of the outflow of W28A2 to
be 3.5 X 10''^^ erg, which requires very high matter density to account for the observed gamma-ray
flux in both pion and bremsstrahlu ng models. If the matter density is as high as ~ 10^ cm"^ as
suggested by iKlaassen et al.l (|2006() . the resulting total energy would be less than a few percent of
the total kinematic energy of the outflow of W28A2.

Cosmic rays that escaped from the SNR W28 in earlier epochs would be another possi-
ble explanation given the energetics. Note that Source S would lie at a projected distance of
~ 10 • (J/2kpc) pc from the southern circular boundary of W28. If we assume that cosmic rays
are uniformly radiated by the SNR, the cosmic-ray flux bombarding the molecular cloud from the
SNR can be scaled by the fraction of the solid angle of the molecular cloud, ~ 8 x 10"^. Using this

14

scaling factor, the total energies in the whole SNR W28 could be calculated from the total energy
observed in this source to be ~ 2 x 10"^^ erg for protons in the tt" model. This is a reasonable
value for an SNR. The energy density of such cosmic rays is enhanced to be at least ~ 2 eV cm"^
at the molecular clouds under the assumption that all of the associated clouds (M ~ 4 x 10"^ Mq
at 2 kpc) are interacting with energetic particles. The molecular clouds can accumulate those
cosmic rays in the past since the diffusion coefficient is expected to be low in dense environ-

ments ( Aharonian & Atovanlll996 ).

While we detected LAT sources spatially coincident with HESS J 1801-233 and HESS J1800-230B,
no strong LAT counterpart can be found at the positions of HESS J1800-230A and HESS J1800-230C.
We found evidence of a spectral break between 10 GeV and 100 GeV for HESS J1800-230A. The
velocity with respect to the local standard of rest of the molecular emission po sitionally coinci-

dent w ith this source peaks at ~16 km s \ corresponding to a distance of ~4 kpc (|Aharonian et al

20081) . However, the velocity difference between this cloud and the one associated with HESS J 1 800-230B
is only ~ 10 km s"^ and can be attributed to the proper motion of the cloud complex that contains
both clouds. If the distance to the cloud associated with HESS J1800-230A is 4 kpc, then it is
unlikely that this source is associated with W28 since the distance to W28 is ~2 kpc. H U regions,
G6. 1-0.6 and/or 06. 225-0. 569 could be a source of accelerated particles. If the cloud associ-
ated with HESS J1800-230A is part of the cloud complex that includes the cloud associated with
HESS J1800-230B, then the distance would be 2 kpc, which is very similar to the clouds associ-
ated with W28. In this case, particles escaped from W28 may be responsible for the gamma-ray
emissions at both HESS J1800-230A and B and further studies of the origin of the difference in
the spectral shapes in these sources may provide good constraints on the particle diffusion process.

5. CONCLUSIONS

We have investigated two LAT sources in the W28 field. Source N (IFGL J1801.3-2322c)
which is located at the northeast boundary of the SNR W28 is positionally coincident with shocked
molecular clouds, and is spatially extended. The spectrum has a break around 1.0 GeV and
smoothly connects to the TeV spectrum, suggesting a physical relationship. Decay of 7r°s pro-
duced by the interaction of an SNR with molecular clouds naturally explains the gamma rays from
Source N based on the spatial correlation between GeV gamma rays and molecular clouds and
the energetics of cosmic rays. Electron bremsstrahlung can not be ruled out completely although
it requires a low density and low magnetic field in contradiction with the association with the
molecular clouds. W28 is the most plausible energy source due to the observational evidence of
interaction with the molecular clouds. The soft spectrum of the gamma rays may be explained
by the time evolution of non-thermal emission from molecular clouds illuminated by cosmic rays

15

from a nearby SNR due to energy-dependent diffusion of cosmic rays.

The Source S (IFGL J1800.5-2359c) was found to have no significant extension and spatially
coincides with the TeV source HESS J1800-240B, molecular clouds, and the ultracompact H II
region W28A2 (G5.89-0.39). The compact H II region W28A2 and the SNR W28 are possible
energy sources, but the W28A2 hypothesis requires extremely dense gas.

While no significant LAT counterpart is found at the positions of HESS J1800-230A and
HESS J1800-230C, the LAT upper limits for HESS J1800-230A coupled with the H.E.S.S. data
points imply a spectral break between 10 GeV and 100 GeV.

The Fermi LAT Collaboration acknowledges generous ongoing support from a number of
agencies and institutes that have supported both the development and the operation of the LAT as
well as scientific data analysis. These include the National Aeronautics and Space Administration
and the Department of Energy in the United States, the Commissariat a I'Energie Atomique and
the Centre National de la Recherche Scientifique / Institut National de Physique Nucleaire et de
Physique des Particules in France, the Agenzia Spaziale Italiana and the Istituto Nazionale di Fisica
Nucleare in Italy, the Ministry of Education, Culture, Sports, Science and Technology (MEXT),
High Energy Accelerator Research Organization (KEK) and Japan Aerospace Exploration Agency
(JAXA) in Japan, and the K. A. Wallenberg Foundation, the Swedish Research Council and the
Swedish National Space Board in Sweden.

Additional support for science analysis during the operations phase is gratefully acknowl-
edged from the Istituto Nazionale di Astrofisica in Italy and the Centre National d'Etudes Spatiales
in France.

REFERENCES

Abdo, A. A., et al. (The Fermi LAT Collaboration) 2009, ApJ, 706, LI

Abdo, A. A., et al. (The Fermi LAT Collaboration) 2010a, Science in press

Abdo, A. A., et al. (The Fermi LAT Collaboration) 2010b, to appear in ApJ, a rXiv: 1002.21981

Abdo, A. A., et al. (The Fermi LAT Collaboration) 2010c, to appear in ApJS. IarXiv:0910. 16081

Abdo, A. A., et al. (The Fermi LAT Collaboration) 2010d, submitted to Ap&SS. EXiv: 1002.22801

Acciari, V. A., et al. 2009, ApJ, 698, LI 33

Acord, J. M., Walmsley, C. M., & Churchwell, E. 1997, ApJ, 475, 693

-16-

Aharonian, F. A., Atoyan, A. M. 1996, A&A, 309, 917

Aharonian, F. A., Drury, L. O., & Volk, H. J. 1994, A&A, 285, 645

Aharonian, F, et al. 2008, A&A, 481, 401

Albert J., et al. (The MAGIC Collaboration) 2007, ApJ, 664, L87

Arikawa, Y, Tatematsu, K., Sekimoto, Y., & Takahashi, T. 1999, PAST, 51, L7

Atoyan A. M., Aharonian A. F, & Volk, H. J. 1995, Phys. Rev. D, 52, 3265

Atwood, W. B., et al. (The Fermi LAT Collaboration) 2009, ApJ, 697, 1071

Blandford, R. D. & Eichler, D. 1987, Phys. Rep., 154, 1

Brogan C. L., Gelfand, J. D., Gaensler, B. M., Kassim, N. E., & Lazio, T. J. W. 2006, ApJ, 639,
L25

Buckley, J. H., et al. 1998, A&A, 329, 639

Claussen, M. J., Frail, D. A., Goss, W. M., & Gaume, R. A. 1997, ApJ, 489, 143

Claussen, M. J., Goss, W. M., Frail, D. A., & Desai, K. 1999, ApJ, 522, 349

Dubner, G. M., Velazquez, P F, & Goss, W. M., 2000, ApJ, 120, 1933

Frail, D. A., Goss, W. M., & Slysh, V. I. 1994, ApJ, 424, LI 11

Gabici, S., & Aharonian, F A. 2007, ApJ, 665, L131

Gomez, Y , Rodriguez, L. F, Garay, G., & Moran, J. M. 1991, ApJ, 377, 519

Goudis, C. 1976, Ap&SS, 40, 91

Handa, T., Sofue, Y, Nakai, N., Hirabayashi, H., & Inoue, M. 1987, PASJ, 39, 709

Hartman, R. C, et al. 1999, ApJS, 123, 79

Harvey, P M., & Forveille, T. 1988, A&A, 197, L19

Lozinskaya, T. A. 1974, Soviet Astronomy, 17, 603

Lozinskaya, T. A., Supemovae and Stellar Wind in the Interstellar Medium (New York: AIP)

Mattox, J. R., et al. , 1996, ApJ, 461, 396

17

Mizuno, A., & Fukui, Y. 2004, Milky Way Surveys: The Structure and Evolution of our Galaxy,
317,59

Mori, M., 2009, Astropart. Phys., 31, 341

Kaspi, V. M., Lyne, A. G., Manchester, R. N., Johnston, S., D'Amico, N., & Shemar, S. L. 1993,
ApJ, 409, L57

Klaassen, P. D., Plume, R., Ouyed, R., von Benda-Beckmann, A. M., & Di Francesco, J. 2006,
ApJ, 648, L1079

Kovalenko, A. V., Pynzar', A. V., & Udal'Tsov, V. A. 1994, Astronomy Reports, 38, 95

Kuchar, T. A., & Clark, F. O. 1997, ApJ, 488, 224

Lockman, F J. 1989, ApJS, 71, 469

Pittori, C, et al. 2009, A&A, 506, 1563

Porter, T., et al. 2008, ApJ, 682, 400

Rando, R. et al..' arXiv:0907.0626l

Reach, W. T., Rho, J., & Jarrett, T. H. 2005, ApJ, 618, 297

Rho, J., & Borkowski, K. J. 2002, ApJ, 575, 201

SoUins, P K., et al. 2004, ApJ, 616, L35

Takeuchi, T. et al. 2010, PAS J, in press

Velazquez, P F, Dubner, G. M., Goss, W. M., & Green, A. J. 2002, AJ, 124, 2145

Wootten, A. 1981, ApJ, 245, 105

This preprint was prepared with the AAS LMj5C macros v5.2.

18

CM

(0

c
o
O

,. X ....... ■

--20:00:00 \

V ^^^^Hl^

PSRJ1 809-2332

\ W28

--24:00:00

\*«1 -

ttAVfei -

MfwV^PS* .

7V%i

-18:20:00

18:00:00 ^^^^ -

1 ^T^ 1

Right Ascension (J2000)

0.6 0.8 1

[counts/pixel]

Fig. 1 . — Fermi LAT 2-10 GeV counts map around the supernova remnant (SNR) W28. The count
map is smoothed by a Gaussian kernel of a = 0.°2, with the pixel size of 0.°025. The white line
from top left to bottom right indicates the Galactic plane.

19

Right Ascension (J2000)

0.5 0.6 0.7

[counts/pixel]

Righit Ascension (J2000)

0.1

0.2

0.3

0.4

0.5 0.6 0.7

[counts/pixel]

Riglit Ascension (J2000)

0.1

0.2

0.3

0.5 0.8 0.7

[counts/pixel]

Fig. 2. — Close-up views of the LAT 2-10 GeV count map around W28. The counts map is
smoothed by a Gaussian kernel of cr = 0.°2, with the pixel size of 0.°025. The inset of each figure
shows the effective LAT PSF for a photon spectral index of 2.5. A green circle in the north of
each figure indicates the best-fit disk size for Source N. A green cross indicate s the position of
Sourc e S. White diam onds indicate H II regions; W28A2 (see text), 06. 1-0.6 (IKuchar & Clark
1997k G6.225-0.569 (lLockmanlll989l) . The diamond on the right is W28A2. Black contours in
(a) sh ow the H.E.S.S . significance map for TeV gamma rays at 20, 40, 60 and 80% of the peak
value (JAharonian et a l. 2008). Bright TeV spots in the south are HESS J 1800-240 A, B and C as
indicated in the figure. Black contours in (b) give CO (J=l-0) line intensity taken by NANTEN at
25, 50, 75 % levels, for the velocity range f rom km s~^ to 20 km s~\ corresponding to k inematic

distances of approximately to ~ 4 kpc (JMizuno & Fukuil 120041: iTakeuchi et al.

2010). Black

contours in (c) indicate the VLA 90cm image at 25, 50, 75 % of the peak intensity (JBrogan Il2006|)
Outer boundaries of SNRs, as determined by the radio images, are drawn as white dashed circles
A white plus sign shows the position of PSR J 180 1-23.

20

lO-" 10

Energy (GeV)

10^ 10

Energy (GeV)

Fig. 3. — (a) Spectral energy distribution of the Fermi LAT source on the northeast boundary of
the supernova remnant W28. The red squares in the GeV regime are the LAT data. Horizontal
bars for data points indicate the energy range used in the fit. Vertical bars show statistical errors
in red and systematic errors in black. Upper limits are obtained at the 90% confidence level in
energy bins where the likelihood ratio va lue is < 9. The black circles represent data points for
HESS J1801-233 jAharonian et allboosl) . (b) Spectral energy distribution of the LAT source to
the south of W28. The red region is the 68% co nfidence range of the LAT spectrum. The black
circles show data points for HESS J1800-240B (lAharonian et al.ll2008k

lO-" 10

Energy (GeV)

10^ 10

Energy (GeV)

Fig. 4. — Upper limits obtained from the LAT data at the 90% confidence level in each energy
bin at the positions of (a) HESS J1800-240A and (b) HESS J1800-240C, on the assumption of
a phot on index of 2 for the power-law function. Black circles are H.E.S.S. spectral measure-
ments ( Aharonian et al.ll2008b .

21

(0

CM

10^ =

E

o

1

>

o

LU

10^

■o

u.

■o

CM

LU

1

10^

1 =-

E (b) Brems dominated

1 r

T ■— I ■— I ' — I ■— I ■— I ■— I ■— I ■— I ■— I ' — I ■— I ■— I ' — I ■— I ■— I ■— r

(a) Pion-decay dominated

H H H H H H H H <■

(c) IC dominated

_i i_i i_i i_i i_i i_i i_i i_i i_z.

H — H — H — H — H — H — H — ' I ^ i I — H — H — H — <H — H — 'I 'M — H — H — ^1 — H*— H — '

10

-1

10"^

10^

Energy (eV)

Fig. 5. — Multi-band spectra of the Fermi LAT source on the northeast boundary of the supernova
remnant W28. The red squares in the GeV regime are the LAT data, where t he red and black

errors on the flux ar e statistical and systematic, respectively. The radio emission (iKovalenko et al

19941 : iDubner et al.l 12000 ) is modeled by synchrotron radiation, while the gamma-ray emission
is modeled by different combinations of 7r°-decay (long-dashed curve), bremsstrahlung (dashed
curve), and inverse Compton (IC) scattering (dotted curve). Details of the models are described in
the text.

22-

(0

CM
I

E
o

>

LU

■o

CM

LU

10 =

1 ^

10

1

= (b) Brems dominated

10

1

-i 1— I 1— I 1 — I 1— I 1— I 1— I 1— I 1— I 1— I 1 — I 1— I 1— I 1 — I 1— I 1— I 1— r-

(a) Pion-decay dominated

A

H H-

-H H H H H H H H H ' I "< I H ^

H — H — H — H — H — H — H — H — H — H-

(c) IC dominated

J i_it i_i i_i i_i i_i i_i ■ I ^ ■ I i_i i_i L

J I 1

10^^ 10^^

Energy (eV)

Fig. 6. — Multi-band spectra of the Fermi LAT source to the south of the supernova remnant W28.
The red squares in the GeV regime are the LAT data, where the red and black errors on the flux

are st atistical and systematic, respectively. The radio emission (JHanda et al.lll987l : iGomez et al

19911) is modeled by synchrotron radiation, while the gamma-ray emission is modeled by differ-
ent combinations of 7r°-decay (long-dashed curve), brems strahlung (dashed curve), and inverse
Compton (IC) scattering (dotted curve). The sum of the three components is shown as a solid
curve. Details of the models are described in the text.

23

Table 1: Model parameters for the Fermi LAT sources around W28.

Model

Kep^

aL^

Pb^

an'

B

nn'

Wp'

wj

(GeV c-i)

(fiG)

(cm-3)

(10^" erg)

(lO^'^erg)

Source N

(a) Pion

0.01

1.7

2

2.7

160

100

1.3

1.9 X 10-2

(b) Bremsstrahlung

1

1.7

1

2.7

4

5

1.9

4.9

(c) Inverse Compton^

1

1.7

5

3.6

0.6

0.02

54

90

Source S

(a) Pion

0.01

1.7

2

2.4

1.2 X 10^

10^

1.5 X 10-2

6.5 X 10-5

(b) Bremsstrahlung

1

1.7

1

2.4

4

10^

1.1 X 10-3

2.1 X 10-3

(c) Inverse Compton^

1

1.7

15

3.2

2

0.05

1.7

2.3

The ratio of electron and proton distribution functions at 1 GeV c ' .

''A momentum distribution of particles is assumed to be a broken power-law, where the indices and the break momentum

are identical for both accelerated protons and electrons. The q;l is a photon index below the momentum break.
"^The pb is a momentum break for particle distribution.

''The photon index for broken power-law functions above the momentum break.
'"Average hydrogen number density of ambient medium.

^The distance is assumed to be 2 kpc. The total energy is calculated for particles > 100 MeV c"'.
*Seed photons for inverse Compton scattering of electrons include the CMB, two infrared (Tjr = 29,4.9 x 10^ K,

t/iR = 0.29,5.3 X 10"^ eV cm"^, respectively), and two optical components (Topt = 3.6 x 10-^, 1.0 x 10"* K, f/opt =

0.37,0.13 eV cm"-', respectively) in the vicinity of W28, assuming a distance of 2 kpc.