Faraday rotation occurs along every line of sight in the Galaxy; Rotation Measure (RM) synthesis allows a three-dimensional representation of the interstellar magnetic field. This study uses data from the Global Magneto-Ionic Medium Survey, a combination of single-antenna spectro-polarimetric studies, including northern sky data from the DRAO 26-m Telescope (1270-1750 MHz) and southern sky data from the Parkes 64-m Telescope (300-480 MHz). From the synthesized Faraday spectral cubes we compute the zeroth, first, and second moments to find the total polarized emission, mean and RM-width of the polarized emission. From DRAO first moments we find a weak vertical field directed from Galactic North to South, but Parkes data reveal fields directed towards the Sun at high latitudes in both hemispheres: the two surveys clearly sample different volumes. DRAO second moments show feature widths in Faraday spectra increasing with decreasing positive latitudes, implying that longer lines of sight encounter more Faraday rotating medium, but this is not seen at negative latitudes. Parkes data show the opposite: at positive latitudes the second moment decreases with decreasing latitude, but not at negative latitudes. Comparing first moments with RMs of pulsars and extragalactic sources and a study of depolarization together confirm that the DRAO survey samples to larger distances than the Parkes data. Emission regions in the DRAO survey are typically 700 to 1000 pc away, slightly beyond the scale-height of the magneto-ionic medium; emission detected in the Parkes survey is entirely within the magneto-ionic disk, less than 500 pc away.

The interstellar medium is a complex 'ecosystem' with gas constituents in the atomic, molecular, and ionized states, dust, magnetic fields, and relativistic particles. The Canadian Galactic Plane Survey has imaged these constituents with angular resolution of the order of arcminutes. This paper presents radio continuum data at 408 MHz over the area 52 degrees < longitude < 193 degrees, -6.5 degrees < latitude < 8.5 degrees, with an extension to latitude = 21 degrees in the range 97 degrees < longitude < 120 degrees, with angular resolution 2.8' x 2.8' cosec(declination). Observations were made with the Synthesis Telescope at the Dominion Radio Astrophysical Observatory as part of the Canadian Galactic Plane Survey. The calibration of the survey using existing radio source catalogs is described. The accuracy of 408-MHz flux densities from the data is 6%. Information on large structures has been incorporated into the data using the single-antenna survey of Haslam (1982). The paper presents the data, describes how it can be accessed electronically, and gives examples of applications of the data to ISM research.

The Fan Region is one of the dominant features in the polarized radio sky, long thought to be a local (distance < 500 pc) synchrotron feature. We present 1.3-1.8 GHz polarized radio continuum observations of the region from the Global Magneto-Ionic Medium Survey (GMIMS) and compare them to maps of Halpha and polarized radio continuum intensity from 0.408-353 GHz. The high-frequency (> 1 GHz) and low-frequency (< 600 MHz) emission have different morphologies, suggesting a different physical origin. Portions of the 1.5 GHz Fan Region emission are depolarized by about 30% by ionized gas structures in the Perseus Arm, indicating that this fraction of the emission originates >2 kpc away. We argue for the same conclusion based on the high polarization fraction at 1.5 GHz (about 40%). The Fan Region is offset with respect to the Galactic plane, covering -5{\deg} < b < +10{\deg}; we attribute this offset to the warp in the outer Galaxy. We discuss origins of the polarized emission, including the spiral Galactic magnetic field. This idea is a plausible contributing factor although no model to date readily reproduces all of the observations. We conclude that models of the Galactic magnetic field should account for the > 1 GHz emission from the Fan Region as a Galactic-scale, not purely local, feature.

