Skip to main content

Full text of "Discovery of very high energy gamma rays from PKS 1424+240 and multiwavelength constraints on its redshift"

See other formats


Discovery of very high energy gamma rays from PKS 1424+240 
and multiwavelength constraints on its redshift 

VERITAS collaboration: V. A. Acciari" 1 , E. Aliu" 2 , T. Ai\en v3 , T. Aune" 4 , 

M. Bautista" 5 , M. Beilicke" 6 , W. Benbow" 1 , M. Bottcher" 7 , D. Boltuch" 2 , 
S. M. Bradbury 1 ' 8 , J. H. Buckley 1 ' 6 , V. Bugaev" 6 , K. Byrum" 9 , A. Cannon 1 ' 10 , 
A. Cesarini" 11 , Y. C. Chow" 3 , L. Ciupik" 12 , P. Cogan" 5 , W. Cui" 13 , C. Duke" 14 , 
A. Falcone" 15 , J. P. Finley" 13 , G. Finnegan" 16 , L. Fortson" 12 , A. Furniss" 4 '*, N. Galante" 1 , 
D. Gall" 13 , G. H. Gillanders" 11 , S. Godambe" 16 , J. Grube" 10 , R. Guenette" 5 , G. Gyuk" 12 , 
D. Hanna" 5 , J. Holder" 2 , C. M. Hui" 16 , T. B. Humensky" 17 , P. Kaaret" 18 , N. Karlsson" 12 , 
M. Kertzman" 19 , D. Kieda" 16 , A. Konopelko" 20 , H. Krawczynski" 6 , F. Krennrich" 21 , 
M. J. Lang" 11 , S. LeBohec" 16 , G. Maier" 5 , S. McArthur" 6 , A. McCann" 5 , M. McCutcheon" 5 
J. Millis" 13 '" 22 , P. Moriarty" 23 , T. Nagai" 21 , R. A. Ong" 3 , A. N. Otte" 4 '*, D. Pandel" 18 , 

J. S. Perkins" 1 , A. Pichel" 24 , M. Pohl" 21 , J. Quinn" 10 , K. Ragan" 5 , L. C. Reyes" 25 , 
P. T. Reynolds" 26 , E. Roache" 1 , H. J. Rose" 8 , M. Schroedter" 21 , G. H. Sembroski" 13 , 
G. Demet Senturk" 27 , A. W. Smith" 9 , D. Steele" 12 , S. P. Swordy" 17 , M. Theiling" 1 , 
S. Thibadeau" 6 , A. Varlotta" 13 , V. V. Vassiliev" 3 , S. Vincent" 16 , R. G. Wagner" 9 , 
S. P. Wakely" 17 , J. E. Ward" 10 , T. C. Weekes" 1 , A. Weinstein" 3 , T. Weisgarber" 17 , 
D. A. Williams" 4 , S. Wissel" 17 , M. Wood" 3 , B. Zitzer" 13 , 
Fermi LAT collaboration: A. A. Abdo 1,2 , M. Ackermann 3 , M. Ajello 3 , L. Baldini 4 , 
J. Ballet 5 , G. Barbiellini 6 ' 7 , D. Bastieri 8 - 9 , B. M. Baughman 10 , K. Bechtol 3 , R. Bellazzini 4 , 

B. Berenji 3 , R. D. Blandford 3 , E. D. Bloom 3 , E. Bonamente 11,12 , A. W. Borgland 3 , 
J. Bregeon 4 , A. Brez 4 , M. Brigida 13 ' 14 , P. Bruel 15 , T. H. Burnett 16 , G. A. Caliandro 13 ' 14 , 

R. A. Cameron 3 , P. A. Caraveo 17 , J. M. Casandjian 5 , E. Cavazzuti 18 , C. Cecchi 11,12 , 
6. Qelik 19 ' 20 ' 21 , A. Chekhtman 1 ' 22 , C. C. Cheung 19 , J. Chiang 3 '*, S. Ciprini 11,12 , R. Claus 3 , 
J. Cohen-Tanugi 23 , J. Conrad 24,25,26 , S. Cutini 18 , C. D. Dermer 1 , A. de Angelis 27 , 
F. de Palma 13,14 , E. do Couto e Silva 3 , P. S. Drell 3 , A. Drlica- Wagner 3 , R. Dubois 3 , 



- 2 - 



D. Dumora 28 - 29 , C. Farmer 23 , C. Favuzzi 13 ' 14 , S. J. Fegan 15 , W. B. Focke 3 , P. Fortin 15 , 
28 M. Frailis 27 , Y. Fukazawa 30 , P. Fusco 13 - 14 , F. Gargano 14 , D. Gasparrini 18 , N. Gehrels 19 - 31 , 

S. Germani 11 ' 12 , B. Giebels 15 , N. Giglietto 13 ' 14 , P. Giommi 18 , F. Giordano 13 - 14 , 
30 T. Glanzman 3 , G. Godfrey 3 , I. A. Grenier 5 , J. E. Grove 1 , L. Guillemot 28 ' 29 , S. Guiriec 32 , 
si Y. Hanabata 30 , E. Hays 19 , R. E. Hughes 10 , M. S. Jackson 24 ' 25 ' 33 , G. Johannesson 3 , 

32 A. S. Johnson 3 , W. N. Johnson 1 , T. Kamae 3 , H. Katagiri 30 , J. Kataoka 34 ' 35 , N. Kawai 34 ' 36 , 

33 M. Kerr 16 , J. Knodlseder 37 , M. L. Kocian 3 , M. Kuss 4 , J. Lande 3 , L. Latronico 4 , 

34 F. Longo 6 ' 7 , F. Loparco 13 ' 14 , B. Lott 28 - 29 , M. N. Lovellette 1 , P. Lubrano 11 ' 12 , 

