A Comparison Study of Intermediate Frequency BL Lacs and High Frequency BL Lacs Detected by VERITAS and Fermi LAT
Total Page:16
File Type:pdf, Size:1020Kb
A Comparison Study of Intermediate Frequency BL Lacs and High Frequency BL Lacs Detected by VERITAS and Fermi LAT Anna Louise Cannon Submitted to University College Dublin for the degree of Ph.D. in the College of Physics Febuary 2011 UCD School of Physics Prof. Lorraine Hanlon Supervised by Dr. John Quinn Abstract This thesis presents the comprehensive analysis and results of an extensive comparison study between two subclasses of blazars : Intermediate Frequency BL Lacs and High Frequency BL Lacs Detected by VERITAS and Fermi LAT. VERITAS (the Very Energetic Radiation Imaging Telescope Array System) is an array of four 12m telescopes, located at the Fred Lawrence Whipple Observatory in Southern Arizona and is sensitive in the energy range of 100 GeV - 30 TeV. Fermi LAT (Large Area Telescope) is a pair production satellite sensitive between 20 MeV - 300 GeV. Since the launch of Fermi LAT, there is for the first time, continuous coverage from MeV to TeV with an overlap between 100 and 300 GeV, crucial for constraining the inverse Compton peak and providing a clearer picture of the sources themselves, their environments and the processing taking place between 20 MeV and 30 TeV. Three HBLs (1ES 0502+675,1ES 0806+524 and 1ES 1959+650) and three IBLs (W Comae, 3C66A and PKS 1424+240) were selected based on previous results and publications. Quasi- contemporaneous Fermi LAT and VERITAS data for each source were analysis and combined to look closely at spectra. The inverse Compton peak was constrained and temporal analysis was performed. Multiwavelength spectral energy distributions were constructed and fit with both a single zone synchrotron model and synchrotron-external Compton combination model. ii Acknowledgements I would first like to express my utmost gratitude to my thesis supervisor Dr. John Quinn for his constant encouragement and patience throughout the course of my PhD. To everybody in the High Energy Astrophysics group in UCD, thank you for your help and efforts in the last few months of my write up, in particular, Anna O'Faolain deBhroithe for proof-reading endless pages of my thesis (I'll return the favour when it's your turn, I promise!). I would like to mention Dr Remo Huglii for the hours of conversation on just about everything from statistics to traveling and Leonard English for the music supply, most of it good. Thank you to Professor Lorraine Hanlon, Bairbre Fox and Marian Hanson who were always so approachable and willing to help out with anything. A special cheers to Dr. John Ward for his friendship over the last four years and was answering my never-ending questions. I wish to thank Dr Trevor Weekes for organizing my observing trips to Arizona and for being supportive. Also, thank you to Dr. Ken Gibbs who will be missed at basecamp. To Dr. Julie McEnery, thank you for the opportunity to work within the Fermi LAT collaboration and for organizing my trips to Goddard. To my good friends Arun Madhavan (co-creator of \tpg") and Dr. Sagar Godambe for making 12 hour observing shifts fun and for their advice and encouragement over the last year. Thank you to Dr. Teddy Cheung for the many, many laughs and for making my Maryland trips so entertaining. I would also like to mention Dr. Luis Reyes, Dr. Jeff Grube, Dr. Tyrel Johnson, William McConville and Sean Griffin who were always willing to answer physics/analysis queries and help out. It was greatly appreciated. Finally, I'd like to thank my my family and friends. To my parents for their unwavering support and encouragement, my sisters Christine (for helping me proof-read) and Mariosa (for keeping me entertained while I was writing up), my uncle Joseph and aunt Caitriona for making sure I was always fed and watered while in Dublin and my uncle John for the hours of science conversations. A special thank you to my closest friends Orla O'Neill and Susan Moore for their friendship and well-needed distractions. iii For Louis iv There are more things in heaven and earth, Horatio Than are dreamt of in your philosophy. - William Shakespeare, Hamlet Contents List of Figures vii List of Tables xix Chapter 1. Introduction 1 1.1. History of γ-ray astronomy 1 1.2. Brief Account of Satellite Experiments 1 1.3. Brief Account of Ground Based Experiments 7 1.4. Thesis Outline 10 1.5. My Contributions 11 Chapter 2. Active Galactic Nuclei 13 2.1. Introduction 13 2.2. Blazars 17 2.3. Variability 20 2.4. Relativistic Jets 21 2.5. Jet Dynamics 25 2.6. Synchrotron Radiation 27 2.7. Leptonic Models of Gamma-Ray Production 30 2.8. Hadronic Models 34 2.9. Motivation 36 Chapter 3. Cherenkov Radiation and its Detection from Extensive Air Showers 47 3.1. Cherenkov Radiation 47 3.2. Gamma-ray interactions with matter 49 3.3. Extensive Air Showers 50 3.4. Detection 51 3.5. Atmospheric Cherenkov Telescopes 52 3.6. Imaging Atmospheric Cherenkov Technique 53 Chapter 4. VERITAS and Fermi LAT Instrument Description 55 4.1. The VERITAS array 55 4.2. The Fermi LAT satellite 65 v CONTENTS vi Chapter 5. Data Analysis 76 5.1. Introduction 76 5.2. VERITAS 76 5.3. Fermi LAT 94 Chapter 6. Results 107 6.1. Introduction 107 6.2. Data Analysis 107 6.3. 1ES 0502+675 109 6.4. 1ES 0806+524 119 6.5. 1ES 1959+650 130 6.6. W Comae 139 6.7. 3C66A 150 6.8. PKS 1424+240 160 6.9. Constraining the IC peak 170 6.10. SED modeling 177 6.11. GeV dip 199 6.12. Discussion 201 Chapter 7. Conclusions and Future work 206 7.1. Thesis results 206 7.2. The Future of γ-ray Blazar Science 207 Bibliography 209 Appendix A. Astronomer's Telegrams 217 List of Figures 1.1 Schematic drawing of the γ-ray detector onboard NASA's Explorer satellite, taken from Kraushaar et al.,1965. 3 1.2 The second COS-B satellite catalog sky map, displaying 25 detections of γ-ray emitting sources above 100 MeV. The filled circles denote fluxes above 1.3 10−6 ph cm−2 cm−1 × and the open circles denote fluxes below this. This flux value is chosen to demonstrate the flux magnitude concentration in the galactic plane. The shaded region indicates unsearched areas - (Swanenberg et al., 1981) 4 1.3 Schematic diagram of EGRET's spark chamber, calorimeter and anti-coincidence shield (Thompson et al., 1993) 5 1.4 The third EGRET catalog γ-ray emission sky map, taken from (Hartman et al., 1999). The black diamonds represent AGN with the exception of Cen A, the black squares represent pulsars, the black triangle represents the Large Magellanic Cloud (LMC), the black circle represent solar flares and the open circles represent unidentified EGRET sources. 6 1.5 The Whipple 10m telescope located in Southern Arizona - (www.astro.wisc.edu/ larson/Webpage/Gamma.html) 9 ∼ 1.6 A HEGRA γ-ray telescope (www.sciencephoto.com/media/322268/enlarge). 10 2.1 Illustration of an AGN, highlighting the main components of AGNs and the classification dependent on orientation (http://www.auger.org/news/PRagn/about AGN.html). 15 2.2 AGN Classification including types 0, 1 and 2 classification (Adapted from Toner, 2008). 16 2.3 How AGN classification is influenced by viewing angle (Toner, 2008). 17 2.4 Multiwavelength SED of the blazar Markarian 421, including XMM-Newton OM (triangles), XMM-Newton EPN (filled circles), MAGIC (squares), Whipple (diamonds) and VERITAS (crosses). SSC model for 2006 data (green) and 2008 data (red) (Acciari et al., 2009). 18 2.5 From top to bottom (referring to radio luminosity) FSRQs (empty boxes), 1 Jy BL Lac sample (filled boxes) and Slew survey BL Lac sample (triangles) - (Fossati et al., 1998). 20 vii LIST OF FIGURES viii 2.6 Radio power at 5 GHz vs. the synchrotron peak frequency for FSRQ (filled points) and BL Lacs (crosses) for the DXRBS sample. The dotted lines denote the two quadrants (top-left and bottom-right) occupied by the sources studied by Fossati et al 1998. (Padovani, 2007). 21 2.7 Whipple 10m lightcurve for Markarian 421 from May 15th 1996, binned in 4.5 minute intervals - (Gaidos et al.,1996) 22 2.8 Illustrations of Toroidal Magnetic Field Acceleration (a) and Magneto-centrifugal acceleration from the accretion disk (b), (Sauty, Tsinganos & Trussoni, 2002) 23 2.9 Illustration of superluminal motion (physicsmadeeasy.wordpress.com/physics-made- easy/cosmology-iii/). 26 2.10The spectrum of the synchrotron radiation of a single electron with linear axes. 28 2.11The spectrum of the synchrotron radiation of a single electron with logarithmic axes. 28 2.12An illustration of Inverse Compton scattering of a low energy photon up to γ-ray energies (www.cv.nrao.edu). 30 2.13Daily average Whipple 10m lightcurve for 1ES 1959+650 for May - July 2002, with rate displayed in Crab units. The dashed line corresponds to the quickest flux change, indicating a variability timescale on the order of 7 hours (Holder et al., 2003) 38 2.14Multiwavelength SED of 1ES 1959+650 using November - December 2007 data, with the addition of 2002 HEGRA data. (Bottacini et al, 2010) 39 2.15Lightcurve in months bins in units of integral flux about 300 GeV from 1ES 0806+524 observations. A Chi-squared probability of being constant yields a value of 0.24, indicating little to no variability (Acciari et al., 2009). 40 2.16Broadband multiwavelength SED of 1ES 0806+524. The Swift data was taken on two separate days: grey points - 8th March 2008 and black points - 12th March 2008.