Measurement of the brightness temperature of extended radio emission demands knowledge of the gain (or aperture efficiency) of the telescope and measurement of the polarized component of the emission requires correction for the conversion of unpolarized emission from sky and ground to apparently polarized signal. Radiation properties of the John A. Galt Telescope at the Dominion Radio Astrophysical Observatory were studied through analysis and measurement in order to provide absolute calibration of a survey of polarized emission from the entire northern sky from 1280 to 1750 MHz, and to understand the polarization performance of the telescope. Electromagnetic simulators CST and GRASP-10 were used to compute radiation patterns of the telescope in all Stokes parameters, and aperture efficiency. Aperture efficiency was also evaluated using geometrical optics and was measured using Cyg A. Measured aperture efficiency varied smoothly with frequency between values of 0.49 and 0.54; GRASP-10 yielded values 6.5% higher but with closely similar variation with frequency. Overall error across the frequency band is 3%, but values at any two frequencies are relatively correct to ~1%. Dominant influences on aperture efficiency are illumination taper of the feed radiation pattern and shadowing by the feed-support struts. A model of ground emission was developed based on measurements and on empirical data from remote sensing of the Earth from satellite-borne telescopes. This model was convolved with the computed antenna response to estimate conversion of ground emission into spurious polarized signal. The computed spurious signal is comparable to measured values, but is not accurate enough to be used to correct observations. A simpler model, in which the ground is considered as an unpolarized emitter with a brightness temperature of ~240 K, is shown to have useful accuracy when compared to measurements.

We present radio continuum and polarization images of the North Polar Spur (NPS) from the Global Magneto-Ionic Medium Survey (GMIMS) conducted with the Dominion Radio Astrophysical Observatory 26-m Telescope. We fit polarization angle versus wavelength squared over 2048 frequency channels from 1280 to 1750 MHz to obtain a Faraday Rotation Measure (RM) map of the NPS. Combining this RM map with a published Faraday depth map of the entire Galaxy in this direction, we derive the Faraday depth introduced by the NPS and the Galactic interstellar medium (ISM) in front of and behind the NPS. The Faraday depth contributed by the NPS is close to zero, indicating that the NPS is an emitting only feature. The Faraday depth caused by the ISM in front of the NPS is consistent with zero at b>50 degree, implying that this part of the NPS is local at a distance of approximately several hundred parsecs. The Faraday depth contributed by the ISM behind the NPS gradually increases with Galactic latitude up to b=44 degree, and decreases at higher Galactic latitudes. This implies that either the part of the NPS at b<44 degree is distant or the NPS is local but there is a sign change of the large-scale magnetic field. If the NPS is local, there is then no evidence for a large-scale anti-symmetry pattern in the Faraday depth of the Milky Way. The Faraday depth introduced by the ISM behind the NPS at latitudes b>50 degree can be explained by including a coherent vertical magnetic field.

Astronomy and astrophysics are witnessing dramatic increases in data volume as detectors, telescopes and computers become ever more powerful. During the last decade, sky surveys across the electromagnetic spectrum have collected hundreds of terabytes of astronomical data for hundreds of millions of sources. Over the next decade, the data volume will enter the petabyte domain, and provide accurate measurements for billions of sources. Astronomy and physics students are not traditionally trained to handle such voluminous and complex data sets. In this paper we describe astroML; an initiative, based on Python and scikit-learn, to develop a compendium of machine learning tools designed to address the statistical needs of the next generation of students and astronomical surveys. We introduce astroML and present a number of example applications that are enabled by this package.

This document specifies a standard format for vocabularies based on the W3C's Resource Description Framework (RDF) and Simple Knowledge Organization System (SKOS). By adopting a standard and simple format, the IVOA will permit different groups to create and maintain their own specialised vocabularies while letting the rest of the astronomical community access, use, and combine them. The use of current, open standards ensures that VO applications will be able to tap into resources of the growing semantic web. The document provides several examples of useful astronomical vocabularies.

The DRAO ST was used for the past 15 years as the primary instrument for the Canadian Galactic Plane Survey. This has been a spectacularly successful project, advancing our understanding of the Milky Way Galaxy through panoramic imaging of the main constituents of the Interstellar Medium. Observations for the CGPS are now complete and the Synthesis Telescope at DRAO has returned to proposal-driven mode. The Dominion Radio Astrophysical Observatory invites astronomers to apply for observing time with the DRAO Synthesis Telescope. The DRAO ST provides radio observations of atomic hydrogen and radio continuum emission, including the polarized signal, with high spatial dynamic range and arcminute resolution. Imaging techniques developed for the CGPS have made the telescope into a front-line instrument for wide-field imaging, particularly of polarized emission. We will discuss telescope characteristics, show examples of data to demonstrate the unique capabilities of the ST, and explain where and how to apply for observing time.