G. M. Madejski 3 , A. Makeev 1 ' 22 , M. N. Mazziotta 14 , J. E. McEnery 19 , C. Meurer 24 ' 25 , 

36 P. F. Michelson 3 , W. Mitthumsiri 3 , T. Mizuno 30 , A. A. Moiseev 20 ' 31 , C. Monte 13 ' 14 , 

37 M. E. Monzani 3 , A. Morselli 38 , I. V. Moskalenko 3 , S. Murgia 3 , P. L. Nolan 3 , J. P. Norris 39 , 

38 E. Nuss 23 , T. Ohsugi 30 , N. Omodei 4 , E. Orlando 40 , J. F. Ormes 39 , D. Paneque 3 , 

39 D. Parent 28 ' 29 , V. Pelassa 23 , M. Pepe 11 - 12 , M. Pesce-Rollins 4 , F. Piron 23 , T. A. Porter 41 , 

40 S. Raino 13 - 14 , R. Rando 8 - 9 , M. Razzano 4 , A. Reimer 42 ' 3 , O. Reimer 42 - 3 , T. Reposeur 28 - 29 , 

A. Y. Rodriguez 43 , M. Roth 16 , F. Ryde 33 ' 25 , H. F.-W. Sadrozinski 41 , D. Sanchez 15 , 

42 A. Sander 10 , P. M. Saz Parkinson 41 , J. D. Scargle 44 , C. Sgro 4 , M. S. Shaw 3 , E. J. Siskind 45 , 

43 P. D. Smith 10 , G. Spandre 4 , P. Spinelli 13 ' 14 , M. S. Strickman 1 , D. J. Suson 46 , H. Tajima 3 , 

H. Takahashi 30 , T. Tanaka 3 , J. B. Thayer 3 , J. G. Thayer 3 , D. J. Thompson 19 , 
L. Tibaldo 8 ' 5 ' 9 , D. F. Torres 47 ' 43 , G. Tosti 11 ' 12 , A. Tramacere 3 ' 48 , Y. Uchiyama 49 ' 3 , 
46 T. L. Usher 3 , V. Vasileiou 19 ' 20 ' 21 , N. Vilchez 37 , V. Vitale 38 ' 50 , A. P. Waite 3 , P. Wang 3 , 

B. L. Winer 10 , K. S. Wood 1 , T. Ylinen 33 ' 51 ' 25 , M. Ziegler 41 , 
48 and S. D. Barber 01 , D. M. Terndrup o2 '° 3 



- 3- 



vl Fred Lawrence Whipple Observatory, Harvard-Smithsonian Center for Astrophysics, 
Amado, AZ 85645, USA 

v2 Department of Physics and Astronomy and the Bartol Research Institute, University of 
Delaware, Newark, DE 19716, USA 

v3 Department of Physics and Astronomy, University of California, Los Angeles, CA 90095, 
USA 

v4 Santa Cruz Institute for Particle Physics and Department of Physics, University of Cal- 
ifornia, Santa Cruz, CA 95064, USA 

v5 Physics Department, McGill University, Montreal, QC H3A 2T8, Canada 
v6 Department of Physics, Washington University, St. Louis, MO 63130, USA 
v7 Astrophysical Institute, Department of Physics and Astronomy, Ohio University, Athens, 
OH 45701 

v8 School of Physics and Astronomy, University of Leeds, Leeds, LS2 9JT, UK 
v9 Argonne National Laboratory, 9700 S. Cass Avenue, Argonne, IL 60439, USA 
vl0 School of Physics, University College Dublin, Belfield, Dublin 4, Ireland 
v11 School of Physics, National University of Ireland, Galway, Ireland 
v12 Astronomy Department, Adler Planetarium and Astronomy Museum, Chicago, IL 60605, 
USA 

vl3 Department of Physics, Purdue University, West Lafayette, IN 47907, USA 
vl4 Department of Physics, Grinnell College, Grinnell, IA 50112-1690, USA 
vl5 Department of Astronomy and Astrophysics, 525 Davey Lab, Pennsylvania State Uni- 
versity, University Park, PA 16802, USA 

vl6 Department of Physics and Astronomy, University of Utah, Salt Lake City, UT 84112, 
USA 

vl7 Enrico Fermi Institute, University of Chicago, Chicago, IL 60637, USA 
vl8 Department of Physics and Astronomy, University of Iowa, Van Allen Hall, Iowa City, 



-4- 



IA 52242, USA 

vl9 Department of Physics and Astronomy, DePauw University, Greencastle, IN 46135-0037, 
USA 

v20 Department of Physics, Pittsburg State University, 1701 South Broadway, Pittsburg, KS 
66762, USA 

v21 Department of Physics and Astronomy, Iowa State University, Ames, IA 50011, USA 
v22 now at Department of Physics, Anderson University, 1100 East 5th Street, Anderson, 
IN 46012 

v23 Department of Life and Physical Sciences, Galway-Mayo Institute of Technology, Dublin 
Road, Galway, Ireland 

v24 Instituto de Astronomia y Fisica del Espacio, Casilla de Correo 67 - Sucursal 28, 
(C1428ZAA) Ciudad Autnoma de Buenos Aires, Argentina 
v25 Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA 
v26 Department of Applied Physics and Instrumentation, Cork Institute of Technology, Bish- 
opstown, Cork, Ireland 
v27 Columbia Astrophysics Laboratory, Columbia University, New York, NY 10027, USA 
1 Space Science Division, Naval Research Laboratory, Washington, DC 20375, USA 
2 National Research Council Research Associate, National Academy of Sciences, Wash- 
ington, DC 20001, USA 

3 W. W. Hansen Experimental Physics Laboratory, Kavli Institute for Particle Astro- 
physics and Cosmology, Department of Physics and SLAC National Accelerator Laboratory, 
Stanford University, Stanford, CA 94305, USA 

