Open Collections

UBC Theses and Dissertations

UBC Theses Logo

UBC Theses and Dissertations

Investigating the Milky Way's magnetism : a statistical analysis of the structural scales in polarized.. Stutz, Robert Adrian 2011-12-31

You don't seem to have a PDF reader installed, try download the pdf

Item Metadata

Download

Media
[if-you-see-this-DO-NOT-CLICK]
ubc_2011_fall_stutz_robert.pdf [ 2.66MB ]
Metadata
JSON: 1.0072378.json
JSON-LD: 1.0072378+ld.json
RDF/XML (Pretty): 1.0072378.xml
RDF/JSON: 1.0072378+rdf.json
Turtle: 1.0072378+rdf-turtle.txt
N-Triples: 1.0072378+rdf-ntriples.txt
Original Record: 1.0072378 +original-record.json
Full Text
1.0072378.txt
Citation
1.0072378.ris

Full Text

Investigating the Milky Way’s Magnetism: A statistical analysis of the structural scales in polarized radio observations of our Galaxy by Robert Adrian Stutz B.Sc. (Hons.), The University of British Columbia, 2009 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF SCIENCE in The College of Graduate Studies (Interdisciplinary Studies) THE UNIVERSITY OF BRITISH COLUMBIA (Okanagan) September 2011 c Robert Adrian Stutz 2011Abstract Angular power spectra are calculated and presented for the entirety of the Canadian Galactic Plane Survey polarization dataset at 1.4 GHz. The data analyzed is a combination of data from the 100-m E elsberg Telescope, the 26-m telescope at the Dominion Radio Astrophysical Observatory, and the Synthesis Telescope at the Dominion Radio Astrophysical Observatory, allowing all relevant scales to be sampled down to arcminute resolution. Calculated angular power spectra cover multipoles from ‘  60 to ‘  104, and display at least three distinct regions within this range. A power-law region at low multipoles, a  attening region at intermediate multipoles, and a taper region at high multipoles necessitate a model  tting routine for analysis. Power-law indices are found to have a mode of 2.3, similar to other results, though slightly steeper. An increasing trend in power-law index is found with increasing latitude. A transition from small-scale to large-scale structure is evident at b = 10 . Flattening of the power spectra is also investigated and found to be likely due to background point sources, though this remains inconclusive. iiTable of Contents Abstract . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ii Table of Contents . . . . . . . . . . . . . . . . . . . . . . . . . . . . iii List of Figures . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . v Acknowledgements . . . . . . . . . . . . . . . . . . . . . . . . . . . ix 1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 1.1 Polarization Studies and Angular Power Spectra . . . . . . . 1 1.2 Canadian Galactic Plane Survey Polarization Maps . . . . . 3 1.2.1 E elsberg Telescope Observations . . . . . . . . . . . 5 1.2.2 26-m Telescope Observations . . . . . . . . . . . . . . 5 1.2.3 The Synthesis Telescope at the DRAO . . . . . . . . 5 1.2.4 Combination of Synthesis and Single Dish Datasets . 6 2 Analysis and Results . . . . . . . . . . . . . . . . . . . . . . . . 10 2.1 Calculating Angular Power Spectra . . . . . . . . . . . . . . 10 2.2 Model Fitting to APS . . . . . . . . . . . . . . . . . . . . . . 12 2.3 Point Source Removal . . . . . . . . . . . . . . . . . . . . . . 13 2.4 Maps of Spectral Indices . . . . . . . . . . . . . . . . . . . . 14 2.5 Linear Fitting to Flat Power . . . . . . . . . . . . . . . . . . 14 2.5.1 Maps of Flattening . . . . . . . . . . . . . . . . . . . 15 2.6 Modeling Flat Power Contributions . . . . . . . . . . . . . . 15 2.6.1 Noise Contributions to APS . . . . . . . . . . . . . . 16 2.6.2 Point Source Contributions to APS . . . . . . . . . . 17 2.7 Investigating the Non-Gaussian Taper . . . . . . . . . . . . . 18 3 Discussion and Conclusion . . . . . . . . . . . . . . . . . . . . 34 3.1 Spectral Indices and Flattening across the CGPS . . . . . . . 34 3.1.1 Model Failures . . . . . . . . . . . . . . . . . . . . . . 34 3.1.2 Spectral Index and Flattening Features . . . . . . . . 35 iiiTable of Contents 3.1.3 Spectral Index Statistics . . . . . . . . . . . . . . . . 36 3.1.4 Spectral Index and Galactic Latitude . . . . . . . . . 40 3.1.5 Flattening Statistics . . . . . . . . . . . . . . . . . . . 41 3.1.6 Flattening and Galactic Latitude . . . . . . . . . . . 41 3.2 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42 Bibliography . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45 Appendices A Synchrotron Emission . . . . . . . . . . . . . . . . . . . . . . . 48 B Polarization . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51 B.1 Stokes Parameters . . . . . . . . . . . . . . . . . . . . . . . . 52 B.2 Faraday Rotation . . . . . . . . . . . . . . . . . . . . . . . . 53 B.3 Depolarization Mechanisms . . . . . . . . . . . . . . . . . . . 54 ivList of Figures 1.1 DRAO synthesis telescope con guration. Antenna separa- tions given in units of L = 4:286m. The separation between antennas 1 and 2 can be modi ed from 3L to 14L by move- ment along the track. . . . . . . . . . . . . . . . . . . . . . . . 7 1.2 Arrangement of continuum and spectroscopy bands around 1420 MHz for the synthesis receiver. Hydrogen spectroscopy is done in the S band; continuum observations occur in bands A,B,C,D. Frequencies below the bands give the intermediate frequencies that the signals are converted to. Protected bands allocated to radio astronomy are also indicated. . . . . . . . . 8 1.3 Fig. 3 from Landecker et al. [2010]. Weighting of the three datasets over the spatial coverage of the survey. Also shown is the apodization function applied to the Synthesis Telescope data at the longest baselines. . . . . . . . . . . . . . . . . . . 9 2.1 Sample power spectrum of a single 480x480 pixel submap centered at l = 120:7 , b =  1:7 . Power-law behaviour is evident at low multipoles, as is the tapering e ect as ‘! 104. Pixel e ects are present beyond the taper, these are discussed in the text. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11 2.2 Sample power spectra of two 480x480 pixel submaps located at l = 150:0 , b =  1:3 and b = 0:0 respectively. The model of best  t is overplotted in red, and its parameters are given as well. The model consistently diverges from the data at the highest multipoles, where the taper appears ill- t by the Gaussian model. . . . . . . . . . . . . . . . . . . . . . . . 19 vList of Figures 2.3 Sample power spectra of two 480x480 pixel submaps located at l = 150:0 , b = 1:3 and b = 2:7 respectively. The model of best  t is overplotted in red, and its parameters are given as well. The model consistently diverges from the data at the highest multipoles, where the taper appears ill- t by the Gaussian model. . . . . . . . . . . . . . . . . . . . . . . . . . 20 2.4 Angular power spectra and best- t model with parameters for a submap centered at l = 174:7 , b =  1:7 . The submap of this region contains signi cant power in a point source. The top power spectrum is computed without point source removal, while the bottom power spectrum has undergone point source removal. Of note are the varying model param- eters and the overall steepening of the power spectrum after point source subtraction. . . . . . . . . . . . . . . . . . . . . . 21 2.5 Angular power spectra and best- t model with parameters for a submap centered at l = 167:3 , b =  1:3 . The submap of this region contains signi cant power in a point source. The top power spectrum is computed without point source removal, while the bottom power spectrum has undergone point source removal. The model  t is superior in the latter  gure. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22 2.6 Angular power spectra and best- t model with parameters for a submap centered at l = 156:7 , b =  1:7 . The submap of this region contains little power in point sources. The top power spectrum is computed without point source removal, while the bottom power spectrum has undergone point source removal. No signi cant changes occur in the model parameters. 23 2.7 Power spectral index ( ) maps created from the APS anal- ysis of supermosaic 6. The top  gure has not undergone point source removal, whereas the bottom  gure has under- gone point source subtraction. Evident in the top map are the sharp transitions in  due to the presence of bright point sources which dominate the power in the di use emission at large multipoles. . . . . . . . . . . . . . . . . . . . . . . . . . 24 2.8 Maps of power spectral index values ( ) in colour and polar- ized intensity in black and white for l = 66 to l = 104 ; all  values from model  tting to power spectra. . . . . . . . . . 25 2.9 Maps of power spectral index values ( ) in colour and polar- ized intensity in black and white for l = 102 to l = 140 . . . 26 viList of Figures 2.10 Maps of power spectral index values ( ) in colour and polar- ized intensity in black and white for l = 138 to l = 176 . . . 27 2.11 Power spectral-index ( ) map for l = 101 to l = 116 , b = 3 to b = 17 . Polarized intensity map of the same region in black and white. This is the high-latitude extension of the CGPS. Note that the colourbar di ers from those of the previous spectral index maps due to a greater range in values in this region. A clear break in  is evident at b = 10 . . . . 28 2.12 Flatness parameter ( ) maps for l = 67 to l = 121 . Higher values (red) indicate less  attening, while lower values (blue) indicate signi cant  attening. . . . . . . . . . . . . . . . . . . 29 2.13 Flatness parameter ( ) maps for l = 121 to l = 175 . Higher values (red) indicate less  attening, while lower values (blue) indicate signi cant  attening. . . . . . . . . . . . . . . . . . . 30 2.14 Flattening ( ) map for l = 101 to l = 116 , b = 3 to b = 17 , the high latitude extension. As seen in the  map of the same region, a sharp transition is evident around b = 10 . 31 2.15 Power spectrum of a single synthesis  eld (K3) without added single-dish information. Of note is the  attening of the power spectrum for ‘ > 1000; this  attening is also evident in the APS of the combined dataset. . . . . . . . . . . . . . . . . . . 32 2.16 Power spectrum for a submap of the Synthesis Telescope ob- servation of  eld K3 in Stokes V . This power spectrum rep- resents the noise power present in this  eld, and by extension an approximation of the noise present in the supermosaics at synthesis multipoles. . . . . . . . . . . . . . . . . . . . . . . . 32 2.17 Angular power spectrum of a model point source with 0.4Jy of  ux density. Map containing point source is 480 pixels to a side, with an image scale equal to that of the supermosaics at 0:00556 /px. . . . . . . . . . . . . . . . . . . . . . . . . . . 33 2.18 Angular power spectrum of a model point source with 0.4Jy of  ux density. The observational e ect of the Synthesis beam is approximated by rotating and elongating the point source into an ellipse. The model does not accurately  t the taper region for this elliptical beam shape. . . . . . . . . . . . . . . 33 3.1 Histogram of power spectral-index values derived from Stokes PI supermosaics. . . . . . . . . . . . . . . . . . . . . . . . . . 36 3.2 Histogram of power spectral-index values derived from Stokes Q supermosaics. . . . . . . . . . . . . . . . . . . . . . . . . . 37 viiList of Figures 3.3 Histogram of power spectral-index values derived from Stokes U supermosaics. . . . . . . . . . . . . . . . . . . . . . . . . . . 38 3.4 Histogram of  values from the high latitude extension of the CGPS in Stokes I. Two peaks are visible due to the sharp increase in  at b = 10 . . . . . . . . . . . . . . . . . . . . . . 39 3.5 Power spectral index  plotted versus galactic latitude b to view averaged trends across the visible galactic disk. . . . . . 40 3.6 Histogram of  values derived from the Stokes I supermosaics. 42 3.7 Histogram of  values derived from the Stokes I map of the high latitude region. . . . . . . . . . . . . . . . . . . . . . . . 43 3.8 Flatness parameter  plotted versus galactic latitude b. . . . 44 viiiAcknowledgements I would like to extend my thanks to my supervisor, Professor Erik Rosolowsky, for his continued support and helpful advice throughout the investigation and writing of this thesis. I would also like to thank Dr. Tom Landecker for providing both the original research direction as well as con- tinual advice and invaluable knowledge of the subject matter. Dr. Roland Kothes and Dr. Marijke Haverkorn also deserve considerable gratitude for their helpful input and illuminating discussions. The sta at the Dominion Radio Astrophysical Observatory in general were exceedingly helpful and supportive, for which I am very grateful. Finally, I must thank my partner Natasha Drew for continuing to support and encourage me in my attempts at completing this thesis, and for bearing the brunt of the fallout associated with this endeavour. ixChapter 1 Introduction 1.1 Polarization Studies and Angular Power Spectra The magnetic  eld of our Milky Way Galaxy is commonly divided into two distinct but equally important components: the large-scale  eld on scales of the spiral arms of the Galaxy, and the small-scale  eld on scales of stellar activity. Both components are locked into the interstellar medium plasma by their interactions with charged particles. The large-scale  eld has a magnetic  eld strength of  2 G at the solar radius, and appear to follow the spiral arms of increased matter density within the Galaxy, stay- ing parallel to the plane of the galactic disk. The magnetic  eld strength increases to  10 G at the galactic center [Landecker, 2011]. The magnetic  eld direction of this large-scale  eld undergoes a reversal inside the solar radius, and other  eld reversals have been proposed but remain a subject of debate [Beck, 2009]. This large-scale  eld is not sampled directly by the observations presented in this thesis, but it remains an important topic in galactic evolution. The small-scale component of the magnetic  eld is a di use quasi-random  eld, exhibiting none of the overlying structure associated with the large- scale  eld. Individual structures in the di use magnetic  eld reach scales on the order of tens of parsecs. Nonetheless, the small-scale component has a similar strength to that of the large-scale  eld, reaching  3 G at the solar radius [Landecker, 2011]. As the magnetic  eld is locked into the interstellar medium, any perturbations of this medium due stellar events also a ect the magnetic  eld; events like supernovae directly alter this small-scale di use magnetic  eld. The turbulence of the di use  eld is greater in the arm regions compared to the interarm regions due to the increased stellar activity within the Galaxy’s spiral arms. This component of the magnetic  eld is indirectly studied in this thesis by analysing the continuum emission associated with the magnetic  eld lines. The continuum emission that is detected by radio telescopes when ob- 11.1. Polarization Studies and Angular Power Spectra serving the plane of our Galaxy is due to the spiraling of cosmic rays around the magnetic  eld lines of the Galaxy. These magnetic  eld lines may be part of the large regular component of the  eld or small-scale irregular com- ponents. Cosmic ray electrons are a ected by this magnetic  eld and begin to spiral along the  eld lines, producing continuum synchrotron emission [Burke, 2010]. This emission may be partially polarized, and as it propa- gates on its way to Earth, it is a ected and altered by the warm magnetized interstellar medium through which it passes. Observations of this emission by radio telescopes on Earth thus hold information both about the origi- nal emission mechanism and strength as well as the medium through which it travels. Unfortunately the study of radio polarization is complicated by the lack of correlation between radio polarization images and total intensity radio images. Faraday depolarization and other depolarization mechanisms alter the original emission to the point where it bears little resemblance to physical objects at the point of emission. Increasing the understanding of polarization observations is essential for understanding evolutionary galactic dynamics and the prevalence of our Galaxy’s cosmic magnetism. The magnetic  eld of our Galaxy is expected to have an impact on many of the processes within the Galaxy from indi- vidual star formation to the density waves in the Galaxy that produce its spiral arms. Increasing the understanding of the cosmic magnetism in our Galaxy will allow us to more completely describe the dynamical processes which shape its interstellar medium and the structures and objects in that interstellar medium which it in uences. Structures visible in polarized intensity are often entirely missing in total intensity maps of the same regions. Interpretation of polarization maps is further confounded by a lack of knowledge on the distance scales involved in depolarization e ects. That is, it is often impossible to determine whether the modulation of polarized emission occurs at a signi cant distance from the point of emission, or whether it occurs at a much greater distance along the path of propagation through the warm interstellar medium. These factors compound to produce polarization images which not only appear completely di erent from total intensity images but are also hard to classify into indi- vidual objects. Distant emission may be modulated many times on its path to Earth, and it is nearly impossible to determine the foreground object responsible for the most signi cant e ects [Landecker et al., 2010]. Polarization observations are subject to the same resolution constraints as their corresponding total intensity observations; thankfully, they also per- mit the use of synthesis radio telescopes to improve their resolution. The fundamental Rayleigh criterion states that the angular resolution of a tele- 21.2. Canadian Galactic Plane Survey Polarization Maps scope is proportional to the observed wavelength and inversely proportional to the diameter of the telescope; a smaller angular resolutions allows a tele- scope to resolve smaller objects. As radio telescopes operate at the longest observable wavelengths they must compensate by increasing their dish sizes. Alternatively, aperture synthesis techniques allow pairs of radio telescopes to act as a single-dish telescope of a diameter equal to their separation. The use of synthesis telescopes to sample the smallest scales in radio observa- tions has its own drawbacks, namely that large-scale structure is lost due to its inability to be sampled by the phyisically separated antenna of such a telescope. Large single-dish telescopes operating in the radio regime sam- ple such large scales accurately, but cannot approach the size of synthesis arrays [Verschuur, 2007]. Only with a combination of single-dish and syn- thesis telescope observations can all scales of the galactic polarized emission be sampled accurately. The analysis of polarized intensity data is often done in terms of the structural scales in the images. Angular power spectra, which convey the power present at speci c angular scales, are used to determine the distribu- tion of power in the observations. Since identi cation and taxonomy of indi- vidual objects in polarized intensity maps is di cult, such broad-approach statistical analyses allow us to describe the observed polarized medium with- out having to deal with classi cation issues. In previous studies the power spectra are generally given as a function of multipole ‘, which is inversely proportional to angular scales in the observations. Speci cally, ‘  180 = , where  corresponds to the angular scales in the image in degrees. Fur- ther discussion will focus on both angular scales for simplicity and their corresponding multipoles for the sake of comparison with previous studies. The radio continuum survey used in this thesis is detailed below in Sec- tion 1.2, and analyzed in Section 2. The results of the analysis are described in Section 3.1. The source of the observed radio emission - synchrotron emission - is discussed in Appendix A. Important polarization de nitions and depolarization e ects are outlined in Appendix B. 1.2 Canadian Galactic Plane Survey Polarization Maps The polarization maps analyzed in this thesis are part of the Canadian Galactic Plane Survey (CGPS) datasets presented in Landecker et al. [2010]. This survey is an international collaboration led by Canadian researchers, with its principal aim being to fully map the total intensity and the polarized 31.2. Canadian Galactic Plane Survey Polarization Maps intensity of the portion of the Milky Way galaxy visible in the northern hemisphere. Continuum polarized emission around 1420 MHz was measured for 1060 square degrees of the nothern Galactic plane down to a resolution of  1 arcminute. The survey extends from a galactic longitude of l = 66 to l = 175 over the latitude range from b =  3 to b = 5 . The survey area is divided into six data maps, termed supermosaics, that each cover the entirety of the observed galactic latitudes and a range of galactic longitudes. In addition to the main survey a high-latitude extension up to b = 17:5 was observed from l = 101 to l = 116 . The CGPS polarization survey is unique both in size and in the fact that the  nal dataset is a combination of observations from three seperate telescopes. Encompassing observations from  1995 to  2005, the data is a combination of observations from the E elsberg 100-m Telescope, the Dominion Radio Astrophysical Observatory (DRAO) 26-m Telescope, and the DRAO Synthesis Telescope. The aim of this combination is to fully sample the polarized sky on all angular scales, especially the largest scales that are missed by other synthesis-only surveys. A brief analysis of the polarization maps is also given in Landecker et al. [2010], with the most salient points discussed here. As seen throughout the literature of polarization studies, the polarized intensity sky and total in- tensity sky have few similarities. This is visually obvious in Figures 5 - 7 in Landecker et al. [2010], and is generally agreed to be due to Faraday ro- tation e ects and depolarization. The polarization fraction throughout the survey is near 10%, signi cantly lower than the theoretical maximum of 70%, providing further evidence for the existence of a signi cant depolarization mechanism. The polarized sky also appears to have features on wildly vary- ing scales. Large-scale polarization features are evident, and are assumed to be connected with the large-scale features of the galactic magnetic  eld itself, being roughly aligned with the spiral arms of the Galaxy. Small-scale variations appear more randomly distributed, and may re ect the irregular component of the magnetic  eld and the turbulence of the ionized gas. Power may be shifted from larger scales to smaller scales by depolarization e ec- tively randomizing areas of the large-scale emission. Landecker et al. [2010] note that at the survey frequency of 1420 MHz the irregular component is obvious but does not dominate the larger regular component. As Faraday rotation increases at longer wavelengths, it is conceivable that a wavelength exists beyond which the regular component of the polarized emission is en- tirely obscured. This thesis concerns itself mainly with the distribution of power between the large-scale regular component of the emission (and thus magnetic  eld) and the small-scale irregular component. 41.2. Canadian Galactic Plane Survey Polarization Maps 1.2.1 E elsberg Telescope Observations Part of the CGPS polarization data is made up of observations from the 100-m Telescope near E elsberg, Germany. The polarization dataset observed at the E elsberg 100-m telecope are part of the results from the E elsberg Medium Latitude Survey, which ran from 1994 to 2005 [Reich et al., 2004]. While various polarimeters were utilized during the survey, the telescope nominally observed at a central frequency of 1400 MHz with a bandwidth of 20 MHz. At this frequency the telescope’s main beam is approximated by a Gaussian with a half-power-beam-width of 9.35’, and the total system noise temperature is around 26 K. While the total intensity is noise-limited around 15 mK, the sensitivity in Stokes Q and U individually is as low as 8 mK. The receiving system is further discussed in Uyaniker et al. [1998], and any further parameters important for the CGPS data are outlined in Landecker et al. [2010]. 1.2.2 26-m Telescope Observations A second component of the full datasets presented in Landecker et al. [2010] is provided by the 26-m Telescope located at the DRAO. This tele- scope is actually a 25.6 m equatorially mounted, axially symmetric dish. At an observing frequency of 1.4 GHz this dish has a beamwidth of 36’, signi - cantly lower than the resolutions for the other two telescopes. However, the data are corrected for ground radiation, allowing the telescope to sample large angular scales. The system temperature of the dish is around 125 K, with a sensitivity of 12 mK in Stokes Q and U. The bandwidth of the po- larization measurements around 1410 MHz is 12 MHz, for which range the polarimeter provides the normal power channels of RR, LL and the cross products RL and LR. Parameters of signi cance for the survey can be found in Landecker et al. [2010], with a full treatment of the receiving system presented in Wolleben et al. [2006]. 1.2.3 The Synthesis Telescope at the DRAO The third component of the  nalized datasets presented in Landecker et al. [2010] is provided by the Synthesis Telescope at the DRAO. The data presented and analyzed in this thesis deals primarily with the combined polarization datasets, but some emphasis is also given to the synthesis data alone. For this reason a more thorough treatment of the synthesis telescope’s parameters is given below. The present data were observed between 1995 51.2. Canadian Galactic Plane Survey Polarization Maps April and 2004 March, so the given parameters are an average description of the telescope array and its characteristics over that period. The synthesis telescope located at the DRAO consists of seven equa- torially mounted parabolic antennas arrayed in an east-west con guration. An array of seven antennas gives only 21 instantaneous baselines, thus the east-west con guration was chosen to allow for Earth’s rotation to vary the baseline vectors. Four of the antennas remain in a  xed position in the array, while three are moveable along a rail track running east to west. Having such movable elements in the array allows for intermediate baselines to be sampled. Five of the antennas have a diameter of 8:53m, the other two have a diameter of 9:14m. The con guration of the antennas as seen in Figure 1.1 has a baseline increment of L = 4:286m, approximately half the antenna diameters of 8:53m. Complete coverage of all baselines from 3L to 141L (plus 144L) can be obtained with 12 di erent array con gurations. The to- tal time required for complete sampling of the entire aperture is then 12 12 hours [Landecker et al., 2000]. The minimum usable baseline corresponds to 3L ’ 12:86m due to me- chanical interference considerations. This corresponds to a maximum size for sampled structures of  40’ at 1420 MHz. The largest baselines sampled by the synthesis telescope are 144L ’ 617:1m, giving a resolution of 58". [Landecker et al., 2000]. The di erence in antenna diameters results in signi cantly di erent weight- ings of u-v samples, but a method of corrections has been devised by Willis [1999]. Antenna noise temperatures also vary from a maximum of 26.0 K for one antenna to a minimum of 17.8 K for another. The synthesis telescope had a total system temperature of 60 K for most of the survey, improv- ing to 45 K in 2003. The telescope’s sensitivity was consequently improved from 0.30 mJy beam 1 to 0.23 mJy beam 1. The frequency coverage of the synthesis receiver bands is displayed in Figure 1.2. Continuum polarization measurements utilize bands A,B,C, and D of the receiver, each of width 7.5 MHz. 1.2.4 Combination of Synthesis and Single Dish Datasets The CGPS dataset presented in Landecker et al. [2010] is unique in that it is the  rst extensive polarization survey that accurately samples informa- tion on all scales by combining synthesis and single-dish telescope data. The data from the DRAO Synthesis Telescope, the 100-m E elsberg Telescope, and the DRAO 26-m Telescope are combined in an attempt to mitigate problems with the individual telescopes. The motivation and results of this 61.2. Canadian Galactic Plane Survey Polarization Maps Figure 1.1: DRAO synthesis telescope con guration. Antenna separations given in units of L = 4:286m. The separation between antennas 1 and 2 can be modi ed from 3L to 14L by movement along the track. combination are brie y outlined below, with full details available in Lan- decker et al. [2010]. The motivation for combining datasets arises from the issues that each telescope faces. While the Synthesis Telescope has an adequate maximum resolution corresponding to a multipole of ‘  7700, it is limited to observing baselines longer than 12.86 m, the minimum phyical separation of antennas. This corresponds to an angular size of  40’, or a multipole of ‘  270. Further complications arise as the ground radiation is detected as polarized emission in the sidelobes of both antennas when spaced so closely. This makes the 12.86 m baselines unusable for polarization imaging, so larger scales and extended structure remain unsampled by the Synthesis Telescope [Landecker et al., 2010]. The 100-m E elsberg Telescope has an adequate resoluion of  9:35’ and can sample larger scales, but it is similarly limited by the data processing. A linear baseline is subtracted from each  eld in order to set the edges to zero, which e ectively removes ground radiation. Unfortunately any real data at scales approaching the  eld sizes would also be removed by this process. The E elsberg  elds range in size from 5 to 16 ; any structures on the size scale of the  eld or larger will be lost. Taking the 5  elds as the limiting case, the E elsberg Telescope nominally samples multipoles between ‘  40 and ‘  1200. The data from the smaller 26-m Telescope has correct zero levels and does not require removal of the ground radiation, allowing the largest scales and smallest multipoles to be sampled. However, the resolution for the 26-m Telescope of  36’, corresponding to a multipole of ‘  300, is signi cantly 71.2. Canadian Galactic Plane Survey Polarization Maps Figure 1.2: Arrangement of continuum and spectroscopy bands around 1420 MHz for the synthesis receiver. Hydrogen spectroscopy is done in the S band; continuum observations occur in bands A,B,C,D. Frequencies below the bands give the intermediate frequencies that the signals are converted to. Protected bands allocated to radio astronomy are also indicated. worse than the other two telescopes. The separations between drift scans of the 26-m Telescope are also often too large to provide a fully sampled coverage of the survey; the sky sampling is only 41% of full Nyquist. The E elsberg dataset is used to  ll in the gaps in the 26-m data. Comparisons of the datasets were undertaken to demonstrate the fea- sibility and accuracy of combining the data. A pixel-by-pixel comparison plot of the same sky region in the 26-m and 100-m datasets does not give a slope of exactly 1 as expected, nor does it pass through the origin. This systematic o set remains unexplained, but Landecker et al. [2010] nonethe- less conclude that the two intensity scales are within  5%. A comparison of the E elsberg and Synthesis datasets shows that at baselines below 15 m the power in the E elsberg data exceeds that from the Synthesis Tele- scope. This is expected due to the inability of the Synthesis Telescope to sample these short baselines. Apart from this expected discrepancy the two 81.2. Canadian Galactic Plane Survey Polarization Maps Figure 1.3: Fig. 3 from Landecker et al. [2010]. Weighting of the three datasets over the spatial coverage of the survey. Also shown is the apodiza- tion function applied to the Synthesis Telescope data at the longest baselines. datasets are found to be within  10%. A combination of the three datasets allows for the 26-m data to sample the largest structures, the E elsberg data to  ll in the gaps in the 26-m data, and the Synthesis telescope data to provide small-scale information with its superior resolution. A simpli ed approximation of the weighting given to each telescope dataset across the baselines of the survey is given in Figure 1.3 taken from Landecker et al. [2010]. A Gaussian apodization function is also applied at the longest baselines strictly in the Synthesis regime, in order to diminish the contributions of the sidelobes to the antenna beam. The applied Gaussian reaches 20% at the longest baseline of 617m. The transitions between telescopes in Figure 1.3, where the weightings reach 0.5, occur at a baseline of 9m for the 26-m to 100-m transition, and a baseline of  33m for the 100-m to Synthesis transition. The correspond- ing multipoles are ‘  110 and ‘  400 respectively. The extent of the apodization function means the Synthesis Telescope samples multipoles up to ‘  7700, which is then also the upper limit of the combined dataset. 9Chapter 2 Analysis and Results 2.1 Calculating Angular Power Spectra To determine the distribution of power across the angular scales of the CGPS polarization data we calculate and analyze angular power spectra (APS) of the data. This calculation of APS can also be used for the single- dish maps. The method for calculating APS closely follows that of Haverkorn et al. [2003], though the analysis di ers signi cantly due to slightly better resolution and the inclusion of larger scales in the CGPS data. The map to be analyzed is Fourier transformed into the spatial frequency domain and then multiplied by its complex conjugate to calculate power: PSX(‘) = F(X) F  (X): (2.1) X may be any of the Stokes parameters Q, U , or PI for which polarization maps exist, though power spectra presented here are generally calculated for polarized intensity PI unless otherwise noted. Due to the nature of the Fourier transform, this is best accomplished by taking square maps or submaps from the supermosaics. The calculated PSX is a two-dimensional representation of the power data, so it is averaged in radial bins to determine the one-dimensional power behaviour. The radial bins are taken outward from the center of the Fourier transformed image in one-pixel increments. The standard deviation within each ring is also calculated and is divided by the square root of the number of data points to give the standard error in the mean. The units of the mean - the reported power - are K2, as the original maps are in units of K; taking the square of the Fourier transform of the brightness temperature maps gives a result in units of brightness temperature squared. As the power spectra are calculated in pixel-based rings, a conversion from pixel scales to angular scales is necessary. The scales covered by any power spectrum depend on the size of the map analyzed, and the Fourier transform essentially inverts the scale so that the angular scale is simply  =  map=dpx.  map is the width of the analyzed map in degrees, and dpx is the distance of the radial bins in pixels. The  calculated in this way 102.1. Calculating Angular Power Spectra are converted to multipole ‘ by the relation ‘  180 = , and any bad data points are discarded. Figure 2.1: Sample power spectrum of a single 480x480 pixel submap cen- tered at l = 120:7 , b =  1:7 . Power-law behaviour is evident at low multipoles, as is the tapering e ect as ‘ ! 104. Pixel e ects are present beyond the taper, these are discussed in the text. A representative power spectrum is given in Figure 2.1. The power decreases approximately linearly from the lowest sampled multipole to ‘   103. On this log-log plot the linear portion is a power law. We consider the trend in this portion of the plot as a simple power law of the form PSX(‘) / ‘  . The spectral index  of polarization power spectra is usually the result of interest, as in eg. Haverkorn et al. [2003], Giardino et al. [2002], Carretti et al. [2010], since it describes the distribution of power. Other features of the power spectra are rarely discussed, though Haverkorn et al. [2003] does mention a  attening at high multipoles due to noise in the image. The applied taper discussed in Section 1.2.4 appears in Figure 2.1 from approximately ‘  3000to  7700. Power beyond this multipole in Figure 2.1 cannot be a sampling of real signals, only pixel noise and other imaging e ects. This is also evident in the plot itself, as the taper is replaced by these e ects around ‘  104. Any further analysis and plots of power spectra has been limited to ‘  7700. Only slightly visually evident in Figure 2.1, but easily found analyti- cally, is a break in the power-law behaviour around ‘  103. At multipoles lower than this value the power law is considerably steeper than at higher 112.2. Model Fitting to APS multipoles. Unfortunately the applied taper confounds the interpretation of this  attening behaviour at high multipoles; nonetheless an investigation and interpretation is attempted in Section 2.6. In general the APS show complex behaviour across the multipole range and are not fully described by a single power-law with some noise  attening at high multipoles; a more detailed analytical model is necessary. 2.2 Model Fitting to APS In order to determine accurate spectral index values while accounting for the  attening and taper e ects at higher multipoles, we apply least-squares  tting to a given model. The  attening  oor is modeled as a constant o set and the taper as a Gaussian apodization function, giving the following model: P =  A+B (l)   exp   l2 2 2  : (2.2) The spectral index  determined by this model should be comparable to those in the literature, which generally  t only multipoles below ‘ = 1000. Sample power spectra for a single galactic longitude are presented in Figures 2.2 and 2.3. Two things of note become apparent in the attempted  ts of the model to a single submap power spectrum as in Figures 2.2 and 2.3. First, it appears that the lowest sampled multipole has less power than the power- law  t expects in certain cases; and second, the Gaussian model taper does not accurately  t the observed taper at high multipoles. Variations in the di erent power spectra can reduce the discrepancy at the lowest multipole, but the non-Gaussian taper e ect appears consistantly. Whether the e ect is due to an incorrect model or a physically non-Gaussian taper is investigated in Section 2.7. Multiple APS are calculated for each supermosaic so that variations in the  tting parameters across the galactic plane can be studied. A submap with 480 pixels to a side is taken from the full supermosaic, and the above analysis is applied to the submap. The submap is shifted by 60 pixels in the x and y directions, and analyzed again. This process is repeated until the entire supermosaic is sampled by submaps of 480 pixels which are themselves oversampled by a factor of 8. The supermosaics have a plate scale of 0:00556 /px, so a 480 pixel submap has  2:7 to a side, and covers an area of  7:1 2. As the submap size limits the range of multipoles on the lower end, a 480 pixel submap allows multipoles down to ‘  67 122.3. Point Source Removal to be accurately sampled. Information at larger scales is still present in the images, but will contribute only to the total power in the image and not individual multipoles. We determine model parameters for each of the 480 pixel submaps, which we then combine in an image to determine the parameter trends across the supermosaics. 2.3 Point Source Removal Point sources on the scale of the Synthesis Telescope beam size are ap- parent in the polarized intensity maps, and are found to add considerable power to the APS in which they are included. In some instances the pres- ence of strong point sources leads to the model  t failing to converge to a solution, though this appears to be due to remaining imaging artifacts that were not removed fully in the cleaning process. Such image artifacts are prevalent especially around the strongest sources. Point sources are primar- ily expected to be extragalactic radio galaxies, showing no dependence on galactic latitude; a small portion of the point sources may be compact super- nova remnants within the Galaxy itself, these being indistinguishable from extragalactic point sources. To study the di use polarized emission within our galactic plane the e ects of such point sources should be mitigated, and this is accomplished in this investigation by their removal prior to the APS analysis. We apply a standard Roberts edge enhancement  lter to the polarization maps and any objects in the maps remaining above a threshold of 0.3K are marked for analysis and possible removal. Such objects are individually  tted with a two-dimensional Gaussian, and if the Gaussian full-widths at half maximum are less than 8 pixels (0:04 ) the Gaussian  t is subtracted from the point source. The baseline of the  t remains, setting the previous point source to a brightness temperature comparable to its surroundings. It should be noted that such removal still leaves a feature in the polarization maps on the size scale of the telescope beam, leaving at least some of the original power in the map. If the attempted Gaussian  tting fails we discard the object as being a non point source and it remains in the polarization maps. Figures 2.4, 2.5, and 2.6 show a comparison of APS before and after point source removal. All three sample regions are taken from the sixth supermosaic, with the  rst two having signi cant point sources in the area, while the third has only minimal point sources and serves as a comparison. The model parameters for the regions containing point sources are altered 132.4. Maps of Spectral Indices by the removal of the point sources: the constant power o set parameter A is reduced while the power spectral index  is increased. Both e ects are expected as the removal of the point sources signi cantly lowers the total power in the region, which has a larger relative impact at higher multipoles where the power is already lower. The reduction in small scale structure by removal of the point sources steepens the power spectra, which should then be representative of the di use emission in the Galaxy. A comparison of the entire point-source-subtracted supermosaic 6 and a non-altered supermosaic 6 is presented in Figure 2.7. We note especially the lighter regions at l = 174:7 and l = 167:3 , where the original supermosaic analysis shows sharp transitions to lower  . A single strong point source will be included in multiple submaps, thus reducing the spectral index in an area equivalent to a submap. The point-source-subtracted supermosaic shows more gradual transitions in these regions, indicating that the removal is e ective. This behaviour is consistent throughout the supermosaics. 2.4 Maps of Spectral Indices The maps of  determined from model- tting to the supermosaic data are presented in Figures 2.8, 2.9, 2.10, and 2.11. The supermosaics are attached for comparison. Some information at the top and bottom edges of the supermosaics is lost due to the size of the submaps at 480 pixels. Smaller submaps would allow a more thorough sampling toward the edges of the maps, but the smaller size would also limit the multipole ranges in each power spectrum. In the maps of spectral indices, darker pixels indicate steeper power spectra, and thus excess power on large scales. Lighter pix- els indicate regions with signi cantly more power on smaller scales. White pixels denote a failure of the model  tting procedure to converge to a reason- able solution. Most  tting failures are attributable to strong sources with residual imaging artifacts. Further discussion of distinct features in these  maps and their overall statistics is given in Section 3.1. 2.5 Linear Fitting to Flat Power The extent of the  attening behaviour varies across the supermosaics, with some submap power spectra showing no  attening and others showing  at power across a wide range of multipoles. We parameterize the extent and  atness of the  attening e ect with the linear slope in log-log space of the power spectra from ‘ = 1000 to ‘ = 3000. Linear-least-absolute- 142.6. Modeling Flat Power Contributions deviation  tting is applied to this multipole range. The best- t slope of this range is de ned as the  atness parameter  . Submaps which show an increased  attening will have a smaller  tted slope in this region, thus a smaller  atness parameter  indicates a greater  attening e ect. 2.5.1 Maps of Flattening The same submap regions used for creating the spectral index maps are also used to create the  atness parameter maps, allowing us to sample the  atness of the power spectra across the galactic plane observed by the CGPS. The results of this are presented in Figures 2.12 and 2.13. As in the power spectral index maps, information at the edges is lost due to the submap size. As the  attening occurs at multipoles sampled exclusively by the Synthesis Telescope it would be possible to use smaller submaps and still accurately measure its e ects. This would allow us to reduce the information lost at the edges while also increasing the resolution of the sampling. This is not attempted here for consistency with the spectral index maps. In the maps of  atness, data with  > 2 indicates a region with a non-  attened power spectrum; for 1 <  < 2, the power spectrum visually shows some  attening; and for  < 1 the  attening is obvious and signi cant. The maps are colourized with white representing the break from non- attening to  attening; red colours indicate little  attening, and blue colours indicate signi cant  attening. Areas where the model  tting has failed are also white as in the spectral index maps. 2.6 Modeling Flat Power Contributions The  attening observed in the angular power spectra is unexplained though not fully unexpected. Haverkorn et al. [2003] observes the  attening in power spectra of Stokes Q maps, and notes that it is noise, truncating the  tting at lower multipoles to avoid including its e ects in the spectral index calcuations. While its appearance in the present data is expected, its extent is both unexpected and problematic. Prevalent for multipoles of ‘ > 103, the  attening dominates the synthesis multipoles which range from ‘  400 to ‘  8000. A power spectrum from a  eld observed with the synthesis telescope alone is shown in Figure 2.15. No data from the single-dish observations has been added, and the  attening is apparent out to the highest multipoles because the taper has not been applied. The extent of the  at portion of the 152.6. Modeling Flat Power Contributions power spectrum is seen to cover most of the signal sampled by the Synthesis Telescope. The  eld plotted in Figure 2.15 is located at l  90, b  +1, an area with signi cant  attening according to the  maps. Noise is present in the data, but it should not overwhelm the signal to the extent demonstrated by the  attening e ects. There are other possible sources of  at power spectra that may have some signi cance, and these are discussed alongside the e ects of noise below. 2.6.1 Noise Contributions to APS The stated sensitivity of the synthesis observations is 76 mK, from which the theoretical noise  oor can be determined. Creating a noise map and taking its power spectrum is discussed in Section 2.7, analytical calculations of the noise power and an observational measure of the noise are given here. The expected contribution of the noise power to the overall power spec- trum can be calculated from the stated sensitivities of 76 mK in the synthesis data. The noise power depends on this sensitivity and the number of beams in the observed  eld: PSnoise = 2 (0:076)2 Nbeams (2.3) The number of beams is approximated by a ratio of the pixels in the map and the number of pixels in a single synthesis beam on that map: Nbeams = w 2 px=   2beam 8 ln(2) (2.4) For a regular submap of 480 pixels to a side (wpx = 480), the noise power is calculated to be 2:5  10 8K2. This noise power is below the  attening observed in the power spectra. An observational measure of the noise is also available as long as the  eld contains few imaging artifacts. The Synthesis Telescope at DRAO does not determine Stokes V directly, but it does create Stokes V maps during the data processing of Stokes Q and U maps that remain empty of real signal. However, the Stokes V maps do contain the noise and imaging artifacts present in Q and U . Calculating the power spectrum of Stokes V of a given Synthesis Telescope observation should give a measure of the noise power in that observation. Such a power spectrum is shown in Figure 2.16 for the observed  eld K3. The majority of the power spectrum is below a value of 10 8, which when compared to the sample power spectrum in Figure 2.1 is found to be approximately an order of magnitude below the power at the  attening. The power spectrum is in fairly good agreement with the 162.6. Modeling Flat Power Contributions above calculated value for the noise power. Noise does not appear to contain enough power to be responsible for the  attening e ect. 2.6.2 Point Source Contributions to APS To further investigate the  attening of the power spectra we created var- ious model maps with point source emission and determined their angular power spectra. Compact unresolved extragalactic radio point sources that were not subtracted by the removal algorithm are an obvious candidate for providing a  at power spectrum due to the nature of their Fourier trans- forms. We consider signal in a single image pixel to be a two dimensional Dirac delta function, and calculate the Fourier transform of a delta function as follows: Ff (x)g = Z 1  1  (x) exp( {2 ux)dx = 1: (2.5) Since the power spectrum is proportional to the square of the Fourier trans- form we conclude that the expected power spectrum will have equal power on all scales. Actual radio observations of point sources do not generate a perfect Dirac delta function in a single pixel, instead smearing the signal across the beam of the telescope. As the  attening occurs entirely within the range of multipoles sampled by the Synthesis Telescope its number of pixels in a beam is determined: Npx = 2  beam 2 8 ln(2)%2 (2.6) The half-power beam width HPBW for the synthesis telescope is approxi- mately 10, or 0:0167 , and % is the image scale in degrees per pixel. The supermosaics have an image scale of 0:00556 /px, so the synthesis telescope beam will cover 10:2 pixels. An observed point source can therefore be approximated with a Gaussian of the given HPBW and pixel coverage. Placing a single source into an otherwise blank map of the same size as used for the submaps in Section 2.1 produces the power spectrum seen in Figure 2.17. The non-delta function behaviour is evident at high multipoles where the Gaussian appears while most of the power spectrum shows the  atness expected from a point source. The  ux density of the point source was set to 0.4 Jy in order to match the observed powers in the  attening. This is equivalent to a peak surface brightness temperature of 6.7K. Polar- ized intensity emission as seen in Figures 5 to 7 in Landecker et al. [2010] is consistently below this brightness temperature, but the sum of multiple 172.7. Investigating the Non-Gaussian Taper point sources in the image may provide su cient  ux to reach such a tem- perature. As expected, the power spectrum is  at at the correct multipoles, and with a  ux density of 0.4 Jy, also at approximately the correct power. 2.7 Investigating the Non-Gaussian Taper The model- tting diverges from the observed power spectra in the taper region. The taper applied to the synthesis data is a Gaussian that reaches 20% of its peak at the longest baseline of 617m, corresponding to a multipole of ‘ = 7700. The other two telescopes contribute zero power at such a high multipole, thus having no e ect on the taper or the accuracy of the  t. As seen in Figure 2.2, at the beginning of the taper the Gaussian  t has de cient power relative to the observed tapering. At the end of the taper, the highest multipoles, the Gaussian  t has excess power relative to the observed tapering. The Fourier transform of a Gaussian is also a Gaussian in the frequency domain, so the  tting should be accurate. A perfectly Gaussian beam is  t well by a Gaussian model as seen in the model point source  tting in Figure 2.17, though this is an idealized situa- tion. The actual beam changes shape as the declination of the observation increases; the synthesis beam has an actual resolution of 58"  58" csc  , where  is the declination of the observed  eld. An elliptical trend in the supermosaics such as the e ects of the synthesis beam losing circularity will create a non-circular Fourier transformed image as well. The averaging in radial bins would no longer correctly sample the Fourier transformed image in such a case, and this may be the cause of the apparent non-Gaussian taper. Figure 2.18 is the same point source as in Figure 2.17; however, the point spread function of the beam has been altered to an ellipse, with the major axis at double the minor axis of the ellipse, and the major axis ro- tated to 45 from the x axis. The model taper fails to correctly  t the point source’s power spectrum, so this changing ellipticity across the supermosaics is the probable cause of the non-Gaussian tapering. 182.7. Investigating the Non-Gaussian Taper Figure 2.2: Sample power spectra of two 480x480 pixel submaps located at l = 150:0 , b =  1:3 and b = 0:0 respectively. The model of best  t is overplotted in red, and its parameters are given as well. The model consistently diverges from the data at the highest multipoles, where the taper appears ill- t by the Gaussian model. 192.7. Investigating the Non-Gaussian Taper Figure 2.3: Sample power spectra of two 480x480 pixel submaps located at l = 150:0 , b = 1:3 and b = 2:7 respectively. The model of best  t is overplotted in red, and its parameters are given as well. The model consistently diverges from the data at the highest multipoles, where the taper appears ill- t by the Gaussian model. 202.7. Investigating the Non-Gaussian Taper Figure 2.4: Angular power spectra and best- t model with parameters for a submap centered at l = 174:7 , b =  1:7 . The submap of this region contains signi cant power in a point source. The top power spectrum is computed without point source removal, while the bottom power spectrum has undergone point source removal. Of note are the varying model param- eters and the overall steepening of the power spectrum after point source subtraction. 212.7. Investigating the Non-Gaussian Taper Figure 2.5: Angular power spectra and best- t model with parameters for a submap centered at l = 167:3 , b =  1:3 . The submap of this region contains signi cant power in a point source. The top power spectrum is computed without point source removal, while the bottom power spectrum has undergone point source removal. The model  t is superior in the latter  gure. 222.7. Investigating the Non-Gaussian Taper Figure 2.6: Angular power spectra and best- t model with parameters for a submap centered at l = 156:7 , b =  1:7 . The submap of this region contains little power in point sources. The top power spectrum is computed without point source removal, while the bottom power spectrum has un- dergone point source removal. No signi cant changes occur in the model parameters. 232.7. Investigating the Non-Gaussian Taper Figure 2.7: Power spectral index ( ) maps created from the APS analysis of supermosaic 6. The top  gure has not undergone point source removal, whereas the bottom  gure has undergone point source subtraction. Evident in the top map are the sharp transitions in  due to the presence of bright point sources which dominate the power in the di use emission at large multipoles. 242.7. Investigating the Non-Gaussian Taper Figure 2.8: Maps of power spectral index values ( ) in colour and polarized intensity in black and white for l = 66 to l = 104 ; all  values from model  tting to power spectra. 252.7. Investigating the Non-Gaussian Taper Figure 2.9: Maps of power spectral index values ( ) in colour and polarized intensity in black and white for l = 102 to l = 140 . 262.7. Investigating the Non-Gaussian Taper Figure 2.10: Maps of power spectral index values ( ) in colour and polarized intensity in black and white for l = 138 to l = 176 . 272.7. Investigating the Non-Gaussian Taper Figure 2.11: Power spectral-index ( ) map for l = 101 to l = 116 , b = 3 to b = 17 . Polarized intensity map of the same region in black and white. This is the high-latitude extension of the CGPS. Note that the colourbar di ers from those of the previous spectral index maps due to a greater range in values in this region. A clear break in  is evident at b = 10 . 282.7. Investigating the Non-Gaussian Taper Figure 2.12: Flatness parameter ( ) maps for l = 67 to l = 121 . Higher values (red) indicate less  attening, while lower values (blue) indicate sig- ni cant  attening. 292.7. Investigating the Non-Gaussian Taper Figure 2.13: Flatness parameter ( ) maps for l = 121 to l = 175 . Higher values (red) indicate less  attening, while lower values (blue) indicate sig- ni cant  attening. 302.7. Investigating the Non-Gaussian Taper Figure 2.14: Flattening ( ) map for l = 101 to l = 116 , b = 3 to b = 17 , the high latitude extension. As seen in the  map of the same region, a sharp transition is evident around b = 10 . 312.7. Investigating the Non-Gaussian Taper Figure 2.15: Power spectrum of a single synthesis  eld (K3) without added single-dish information. Of note is the  attening of the power spectrum for ‘ > 1000; this  attening is also evident in the APS of the combined dataset. Figure 2.16: Power spectrum for a submap of the Synthesis Telescope ob- servation of  eld K3 in Stokes V . This power spectrum represents the noise power present in this  eld, and by extension an approximation of the noise present in the supermosaics at synthesis multipoles. 322.7. Investigating the Non-Gaussian Taper Figure 2.17: Angular power spectrum of a model point source with 0.4Jy of  ux density. Map containing point source is 480 pixels to a side, with an image scale equal to that of the supermosaics at 0:00556 /px. Figure 2.18: Angular power spectrum of a model point source with 0.4Jy of  ux density. The observational e ect of the Synthesis beam is approximated by rotating and elongating the point source into an ellipse. The model does not accurately  t the taper region for this elliptical beam shape. 33Chapter 3 Discussion and Conclusion 3.1 Spectral Indices and Flattening across the CGPS 3.1.1 Model Failures The model  tting procedure is seen to fail in four distinct major regions which generally correspond to extremely bright point sources in the polarized intensity maps of Landecker et al. [2010]. In order of increasing galactic longitude, these regions are outlined below. The  rst, visible at l  69 , b   1 in Figure 2.8 does not correspond to a single bright source, but a collection of very bright non-point-source emission features. This model failing is unique in this, as all others are clearly attributable to a single bright source as indicated by their square (submap-sized) blank regions. The second model  t failure occurs at l  120 , b  0:4 in both Figures 2.8 and 2.9 due to its proximity to the edges of the relevant supermosaics. This appears to be caused by strong polarized emission from Tycho’s su- pernova remnant. Visible in the polarized intensity maps as a very strong source of emission (T  0:6K), the peak emission appears larger than a sin- gle beam of the Synthesis telescope and would thus escape the point-source removal algorithm. The addition of large amounts of power to the APS then causes the model  tting failure. The third model  t failure occurs at l  131 , b  3:1 . The center of this failure region corresponds to object 3C58, a strong pulsar with an associated synchrotron nebula. The pulsar is not clearly visible in the polar- ized intensity map due to its location in a region of generally high polarized emission, but there are rings of image artifacts surrounding its location. It is more clearly evident in the total intensity map of the region. In this case the artifacts present in the polarization map would cause a model failure by introducing signi cant power on the size scales of the rings. The fourth failure of the model to converge to a solution occurs at l  160 , b  0:0 . The radio galaxy 3C129 is a strong extragalactic source near this location, and it appears in the total intensity map as an extended 343.1. Spectral Indices and Flattening across the CGPS object. Though a strong point-like object that signi cantly alters the power spectra, 3C129 was not removed by the point source removal algorithm due to its extended nature in the map. 3.1.2 Spectral Index and Flattening Features The regular component of the Galaxy’s magnetic  eld is in the direction along l = 85 . This region is expected to have higher than average Faraday rotation and its resultant depolarization e ects. However, the region also contains large-scale features that ensure high values of  . As the spectral index is measuring the balance of power between large and small scales, an increase in power in the small scale can be balanced by an increase in the large scales. The region around l = 85 does however show an increase in  attening compared to its surroundings. The Cygnus region peaks at l = 80 , and is visible in both the  and  maps. The spectral index increases slightly in the region, and is accompanied by a decrease in the  attening parameter. The Cygnus region is located in the middle of the longitudes that cover local features, which stretch from l = 66 to l = 100 [Landecker et al., 2010]. This region of local features has signi cant variation in both  and  , but is punctuated by a few distinct objects that alter the power spectra. The supernova remnant CTB80 at l = 69 , b = 2:7 causes a strong decrease in  attening similar to the decrease seen around Cygnus; it also signi cantly steepens the power spectra as seen in the  map of the region. The supernova remnant CTB104A at l = 94 , b =  0:2 is not evident in either the  or  maps, so the presence of a supernova remnant does not guarantee an increase in  nor a decrese in  . However, CTB104A has visually less structure in the polarized intensity maps than CTB80, so its lack of e ect is expected. The power spectra around Cassiopeia A at l = 112 , b =  2 are very steep, with some of the highest  values in any of the maps. Unfortunately these power spectra are not indicative of the di use emission in the region due to the removal of Cassiopeia A in the polarized intensity maps and the presence of signi cant imaging artifacts [Landecker et al., 2010]. The superbubble GSH166-01-17 identi ed in Landecker et al. [2010] is not distinct from its surroundings in either the  or  maps. In general it is di cult to attribute changes in  or  to speci c objects in the polarization maps. This is due to the ambiguity in whether the observed polarization features are due to intrinsically polarized emission in the background or due to the depolarization and vector-averaging as this emission passes through 353.1. Spectral Indices and Flattening across the CGPS foreground features. Examining the overall statistics and broader galactic trends of the results is more informative than attempting to identify indi- vidual features, as is often the case with polarization images due to their lack of resemblance to the total intensity maps in which objects are easily identi ed. 3.1.3 Spectral Index Statistics The histogram in Figure 3.1 depicts the spread of spectral index values calculated from the Stokes PI supermosaics. Figures 3.2 and 3.3 give the histograms of the  values calculated from the Stokes Q and U maps re- spectively. These histograms do not contain values from the high latitude extension. The Stokes PI histogram peak is at a spectral index of  = 2:3; Stokes Q derived values peak at  = 2:2; Stokes U derived values peak at  = 2:5. Stokes Q and U results are presented for comparison, but all dis- cussion of spectral index features and large-scale trends are for Stokes PI derived spectral index value maps. Figure 3.1: Histogram of power spectral-index values derived from Stokes PI supermosaics. The result for the mode of  is similar to that observed in the Horologium and Auriga regions by Haverkorn et al. [2003], where the spectral indices are found to be 2:34  0:31 and 2:20  0:24 respectively. The Haverkorn et al. [2003] values come from polarization measurements at 350 MHz with a synthesis telescope su ering from missing short spacings, and at a much lower frequency where Faraday rotation e ects are expected to be more 363.1. Spectral Indices and Flattening across the CGPS Figure 3.2: Histogram of power spectral-index values derived from Stokes Q supermosaics. signi cant; the  elds are also at galactic latitudes of 7 and 16 respectively, so direct comparisons may not be accurate. Observations around 2.4 GHz by Giardino et al. [2001] with a single- dish telescope able to sample up to multipole ‘ = 1000 give similar spectral indices as well, with a mean value of  = 2:43  0:05. Similar latitudes are covered as in the present study, but the higher frequency is expected to be less a ected by Faraday rotation, resulting in steeper power spectra as more power remains at larger scales. A later set of observations at 2.4 GHz by Carretti et al. [2010] over galactic latitudes from the South Galactic Pole to the galactic meridian found similarly steep angular power spectra at latitudes of b > 40 . However, for b < 20 , the regions most applicable to the results presented in this thesis, Carretti et al. [2010]  nd an average spectral index of  = 1:8. The above comparisons su er from being at di erent frequencies where the e ects of Faraday rotation may be di erent. Observations at 1.4 GHz have been less complete at sampling galactic latitudes, but some studies exist. Carretti et al. [2005] studies regions at b = 63 and  40 . Their results are lower than those here, with average spectral indices of   1:7 and 2:0 respectively. Carretti et al. [2010] notes that there is a transition from the galactic disk to the halo between b = 20 and 40 , so these spectral indices are sampling visibly di erent polarization skies. Directly comparable to the results presented in this thesis is the analysis 373.1. Spectral Indices and Flattening across the CGPS Figure 3.3: Histogram of power spectral-index values derived from Stokes U supermosaics. by Tucci et al. [2002] of the Southern Galactic Plane Survey test  eld of 28 deg2 located on the galactic plane, observed at 1.4 GHz. They determine a spectral index of   1:7, lower than the mode of the current results but not beyond the bounds of the histogram. Similar regions exist in the data analyzed in this thesis, so this lower value for  may not be representative of all southern galactic plane regions. La Porta and Burigana [2006] use single-dish data to determine spectral indices at 1411 MHz. They report spectral indices of   2  3 for a broad latitude coverage. It should also be noted that the point-source removal is currently unique to this thesis, the e ect of which would be an increase of  in regions where point sources were previously signi cant. The high latitude region has two populations of  values, with a clear break at a galactic latitude of b = 10 . Figure 3.4 shows the histogram of the  values in the high latitude region. For b > 10 , the most common spectral index is  = 3:8. For latitudes closer to the galactic meridian the spectral indices are considerably lower with a mode of  = 1:9. Landecker et al. [2010] notes that the polarization map of the high latitude extension has a visual separation between b = 8 and b = 10 ; small-scale structure is visible for b < 10 , and is absent for b > 10 . This visual separation is analytically veri ed by our results for  which suggest the presence of local features for b < 10 and a regular component of the magnetic  eld that remains present above this latitude. The latitude extension results are 383.1. Spectral Indices and Flattening across the CGPS comparable to those of Carretti et al. [2010], who  nd a transition from disk to halo in the region 20 < b < 40 . Both our results and those of Carretti et al. [2010] are limited by the ranges of galactic longitudes covered. Local features not representative of the overall galactic trend may be present in either set of observations. A full observation of the galactic plane at larger galactic latitudes would indicate where the true transition lies. Covariances between the model parameters are also determined in each  t, and their averages across all calculated power spectra are given here. The strongest covariance is between  and B, at an average of 0.995.  and A have an average covariance of -0.927, while the other covariances are insigni cant and vary between 0.5 and 0.7. Figure 3.4: Histogram of  values from the high latitude extension of the CGPS in Stokes I. Two peaks are visible due to the sharp increase in  at b = 10 . The values determined in this investigation are within the accepted ranges of spectral index values in the literature. Though there is no clear accepted value for the spectral indices at 1.4 GHz for low galactic latitudes, it appears that the presented values are on the steeper end when compared to higher latitude regions. It should be noted that there is no consensus on the appropriate multipole range to be used to calculate spectral index values, and most authors arbitrarily choose the range such that non-linear portions of the power spectra are entirely excluded. Our work has attempted to include the full range of the power spectra by utilizing a model of best  t, and thus our results are expected to di er from those in the literature. 393.1. Spectral Indices and Flattening across the CGPS 3.1.4 Spectral Index and Galactic Latitude Figure 3.5 shows how the spectral index changes with galactic latitude in steps of b = 0:1 . Median values for the spectral index are given at each latitude to avoid outlier e ects of bad model  ts. The spectral index values trend upwards as galactic latitude increases. This indicates a steepening of power spectra as observations move away from the galactic equator; ob- servationally this describes a shifting of power from smaller to larger scales in the polarization maps. This trend may vary signi cantly with speci c subsections of longitudes due to locally varying features; the average over the entirety of the longitudinal range nonetheless appears signi cant. Figure 3.5: Power spectral index  plotted versus galactic latitude b to view averaged trends across the visible galactic disk. Carretti et al. [2010] and Haverkorn et al. [2003] are observations with a broad enough range of latitudes to determine a trend. Carretti et al. [2010]  nds a general increase of the spectral index  from b = 0 to b   60 , and a decrease of the spectral index from b   60 to b   90 . Haverkorn et al. [2003]  nds a decrease in the spectral index with increasing galactic latitude up to b = 25 , however the latitudes below b = 10 are not well sampled, and there is no change from b = 5 to 10 . The results of the present investigation do not disagree with either of these previous studies, though they themselves disagree for b = 10 to 25 . As the spectral indices seem to respond strongly to individual features that may be physically close, trends with galactic latitude may be overwhelmed easily by local features. If local features dominate then our investigation may be the most accurate 403.1. Spectral Indices and Flattening across the CGPS measure of the broader latitudinal trend due to the unprecedented range of longitudinal coverage. 3.1.5 Flattening Statistics A histogram of the  values is included in Figure 3.6, and is seen to peak at  = 1:3. The histogram also indicates that the majority of the power spectra analyzed did not have a signi cant  attening trend. Figure 3.7 is the histogram of  values for the high latitude extension; similar to the  histogram of this region, this histogram appears to contain two distinct populations, one with very signi cant  attening and one with generally in- signi cant  attening. As seen previously in Figure 2.14 the  attening trend is signi cant for b > 10 . The obvious anticorrelation of  and  in the high latitude region is not supported by the results of the supermosaic analyses. 3.1.6 Flattening and Galactic Latitude Figure 3.8 depicts the  attening trend as a function of galactic latitude for the supermosaics. Median values were found in 0:1 latitude intervals. As an increase in  atness parameter  signi es a decrease in the amount of  attening, it appears that the strength of the  attening e ect increases as observations move away from the galactic equator. The sudden drop in  at b  1:9 comes shortly after a sudden increase in  at b  1:6 . Some correlation is expected, but the signi cance of this transition is unknown. No similar measurements exist in the literature, so these conclusions are presented without veri cation or support. If the  attening is due to remaining point sources in the maps we would not expect a dependence on galactic latitude, as the point sources are extra- galactic and should be uniformly distributed across the sky. If the  attening is due to noise in the maps, there is no reason to expect an increase with galactic latiude either, so the observed trend does not support a speci c cause of the  attening. Both the noise and background point sources would produce a noise-like power spectrum at a constant power as discussed in Section 2.6.2. The total power in polarized emission increases towards the galactic equator, so the constant power from noise or point sources would have proportionally less e ect at lower latitudes, possibly falling below the actual power of the Faraday-rotated structure in the maps. This could ex- plain the latitude dependence of the  attening while allowing it to keep a constant power due to either noise or remaining point sources. Unfortu- nately little correlation was found between the model parameter A and the 413.2. Conclusion Figure 3.6: Histogram of  values derived from the Stokes I supermosaics.  attening parameter  . 3.2 Conclusion Power spectral indices were calculated for the entirety of the CGPS po- larization dataset presented by Landecker et al. [2010]. Unparalleled in its galactic coverage and resolution, the survey is truly set apart by its inclu- sion of all angular scales. By combining single-dish and synthesis telescope data the CGPS can accurately measure polarization that would otherwise be missing the largest of structures. The calculated spectral indices con- stitute the largest such analysis, and the  rst mapping of spectral indices. The results presented here determine a dependence of the spectral indices on galactic latitude, with a shift of power to larger scales as observations move away from the galactic equator. The high latitude extension of the CGPS clearly shows a transition in power from smaller scales to larger scales at the galactic latitude of b = 10 . Instrumental and observational e ects of the data were also examined, such as the non-Gaussian taper e ect that may need to be corrected for future studies of the CGPS maps. A  attening of power spectra was also observed for l > 1000, covering most of the multipoles sampled by the Synthesis Telescope. Investigation of this  attening indicates an increase in  attening as galactic latitude increases, but its cause remains uncertain. It is unlikely to be noise or extragalactic point sources if the latitude dependence 423.2. Conclusion Figure 3.7: Histogram of  values derived from the Stokes I map of the high latitude region. is valid, and quantifying the e ects of Faraday depolarization in causing the  attening are di cult. The polarized sky still holds much to be discovered, and a full classi-  cation of the features in the polarized sky may soon be possible as large surveys such as the CGPS reveal polarization structure on all scales. De- coupling small-scale magnetic  eld variances and electron densities from the regular component of the magnetic  eld remains di cult at 1.4 GHz, and is not attempted in this thesis; nonetheless, signi cant statistical results regarding the polarized sky have been determined and presented. 433.2. Conclusion Figure 3.8: Flatness parameter  plotted versus galactic latitude b. 44Bibliography R. Beck. Galactic and extragalactic magnetic  elds - a concise review. As- trophysics and Space Sciences Transactions, 5:43{47, October 2009. Bernard F. Burke. An Introduction to Radio Astronomy. Cambridge Uni- versity Press, New York, New York, 3rd edition, 2010. E. Carretti, G. Bernardi, R. J. Sault, S. Cortiglioni, and S. Poppi. High Galactic latitude polarized emission at 1.4 GHz and implications for cosmic microwave background observations. Monthly Notices of the Royal Astronomical Society, 358:1{12, March 2005. doi: 10.1111/j.1365- 2966.2005.08761.x. E. Carretti, M. Haverkorn, D. McConnell, G. Bernardi, N. M. McClure- Gri ths, S. Cortiglioni, and S. Poppi. The Parkes Galactic Meridian Survey: observations and CMB polarization foreground analysis. Monthly Notices of the Royal Astronomical Society, 405:1670{1689, July 2010. doi: 10.1111/j.1365-2966.2010.16608.x. J.J. Condon and S.M. Ransom. Nrao essential radio astronomy course, 04 2011. URL http://www.cv.nrao.edu/course/astr534/ERA.shtml. G. Giardino, A. J. Banday, P. Fosalba, K. M. G orski, J. L. Jonas, W. O’Mullane, and J. Tauber. The angular power spectrum of radio emission at 2.3 GHz. Astronomy and Astrophysics, 371:708{717, May 2001. doi: 10.1051/0004-6361:20010384. G. Giardino, A. J. Banday, K. M. G orski, K. Bennett, J. L. Jonas, and J. Tauber. Towards a model of full-sky Galactic synchrotron intensity and linear polarisation: A re-analysis of the Parkes data. Astronomy and Astrophysics, 387:82{97, May 2002. doi: 10.1051/0004-6361:20020285. M. Haverkorn, P. Katgert, and A. G. de Bruyn. Characteristics of the structure in the Galactic polarized radio background at 350 MHz. As- tronomy and Astrophysics, 403:1045{1057, June 2003. doi: 10.1051/0004- 6361:20030464. 45Bibliography M. Haverkorn, P. Katgert, and A. G. de Bruyn. Structure in the polarized Galactic synchrotron emission, in particular \depolarization canals". Astronomy and Astrophysics, 427:549{559, November 2004. doi: 10.1051/0004-6361:200400051. R. Kothes and T. L. Landecker. Polarization Structures in the Canadian Galactic Plane Survey: A Multispectral View. In B. Uyaniker, W. Reich, & R. Wielebinski, editor, The Magnetized Interstellar Medium, pages 33{ 38, February 2004. J. D. Kraus. Radio Astronomy. Cygnus-Quasar Books, 2nd 1986. L. La Porta and C. Burigana. A multifrequency angular power spectrum analysis of the Leiden polarization surveys. Astronomy and Astrophysics, 457:1{14, October 2006. doi: 10.1051/0004-6361:20054321. T. L. Landecker. The Role of Magnetic Fields in the Interstellar Medium of the Milky Way. , pages 80{+, July 2011. T. L. Landecker, P. E. Dewdney, T. A. Burgess, A. D. Gray, L. A. Higgs, A. P. Ho mann, G. J. Hovey, D. R. Karpa, J. D. Lacey, N. Prowse, C. R. Purton, R. S. Roger, A. G. Willis, W. Wyslouzil, D. Routledge, and J. F. Vaneldik. The synthesis telescope at the Dominion Radio Astrophysi- cal Observatory. Astronomy and Astrophysics Supplement, 145:509{524, September 2000. doi: 10.1051/aas:2000257. T. L. Landecker, W. Reich, R. I. Reid, P. Reich, M. Wolleben, R. Kothes, B. Uyaniker, A. D. Gray, D. Del Rizzo, E. Furst, A. R. Taylor, and R. Wielebinski. A Survey of the Polarized Emission from the Galactic Plane at 1420 MHz with Arcminute Angular Resolution. ArXiv e-prints, April 2010. W. Reich, E. F urst, P. Reich, B. Uyaniker, R. Wielebinski, and M. Wolleben. The E elsberg 1.4 GHz Medium Galactic Latitude Survey (EMLS). In B. Uyaniker, W. Reich, & R. Wielebinski, editor, The Magnetized Inter- stellar Medium, pages 45{50, February 2004. K. Rohlfs and T.L. Wilson. Tools of Radio Astronomy. Springer-Verlag Telos, 2nd edition, 1996. George B. Rybicki and Alan P. Lightman. Radiative Processes in Astro- physics. Wiley-VCH, 1985. 46Bibliography D. D. Sokolo , A. A. Bykov, A. Shukurov, E. M. Berkhuijsen, R. Beck, and A. D. Poezd. Depolarization and Faraday e ects in galaxies. Monthly Notices of the Royal Astronomical Society, 299:189{206, August 1998. doi: 10.1046/j.1365-8711.1998.01782.x. M. Tucci, E. Carretti, S. Cecchini, L. Nicastro, R. Fabbri, B. M. Gaensler, J. M. Dickey, and N. M. McClure-Gri ths. Polarization Angular Spectra of Galactic Synchrotron Emission on Arcminute Scales. The Astrophysical Journal, 579:607{615, November 2002. doi: 10.1086/342793. B. Uyaniker, E. Fuerst, W. Reich, P. Reich, and R. Wielebinski. A 1.4 GHz radio continuum and polarization survey at medium Galactic latitudes. I. Observation and reduction technique. Astronomy and Astrophysics Sup- plement, 132:401{411, November 1998. doi: 10.1051/aas:1998449. Gerrit L. Verschuur. The Invisible Universe. Springer, New York, New York, 2nd edition, 2007. A.G. Willis. Some Synthesis Telescope imaging algorithms to remove non- isoplanatic and other nasty artifacts. Astronomy and Astrophysics Sup- plement Series, 136:603{614, May 1999. doi: 10.1051/aas:1999237. M. Wolleben, T. L. Landecker, W. Reich, and R. Wielebinski. An absolutely calibrated survey of polarized emission from the northern sky at 1.4 GHz. Observations and data reduction. Astronomy and Astrophysics, 448:411{ 424, March 2006. doi: 10.1051/0004-6361:20053851. 47Appendix A Synchrotron Emission As previously noted, the main contribution to the observed polarized structures is provided by Faraday rotation of Galactic synchrotron emis- sion. In fact, synchrotron emission dominates emission at all frequencies below 30 GHz for star-forming galaxies such as our own [Condon and Ran- som, 2011]. This synchrotron emission, or magnetobrehmsstrahlung, occurs when charged particles such as electrons are accelerated in a magnetic  eld. By Larmor’s formula, acceleration of a charged particle will produce elec- tromagnetic radiation. In the case of synchrotron radiation, the charged particles are relativistic electrons that undergo a uniform circular motion with a constant orbital velocity, resulting in a helical path around the mag- netic  eld lines. Applying the relativistic corrections to the gyration frequency expected from such electrons gives a frequency in Gaussian units of !B = eB  mec : (A.1) Here  is the usual Lorentz factor. In order to calculate the power emitted, the Thomson cross section of the electrons for scattering of electromagnetic radiation is determined as  T = 8 3  e2 mec2  2 : (A.2) Finally, the e ects of the magnetic  eld are applied as a magnetic energy density UB = B 2 8 , giving an emitted power per electron from Larmor’s for- mula of P = 2 T 2 2cUB sin 2  : (A.3) Here  = vc , and  is the pitch angle between the direction of motion of the electron and the direction of the magnetic  eld. By averaging over a large number of electrons with varying pitch angles but identical velocities the average power can be determined: hP i = 4 3  T 2 2cUB: (A.4) 48Appendix A. Synchrotron Emission This emitted power is much greater than classically expected from a gyrating electron; it is the relativistic corrections that allow synchrotron emission to be the dominant process at our frequencies of interest [Rybicki and Lightman, 1985]. Aside from the overall power emitted by synchrotron radiation, the de- pendence of this power on frequency is worth considering. The radiation emitted from such a gyration e ect is radiated in the instantaneous direc- tion of motion of the emitting electron. Due to the helical path of the electrons, the aggregate emission is strongly beamed in a cone of emission with an opening angle identical to the pitch angle and its axis pointed along the magnetic  eld. Relativistic e ects cause this cone of emission to narrow as the velocity approaches the speed of light; the emission cone opening is e ectively dictated by the electron energies. Furthermore, the emitted syn- chrotron radiation is only detectable when the emission cone aligns with the observational beam, which occurs very brie y but periodically [Rohlfs and Wilson, 1996]. Each synchrotron pulse is extremely narrow compared to the time be- tween pulses, with a typical pulse duration on the order of 10 10 seconds and inter-pulse spacings on the order of 103 seconds. The Fourier transform of the pulse train gives a broad spectrum - the observed power spectrum of the emission; this transform consists of a series of peaks so closely spaced in the frequency domain that they appear essentially continuous. Furthermore, as noted in Condon and Ransom [2011], any  uctuations in the electron energy or magnetic  eld strength causes frequency shifts in the observed spectrum larger than the separations between these closely spaced peaks. The emission from a single electron appears nearly continuous. Applying a cosmic-ray electron energy distribution as observed in our galaxy smears out the spectrum and predicts an accurate measure of the slope of the syn- chrotron emission in the frequency domain. As calculated in Condon and Ransom [2011], the emission coe cient within our galaxy,   , is   / B 1:7 0:7: (A.5) Of primary importance here is the negative spectral index of 0.7 that de- termines the shape of the radio spectrum below 30 GHz. It should also be noted that the spectrum at frequencies below 1 GHz can be  attened con- siderably by free-free absorption, though this e ect is unimportant at the frequencies considered in this thesis. With a spectral index of 0.7 the  28 MHz spread in centre frequencies for the bands of the Synthesis Telescope (discussed in Section 1.2.3) amounts to a di erence of  4% in amplitude scales. 49Appendix A. Synchrotron Emission Synchrotron emission also has an intrinsic polarization due to the align- ment of the spiraling electrons. The radiation from a single charge is ellip- tically polarized, but a distribution of charges will cancel out this elliptical component, leaving the radiation as partially linearly polarized perpendic- ular to the spiral direction. Comparing the powers emitted parallel and perpendicular to the magnetic  eld (components of hP i above) the maxi- mum fraction of polarization is determined to be between 70 and 75 %. The variation in fractional polarization is due to the slight uncertainty in the power law distribution of electron energies [Rybicki and Lightman, 1985]. As discussed in Section 1.2, the observed polarization fractions are closer to 10%, but this may be due to signi cant depolarization e ects. 50Appendix B Polarization From Maxwell’s equations of electromagnetism we determine that any electromagnetic radiation consists of transverse waves of oscillating electric and magnetic  elds. In free space the wave equation for the propagation of the electric  eld is r2" = 1 c2 @2" @t2 : (B.1) The magnetic  eld obeys a similar equation, but is found perpendicular to the electric  eld with the cross product of the two  elds determining the direction of propagation of the wave. The electric  eld is considered here as it is chosen by convention to represent the polarization of the wave. At a signi cant distance from the source of the radiation the wavefronts are e ectively plane-parallel, and a general solution of the wave equation for the electric  eld takes the following form: "(z; t) = (x^"x + y^"y) exp[i2 ( t kz)]: (B.2) In this solution the electric  eld has an inverse wavelength of k = 1= and a frequency of  . The direction of wave propagation is z, and the preceding multiplicative factor expresses the wave amplitude in terms of separate complex amplitudes "x and "y [Burke, 2010]. The polarization of the wave is thus described by the relative amplitudes of "x and "y and their relative phase   =  1   2. In this formulation the x and y axes specify a plane that lies perpendicular to the direction of propagation. If the orientation of the electric  eld remains constant as the wave propagates, then the wave is said to be linearly polarized. If instead the amplitudes are equal but the phases of "x and "y di er by 90 degrees the wave is circularly polarized, and the combined electric wave vector rotates in the x,y plane. These are the two limiting cases of the general result in which the projection of " on the x,y plane describes an ellipse in time - elliptical polarization. The wave vector may rotate in either direction, so right-hand elliptical polarization is de ned as the vector rotating clockwise when viewed along the direction of propagation; left-hand polarization is the counter-clockwise case [Burke, 2010]. 51B.1. Stokes Parameters In general the axes of the ellipse need not be aligned with the obser- vational x and y axes, being rotated instead by some angle  . For this reason the electric  eld vector may also be described by separating the sin- gle amplitude of "0 into "a and "b, which represent the electric  eld vector components along the major and minor axes of the ellipse respectively: "a = "0 cos ; "b =  "0 sin : (B.3) The major and minor axes of the ellipse are separated by   =2, and the parameter  describes the classi cation of the polarization. For  = 0 the light is plane-polarized, and for  =   =4 the light is circularly polarized. Thus for the plane waves considered here the polarization is fully speci ed by the parameter  , the orientation angle of the ellipse  , and the total amplitude "0 [Kraus, 1986]. B.1 Stokes Parameters The polarization state of an electromagnetic wave can be fully and suc- cinctly represented by the Stokes parameters I, Q, U, and V as de ned by George Gabriel Stokes in 1852 [Burke, 2010]. In terms of the elliptical parameters de ned at the end of Section B, the Stokes parameters for a monochromatic polarized wave are: I  "20; Q  " 2 0 cos(2 ) cos(2 ); U  " 2 0 cos(2 ) sin(2 ); V  " 2 0 sin(2 ) The Stokes parameters are a particularly useful characterization of the po- larization as they are directly related to observable powers. Stokes I is proportional to the total energy  ux; Stokes Q and U are determined from LHCP and RHCP power detection; Stokes V is not determined in the cur- rent dataset. For Stokes V equal to zero the polarization is entirely linear; conversely, for V="20 = +1 the polarization is entirely (right-hand) circular. For fully polarized waves Stokes I is related directly to the polarized Stokes parameters: I2 = Q2 + U2 + V 2 (B.4) As the radiation  eld is described by the three parameters "0,  , and  , the four Stokes parameters completely describe the radiation  eld. In the case of completely unpolarized radiation I is the speci c intensity of the radiation, with Q, U, and V all zero, and B.4 no longer applies. Instead, a degree of polarization PI can be de ned as the ratio of the quadrature sum 52B.2. Faraday Rotation of the polarization parameters Q, U, V and the speci c intensity: PI = (Q2 + U2 + V 2)1=2 I : (B.5) The observations analyzed in this thesis are made in two orthogonal po- larizations, which allows us to specify Q and U, while leaving Stokes V unknown. In this case the degree of linear polarization can be determined by the above equation with V = 0. For the present data the expected de- gree of linear polarization is less than 70 %, the theoretical maximum for synchrotron emission [Sokolo et al., 1998]. Finally, the polarization angle  can be determined by rearranging Stokes Q and U:  = 1 2 arctan  U Q  : (B.6) It should also be noted that in reality we do not measure single monochro- matic waves, but rather a time-average of many superpositions over some  nite bandwidth of the receiver, so the above equations should be consid- ered as time- and bandwidth-averaged values determined from the power measurements of the telescope receiver [Burke, 2010]. B.2 Faraday Rotation Radiation emitted within the Galaxy does not travel to the telescopes on Earth unhindered or una ected by the media it traverses. Interactions with the interstellar medium may cause the radiation to be absorbed, scattered, or otherwise fundamentally altered. Of importance for the observations presented in this thesis, and generally for any radio polarization studies, are the e ects of Faraday rotation. As an electromagnetic wave passes through a plasma the electrons in the plasma respond to the changing electric  eld by oscillating and re-radiating at the wave frequency, giving the plasma a certain refractive index [Burke, 2010]. The presence of a static magnetic  eld forces the moving electrons into curved paths reminiscient of the synchrotron mechanism responsible for the original emission. However, the response of the polarized wave components to this refractive index di ers depending on the direction of polarization. Considering the linearly polarized emission to be a combination of left and right circular polarizations, the e ect causes the circular polarizations to move at di erent speeds through the medium. This relative phase shift of the circular polarizations amounts to a rotation of the linear polarization [Burke, 2010]. 53B.3. Depolarization Mechanisms In essence the plasma through which the electromagnetic wave is prop- agating is birefringent, with a refractive index for the two circular polariza- tions of n2 = 1 !2p !(!  !B) : (B.7) Here the gyration frequency is given by the nonrelativistic form: !B = eB mec : (B.8) As the electrons are required to gyrate perpendicular to the wave direction in order to have a signi cant e ect, the Faraday rotation e ect is dominated by the component of the magnetic  eld aligned with the direction of wave propagation. In the radio regime, the angle of total rotation along a ray path is given by  R:  R = RM 2 (B.9) where RM is the rotation measure of the ray path in rad m 2. For distances in parsecs, density in cm 3, and magnetic  eld in microgauss, this rotation measure is RM = 0:81 Z BN cos  dl: (B.10) In studies of the rotation measure along various sight-lines the value of RM can be determined from observations at multiple frequencies by compar- ing the angles of rotation. As this thesis deals with single-frequency mea- surements alone, no calculations of rotation measure are made. It should be noted that a rough analysis of rotation measure is possible due to the four frequency bands utilized in the synthesis telescope, but the sensitiv- ity required is lacking. Nonetheless, Faraday rotation causes much of the small-scale variation seen in the polarization angle [Landecker et al., 2000]. Faraday rotation may also be responsible for signi cant depolarization as discussed in Section B.3. B.3 Depolarization Mechanisms As this thesis and the observations considered herein are polarization studies, any e ects that can cause erroneous polarization features need to be determined in detail. The various e ects that can cause depolarization or a loss of polarization in the measurements are discussed here. 54B.3. Depolarization Mechanisms An important source of depolarization is di erential Faraday rotation or depth depolarization. Faraday rotation in isolation does not cause de- polarization since it merely changes the polarization angle. However, if a medium contains both thermal and relativistic electrons and some semi- ordered magnetic  eld then the synchrotron emission at varying depths in the medium will undergo di erent amounts of Faraday rotation. With a signi cantly broad distribution of Faraday rotations these components can begin to cancel and e ectively reduce the polarized intensity. This results in a reduced polarization fraction of the observed emission dependent on the rotation measure along the line of sight [Sokolo et al., 1998]. In such circumstances the observed polarized intensity is given by I = I0     sin(2RM 2) 2RM 2     (B.11) [Haverkorn et al., 2004]. Depth depolarization also contains a contribution from variations of the emitted polarization angles along a line of sight; that is, the polarization of the synchrotron emission along the line of sight may vary, and the superposition of these vectors may lead to partial cancellation. The other signi cant source of depolarization occurs across the beam of the receiver, and is termed beam depolarization. This depolarization is sim- ply the result of vector addition of di ering polarizations for di erent lines of sight within a single telescope beam, and has no dependence on wavelength. This e ect strongly contributes to depolarization canals in polarized inten- sity maps [Haverkorn et al., 2004]. Such canals appear as  lamentary struc- ture of one beam width with lower polarized intensity than the surroundings. The polarization angle also changes by 90 across each canal. Haverkorn et al. [2004] argue that these canals exist between two regions which have a polarization angle di erence of   = (n + 1=2)180 (n = 0; 1; 2; :::), which causes the vector addition across the beam to return almost zero polarized intensity. If this rapid transition in polarization angle occurs within a beam then the resulting region of depolarization will clearly have a width of a single beam. The analysis by Landecker et al. [2010] of the CGPS maps used in this thesis indicates that along directions close to a galactic longitude of 85 there is signi cant Faraday rotation and depolarization e ects. The ordered component of the local magnetic  eld is in the direction of ‘ = 85 , allowing for greater Faraday rotation. Shifting the line of sight closer to the anticenter of the Galaxy shows decreased Faraday rotation as the line of sight becomes perpendicular to the regular component of the magnetic  eld. A strong decrease in depolarization is seen with an increase in polarization fraction. 55B.3. Depolarization Mechanisms The synchrotron emission in the direction of the anticenter also decreases as evidenced by the reduction in total intensity. Depolarization is also strong in nearby HII regions where Faraday rota- tion is produced by the abundance of free electrons; this e ect is seen in the polarization images as a sudden drop in polarized intensity at the edges of HII regions. (HII regions can be viewed in total intensity or H maps.) The HII regions depolarize more distant emission, and may have little intrin- sic synchrotron emission on the near side to our observations. Less drastic depolarization occurs when distant emission from the Perseus Arm passes through HII regions in the closer Local Arm. The polarized emission in these directions is fairly smooth, probably due to synchrotron emission on the near side of the HII regions which then undergoes little further Faraday rotation. Further discussion of speci c depolarization regions can be found in Wolleben et al. [2006], Kothes and Landecker [2004], Landecker et al. [2010]. There exists a distance beyond which we are unable to measure polariza- tion due to the combined e ects of di erential Faraday rotation and beam depolarization. Beyond this distance all polarized emission has been e ec- tively depolarized, thus it is termed the polarization horizon. The exact distance depends on the frequency of the measurements, beamwidth, and the direction of the line of sight. Between ‘  66 and ‘  100 the ob- served polarization features are predominantly local - within 2 kpc [Kothes and Landecker, 2004]; at greater longitudes the features are a combination of large, smooth local features and seemingly random small-scale structure from the Perseus arm. Near the galactic anticenter the polarization horizon is beyond the Perseus arm, and possibly beyond the extent of the Galaxy [Landecker et al., 2010]. 56

Cite

Citation Scheme:

        

Citations by CSL (citeproc-js)

Usage Statistics

Country Views Downloads
United States 21 8
China 8 0
Russia 4 0
Germany 2 10
Brazil 1 7
City Views Downloads
Ashburn 16 0
Shenzhen 7 0
Unknown 6 17
Wilmington 2 0
Washington 2 0
Sunnyvale 1 0
Saint Petersburg 1 0
Beijing 1 0

{[{ mDataHeader[type] }]} {[{ month[type] }]} {[{ tData[type] }]}
Download Stats

Share

Embed

Customize your widget with the following options, then copy and paste the code below into the HTML of your page to embed this item in your website.
                        
                            <div id="ubcOpenCollectionsWidgetDisplay">
                            <script id="ubcOpenCollectionsWidget"
                            src="{[{embed.src}]}"
                            data-item="{[{embed.item}]}"
                            data-collection="{[{embed.collection}]}"
                            data-metadata="{[{embed.showMetadata}]}"
                            data-width="{[{embed.width}]}"
                            async >
                            </script>
                            </div>
                        
                    
IIIF logo Our image viewer uses the IIIF 2.0 standard. To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.24.1-0072378/manifest

Comment

Related Items