Context: Observations of polarized emission are a significant source of information on the magnetic field that pervades the Interstellar Medium of the Galaxy. Despite the acknowledged importance of the magnetic field in interstellar processes, our knowledge of field configurations on all scales is seriously limited. Aims: This paper describes an extensive survey of polarized Galactic emission at 1.4 GHz that provides data with arcminute resolution and complete coverage of all structures from the broadest angular scales to the resolution limit, giving information on the magneto-ionic medium over a wide range of interstellar environments. Methods: Data from the DRAO Synthesis Telescope, the Effelsberg 100-m Telescope, and the DRAO 26-m Telescope have been combined. Angular resolution is ~1' and the survey extends from l = 66 deg to l = 175 deg over a range -3 deg < b < 5 deg along the northern Galactic plane, with a high-latitude extension from l = 101 deg to l = 116 deg up to b = 17.5 deg. This is the first extensive polarization survey to present aperture-synthesis data combined with data from single antennas, and the techniques developed to achieve this combination are described. Results: The appearance of the extended polarized emission at 1.4 GHz is dominated by Faraday rotation along the propagation path, and the diffuse polarized sky bears little resemblance to the total-intensity sky. There is extensive depolarization, arising from vector averaging on long lines of sight, from HII regions, and from diffuse ionized gas seen in H-alpha images. Preliminary interpretation is presented of selected polarization features on scales from parsecs (the planetary nebula Sh 2-216) to hundreds of parsecs (a supperbubble GSH 166-01-17) to kiloparsecs (polarized emission in the direction of Cygnus X).

We explore the multidimensional, multiwavelength selection of quasars from mid-IR (MIR) plus optical data, specifically from Spitzer-IRAC and the Sloan Digital Sky Survey (SDSS). We apply modern statistical techniques to combined Spitzer MIR and SDSS optical data, allowing up to 8-D color selection of quasars. Using a Bayesian selection method, we catalog 5546 quasar candidates to an 8.0 um depth of 56 uJy over an area of ~24 sq. deg; ~70% of these candidates are not identified by applying the same Bayesian algorithm to 4-color SDSS optical data alone. Our selection recovers 97.7% of known type 1 quasars in this area and greatly improves the effectiveness of identifying 3.5<z<5 quasars. Even using only the two shortest wavelength IRAC bandpasses, it is possible to use our Bayesian techniques to select quasars with 97% completeness and as little as 10% contamination. This sample has a photometric redshift accuracy of 93.6% (Delta Z +/-0.3), remaining roughly constant when the two reddest MIR bands are excluded. While our methods are designed to find type 1 (unobscured) quasars, as many as 1200 of the objects are type 2 (obscured) quasar candidates. Coupling deep optical imaging data with deep mid-IR data could enable selection of quasars in significant numbers past the peak of the quasar luminosity function (QLF) to at least z~4. Such a sample would constrain the shape of the QLF and enable quasar clustering studies over the largest range of redshift and luminosity to date, yielding significant gains in our understanding of quasars and the evolution of galaxies.

We present a catalog of 1,172,157 quasar candidates selected from the photometric imaging data of the Sloan Digital Sky Survey (SDSS). The objects are all point sources to a limiting magnitude of i=21.3 from 8417 sq. deg. of imaging from SDSS Data Release 6 (DR6). This sample extends our previous catalog by using the latest SDSS public release data and probing both UV-excess and high-redshift quasars. While the addition of high-redshift candidates reduces the overall efficiency (quasars:quasar candidates) of the catalog to ~80%, it is expected to contain no fewer than 850,000 bona fide quasars -- ~8 times the number of our previous sample, and ~10 times the size of the largest spectroscopic quasar catalog. Cross-matching between our photometric catalog and spectroscopic quasar catalogs from both the SDSS and 2dF Surveys, yields 88,879 spectroscopically confirmed quasars. For judicious selection of the most robust UV-excess sources (~500,000 objects in all), the efficiency is nearly 97% -- more than sufficient for detailed statistical analyses. The catalog's completeness to type 1 (broad-line) quasars is expected to be no worse than 70%, with most missing objects occurring at z<0.7 and 2.5<z<3.0. In addition to classification information, we provide photometric redshift estimates (typically good to Delta z +/- 0.3 [2 sigma]) and cross-matching with radio, X-ray, and proper motion catalogs. Finally, we consider the catalog's utility for determining the optical luminosity function of quasars and are able to confirm the flattening of the bright-end slope of the quasar luminosity function at z~4 as compared to z~2.