4 Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, 1-56127 Pisa, Italy 
5 Laboratoire AIM, CEA-IRFU/CNRS/Universite Paris Diderot, Service d'Astrophysique, 
CEA Saclay, 91191 Gif sur Yvette, France 

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



- 5 - 



7 Dipartimento di Fisica, Universita di Trieste, 1-34127 Trieste, Italy 
8 Istituto Nazionale di Fisica Nucleare, Sezione di Padova, 1-35131 Padova, Italy 
9 Dipartimento di Fisica "G. Galilei", Universita di Padova, 1-35131 Padova, Italy 
10 Department of Physics, Center for Cosmology and Astro-Particle Physics, The Ohio 
State University, Columbus, OH 43210, USA 

11 Istituto Nazionale di Fisica Nucleare, Sezione di Perugia, 1-06123 Perugia, Italy 
12 Dipartimento di Fisica, Universita degli Studi di Perugia, 1-06123 Perugia, Italy 
13 Dipartimento di Fisica "M. Merlin" dell' Universita e del Politecnico di Bari, 1-70126 
Bari, Italy 

14 Istituto Nazionale di Fisica Nucleare, Sezione di Bari, 70126 Bari, Italy 
15 Laboratoire Leprince-Ringuet, Ecole polytechnique, CNRS/IN2P3, Palaiseau, France 
16 Department of Physics, University of Washington, Seattle, WA 98195-1560, USA 
17 INAF-Istituto di Astrofisica Spaziale e Fisica Cosmica, 1-20133 Milano, Italy 
18 Agenzia Spaziale Italiana (ASI) Science Data Center, 1-00044 Frascati (Roma), Italy 
19 NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 

20 Center for Research and Exploration in Space Science and Technology (CRESST), NASA 
Goddard Space Flight Center, Greenbelt, MD 20771, USA 

21 University of Maryland, Baltimore County, Baltimore, MD 21250, USA 
22 George Mason University, Fairfax, VA 22030, USA 

23 Laboratoire de Physique Theorique et Astroparticules, Universite Montpellier 2, 

CNRS/IN2P3, Montpellier, France 

24 Department of Physics, Stockholm University, AlbaNova, SE-106 91 Stockholm, Sweden 
25 The Oskar Klein Centre for Cosmoparticle Physics, AlbaNova, SE-106 91 Stockholm, 

Sweden 

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



- 6- 



Wallenberg Foundation 

27 Dipartimento di Fisica, Universita di Udine and Istituto Nazionale di Fisica Nucleare, 
Sezione di Trieste, Gruppo Collegato di Udine, 1-33100 Udine, Italy 

28 Universite de Bordeaux, Centre d'Etudes Nucleaires Bordeaux Gradignan, UMR 5797, 
Gradignan, 33175, France 

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

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

31 University of Maryland, College Park, MD 20742, USA 
32 University of Alabama in Huntsville, Huntsville, AL 35899, USA 

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

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

35 Waseda University, 1-104 Totsukamachi, Shinjuku-ku, Tokyo, 169-8050, Japan 
36 Cosmic Radiation Laboratory, Institute of Physical and Chemical Research (RIKEN), 

Wako, Saitama 351-0198, Japan 

37 Centre d'Etude Spatiale des Rayonnements, CNRS/UPS, BP 44346, F-30128 Toulouse 

Cedex 4, France 

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

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

40 Max-Planck Institut fur extraterrestrische Physik, 85748 Garching, Germany 

41 Santa Cruz Institute for Particle Physics, Department of Physics and Department of 



Received 



accepted 



Astronomy and Astrophysics, University of California at Santa Cruz, Santa Cruz, CA 95064, 
USA 

42 Institut fiir Astro- und Teilchenphysik and Institut fur Theoretische Physik, Leopold- 

Franzens-Universitat Innsbruck, A-6020 Innsbruck, Austria 

43 Institut de Ciencies de l'Espai (IEEC-CSIC), Campus UAB, 08193 Barcelona, Spain 
44 Space Sciences Division, NASA Ames Research Center, Moffett Field, CA 94035-1000, 

USA 

45 NYCB Real-Time Computing Inc., Lattingtown, NY 11560-1025, USA 
46 Department of Chemistry and Physics, Purdue University Calumet, Hammond, IN 
46323-2094, USA 

47 Instituci6 Catalana de Recerca i Estudis Avangats, Barcelona, Spain 
48 Consorzio Interuniversitario per la Fisica Spaziale (CIFS), 1-10133 Torino, Italy 
49 Institute of Space and Astronautical Science, JAXA, 3-1-1 Yoshinodai, Sagamihara, 

Kanagawa 229-8510, Japan 

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

51 School of Pure and Applied Natural Sciences, University of Kalmar, SE-391 82 Kalmar, 

Sweden 

ol Homer L. Dodge Department of Physics and Astronomy, The University of Oklahoma, 
440 W. Brooks St., Norman, OK 73019, USA 

o2 Department of Astronomy, The Ohio State University, 140 West 18th Avenue, Columbus, 
OH 43210, USA 

o3 National Science Foundation, 4201 Wilson Boulevard, Arlington, VA 22230, USA 
* Corresponding author, nepomuk.otte@gmail.com, amy.furniss@gmail.com, jchi- 
ang@ slac . st anf or d . edu 



- 8- 



so ABSTRACT 

51 We report the first detection of very-high-energy 1 (VHE) gamma-ray emission 
above 140 GeV from PKS 1424+240, a BL Lac object with an unknown redshift. 
The photon spectrum above 140 GeV measured by VERITAS is well described by 
a power law with a photon index of 3.8 ± 0.5 sta t ± 0.3 syst and a flux normalization 
at 200 GeV of (5.1 ± 0.9 stat ± 0.5 syst ) x 10~ n TeV^cm-V 1 , where stat and syst 
denote the statistical and systematical uncertainty, respectively. The VHE flux is 
steady over the observation period between MJD 54881 and 55003 (2009 Febru- 
ary 19 to June 21). Flux variability is also not observed in contemporaneous high 
energy observations with the Fermi Large Area Telescope (LAT). Contempora- 
neous X-ray and optical data were also obtained from the Swift XRT and MDM 
observatory, respectively. The broadband spectral energy distribution (SED) is 
well described by a one-zone synchrotron self-Compton (SSC) model favoring a 
redshift of less than 0.1. Using the photon index measured with Fermi in com- 
bination with recent extragalactic background light (EBL) absorption models it 
can be concluded from the VERITAS data that the redshift of PKS 1424+240 is 
less than 0.66. 