Polarimetric observations are affected by leakage of unpolarized light into the polarization channels, in a way that varies with the angular position of the source relative to the optical axis. The off-axis part of the leakage is often corrected by subtracting from each polarization image the product of the unpolarized map and a leakage map, but it is seldom realized that heterogeneities in the array shift the loci of the leaked radiation in a baseline-dependent fashion. We present here a method to measure and remove the wide-field polarization leakage of a heterogeneous array. The process also maps the complex voltage patterns of each antenna, which can be used to correct all Stokes parameters for imaging errors due to the primary beams.

We present deep polarimetric observations at 1420 MHz of the European Large Area ISO Survey North 1 region (ELAIS N1) as part of the Dominion Radio Astrophysical Observatory Planck Deep Fields project. By combining closely spaced aperture synthesis fields, we image a region of 7.43 square degrees to a maximum sensitivity in Stokes Q and U of 78 microJy/beam, and detect 786 compact sources in Stokes I. Of these, 83 exhibit polarized emission. We find that the differential source counts (log N - log p) for polarized sources are nearly constant down to p > 500 microJy, and that these faint polarized radio sources are more highly polarized than the strong source population. The median fractional polarization is (4.8 +/- 0.7)% for polarized sources with Stokes I flux density between 1 and 30 mJy; approximately three times larger than sources with I > 100 mJy. The majority of the polarized sources have been identified with galaxies in the Spitzer Wide Area Infrared Extragalactic Survey (SWIRE) image of ELAIS N1. Most of the galaxies occupy regions in the IRAC 5.8/3.6 micron vs. 8.0/4.5 micron color-color diagram associated with dusty AGNs, or with ellipticals with an aging stellar population. A few host galaxies have colors that suggests significant PAH emission in the near-infrared. A small fraction, 12%, of the polarized sources are not detected in the SWIRE data. None of the polarized sources in our sample appears to be associated with an actively star-forming galaxy.

We present measurements of the redshift-space three-point correlation function of 50,967 Luminous Red Galaxies (LRGs) from Data Release 3 (DR3) of the Sloan Digital Sky Survey (SDSS). We have studied the shape dependence of the reduced three-point correlation function (Qz(s,q,theta)) on three different scales, s=4, 7 and 10 h-1 Mpc, and over the range of 1 < q < 3 and 0 < theta < 180. On small scales (s=4 h-1 Mpc), Qz is nearly constant, with little change as a function of q and theta. However, there is evidence for a shallow U-shaped behaviour (with theta) which is expected from theoretical modeling of Qz . On larger scales (s=7 and 10 h-1 Mpc), the U-shaped anisotropy in Qz (with theta) is more clearly detected. We compare this shape-dependence in Qz(s,q,theta) with that seen in mock galaxy catalogues which were generated by populating the dark matter halos in large N-body simulations with mock galaxies using various Halo Occupation Distributions (HOD). We find that the combination of the observed number density of LRGs, the (redshift-space) two-point correlation function and Qz provides a strong constraint on the allowed HOD parameters (M_min, M_1, alpha) and breaks key degeneracies between these parameters. For example, our observed Qz disfavors mock catalogues that overpopulate massive dark matter halos with many LRG satellites. We also estimate the linear bias of LRGs to be b=1.87+/-0.07 in excellent agreement with other measurements.

We present evidence of a large angle correlation between the cosmic microwave background measured by WMAP and a catalog of photometrically detected quasars from the SDSS. The observed cross correlation is (0.30 +- 0.14) microK at zero lag, with a shape consistent with that expected for correlations arising from the integrated Sachs-Wolfe effect. The photometric redshifts of the quasars are centered at z ~ 1.5, making this the deepest survey in which such a correlation has been observed. Assuming this correlation is due to the ISW effect, this constitutes the earliest evidence yet for dark energy and it can be used to constrain exotic dark energy models.

We present measurements of the normalised redshift-space three-point correlation function (Q_z) of galaxies from the Sloan Digital Sky Survey (SDSS) main galaxy sample. We have applied our "npt" algorithm to both a volume-limited (36738 galaxies) and magnitude-limited sample (134741 galaxies) of SDSS galaxies, and find consistent results between the two samples, thus confirming the weak luminosity dependence of Q_z recently seen by other authors. We compare our results to other Q_z measurements in the literature and find it to be consistent within the full jack-knife error estimates. However, we find these errors are significantly increased by the presence of the ``Sloan Great Wall'' (at z ~ 0.08) within these two SDSS datasets, which changes the 3-point correlation function (3PCF) by 70% on large scales (s>=10h^-1 Mpc). If we exclude this supercluster, our observed Q_z is in better agreement with that obtained from the 2dFGRS by other authors, thus demonstrating the sensitivity of these higher-order correlation functions to large-scale structures in the Universe. This analysis highlights that the SDSS datasets used here are not ``fair samples'' of the Universe for the estimation of higher-order clustering statistics and larger volumes are required. We study the shape-dependence of Q_z(s,q,theta) as one expects this measurement to depend on scale if the large scale structure in the Universe has grown via gravitational instability from Gaussian initial conditions. On small scales (s <= 6h^-1 Mpc), we see some evidence for shape-dependence in Q_z, but at present our measurements are consistent with a constant within the errors (Q_z ~ 0.75 +/- 0.05). On scales >10h^-1 Mpc, we see considerable shape-dependence in Q_z.

We outline our first steps towards marrying two new and emerging technologies; the Virtual Observatory (e.g, AstroGrid) and the computational grid. We discuss the construction of VOTechBroker, which is a modular software tool designed to abstract the tasks of submission and management of a large number of computational jobs to a distributed computer system. The broker will also interact with the AstroGrid workflow and MySpace environments. We present our planned usage of the VOTechBroker in computing a huge number of n-point correlation functions from the SDSS, as well as fitting over a million CMBfast models to the WMAP data.

There is a growing need for massive computational resources for the analysis of new astronomical datasets. To tackle this problem, we present here our first steps towards marrying two new and emerging technologies; the Virtual Observatory (e.g, AstroGrid) and the computational grid (e.g. TeraGrid, COSMOS etc.). We discuss the construction of VOTechBroker, which is a modular software tool designed to abstract the tasks of submission and management of a large number of computational jobs to a distributed computer system. The broker will also interact with the AstroGrid workflow and MySpace environments. We discuss our planned usages of the VOTechBroker in computing a huge number of n-point correlation functions from the SDSS data and massive model-fitting of millions of CMBfast models to WMAP data. We also discuss other applications including the determination of the XMM Cluster Survey selection function and the construction of new WMAP maps.

We present new measurements of the quasar autocorrelation from a sample of \~80,000 photometrically-classified quasars taken from SDSS DR1. We find a best-fit model of $\omega(\theta) = (0.066\pm^{0.026}_{0.024})\theta^{-(0.98\pm0.15)}$ for the angular autocorrelation, consistent with estimates from spectroscopic quasar surveys. We show that only models with little or no evolution in the clustering of quasars in comoving coordinates since z~1.4 can recover a scale-length consistent with local galaxies and Active Galactic Nuclei (AGNs). A model with little evolution of quasar clustering in comoving coordinates is best explained in the current cosmological paradigm by rapid evolution in quasar bias. We show that quasar biasing must have changed from b_Q~3 at a (photometric) redshift of z=2.2 to b_Q~1.2-1.3 by z=0.75. Such a rapid increase with redshift in biasing implies that quasars at z~2 cannot be the progenitors of modern L* objects, rather they must now reside in dense environments, such as clusters. Similarly, the duration of the UVX quasar phase must be short enough to explain why local UVX quasars reside in essentially unbiased structures. Our estimates of b_Q are in good agreement with recent spectroscopic results, which demonstrate the implied evolution in b_Q is consistent with quasars inhabiting halos of similar mass at every redshift. Treating quasar clustering as a function of both redshift and luminosity, we find no evidence for luminosity dependence in quasar clustering, and that redshift evolution thus affects quasar clustering more than changes in quasars' luminosity. We provide a new method for quantifying stellar contamination in photometrically-classified quasar catalogs via the correlation function.