52 Subject headings: BL Lacertae objects: individual ( PKS 1424+240 = 

53 VER J1427+237); gamma rays: observations 



x 7-ray emission above 100 GeV 



- 9- 



54 1. Introduction 

55 PKS 1424+240 was detected as a radio source by Condon et al. (1977). It was classified 

56 as a blazar by Impey & Tapia (1988) from optical polarization studies. Fleming et al. 

57 (1993) verified the polarization results and also reported non-thermal X-ray radiation, 

58 further strengthening the classification. 

59 Blazar emission is dominated by non-thermal radiation, which is thought to be related 
eo to charged particle acceleration near a massive compact object in the center of the host 

ei galaxy, or in outflowing relativistic jets. The SED is characterized by two peaks. The lower 

62 peak is widely accepted to be synchrotron radiation from relativistic electrons and occurs 

63 in the IR to X-ray bands. The higher energy peak is in the gamma-ray band, sometimes at 

64 energies as high as a few TeV, and can be created via either inverse- Compton scattering by 
es relativistic electrons or hadronic interactions (for a review see Bottcher 2007, and references 
ee therein). The position of the synchrotron peak of PKS 1424+240 has not been measured, 
67 but it can be constrained from optical and X-ray data to be between 10 15 Hz and 10 17 Hz. 
es Depending on the definition used, PKS 1424+240 is either an intermediate-frequency-peaked 

69 BL Lac (IBL) (Nieppola et al. 2006) or a high-frequency-peaked BL Lac (HBL) (Padovani 

70 & Giommi 1996; Abdo et al. 2009a). 

71 Gamma-ray emission from PKS 1424+240 was not detected by EGRET (Fichtel et al. 

72 1994), but was recently observed with the Fermi LAT pair-conversion telescope (Abdo et 

73 al. 2009b,c). The reported flux above 100 MeV of (6.2 ± 0.8) x lO^cm" 2 s" 1 and hard 

74 spectral index T = 1.80 ± 0.07 (dN/dE oc E~ r ) triggered VERITAS observations. 

75 The redshift of PKS 1424+240 is not known. Scarpa & Falomo (1995) have derived 

76 a lower limit on the redshift of z > 0.06 and Sbarufatti et al. (2005) a limit of z > 0.67, 

77 both assuming a minimum luminosity of the host galaxy. The latter authors also reported 

78 evidence that the ratio of the nucleus to host luminosity is much larger than 100, which is 



-10- 



79 typical for BL Lac objects but complicates photometric determination of the redshift. 

so We report the detection of PKS 1424+240 in VHE gamma rays and contemporaneous 

si observations with Fermi, Swift, and the MDM observatory. Shortly after the VHE discovery 

82 (Ong 2009), it was confirmed by the MAGIC collaboration (Teshima 2009). This discovery 

83 marks the first Fermi-motivated VHE discovery. 

84 2. Observations and Analysis of VERITAS Data 

as The VERITAS observatory, located in southern Arizona at 1.3 km a.s.L, is described in 

8 6 detail in Weekes et al. (2002) and Holder et al. (2006). 

87 PKS 1424+240 was observed with VERITAS between 2009 February 19 and June 21 at 
ss zenith angles between 7° and 30°. The observations were performed in wobble mode (Fomin 

89 et al. 1994) with a 0.5° offset, enabling simultaneous background estimation. About one 

90 third of the data were taken during low levels of moonlight. About 65% of the observations 

91 were conducted with only three telescopes due to the relocation of one telescope, which 

92 began in May and was completed in August 2009. Of the 37.3 hours of data, 28.5 hours 

93 survive standard data quality selection. 

94 Events are reconstructed following the procedure in Acciari et al. (2008). The recorded 

95 shower images are parameterized by their principal moments, giving an efficient suppression 

96 of the far more abundant cosmic-ray background. Two separate sets of cuts are applied 

97 to reject background events, hereafter called soft and medium. These cuts are applied 

98 to the parameters mean scaled width (MSW), and mean scaled length (MSL), apparent 

99 altitude of the maximum Cherenkov emission (shower maximum), and 9 2 , the squared 

100 direction between the position of PKS 1424+240 and the reconstructed origin of the event. 

101 Studies on independent data sets show that a shower-maximum cut significantly improves 



- 11 - 



102 the low energy sensitivity. Soft cuts have a higher sensitivity for sources with soft photon 

103 spectra because of a lower energy threshold resulting from a minimum size cut of 50 

104 photoelectrons. In the medium cuts a minimum size cut of 100 photoelectrons is applied. 

105 Size is a measure of the recorded photoelectrons from a shower and a good indicator of 
we the energy of the primary. For the soft-cuts analysis the remaining cuts are MSW < 1.06, 
107 MSL < 1.30, shower maximum > 7km, and 9 2 < (0.14 ) 2 , and MSW < 1.04, MSL < 1.28, 
los shower maximum > 5 km, and 6 2 < (0.1°) 2 for the medium cuts. The cuts have been 

109 optimized a priori to yield the highest sensitivity for a source with 5% of the Crab Nebula 

no gamma-ray flux. The results are independently reproduced with two different analysis 

in packages explained in Cogan (2008) and Daniel (2008). 

n2 In the soft-cuts analysis, 1907 on-source events remain out of 1.25 x 10 7 triggered 