Cost effective, new antenna technology is required to build the large collecting area being planned for the next generation of radio telescopes. The Large Adaptive Reflector (LAR) is a novel concept being developed in Canada to meet this challenge. A prototype with a 200-350m diameter reflector, operating up to 2 GHz with a bandwidth of 750 MHz is planned. With a collecting area up to ~10% of the planned SKA, and an array feed capable of imaging a 0.3 deg^2 field-of-view, the prototype would address a number of the most compelling questions in modern astrophysics.

We present an 8 sigma detection of cosmic magnification measured by the variation of quasar density due to gravitational lensing by foreground large scale structure. To make this measurement we used 3800 square degrees of photometric observations from the Sloan Digital Sky Survey (SDSS) containing \~200,000 quasars and 13 million galaxies. Our measurement of the galaxy-quasar cross-correlation function exhibits the amplitude, angular dependence and change in sign as a function of the slope of the observed quasar number counts that is expected from magnification bias due to weak gravitational lensing. We show that observational uncertainties (stellar contamination, Galactic dust extinction, seeing variations and errors in the photometric redshifts) are well controlled and do not significantly affect the lensing signal. By weighting the quasars with the number count slope, we combine the cross-correlation of quasars for our full magnitude range and detect the lensing signal at >4 sigma in all five SDSS filters. Our measurements of cosmic magnification probe scales ranging from 60 kpc/h to 10 Mpc/h and are in good agreement with theoretical predictions based on the WMAP concordance cosmology. As with galaxy-galaxy lensing, future measurements of cosmic magnification will provide useful constraints on the galaxy-mass power spectrum.

We present a catalog of 100,563 unresolved, UV-excess (UVX) quasar candidates to g=21 from 2099 deg^2 of the Sloan Digital Sky Survey (SDSS) Data Release One (DR1) imaging data. Existing spectra of 22,737 sources reveals that 22,191 (97.6%) are quasars; accounting for the magnitude dependence of this efficiency, we estimate that 95,502 (95.0%) of the objects in the catalog are quasars. Such a high efficiency is unprecedented in broad-band surveys of quasars. This ``proof-of-concept'' sample is designed to be maximally efficient, but still has 94.7% completeness to unresolved, g<~19.5, UVX quasars from the DR1 quasar catalog. This efficient and complete selection is the result of our application of a probability density type analysis to training sets that describe the 4-D color distribution of stars and spectroscopically confirmed quasars in the SDSS. Specifically, we use a non-parametric Bayesian classification, based on kernel density estimation, to parameterize the color distribution of astronomical sources -- allowing for fast and robust classification. We further supplement the catalog by providing photometric redshifts and matches to FIRST/VLA, ROSAT, and USNO-B sources. Future work needed to extend the this selection algorithm to larger redshifts, fainter magnitudes, and resolved sources is discussed. Finally, we examine some science applications of the catalog, particularly a tentative quasar number counts distribution covering the largest range in magnitude (14.2<g<21.0) ever made within the framework of a single quasar survey.

We present 21-cm HI line and 408 MHz and 1.4 GHz continuum observations of the dwarf galaxy UGCA 86 with the DRAO synthesis telescope. UGCA 86 is detected in the continuum at 408 MHz (S_{408} = 120 +/- 30 mJy) and 1.4 GHz (S_{1400} = 79 +/- 3 mJy). The HI structure of UGCA 86 is complex, with two separate components: a rotating disk and a highly elongated spur that is kinematically disjunct from the disk. The HI disk is centered on the optical galaxy with similar axial ratio and orientation of the major axis. An area of the disk with a peculiar velocity of 25 km/s relative to the regular rotation of the disk is found on the southern side, where most of the star formation is concentrated. The spur is seen along the minor axis of UGCA 86 and overlaps in part with the disk. Towards the optical center of UGCA 86, the velocity difference between the spur and the disk is 40 km/s, about one third of the rotation velocity of the HI disk at 6 kpc from the center. This implies a large radial component of the orbital velocity of the spur, and therefore a significantly non-circular orbit. The velocity dispersion of the disk is 8.8 km/s, whereas the velocity dispersion of the spur varies from 10 km/s to 30 km/s. A possible tidal origin of the spur is considered in view of the proximity of the large Scd galaxy IC 342. However, evidence that the spur is located far outside the plane of the HI disk, and the absence of evidence for a warp in the outer HI disk, poses a problem for the interpretation of the spur as a tidal tail induced by IC 342. Detailed modelling of the IC 342/UGCA 86 system will be required before a tidal origin of the spur can be dismissed conclusively. The possibility that the spur is part of the nascent cloud of UGCA 86 or the remains of a small dwarf galaxy is presented as an alternative interpretation.

We present the two--point correlation function (2PCF) of narrow-line active galactic nuclei (AGN) selected within the First Data Release of the Sloan Digital Sky Survey. Using a sample of 13605 AGN in the redshift range 0.055 < z < 0.2, we find that the AGN auto--correlation function is consistent with the observed galaxy auto--correlation function on scales 0.2h^{-1}Mpc to >100h^{-1}Mpc. The AGN hosts trace an intermediate population of galaxies and are not detected in either the bluest (youngest) disk--dominated galaxies or many of the reddest (oldest) galaxies. We show that the AGN 2PCF is dependent on the luminosity of the narrow [OIII] emission line (L_{[OIII]}), with low L_{[OIII]} AGN having a higher clustering amplitude than high L_{[OIII]} AGN. This is consistent with lower activity AGN residing in more massive galaxies than higher activity AGN, and L_{[OIII]} providing a good indicator of the fueling rate. Using a model relating halo mass to black hole mass in cosmological simulations, we show that AGN hosted by ~ 10^{12} M_{odot} dark matter halos have a 2PCF that matches that of the observed sample. This mass scale implies a mean black hole mass for the sample of M_{BH} ~ 10^8 M_{odot}.

Many fundamental statistical methods have become critical tools for scientific data analysis yet do not scale tractably to modern large datasets. This paper will describe very recent algorithms based on computational geometry which have dramatically reduced the computational complexity of 1) kernel density estimation (which also extends to nonparametric regression, classification, and clustering), and 2) the n-point correlation function for arbitrary n. These new multi-tree methods typically yield orders of magnitude in speedup over the previous state of the art for similar accuracy, making millions of data points tractable on desktop workstations for the first time.

We analyse the observed correlation between galaxy environment and H-alpha emission line strength, using volume-limited samples and group catalogues of 24968 galaxies drawn from the 2dF Galaxy Redshift Survey (Mb<-19.5) and the Sloan Digital Sky Survey (Mr<-20.6). We characterise the environment by 1) Sigma_5, the surface number density of galaxies determined by the projected distance to the 5th nearest neighbour; and 2) rho1.1 and rho5.5, three-dimensional density estimates obtained by convolving the galaxy distribution with Gaussian kernels of dispersion 1.1 Mpc and 5.5 Mpc, respectively. We find that star-forming and quiescent galaxies form two distinct populations, as characterised by their H-alpha equivalent width, EW(Ha). The relative numbers of star-forming and quiescent galaxies varies strongly and continuously with local density. However, the distribution of EW(Ha) amongst the star-forming population is independent of environment. The fraction of star-forming galaxies shows strong sensitivity to the density on large scales, rho5.5, which is likely independent of the trend with local density, rho1.1. We use two differently-selected group catalogues to demonstrate that the correlation with galaxy density is approximately independent of group velocity dispersion, for sigma=200-1000 km/s. Even in the lowest density environments, no more than ~70 per cent of galaxies show significant H-alpha emission. Based on these results, we conclude that the present-day correlation between star formation rate and environment is a result of short-timescale mechanisms that take place preferentially at high redshift, such as starbursts induced by galaxy-galaxy interactions.