in events. The background calculated with the reflected-region method (Berge et al. 2007) is 

H4 1537 events, which leaves an excess of 370 events. Figure 1 shows the corresponding 6 2 

H5 distribution. The statistical significance of the observed excess is 8.5 standard deviations, 

lie a, calculated with Equation 17 of Li & Ma (1983), and including a trials factor of two 

ii7 for the two sets of cuts. In the medium-cuts analysis the post-trials significance is 4.8 a 

us (329 on-source events with an estimated background of 244). The angular distribution of 

n9 the excess events is consistent with a point source. The center of gravity of the excess is 

120 14 h 27 m , s ± 7* tat , 23° 47' 40" ± 2' stat coinciding with the position of PKS 1424+240 in radio 

121 (Fey et al. 2004). The VERITAS source name is VER J1427+237. 

122 Figure 2 shows the light curve of PKS 1424+240 in different energy bands for the time 

123 period overlapping the VERITAS observations. The flux measured by VERITAS above 

124 140 GeV is ~ 5% of the Crab Nebula flux. The VERITAS data from each dark period 2 are 

125 combined into a single bin to produce a light curve, which is consistent with a constant flux, 



2 The ~ 3 week observing period between full moons 



- 12 - 



126 x 2 =0-3 for 3 degrees of freedom (d.o.f.). However, even a doubling in flux would have been 

127 difficult to detect. There is no evidence for strong flaring episodes on shorter timescales. 

128 Figure 3 shows the differential photon spectra derived with the soft-cuts and 

129 medium-cuts analyses, with one overlapping flux point at 260 GeV. The fraction of events 

130 that are used both in the last bin in the soft-cuts analysis and in the second bin in the 

131 medium-cuts analysis is about 2%, small enough to allow a combined fit of the flux points 

132 from the two analyses, with the more significant soft-cuts result at 260 GeV used in the 

133 fit. The combined spectrum is well parameterized (x 2 =2.2 for 4 d.o.f.) by a power law 

134 dN/dE = F ■ (E/E )~ r , where the photon index T is 3.8 ± 0.5 sta t ± 0.3 syst and F is 

13 5 (5.1 ± 0.9 stat ± 0.5 syst ) x HT 11 TeV^cm-V 1 for E = 200 GeV. The combined spectrum is 
we consistent with the fit of the soft-cuts points alone, albeit with half the uncertainty on the 

137 photon index. 

138 3. Multiwavelength Observations 

139 Gamma-ray observations with Fermi-LAT (100 MeV to 300 GeV), X-ray and 

wo optical observations with Swift XRT (0.2-10 keV) and UVOT (170-650 nm), and optical 

141 observations in the R, V and 1 bands at the MDM observatory were obtained simultaneously 

142 or quasi-simultaneously with the VERITAS observations. 

143 The LAT pair-conversion telescope on board the Fermi Gamma-ray Space Telescope 

144 continuously monitors the entire sky between 100 MeV and several hundred GeV (Atwood 

145 et al. 2009). The LAT data overlapping with the VERITAS observations were analyzed by 
we selecting "diffuse" class events that have the highest probability of being photons. Further 
147 event selection was done by only accepting events that come within a 15° radius from 

us PKS 1424+240 and have energies between 0.1 and 300 GeV. Events with zenith angles above 



-13- 



149 105° were excluded to limit contamination by gamma rays coming from the Earth's albedo. 

150 The analysis of the photon spectrum and light curve were performed with the standard 

151 likelihood analysis tools available from HEASARC ScienceTools v9rl5p2. Accidental 

152 coincidences with charged cosmic rays in the detector were accounted for using instrument 

153 response functions P6_V3_DIFFUSE. The background model used to extract the gamma-ray 

154 signal includes a Galactic diffuse emission component and an isotropic component 3 . The 

155 isotropic component includes contributions from the extragalactic diffuse emission as well as 

156 from residual charged particle backgrounds. The spectral shape of the isotropic component 

157 was derived from residual high latitude events after the Galactic contribution had been 

158 modeled. The background model also takes into account unresolved gamma-ray sources in 

159 the region of interest, thus avoiding a bias in the spectral reconstruction. To further reduce 

160 systematic uncertainties in the analysis, the normalization and spectral parameters in the 
lei background model were allowed to vary freely during the spectral point fitting. 

162 The Fermi-LAT flux measurements are shown in the broadband SED in Figure 4. The 

163 flux values are unfolded by assuming an underlying power-law, giving an integrated flux 

164 over the 0.1-300 GeV band (7.04 ± 0.96 stat ± 0.38 sys ) x lO^cm" 2 s -1 , and a differential 

165 photon spectral index IYat = 1-73 ± 0.07 sta t ± 0.05 sys . The light curve of the integral flux 
lee above 100 MeV is plotted with 10-day bins in Figure 2. A fit with a constant yields a 

167 x 2 — H-5 for 11 d.o.f., suggesting no variability. 

168 Target of opportunity observations of nearly 16ksec, distributed over ten observing 

169 periods, were obtained with Swift (Gehrels et al. 2004) following the detection of VHE 

170 emission from PKS 1424+240. The data reduction and calibration of the XRT data 

171 were completed with the HEASoft v6.6.3 standard tools. The XRT data were taken in 



3 http:/ /fermi. gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels. html 



-14- 



172 photon-counting mode and contained modest pile up for nine of the observations, which 

173 was taken into account by masking a region with 3-6 pixel radius around the source. The 

174 outer radius chosen for the signal region was 20 pixels and a background region of similar 

175 size was chosen about 5 arcminutes off source. 

176 X-ray energy spectra could be extracted from all observing periods and are well 

177 described by an absorbed power law using the fixed Galactic column density of neutral 

178 hydrogen from Dickey & Lockman (1990) (iVh = 0.264 x 10 21 cm~ 2 ). The fit spectral index 