We present 1420 MHz (lambda=21cm) observations of polarized emission from an area of 117 degree square in the Galactic plane in Cygnus, covering 82 < l < 95, -3.5 < b < +5.5, a complex region where the line of sight is directed nearly along the Local spiral arm. The angular resolution is ~ 1', and structures as large as 45' are fully represented in the images. Polarization features bear little resemblance to features detected in total power: while the polarized signal arises in diffuse Galactic synchrotron emission regions, the appearance of the polarized sky is dominated by Faraday rotation occurring in small-scale structure in the intervening Warm Ionized Medium. There is no concentration of polarization structure towards the Galactic plane, indicating that both the emission and Faraday rotation occur nearby. We develop a conceptual framework for interpretation of the observations. We can detect only that polarized emission which has its origin closer than the polarization horizon, at a distance d_ph; more distant polarized emission is undetectable because of depth depolarization (differential Faraday rotation) and/or beam depolarization (due to internal and external Faraday dispersion). d_ph depends on the instrument used (frequency and beamwidth) as well as the direction being studied. In our data we find that d_ph ~ 2 kpc, consistent with the polarization features originating in the Local arm.

We present here a new algorithm for the fast computation of N-point correlation functions in large astronomical data sets. The algorithm is based on kdtrees which are decorated with cached sufficient statistics thus allowing for orders of magnitude speed-ups over the naive non-tree-based implementation of correlation functions. We further discuss the use of controlled approximations within the computation which allows for further acceleration. In summary, our algorithm now makes it possible to compute exact, all-pairs, measurements of the 2, 3 and 4-point correlation functions for cosmological data sets like the Sloan Digital Sky Survey (SDSS; York et al. 2000) and the next generation of Cosmic Microwave Background experiments (see Szapudi et al. 2000).

We describe an aperture synthesis radio telescope optimized for studies of the Galactic interstellar medium (ISM), providing the ability to image extended structures with high angular resolution over wide fields. The telescope produces images of atomic hydrogen emission using the 21-cm HI spectral line, and, simultaneously, continuum emission in two bands centred at 1420 MHz and 408 MHz, including linearly polarized emission at 1420 MHz, with synthesized beams of 1 degree and 3.4 degrees at the respective frequencies.

We have used polarimetric imaging to study the magneto-ionic medium of the Galaxy, obtaining 1420 MHz images with an angular resolution of 1' over more than 40 square-degrees of sky around the W3/W4/W5/HB3 HII region/SNR complex in the Perseus Arm. Features detected in polarization angle are imposed on the linearly polarized Galactic synchrotron background emission by Faraday rotation arising in foreground ionized gas having an emission measure as low as 1 cm^{-6} pc. Several new remarkable phenomena have been identified, including: mottled polarization arising from random fluctuations in a magneto-ionic screen that we identify with a medium in the Perseus Arm, probably in the vicinity of the HII regions themselves; depolarization arising from very high rotation measures (several times 10^3 rad m^{-2}) and rotation measure gradients due to the dense, turbulent environs of the HII regions; highly ordered features spanning up to several degrees; and an extended influence of the HII regions beyond the boundaries defined by earlier observations. In particular, the effects of an extended, low-density ionized halo around the HII region W4 are evident, probably an example of the extended HII envelopes postulated as the origin of weak recombination-line emission detected from the Galactic ridge. Our polarization observations can be understood if the uniform magnetic field component in this envelope scales with the square-root of electron density and is 20 microG at the edge of the depolarized region around W4, although this is probably an over-estimate since the random field component will have a significant effect.

The analysis and an efficient scientific exploration of the Digital Palomar Observatory Sky Survey (DPOSS) represents a major technical challenge. The input data set consists of 3 Terabytes of pixel information, and contains a few billion sources. We describe some of the specific scientific problems posed by the data, including searches for distant quasars and clusters of galaxies, and the data-mining techniques we are exploring in addressing them. Machine-assisted discovery methods may become essential for the analysis of such multi-Terabyte data sets. New and future approaches involve unsupervised classification and clustering analysis in the Giga-object data space, including various Bayesian techniques. In addition to the searches for known types of objects in this data base, these techniques may also offer the possibility of discovering previously unknown, rare types of astronomical objects.