179 varies between 2.1 and 2.9 (photon index between 3.1 and 3.9) with a typical statistical 

180 uncertainty of 0.1, while the normalization changes between 1.40 x 10 -2 and 0.74 x 10~ 2 

181 photons keV _1 cm _2 s _1 at 1 keV with a typical uncertainty of 0.07 x 10~ 2 keV _1 cm -2 s _1 . 

182 For the modeling of the SED we use the average spectrum shown in Figure 4. The light 

183 curve shows that the X-ray flux is variable over the ten days of observation. A fit to a 

184 constant flux yields a x 2 of 60 for 9 d.o.f.. UVOT observations were taken in the six V, 

185 B U, Wl, M2 and W2 bands and were calibrated using standard techniques (Poole et al. 

186 2008). The reddening has been accounted for by interpolating the absorption values from 

187 Schlegel et al. (1998) with a galactic spectral extinction model (Fitzpatrick 1999) obtaining 

188 0.663, 0.968, 0.922 mag for the three UV bands Wl, M2, and W2 and an assumed redshift 

189 of z=0. The corresponding light curves are shown in Figure 2. 

190 Data in the optical bands were also obtained with the 1.3 m telescope and 4K 

191 imager of the MDM observatory located on the west ridge of Kitt Peak near Tucson, 

192 Arizona. The CCD was operated in unbinned mode, which produces an image scale of 

193 0.315 arcseconds/pixel. 2-4 images were obtained in the V, R, and I bands during each 

194 observation. Physical magnitudes were computed from differences in the instrumental 

195 magnitudes from the three standard stars in Fiorucci & Tosti (1996), assuming that the 
we magnitudes quoted in that paper are exact. The magnitudes were then corrected for 



- 15 - 



W7 Galactic extinction using extinction coefficients calculated following Schlegel et al. (1998), 

W8 taken from NED 4 , and were then converted into vF v fluxes. During the 14-day span of the 

199 optical photometry, the visual brightness increased by 14% and colors became slightly bluer. 

200 4. Redshift Upper Limit 

201 The observed gamma-ray spectrum above 100 GeV is affected by the absorption of 

202 gamma rays via pair conversion with EBL photons (Nikishov A. I. 1962; Gould & Schreder 

203 1967). Depending on the redshift, this effect can result in a significant softening of the 

204 spectrum. We estimate an upper limit of the redshift of PKS 1424+240 by assuming an 

205 intrinsic VHE spectrum and making use of the recent advances in extragalactic background 

206 light (EBL) modeling. 

207 We assume that the intrinsic spectrum above 140 GeV can be described by a power law. 

208 The hardest photon index that we consider is 1.7, which is the value from the simultaneous 

209 Fermi observations. The use of Fermi observations allows a model independent estimate 

210 of the hardest possible intrinsic spectrum (see also Abdo et al. 2009d). The power law 
2n with an index of 1.7 is absorbed using recent EBL models from Franceschini et al. (2008), 

212 Gilmore et al. (2009), and Finke et al. (2009). After absorption the shape of the spectrum 

213 is fit to the VERITAS spectrum with the normalization as a free parameter, and the best 

214 estimate of the redshift is determined by minimizing \ 2 ■ For an intrinsic index of 1.7 this 

215 best fit redshift is z= 0.5 ± 0.1 stat ± 0.1 syst with a x 2= 4 and 5 d.o.f. . The systematic 

216 uncertainty is estimated from the differences in the EBL models. 

217 Instead of assuming no break in the photon spectrum, a more likely scenario is that the 

218 intrinsic spectrum softens with increasing energy. In this case an index of 1.7 is an upper 



4 http: / / nedwww.ipac.caltech.edu 



-16- 



219 limit of the true photon index and the corresponding upper limit on the redshift is z < 0.66 

220 with a 95% confidence level. 

221 5. Spectral Modeling 

222 The spectral energy distribution, comprising data from all of the observations, is shown 

223 in Figure 4. We model the SED using an improved version of the leptonic one-zone jet 

224 model of Bottcher & Chiang (2002). These calculations include time-dependent particle 

225 injection and evolution, and they allow for quasi-equilibrium solutions in which a slowly 

226 varying broken power-law electron distribution arises from a single power-law injection 

227 function, dn^j/d'y oc r y~ q with a low- and high-energy cutoff 71 and 72, respectively. All 

228 model fits presented here are in the fast-cooling regime, with the cooling break at 71. We 

229 define the magnetic-field equipartition e# as eg = Lb/L £ with Lb the Poynting flux derived 

230 from the magnetic energy density and L e the energy flux of the electrons propagating 

231 along the jet. The corresponding partition fraction for an electron-proton plasma assuming 

232 L p = 10 x L e of cold protons would be one order of magnitude lower. For an in-depth 

233 description of this quasi-equilibrium jet model, see Acciari et al. (2009). 

234 There are few observational constraints on the model parameters for PKS 1424+240 

235 and the redshift is unknown. No superluminal motion has been resolved in this object, and 

236 it has not been monitored well enough to firmly establish a minimum variability timescale 

237 to constrain the size of the emitting region. The different sizes of the emission region R B 

238 assumed here are compatible with the X-ray variability timescale of about a day. We 

239 therefore consider a range of plausible redshifts and adopt model parameters which were 

240 typically adequate for modeling other VHE blazars. The redshifts investigated range from 

241 z = 0.05, similar to the redshift of the nearby HBLs Mrk421 and Mrk501, to z = 0.7. This 

242 covers the redshift range determined in the previous section and is just above the lower 



-17- 



243 limit set by Sbarufatti et al. (2005), z > 0.67. 

244 The shape of the high-energy part of the electron spectrum is well constrained by the 

245 rather steep slope of the X-ray spectrum, which has an average photon index Tx-ray ~ 3.7. 

246 In all fits, the relativistic electrons are injected into the emission region with a fixed q = 5.1. 

247 Lacking direct constraints on the viewing angle 9 ohs , it was chosen such that the Doppler 

248 factor D = (T[l — /3 r cos^obs])' 1 = I\ where T is the bulk Lorentz factor of the emitting 

249 material, and /3pc is the velocity. The model parameters that were varied are shown in Table 

250 1. Figure 4 shows the fits, after EBL absorption using the model of Gilmore et al. (2009). 

251 The SED modeling shows that a reasonable fit can in principle be achieved for any 

252 redshift in the considered range. However, the inset in Figure 4 illustrates that above 

253 z ~ 0.2, the model VHE gamma-ray spectrum becomes increasingly too steep compared 

254 with the observed VERITAS spectrum. Furthermore, for redshifts z > 0.4 the models 

255 require unreasonably large Doppler factors of D > 50. We note that in particular for the 

256 lowest redshift considered, z = 0.05, a good fit can be achieved with almost equipartition 

257 between magnetic-field and electron energy densities. 

258 An attempt to improve the fit in the gamma-ray bands, by including an external 

259 Compton component, results in a steeper VHE gamma-ray spectrum. This is in conflict 

260 with the VERITAS spectrum and a worse representation of the Fermi spectrum. We 

261 therefore conclude that a leptonic fit to the SED of PKS 1424+240 during the VERITAS 

262 observation is possible with a pure SSC model very close to equipartition, in particular if 

263 the redshift of the source is z < 0.1. 



264 



6. Summary 



265 We report the detection of PKS 1424+240 in VHE gamma-rays. The observation with 

266 VERITAS was motivated by the release of the first Fermi source lists (Abdo et al. 2009b, c) 

267 and this is the first time that Fermi observations have led to the discovery of a new source 

268 in the adjacent VHE band. 

269 The VHE spectrum of PKS 1424+240 has a photon index of 3.8±0.5 stat ±0.3 sys , whereas 

270 the spectrum in the Fermi energy range has a photon index of 1.73 ± 0.07 stat ± 0.05 sys , 

271 indicating a break in the spectrum at several tens of GeV. The break can be explained 

272 by a one-zone SSC model assuming a wide range of redshifts or could result from EBL 

273 absorption if the redshift is about 0.5 and the intrinsic photon index is 1.7, from which a 

274 redshift upper limit of 0.66 is inferred. The modeling favors a lower redshift but cannot 

275 exclude that PKS 1424+240 is among the most distant sources detected in the VHE regime. 

276 PKS 1424+240 is the third extragalactic source detected in the VHE regime with 

277 an unknown or uncertain redshift. It is evident that increased efforts are needed to 

278 determine the redshifts of VHE detected blazars. A redshift measurement will allow a 

279 better understanding of the source-intrinsic mechanisms and the absorption effects which 

280 go along with the gamma-ray propagation. 

281 VERITAS is supported by grants from the US Department of Energy, the US National 

282 Science Foundation, and the Smithsonian Institution, by NSERC in Canada, by Science 

283 Foundation Ireland, and by STFC in the UK. We acknowledge the excellent work of the 

284 technical support staff at the FLWO and the collaborating institutions in the construction 

285 and operation of the instrument. N. O. acknowledges the receipt of a Feodor Lynen 

286 fellowship of the Alexander von Humboldt Foundation. 

287 The Fermi LAT Collaboration acknowledges support from a number of agencies 



-19- 



288 and institutes for both development and the operation of the LAT as well as scientific 

289 data analysis. These include NASA and DOE in the United States, CEA/Irfu and 

290 IN2P3/CNRS in France, ASI and INFN in Italy, MEXT, KEK, and JAXA in Japan, and 

291 the K. A. Wallenberg Foundation, the Swedish Research Council and the National Space 

292 Board in Sweden. Additional support from INAF in Italy and CNES in France for science 

293 analysis during the operations phase is also gratefully acknowledged. 

294 This research has made use of the SIMBAD database, operated at CDS, Strasbourg, 

295 France. 

296 Facilities: VERITAS, Swift, Fermi. 



-20- 

REFERENCES 

298 Abdo, A., et al. 2009a, ApJ, in preparation 

299 Abdo, A. A., et al. 2009b, ApJS, 183, 46 
soo Abdo, A. A., et al. 2009c, ApJ, 700, 597 

301 Abdo, A., et al. 2009d, ApJ, submitted 

302 Acciari, V. A., et al. 2008, ApJ, 679, 1427 

303 Acciari V. A., et al. 2009, ApJ, 707, 612 

304 Atwood, W. B., et al. 2009, ApJ, 697, 1071 

30 5 Berge, D., Funk, S., k Hinton, J. 2007, A&A, 466, 1219 

306 Bdttcher, M., & Chiang, J. 2002, ApJ, 581, 127 

307 BSttcher, M. 2007, Ap&SS, 309, 95 

308 Cogan, P. 2008, International Cosmic Ray Conference, 3, 1385 

309 Condon, J. J., Hicks, P. D., & Jauncey, D. L. 1977, AJ, 82, 692 

310 Daniel, M. K. 2008, 30th International Cosmic Ray Conference, 3, 1325 

311 Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215 

312 Fey, A. L., et al. 2004, AJ, 127, 3587 

313 Fichtel, C. E., et al. 1994, ApJS, 94, 551 

314 Finke, J. D., Razzaque, S., & Dermer, C. D. 2009, arXiv:0905.1115 

31 5 Fiorucci, M., & Tosti, G. 1996, A&AS, 116, 403 



316 Fitzpatrick, E. L. 1999, PASP, 111, 63 

317 Fleming, T. A., Green, R. F., Jannuzi, B. T., Liebert, J., Smith, P. S., & Fink, H. 1993, 

318 AJ, 106, 1729 

3w Fomin, V. P., Stepanian, A. A., Lamb, R. C, Lewis, D. A., Punch, M., & Weekes, T. C. 

320 1994, Astroparticle Physics, 2, 137 

321 Franceschini, A., Rodighiero, G., & Vaccari, M. 2008, A&A, 487, 837 

322 Gehrels, N., et al. 2004, ApJ, 611, 1005 

323 Gilmore, R. C, Madau, P., Primack, J. R., Somerville, R. S., & Haardt, F. 2009, 

324 arXiv:0905.1144 

325 Gould, R. J., & Schreder, G. P. 1967, Physical Review , 155, 1404 

326 Hillas, A. M. 1985, International Cosmic Ray Conference, 3, 445 

327 Holder, J. et al. 2006, Astroparticle Physics, 25, 391 

328 Impey, C. D., & Tapia, S. 1988, ApJ, 333, 666 

329 Li, T.-P, & Ma, Y.-Q. 1983, ApJ, 272, 317 

330 Maraschi, L., Fossati, G., Tagliaferri, G., & Treves, A. 1995, ApJ, 443, 578 

331 Nieppola, E., Tornikoski, M., & Valtaoja, E. 2006, A&A, 445, 441 

332 Nikishov A. I. 1962 JETP, 14, 393 

333 Ong, R. A. 2009, The Astronomer's Telegram, 2084, 1 

334 Padovani, P., & Giommi, P. 1996, MNRAS, 279, 526 

335 Poole, T. S., et al. 2008, MNRAS, 383, 627 



-22 - 

sse Sbarufatti, B., Treves, A., & Falomo, R. 2005, ApJ, 635, 173 

337 Scarpa, R., & Falomo, R. 1995, A&A, 303, 656 

338 Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 

339 Siebert, J., Brinkmann, W., Laurent-Muehleisen, S. A., & Matsuoka, M. 2000, Advances in 

340 Space Research, 25, 729 

341 Teshima, M. 2009, The Astronomer's Telegram, 2098, 1 

342 Weekes, T. et al. 2002, Astroparticle Physics, 17, 221 



This manuscript was prepared with the AAS IATj^X macros v5.2. 



- 23 - 



W 

c 550 

0) 

> 

HI 

500 
450 
400 
350 

0.02 0.04 0.06 0.08 0.1 0.12 

9 2 [deg 2 ] 

Fig. 1. — Distribution of 9 2 for VERITAS events selected with soft cuts. The points with 
error bars denote the on-source events. The background is shown by the shaded histogram. 
The dashed vertical line shows the applied 6* 2 -cut. The expected distribution for a point 
source is given by the dotted line. 




-24- 



15 
10 
5 

15 
10 

5 



X10" 



12 



VERITAS (140-600 GeV) 



T 



X10' 



Fermi LAT (0.1 - 300 GeV) 



54880 54900 54920 54940 54960 54980 55000 

MJD 



CM 
■ 

E 

o 



3 
2 



1 — 



2 
1.5 
1 



2.5 
2 

1.5 bz, 



X10" 



* * 



XRT (2-10KeV) 



4 



t!!|ttlM| 

MM 



UVOT 



• 


u 




V 


• 


B 




W1 


• 


M2 




W2 



J I I I I I I I I I I I I I I I I I I I I I L 



MDM Observatory 



© © 

J I I I I L 



54992 54994 



© 


© 


© 


© 


© 


I 




i 




i 


54996 




54998 




55000 



J I I I I L 



O I 
* R 
+ V 



55002 55004 
MJD 



Fig. 2. — Light curves of PKS 1424+240 in VHE gamma rays (VERITAS), HE gamma rays 
(Fermi-L AT), X-rays (Swift XRT), UV (Swift UVOT) and optical (Swift UVOT, MDM). 
The X-ray, UV and optical light curves cover the time period indicated in the upper two 
light curves by the shaded region. The horizontal bars in the VHE and HE light curves 



-25 - 




Energy [TeV] 

Fig. 3. — The time averaged differential photon spectrum of PKS 1424+240 measured by 
VERITAS between February 19 and June 21, 2009. The triangles are from the soft-cuts 
analysis and the squares from the medium-cuts analysis. The flux point at 260 GeV is 
reconstructed in both analysis. The solid lines shows the fit with a power law. The shaded 
area shows the systematic uncertainty of the fit, which is dominated by a 20% uncertainty 
on the energy scale. 



-26- 



Table 1: SSC fit parameters for PKS 1424+240 as a function of assumed redshift. 



Parameter 


z = 0.05 


z = 0.10 


z = 0.2 


z = 0.3 


z = 0.4 


z = 0.5 


z = 0.7 


L e [10 43 erg s- 1 ] 


1.60 


4.12 


10.7 


18.9 


29.2 


47.1 


88.8 


L B [10 43 erg s- 1 ] 


1.66 


5.47 


16.9 


31.1 


45.9 


49.8 


66.2 


7i [10 4 ] 


3.7 


3.7 


3.6 


3.4 


3.2 


3.6 


3.7 


72 [10 5 ] 


4.0 


4.0 


4.0 


4.0 


4.5 


4.0 


4.0 


D 


15 


18 


25 


30 


35 


45 


60 


B[G] 


0.37 


0.31 


0.25 


0.24 


0.25 


0.18 


0.14 




1.04 


1.33 


1.59 


1.65 


1.57 


1.06 


0.75 


R B [10 16 cm] 


1.2 


2.2 


3.4 


4.0 


4.0 


4.5 


5.0 



-27- 




10 10 10 12 10 14 10 16 10 18 10 20 10 22 10 24 10 26 

v[Hz] 



Fig. 4. — SED of PKS 1424+240. The lines show SSC-model fits assuming different redshifts. 
The inset shows a zoom of the SED on the VERITAS data in a z/ 2 F^ representation. The 
Fermi data are presented together with their corresponding power-law fit and one standard 
deviation uncertainty. The upper limits correspond to 95% confidence level.