The standard explanation for galaxy spin starts with the tidal-torque theory (TTT), in which an ellipsoidal dark-matter protohalo, which comes to host the galaxy, is torqued up by the tidal gravitational field around it. We discuss a complementary picture, using the relatively familiar velocity field, instead of the tidal field, whose intuitive connection to the surrounding, possibly faraway matter arrangement is more obscure. In this 'spin from primordial inner motions' (SPIM) concept, implicit in TTT derivations but not previously emphasized, the angular momentum from the gravity-sourced velocity field inside a protohalo largely cancels out, but has some excess from the aspherical outskirts. At first, the net spin scales according to linear theory, a sort of comoving conservation of familiar angular momentum. Then, at collapse, it is conserved in physical coordinates. Small haloes are then typically subject to secondary exchanges of angular momentum. The TTT is useful for analytic estimates. But a literal interpretation of the TTT is inaccurate in detail, without some implicit concepts about smoothing of the velocity and tidal fields. This could lead to misconceptions, for those first learning about how galaxies come to spin. Protohaloes are not perfectly ellipsoidal and do not uniformly torque up, as in a naive interpretation of the TTT; their inner velocity fields retain substantial dispersion. Furthermore, quantitatively, given initial conditions and protohalo boundaries, SPIM is more direct and accurate than the TTT to predict halo spins. We also discuss how SPIM applies to rotating filaments, and the relation between halo mass and spin, in which the total spin of a halo can be thought of as a sum of random contributions.

Here we propose a new scheme to reconstruct the baryon acoustic oscillations (BAO) signal, with key cosmological information, based on deep convolutional neural networks. After training the network with almost no fine-tuning, in the test set, the network recovers large-scale modes accurately: the correlation coefficient between the ground truth and recovered initial conditions still reach $90\%$ at $k \leq 0.2~ h\mathrm{Mpc}^{-1}$, which significantly improves the BAO signal-to-noise ratio until the scale $k=0.4~ h\mathrm{Mpc}^{-1}$. Furthermore, our scheme is independent of the survey boundary since it reconstructs initial condition based on local density distribution in configuration space, which means that we can gain more information from the whole survey space. Finally, we found our trained network is not sensitive to the cosmological parameters and works very well in those cosmologies close to that of our training set. This new scheme will possibly help us dig out more information from the current, on-going and future galaxy surveys.

We present SciServer, a science platform built and supported by the Institute for Data Intensive Engineering and Science at the Johns Hopkins University. SciServer builds upon and extends the SkyServer system of server-side tools that introduced the astronomical community to SQL (Structured Query Language) and has been serving the Sloan Digital Sky Survey catalog data to the public. SciServer uses a Docker/VM based architecture to provide interactive and batch mode server-side analysis with scripting languages like Python and R in various environments including Jupyter (notebooks), RStudio and command-line in addition to traditional SQL-based data analysis. Users have access to private file storage as well as personal SQL database space. A flexible resource access control system allows users to share their resources with collaborators, a feature that has also been very useful in classroom environments. All these services, wrapped in a layer of REST APIs, constitute a scalable collaborative data-driven science platform that is attractive to science disciplines beyond astronomy.

Cosmological N-body simulations are crucial for understanding how the Universe evolves. Studying large-scale distributions of matter in these simulations and comparing them to observations usually involves detecting dense clusters of particles called "halos,'' which are gravitationally bound and expected to form galaxies. However, traditional cluster finders are computationally expensive and use massive amounts of memory. Recent work by Liu et al (Liu et al. (2015)) showed the connection between cluster detection and memory-efficient streaming algorithms and presented a halo finder based on heavy hitter algorithm. Later, Ivkin et al. (Ivkin et al. (2018)) improved the scalability of suggested streaming halo finder with efficient GPU implementation. Both works map particles' positions onto a discrete grid, and therefore lose the rest of the information, such as their velocities. Therefore, two halos travelling through each other are indistinguishable in positional space, while the velocity distribution of those halos can help to identify this process which is worth further studying. In this project we analyze data from the Millennium Simulation Project (Springel et al. (2005)) to motivate the inclusion of the velocity into streaming method we introduce. We then demonstrate a use of suggested method, which allows one to find the same halos as before, while also detecting those which were indistinguishable in prior methods.

A Kavli foundation sponsored workshop on the theme \emph{Petabytes to Science} was held 12$^{th}$ to 14$^{th}$ of February 2019 in Las Vegas. The aim of the this workshop was to discuss important trends and technologies which may support astronomy. We also tackled how to better shape the workforce for the new trends and how we should approach education and public outreach. This document was coauthored during the workshop and edited in the weeks after. It comprises the discussions and highlights many recommendations which came out of the workshop. We shall distill parts of this document and formulate potential white papers for the decadal survey.

We present the multi-GPU realization of the StePS (Stereographically Projected Cosmological Simulations) algorithm with MPI-OpenMP-CUDA hybrid parallelization and nearly ideal scale-out to multiple compute nodes. Our new zoom-in cosmological direct N-body simulation method simulates the infinite universe with unprecedented dynamic range for a given amount of memory and, in contrast to traditional periodic simulations, its fundamental geometry and topology match observations. By using a spherical geometry instead of periodic boundary conditions, and gradually decreasing the mass resolution with radius, our code is capable of running simulations with a few gigaparsecs in diameter and with a mass resolution of $\sim 10^{9}M_{\odot}$ in the center in four days on three compute nodes with four GTX 1080Ti GPUs in each. The code can also be used to run extremely fast simulations with reasonable resolution for fitting cosmological parameters. These simulations are useful for prediction needs of large surveys. The StePS code is publicly available for the research community.

The Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production.

This paper describes the organization of the database and the catalog data products from the Pan-STARRS1 $3\pi$ Steradian Survey. The catalog data products are available in the form of an SQL-based relational database from MAST, the Mikulski Archive for Space Telescopes at STScI. The database is described in detail, including the construction of the database, the provenance of the data, the schema, and how the database tables are related. Examples of queries for a range of science goals are included. The catalog data products are available in the form of an SQL-based relational database from MAST, the Mikulski Archive for Space Telescopes at STScI.

Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3\pi$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$\sigma$ point source limiting sensitivities in the stacked 3$\pi$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ \Delta ra, \Delta dec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

The most common statistic used to analyze large-scale structure surveys is the correlation function, or power spectrum. Here, we show how `slicing' the correlation function on local density brings sensitivity to interesting non-Gaussian features in the large-scale structure, such as the expansion or contraction of baryon acoustic oscillations (BAO) according to the local density. The sliced correlation function measures the large-scale flows that smear out the BAO, instead of just correcting them as reconstruction algorithms do. Thus, we expect the sliced correlation function to be useful in constraining the growth factor, and modified gravity theories that involve the local density. Out of the studied cases, we find that the run of the BAO peak location with density is best revealed when slicing on a $\sim 40$ Mpc/$h$ filtered density. But slicing on a $\sim100$ Mpc/$h$ filtered density may be most useful in distinguishing between underdense and overdense regions, whose BAO peaks are separated by a substantial $\sim 5$ Mpc/$h$ at $z=0$. We also introduce `curtain plots' showing how local densities drive particle motions toward or away from each other over the course of an $N$-body simulation.

Cosmological $N$-body simulations play a vital role in studying models for the evolution of the Universe. To compare to observations and make a scientific inference, statistic analysis on large simulation datasets, e.g., finding halos, obtaining multi-point correlation functions, is crucial. However, traditional in-memory methods for these tasks do not scale to the datasets that are forbiddingly large in modern simulations. Our prior paper proposes memory-efficient streaming algorithms that can find the largest halos in a simulation with up to $10^9$ particles on a small server or desktop. However, this approach fails when directly scaling to larger datasets. This paper presents a robust streaming tool that leverages state-of-the-art techniques on GPU boosting, sampling, and parallel I/O, to significantly improve performance and scalability. Our rigorous analysis of the sketch parameters improves the previous results from finding the centers of the $10^3$ largest halos to $\sim 10^4-10^5$, and reveals the trade-offs between memory, running time and number of halos. Our experiments show that our tool can scale to datasets with up to $\sim 10^{12}$ particles while using less than an hour of running time on a single GPU Nvidia GTX 1080.

In view of future high-precision large-scale structure surveys, it is important to quantify the percent and subpercent level effects in cosmological $N$-body simulations from which theoretical predictions are drawn. One such effect involves deciding whether to zero all modes above the one-dimensional Nyquist frequency, the so-called "corner" modes, in the initial conditions. We investigate this effect by comparing power spectra, density distribution functions, halo mass functions, and halo profiles in simulations with and without these modes. For a simulation with a mass resolution of $m_p \sim 10^{11}\,h^{-1}\,M_{\odot}$, we find that at $z>6$, the difference in the matter power spectrum is large at wavenumbers above $\sim 80$\% of $k_{\rm{Ny}}$, reducing to below 2\% at all scales by $z\sim 3$. Including corner modes results in a better match between low- and high-resolution simulations at wavenumbers around the Nyquist frequency of the low-resolution simulation, but the effect of the corner modes is smaller than the effect of particle discreteness. The differences in mass functions are 3\% for the smallest halos at $z=6$ for the $m_p \sim 10^{11}\,h^{-1}\,M_{\odot}$ simulation, but we find no significant difference in the stacked profiles of well-resolved halos at $z \leq 6$. Thus removing power at $|\mathbf{k}|>k_{\rm{Ny}}$ in the initial conditions of cosmological simulations has a small effect on small scales and high redshifts, typically below a few percent.

We present a flexible template-based photometric redshift estimation framework, implemented in C#, that can be seamlessly integrated into a SQL database (or DB) server and executed on-demand in SQL. The DB integration eliminates the need to move large photometric datasets outside a database for redshift estimation, and utilizes the computational capabilities of DB hardware. The code is able to perform both maximum likelihood and Bayesian estimation, and can handle inputs of variable photometric filter sets and corresponding broad-band magnitudes. It is possible to take into account the full covariance matrix between filters, and filter zero points can be empirically calibrated using measurements with given redshifts. The list of spectral templates and the prior can be specified flexibly, and the expensive synthetic magnitude computations are done via lazy evaluation, coupled with a caching of results. Parallel execution is fully supported. For large upcoming photometric surveys such as the LSST, the ability to perform in-place photo-z calculation would be a significant advantage. Also, the efficient handling of variable filter sets is a necessity for heterogeneous databases, for example the Hubble Source Catalog, and for cross-match services such as SkyQuery. We illustrate the performance of our code on two reference photo-z estimation testing datasets, and provide an analysis of execution time and scalability with respect to different configurations. The code is available for download at https://github.com/beckrob/Photo-z-SQL.

Observational astronomy in the time-domain era faces several new challenges. One of them is the efficient use of observations obtained at multiple epochs. The work presented here addresses faint object detection with multi-epoch data, and describes an incremental strategy for separating real objects from artifacts in ongoing surveys, in situations where the single-epoch data are summaries of the full image data, such as single-epoch catalogs of flux and direction estimates for candidate sources. The basic idea is to produce low-threshold single-epoch catalogs, and use a probabilistic approach to accumulate catalog information across epochs; this is in contrast to more conventional strategies based on co-added or stacked image data across all epochs. We adopt a Bayesian approach, addressing object detection by calculating the marginal likelihoods for hypotheses asserting there is no object, or one object, in a small image patch containing at most one cataloged source at each epoch. The object-present hypothesis interprets the sources in a patch at different epochs as arising from a genuine object; the no-object (noise) hypothesis interprets candidate sources as spurious, arising from noise peaks. We study the detection probability for constant-flux objects in a simplified Gaussian noise setting, comparing results based on single exposures and stacked exposures to results based on a series of single-epoch catalog summaries. Computing the detection probability based on catalog data amounts to generalized cross-matching: it is the product of a factor accounting for matching of the estimated fluxes of candidate sources, and a factor accounting for matching of their estimated directions. We find that probabilistic fusion of multi-epoch catalog information can detect sources with only modest sacrifice in sensitivity and selectivity compared to stacking.

This document describes the structures making up the VOTable standard. The main part of this document describes the adopted part of the VOTable standard; it is followed by appendices presenting extensions which have been proposed and/or discussed, but which are not part of the standard.

We present the methodology and data behind the photometric redshift database of the Sloan Digital Sky Survey Data Release 12 (SDSS DR12). We adopt a hybrid technique, empirically estimating the redshift via local regression on a spectroscopic training set, then fitting a spectrum template to obtain K-corrections and absolute magnitudes. The SDSS spectroscopic catalog was augmented with data from other, publicly available spectroscopic surveys to mitigate target selection effects. The training set is comprised of $1,976,978$ galaxies, and extends up to redshift $z\approx 0.8$, with a useful coverage of up to $z\approx 0.6$. We provide photometric redshifts and realistic error estimates for the $208,474,076$ galaxies of the SDSS primary photometric catalog. We achieve an average bias of $\overline{\Delta z_{\mathrm{norm}}} = 5.84 \times 10^{-5}$, a standard deviation of $\sigma \left(\Delta z_{\mathrm{norm}}\right)=0.0205$, and a $3\sigma$ outlier rate of $P_o=4.11\%$ when cross-validating on our training set. The published redshift error estimates and photometric error classes enable the selection of galaxies with high quality photometric redshifts. We also provide a supplementary error map that allows additional, sophisticated filtering of the data.

We simulate the tidal disruption of a collisionless N-body globular star cluster in a total of 300 different orbits selected to have galactocentric radii between 10 and 30 kpc in four dark matter halos: (a) a spherical halo with no subhalos, (b) a spherical halo with subhalos, (c) a realistic halo with no subhalos, and (d) a realistic halo with subhalos. This allows us to isolate and study how the halo's (lack of) dynamical symmetry and substructures affect the dispersal of tidal debris. The realistic halos are constructed from the snapshot of the Via Lactea II simulation at redshift zero. We find that the overall halo's lack of dynamical symmetry disperses tidal debris to make the streams fluffier, consistent with previous studies of tidal debris of dwarf galaxies in larger orbits than ours in this study. On the other hand, subhalos in realistic potentials can locally enhance the densities along streams, making streams denser than their counterparts in smooth potentials. We show that many long and thin streams can survive in a realistic and lumpy halo for a Hubble time. This suggests that upcoming stellar surveys will likely uncover more thin streams which may contain density gaps that have been shown to be promising probes for dark matter substructures.

We study the parity-odd part (that we shall call Doppler term) of the linear galaxy two-point correlation function that arises from wide-angle, velocity, Doppler lensing and cosmic acceleration effects. As it is important at low redshift and at large angular separations, the Doppler term is usually neglected in the current generation of galaxy surveys. For future wide-angle galaxy surveys such as Euclid, SPHEREx and SKA, however, we show that the Doppler term must be included. The effect of these terms is dominated by the magnification due to relativistic aberration effects and the slope of the galaxy redshift distribution and it generally mimics the effect of the local type primordial non-Gaussianity with the effective nonlinearity parameter $f_{\rm NL}^{\rm eff}$ of a few, we show that this would affect forecasts on measurements of $f_{\rm NL}$ at low-redshift. Our results show that a survey at low redshift with large number density over a wide area of the sky could detect the Doppler term with a signal-to-noise ratio of $\sim 1-20$, depending on survey specifications.

We analyse the correlations between continuum properties and emission line equivalent widths of star-forming and active galaxies from the Sloan Digital Sky Survey. Since upcoming large sky surveys will make broad-band observations only, including strong emission lines into theoretical modelling of spectra will be essential to estimate physical properties of photometric galaxies. We show that emission line equivalent widths can be fairly well reconstructed from the stellar continuum using local multiple linear regression in the continuum principal component analysis (PCA) space. Line reconstruction is good for star-forming galaxies and reasonable for galaxies with active nuclei. We propose a practical method to combine stellar population synthesis models with empirical modelling of emission lines. The technique will help generate more accurate model spectra and mock catalogues of galaxies to fit observations of the new surveys. More accurate modelling of emission lines is also expected to improve template-based photometric redshift estimation methods. We also show that, by combining PCA coefficients from the pure continuum and the emission lines, automatic distinction between hosts of weak active galactic nuclei (AGNs) and quiescent star-forming galaxies can be made. The classification method is based on a training set consisting of high-confidence starburst galaxies and AGNs, and allows for the similar separation of active and star-forming galaxies as the empirical curve found by Kauffmann et al. We demonstrate the use of three important machine learning algorithms in the paper: k-nearest neighbour finding, k-means clustering and support vector machines.

The Lagrangian dynamics of a single fluid element within a self-gravitational matter field is intrinsically non-local due to the presence of the tidal force. This complicates the theoretical investigation of the non-linear evolution of various cosmic objects, e.g. dark matter halos, in the context of Lagrangian fluid dynamics, since a fluid parcel with given initial density and shape may evolve differently depending on their environments. In this paper, we provide a statistical solution that could decouple this environmental dependence. After deriving the probability distribution evolution equation of the matter field, our method produces a set of closed ordinary differential equations whose solution is uniquely determined by the initial condition of the fluid element. Mathematically, it corresponds to the projected characteristic curve of the transport equation of the density-weighted probability density function (PDF). Consequently it is guaranteed that the one-point PDF would be preserved by evolving these local, yet non-linear, curves with the same set of initial data as the real system. Physically, these trajectories describe the mean evolution averaged over all environments by substituting the tidal tensor with its conditional average. For Gaussian distributed dynamical variables, this mean tidal tensor is simply proportional to the velocity shear tensor, and the dynamical system would recover the prediction of Zel'dovich approximation (ZA) with the further assumption of the linearized continuity equation. For Weakly non-Gaussian field, the averaged tidal tensor could be expanded perturbatively as a function of all relevant dynamical variables whose coefficients are determined by the statistics of the field.

Accurate modeling of nonlinearities in the galaxy bispectrum, the Fourier transform of the galaxy three-point correlation function, is essential to fully exploit it as a cosmological probe. In this paper, we present numerical and theoretical challenges in modeling the nonlinear bispectrum. First, we test the robustness of the matter bispectrum measured from N-body simulations using different initial conditions generators. We run a suite of N-body simulations using the Zel'dovich approximation and second-order Lagrangian perturbation theory (2LPT) at different starting redshifts, and find that transients from initial decaying modes systematically reduce the nonlinearities in the matter bispectrum. To achieve 1% accuracy in the matter bispectrum for $z\le3$ on scales $k<1$ $h$/Mpc, 2LPT initial conditions generator with initial redshift of $z\gtrsim100$ is required. We then compare various analytical formulas and empirical fitting functions for modeling the nonlinear matter bispectrum, and discuss the regimes for which each is valid. We find that the next-to-leading order (one-loop) correction from standard perturbation theory matches with N-body results on quasi-linear scales for $z\ge1$. We find that the fitting formula in Gil-Mar\'{\i}n et al. (2012) accurately predicts the matter bispectrum for $z\le1$ on a wide range of scales, but at higher redshifts, the fitting formula given in Scoccimarro & Couchman (2001) gives the best agreement with measurements from N-body simulations.

We simulate tidal streams in the presence and absence of substructures inside the zero redshift snapshot of the Via Lactea II (VL-2) simulation. A halo finder is used to remove and isolate the subhalos found inside the high resolution dark matter halo of VL-2, and the potentials for both the main halo and all the subhalos are constructed individually using the self-consistent field (SCF) method. This allows us to make direct comparison of tidal streams between a smooth halo and a lumpy halo without assuming idealized profiles or triaxial fits. We simulate the kinematics of a star cluster starting with the same orbital position but two different velocities. Although these two orbits are only moderately eccentric and have similar apo- and pericentric distances, we find that the two streams have very different morphologies. We conclude that our model of the potential of VL-2 can provide insights about tidal streams that have not been explored by previous studies using idealized or axisymmetric models.

We simulate deep images from the Hubble Space Telescope (HST) using semi-empirical models of galaxy formation with only a few basic assumptions and parameters. We project our simulations all the way to the observational domain, adding cosmological and instrumental effects to the images, and analyze them in the same way as real HST images ("forward modeling"). This is a powerful tool for testing and comparing galaxy evolution models, since it allows us to make unbiased comparisons between the predicted and observed distributions of galaxy properties, while automatically taking into account all relevant selection effects. Our semi-empirical models populate each dark matter halo with a galaxy of determined stellar mass and scale radius. We compute the luminosity and spectrum of each simulated galaxy from its evolving stellar mass using stellar population synthesis models. We calculate the intrinsic scatter in the stellar mass-halo mass relation that naturally results from enforcing a monotonically increasing stellar mass along the merger history of each halo. The simulated galaxy images are drawn from cutouts of real galaxies from the Sloan Digital Sky Survey, with sizes and fluxes rescaled to match those of the model galaxies. The distributions of galaxy luminosities, sizes, and surface brightnesses depend on the adjustable parameters in the models, and they agree well with observations for reasonable values of those parameters. Measured galaxy magnitudes and sizes have significant magnitude-dependent biases, with both being underestimated near the magnitude detection limit. The fraction of galaxies detected and fraction of light detected also depend sensitively on the details of the model.

The SKA will be capable of producing a stream of science data products that are Exa-scale in terms of their storage and processing requirements. This Google-scale enterprise is attracting considerable international interest and excitement from within the industrial and academic communities. In this chapter we examine the data flow, storage and processing requirements of a number of key SKA survey science projects to be executed on the baseline SKA1 configuration. Based on a set of conservative assumptions about trends for HPC and storage costs, and the data flow process within the SKA Observatory, it is apparent that survey projects of the scale proposed will potentially drive construction and operations costs beyond the current anticipated SKA1 budget. This implies a sharing of the resources and costs to deliver SKA science between the community and what is contained within the SKA Observatory. A similar situation was apparent to the designers of the LHC more than 10 years ago. We propose that it is time for the SKA project and community to consider the effort and process needed to design and implement a distributed SKA science data system that leans on the lessons of other projects and looks to recent developments in Cloud technologies to ensure an affordable, effective and global achievement of SKA science goals.

We investigate the nonlinear evolution of cosmic morphologies of the large-scale structure by examining the Lagrangian dynamics of various tensors of a cosmic fluid element, including the velocity gradient tensor, the Hessian matrix of the gravitational potential as well as the deformation tensor. Instead of the eigenvalue representation, the first two tensors, which associate with the "kinematic" and "dynamical" cosmic web classification algorithm respectively, are studied in a more convenient parameter space. These parameters are defined as the rotational invariant coefficients of the characteristic equation of the tensor. In the nonlinear local model (NLM) where the magnetic part of Weyl tensor vanishes, these invariants are fully capable of characterizing the dynamics. Unlike the Zeldovich approximation (ZA), where various morphologies do not change before approaching a one-dimensional singularity, the sheets in NLM are unstable for both overdense and underdense perturbations. While it has long been known that the coupling between tidal tensor and velocity shear would cause a filamentary final configuration of a collapsing region, we show that the underdense perturbation are more subtle, as the balance between the shear rate (tidal force) and the divergence (density) could lead to different morphologies. Interestingly, this instability also sets the basis for understanding some distinctions of the cosmic web identified dynamically and kinematically. We show that the sheets with negative density perturbation in the potential based algorithm would turn to filaments faster than in the kinematic method, which could explain the distorted dynamical filamentary structure observed in the simulation.

Baryon acoustic oscillations (BAO) are a powerful probe of the expansion history of the universe, which can tell us about the nature of dark energy. In order to accurately characterize the dark energy equation of state using BAO, we must understand the effects of both nonlinearities and redshift space distortions on the location and shape of the acoustic peak. In a previous paper we introduced a novel approach to 2nd order perturbation theory in configuration space using the Zel'dovich approximation, and presented a simple result for the first nonlinear term of the correlation function. In this paper, we extend this approach to redshift space. We show how perform the computation, and present the analytic result for the first nonlinear term in the correlation function. Finally, we validate our result through comparison to numerical simulations.

We present photometric redshifts and associated probability distributions for all detected sources in the Extended Chandra Deep Field South (ECDFS). The work makes use of the most up-to-date data from the Cosmic Assembly Near-IR Deep Legacy Survey (CANDELS) and the Taiwan ECDFS Near-Infrared Survey (TENIS) in addition to other data. We also revisit multi-wavelength counterparts for published X-ray sources from the 4Ms-CDFS and 250ks-ECDFS surveys, finding reliable counterparts for 1207 out of 1259 sources ($\sim 96\%$). Data used for photometric redshifts include intermediate-band photometry deblended using the TFIT method, which is used for the first time in this work. Photometric redshifts for X-ray source counterparts are based on a new library of AGN/galaxy hybrid templates appropriate for the faint X-ray population in the CDFS. Photometric redshift accuracy for normal galaxies is 0.010 and for X-ray sources is 0.014, and outlier fractions are $4\%$ and $5.4\%$ respectively. The results within the CANDELS coverage area are even better as demonstrated both by spectroscopic comparison and by galaxy-pair statistics. Intermediate-band photometry, even if shallow, is valuable when combined with deep broad-band photometry. For best accuracy, templates must include emission lines.

Future galaxy surveys will measure baryon acoustic oscillations (BAOs) with high significance, and a complete understanding of the anisotropies of BAOs in redshift space will be important to exploit the cosmological information in BAOs. Here we describe the anisotropies that arise in the redshift-space galaxy two-point correlation function (2PCF) and elucidate the origin of features that arise in the dependence of the BAOs on the angle between the orientation of the galaxy pair and the line of sight. We do so with a derivation of the configuration-space 2PCF using streaming model. We find that, contrary to common belief, the locations of BAO peaks in the redshift-space 2PCF are anisotropic even in the linear theory. Anisotropies in BAO depend strongly on the method of extracting the peak, showing maximum 3 % angular variation. We also find that extracting the BAO peak of $r^2\xi(r,\mu)$ significantly reduces the anisotropy to sub-percent level angular variation. When subtracting the tilt due to the broadband behavior of the 2PCF, the BAO bump is enhanced along the line of sight because of local infall velocities toward the BAO bump. Precise measurement of the angular dependence of the redshift-space 2PCF will allow new geometrical tests of dark energy beyond the BAO.

As an alternative way of describing the cosmological velocity field, we discuss the evolution of rotational invariants constructed from the velocity gradient tensor. Compared with the traditional divergence-vorticity decomposition, these invariants, defined as coefficients of characteristic equation of the velocity gradient tensor, enable a complete classification of all possible flow patterns in the dark-matter comoving frame, including both potential and vortical flows. We show that this tool, first introduced in turbulence two decades ago, proves to be very useful in understanding the evolution of the cosmic web structure, and in classifying its morphology. Before shell-crossing, different categories of potential flow are highly associated with cosmic web structure, because of the coherent evolution of density and velocity. This correspondence is even preserved at some level when vorticity is generated after shell-crossing. The evolution from the potential to vortical flow can be traced continuously by these invariants. With the help of this tool, we show that the vorticity is generated in a particular way that is highly correlated with the large-scale structure. This includes a distinct spatial distribution and different types of alignment between cosmic web and vorticity direction for various vortical flows. Incorporating shell-crossing into closed dynamical systems is highly non-trivial, but we propose a possible statistical explanation for some of these phenomena relating to the internal structure of the three-dimensional invariants space.

Observations of diffuse Galactic gamma ray emission (DGE) by the Fermi Large Area Telescope (LAT) allow a detailed study of cosmic rays and the interstellar medium. However, diffuse emission models of the inner Galaxy underpredict the Fermi-LAT data at energies above a few GeV and hint at possible non-astrophysical sources including dark matter (DM) annihilations or decays. We present a study of the possible emission components from DM using the high-resolution Via Lactea II N-body simulation of a Milky Way-sized DM halo. We generate full-sky maps of DM annihilation and decay signals that include modeling of the adiabatic contraction of the host density profile, Sommerfeld enhanced DM annihilations, $p$-wave annihilations, and decaying DM. We compare our results with the DGE models produced by the Fermi-LAT team over different sky regions, including the Galactic center, high Galactic latitudes, and the Galactic anti-center. This work provides possible templates to fit the observational data that includes the contribution of the subhalo population to DM gamma-ray emission, with the significance depending on the annihilation/decay channels and the Galactic regions being considered.

Not only source catalogs are extracted from astronomy observations. Their sky coverage is always carefully recorded and used in statistical analyses, such as correlation and luminosity function studies. Here we present a novel method for catalog matching, which inherently builds on the coverage information for better performance and completeness. A modified version of the Zones Algorithm is introduced for matching partially overlapping observations, where irrelevant parts of the data are excluded up front for efficiency. Our design enables searches to focus on specific areas on the sky to further speed up the process. Another important advantage of the new method over traditional techniques is its ability to quickly detect dropouts, i.e., the missing components that are in the observed regions of the celestial sphere but did not reach the detection limit in some observations. These often provide invaluable insight into the spectral energy distribution of the matched sources but rarely available in traditional associations.

Understanding the diversity in spectra is the key to determining the physical parameters of galaxies. The optical spectra of galaxies are highly convoluted with continuum and lines which are potentially sensitive to different physical parameters. Defining the wavelength regions of interest is therefore an important question. In this work, we identify informative wavelength regions in a single-burst stellar populations model by using the CUR Matrix Decomposition. Simulating the Lick/IDS spectrograph configuration, we recover the widely used Dn(4000), Hbeta, and HdeltaA to be most informative. Simulating the SDSS spectrograph configuration with a wavelength range 3450-8350 Angstrom and a model-limited spectral resolution of 3 Angstrom, the most informative regions are: first region-the 4000 Angstrom break and the Hdelta line; second region-the Fe-like indices; third region-the Hbeta line; fourth region-the G band and the Hgamma line. A Principal Component Analysis on the first region shows that the first eigenspectrum tells primarily the stellar age, the second eigenspectrum is related to the age-metallicity degeneracy, and the third eigenspectrum shows an anti-correlation between the strengths of the Balmer and the Ca K and H absorptions. The regions can be used to determine the stellar age and metallicity in early-type galaxies which have solar abundance ratios, no dust, and a single-burst star formation history. The region identification method can be applied to any set of spectra of the user's interest, so that we eliminate the need for a common, fixed-resolution index system. We discuss future directions in extending the current analysis to late-type galaxies.

We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

This document presents (off-line) computing requrements and challenges for Cosmic Frontier science, covering the areas of data management, analysis, and simulations. We invite contributions to extend the range of covered topics and to enhance the current descriptions.

Structures like galaxies and filaments of galaxies in the Universe come about from the origami-like folding of an initially flat three-dimensional manifold in 6D phase space. The ORIGAMI method identifies these structures in a cosmological simulation, delineating the structures according to their outer folds. Structure identification is a crucial step in comparing cosmological simulations to observed maps of the Universe. The ORIGAMI definition is objective, dynamical and geometric: filament, wall and void particles are classified according to the number of orthogonal axes along which dark-matter streams have crossed. Here, we briefly review these ideas, and speculate on how ORIGAMI might be useful to find cosmic voids.

Citizen science, in which volunteers work with professional scientists to conduct research, is expanding due to large online datasets. To plan projects, it is important to understand volunteers' motivations for participating. This paper analyzes results from an online survey of nearly 11,000 volunteers in Galaxy Zoo, an astronomy citizen science project. Results show that volunteers' primary motivation is a desire to contribute to scientific research. We encourage other citizen science projects to study the motivations of their volunteers, to see whether and how these results may be generalized to inform the field of citizen science.

Large-scale bulk flows in the Universe distort the initial density field, broadening the baryon-acoustic-oscillation (BAO) feature that was imprinted when baryons were strongly coupled to photons. Additionally, there is a small shift inward in the peak of the conventional overdensity correlation function, a mass-weighted statistic. This shift occurs when high density peaks move toward each other. We explore whether this shift can be removed by applying to the density field a transform (such as a logarithm) that gives fairer statistical weight to fluctuations in underdense regions. Using configuration-space perturbation theory in the Zel'dovich approximation, we find that the log-density correlation function shows a much smaller inward shift in the position of the BAO peak at low redshift than is seen in the overdensity correlation function. We also show that if the initial, Lagrangian density of matter parcels could be estimated at their Eulerian positions, giving a displaced-initial-density field, its peak shift would be even smaller. In fact, a transformed field that accentuates underdensities, such as the reciprocal of the density, pushes the peak the other way, outward. In our model, these shifts in the peak position can be attributed to shift terms, involving the derivative of the linear correlation function, that entirely vanish in this displaced-initial-density field.

This document summarizes the results of a community-based discussion of the potential science impact of the Mayall+BigBOSS highly multiplexed multi-object spectroscopic capability. The KPNO Mayall 4m telescope equipped with the DOE- and internationally-funded BigBOSS spectrograph offers one of the most cost-efficient ways of accomplishing many of the pressing scientific goals identified for this decade by the "New Worlds, New Horizons" report. The BigBOSS Key Project will place unprecedented constraints on cosmological parameters related to the expansion history of the universe. With the addition of an open (publicly funded) community access component, the scientific impact of BigBOSS can be extended to many important astrophysical questions related to the origin and evolution of galaxies, stars, and the IGM. Massive spectroscopy is the critical missing ingredient in numerous ongoing and planned ground- and space-based surveys, and BigBOSS is unique in its ability to provide this to the US community. BigBOSS data from community-led projects will play a vital role in the education and training of students and in maintaining US leadership in these fields of astrophysics. We urge the NSF-AST division to support community science with the BigBOSS multi-object spectrograph through the period of the BigBOSS survey in order to ensure public access to the extraordinary spectroscopic capability.

Extreme value statistics (EVS) is applied to the distribution of galaxy luminosities in the Sloan Digital Sky Survey (SDSS). We analyze the DR8 Main Galaxy Sample (MGS), as well as the Luminous Red Galaxies (LRG). Maximal luminosities are sampled from batches consisting of elongated pencil beams in the radial direction of sight. For the MGS, results suggest a small and positive tail index $\xi$, effectively ruling out the possibility of having a finite maximum cutoff luminosity, and implying that the luminosity distribution function may decay as a power law at the high luminosity end. Assuming, however, $\xi=0$, a non-parametric comparison of the maximal luminosities with the Fisher-Tippett-Gumbel distribution (limit distribution for variables distributed by the Schechter fit) indicates a good agreement provided uncertainties arising both from the finite batch size and from the batch size distribution are accounted for. For a volume limited sample of LRGs, results show that they can be described as being the extremes of a luminosity distribution with an exponentially decaying tail, provided the uncertainties related to batch-size distribution are taken care of.

We present a method to extract the redshift-space distortions beta parameter in configuration space with a minimal set of cosmological assumptions. We show that a novel combination of the observed monopole and quadrupole correlation functions can remove efficiently the impact of mild non linearities and redshift errors. The method offers a series of convenient properties: it does not depend on the theoretical linear correlation function, the mean galaxy density is irrelevant, only convolutions are used, there is no explicit dependence on linear bias. Analyses based on dark matter N-body simulations and Fisher matrix demonstrate that errors of a few percent on beta are possible with a full sky, 1(Gpc/h)^3 survey centered at a redshift of unity and with negligible shot noise. We also find a baryonic feature in the normalized quadrupole in configuration space that should complicate the extraction of the growth parameter from the linear theory asymptote, but that does not have a major impact with our method.

We propose to describe the variety of galaxies from SDSS by using only one affine parameter. To this aim, we build the Principal Curve (P-curve) passing through the spine of the data point cloud, considering the eigenspace derived from Principal Component Analysis of morphological, physical and photometric galaxy properties. Thus, galaxies can be labeled, ranked and classified by a single arc length value of the curve, measured at the unique closest projection of the data points on the P-curve. We find that the P-curve has a "W" letter shape with 3 turning points, defining 4 branches that represent distinct galaxy populations. This behavior is controlled mainly by 2 properties, namely u-r and SFR. We further present the variations of several galaxy properties as a function of arc length. Luminosity functions variate from steep Schechter fits at low arc length, to double power law and ending in Log-normal fits at high arc length. Galaxy clustering shows increasing autocorrelation power at large scales as arc length increases. PCA analysis allowed to find peculiar galaxy populations located apart from the main cloud of data points, such as small red galaxies dominated by a disk, of relatively high stellar mass-to-light ratio and surface mass density. The P-curve allows not only dimensionality reduction, but also provides supporting evidence for relevant physical models and scenarios in extragalactic astronomy: 1) Evidence for the hierarchical merging scenario in the formation of a selected group of red massive galaxies. These galaxies present a log-normal r-band luminosity function, which might arise from multiplicative processes involved in this scenario. 2) Connection between the onset of AGN activity and star formation quenching, which appears in green galaxies when transitioning from blue to red populations. (Full abstract in downloadable version)

Multi-wavelength astronomical studies require cross-identification of detections of the same celestial objects in multiple catalogs based on spherical coordinates and other properties. Because of the large data volumes and spherical geometry, the symmetric N-way association of astronomical detections is a computationally intensive problem, even when sophisticated indexing schemes are used to exclude obviously false candidates. Legacy astronomical catalogs already contain detections of more than a hundred million objects while the ongoing and future surveys will produce catalogs of billions of objects with multiple detections of each at different times. The varying statistical error of position measurements, moving and extended objects, and other physical properties make it necessary to perform the cross-identification using a mathematically correct, proper Bayesian probabilistic algorithm, capable of including various priors. One time, pair-wise cross-identification of these large catalogs is not sufficient for many astronomical scenarios. Consequently, a novel system is necessary that can cross-identify multiple catalogs on-demand, efficiently and reliably. In this paper, we present our solution based on a cluster of commodity servers and ordinary relational databases. The cross-identification problems are formulated in a language based on SQL, but extended with special clauses. These special queries are partitioned spatially by coordinate ranges and compiled into a complex workflow of ordinary SQL queries. Workflows are then executed in a parallel framework using a cluster of servers hosting identical mirrors of the same data sets.

We present the ORIGAMI method of identifying structures, particularly halos, in cosmological N-body simulations. Structure formation can be thought of as the folding of an initially flat three-dimensional manifold in six-dimensional phase space. ORIGAMI finds the outer folds that delineate these structures. Halo particles are identified as those that have undergone shell-crossing along 3 orthogonal axes, providing a dynamical definition of halo regions that is independent of density. ORIGAMI also identifies other morphological structures: particles that have undergone shell-crossing along 2, 1, or 0 orthogonal axes correspond to filaments, walls, and voids respectively. We compare this method to a standard Friends-of-Friends halo-finding algorithm and find that ORIGAMI halos are somewhat larger, more diffuse, and less spherical, though the global properties of ORIGAMI halos are in good agreement with other modern halo-finding algorithms.

The U.S. Virtual Astronomical Observatory (VAO; http://www.us-vao.org/) has been in operation since May 2010. Its goal is to enable new science through efficient integration of distributed multi-wavelength data. This paper describes the management and organization of the VAO, and emphasizes the techniques used to ensure efficiency in a distributed organization. Management methods include using an annual program plan as the basis for establishing contracts with member organizations, regular communication, and monitoring of processes.

Baryon acoustic oscillations are an excellent technique to constrain the properties of dark energy in the Universe. In order to accurately characterize the dark energy equation of state, we must understand the effects of both the nonlinearities and redshift space distortions on the location and shape of the acoustic peak. In this paper, we consider these effects using the Zel'dovich approximation and a novel approach to 2nd order perturbation theory. The second order term of the Zel'dovich power spectrum is built from convolutions of the linear power spectrum with polynomial kernels in Fourier space, suggesting that the corresponding term of the the Zel'dovich correlation function can be written as a sum of quadratic products of a broader class of correlation functions, expressed through simple spherical Bessel transforms of the linear power spectrum. We show how to systematically perform such a computation. We explicitly prove that our result is the Fourier transform of the Zel'dovich power spectrum, and compare our expressions to numerical simulations. Finally, we highlight the advantages of writing the nonlinear expansion in configuration space, as this calculation is easily extended to redshift space, and the higher order terms are mathematically simpler than their Fourier counterparts.

The relationship between observed tracers such as galaxies and the underlying dark matter distribution is crucial in extracting cosmological information. As the linear bias model breaks down at quasi-linear scales, the standard perturbative approach of the nonlinear Eulerian bias model (EBM) is not accurate enough in describing galaxy clustering. In this paper, we discuss such a model in the context of resummed perturbation theory, and further generalize it to incorporate the subsequent gravitational evolution by combining with a Lagrangian description of galaxies' motion. The multipoint propagators we constructed for such model also exhibit exponential damping similar to their dark matter counterparts, therefore the convergence property of statistics built upon these quantities is improved. This is achieved by applying both Eulerian and Lagrangian resummation techniques of dark matter field developed in recent years. As inherited from the Lagrangian description of galaxy density evolution, our approach automatically incorporates the non-locality induced by gravitational evolution after the formation of the tracer, and also allows us to include a continuous galaxy formation history by temporally weighted-averaging relevant quantities with the galaxy formation rate.

Contrary to the common view voids have very complex internal structure and dynamics. Here we show how the hierarchy of structures in the density field inside voids is reflected by a similar hierarchy of structures in the velocity field. Voids defined by dense filaments and clusters can de described as simple expanding domains with coherent flows everywhere except at their boundaries. At scales smaller that the void radius the velocity field breaks into expanding sub-domains corresponding to sub- voids. These sub-domains break into even smaller sub-sub domains at smaller scales resulting in a nesting hierarchy of locally expanding domains. The ratio between the magnitude of the velocity field responsible for the expansion of the void and the velocity field defining the sub voids is approximately one order of magnitude. The small-scale components of the velocity field play a minor role in the shaping of the voids but they define the local dynamics directly affecting the processes of galaxy formation and evolution. The super-Hubble expansion inside voids makes them cosmic magnifiers by stretching their internal primordial density fluctuations allowing us to probe the small scales in the primordial density field. Voids also act like time machines by "freezing" the development of the medium-scale density fluctuations responsible for the formation of the tenuous web of structures seen connecting proto galaxies in computer simulations. As a result of this freezing haloes in voids can remain "connected" to this tenuous web until the present time. This may have an important effect in the formation and evolution of galaxies in voids by providing an efficient gas accretion mechanism via coherent low-velocity streams that can keep a steady inflow of matter for extended periods of time.

We investigate the use of a logarithmic density variable in estimating the Lagrangian displacement field, motivated by the success of a logarithmic transformation in restoring information to the matter power spectrum. The logarithmic relation is an extension of the linear relation, motivated by the continuity equation, in which the density field is assumed to be proportional to the divergence of the displacement field; we compare the linear and logarithmic relations by measuring both of these fields directly in a cosmological N-body simulation. The relative success of the logarithmic and linear relations depends on the scale at which the density field is smoothed. Thus we explore several ways of measuring the density field, including Cloud-In-Cell smoothing, adaptive smoothing, and the (scale-independent) Delaunay tessellation, and we use both a Fourier space and a geometrical tessellation approach to measuring the divergence. We find that the relation between the divergence of the displacement field and the density is significantly tighter with a logarithmic density variable, especially at low redshifts and for very small (~2 Mpc/h) smoothing scales. We find that the grid-based methods are more reliable than the tessellation-based method of calculating both the density and the divergence fields, though in both cases the logarithmic relation works better in the appropriate regime, which corresponds to nonlinear scales for the grid-based methods and low densities for the tessellation-based method.

We present a study of local post-starburst galaxies (PSGs) using the photometric and spectroscopic observations from the Sloan Digital Sky Survey (SDSS) and the results from the Galaxy Zoo project. We find that the majority of our local PSG population have neither early- nor late- type morphologies but occupy a well-defined space within the colour-stellar mass diagram, most notably, the low-mass end of the "green valley" below the transition mass thought to be the mass division between low-mass star-forming galaxies and high-mass passively-evolving bulge-dominated galaxies. Our analysis suggests that it is likely that a local PSG will quickly transform into "red", low-mass early-type galaxies as the stellar morphologies of the "green" PSGs largely resemble that of the early-type galaxies within the same mass range. We propose that the current population of PSGs represents a population of galaxies which is rapidly transitioning between the star-forming and the passively-evolving phases. Subsequently, these PSGs will contribute towards the build-up of the low-mass end of the "red sequence" once the current population of young stars fade and stars are no longer being formed. These results are consistent with the idea of "downsizing" where the build-up of smaller galaxies occurs at later epochs.

In this work we present an atlas of composite spectra of galaxies based on the data of the Sloan Digital Sky Survey Data Release 7 (SDSS DR7). Galaxies are classified by colour, nuclear activity and star-formation activity to calculate average spectra of high signal-to-noise ratio and resolution (S/N = 132 - 4760 at {Dlambda = 1 A), using an algorithm that is robust against outliers. Besides composite spectra, we also compute the first five principal components of the distributions in each galaxy class to characterize the nature of variations of individual spectra around the averages. The continua of the composite spectra are fitted with BC03 stellar population synthesis models to extend the wavelength coverage beyond the coverage of the SDSS spectrographs. Common derived parameters of the composites are also calculated: integrated colours in the most popular filter systems, line strength measurements, and continuum absorption indices (including Lick indices). These derived parameters are compared with the distributions of parameters of individual galaxies and it is shown on many examples that the composites of the atlas cover much of the parameter space spanned by SDSS galaxies. By co-adding thousands of spectra, a total integration time of several months can be reached, which results in extremely low noise composites. The variations in redshift not only allow for extending the spectral coverage bluewards to the original wavelength limit of the SDSS spectrographs, but also make higher spectral resolution achievable. The composite spectrum atlas is available online at http://www.vo.elte.hu/compositeatlas.

We present a statistical study of the correlations and dimensionality of emission lines carried out on a sample of over 40,000 Sloan Digital Sky Survey (SDSS) galaxies. Using principal component analysis, we found that the equivalent widths of the 11 strongest lines can be well represented using three parameters. We also explore correlations of the emission pattern with the eigenspace representation of the continuum spectrum. The observed relations are used to provide an empirical prescription for expectation values and variances of emission-line strengths as a function of spectral shape. We show that this estimation of emission lines has a sufficient accuracy to make it suitable for photometric applications. The method has already proved useful in SDSS photometric redshift estimation.

Star formation is one of the most important processes in galaxy formation. The luminosity of Halpha recombination and [OII] forbidden emissions remain to be most used in measuring formation rate of massive stars in galaxies. Here we report the inclination dependency of continuum-subtracted and aperture-corrected nebular luminosities, including Halpha, Hbeta, Hgamma, [OII], [NII], of disk-dominated galaxies in the local universe. Their luminosities decrease by a factor of three from face-on to edge-on (axis ratio limit = 0.17) orientations. This dependence is deduced to be caused by extinction due to diffuse dust within the disks with an amplitude of 1.2 mag. The line-luminosity--inclination relation provides a novel way to remove extinction in emission lines and present star formation rate of disk galaxies out to redshift of 1.6.

This document describes the Astronomical Data Query Language (ADQL). ADQL has been developed based on SQL92. This document describes the subset of the SQL grammar supported by ADQL. Special restrictions and extensions to SQL92 have been defined in order to support generic and astronomy specific operations.

This specification defines a simple query protocol for retrieving records from a catalog of astronomical sources. The query describes sky position and an angular distance, defining a cone on the sky. The response returns a list of astronomical sources from the catalog whose positions lie within the cone, formatted as a VOTable. This version of the specification is essentially a transcription of the original Cone Search specification in order to move it into the IVOA standardization process.

BigBOSS is a Stage IV ground-based dark energy experiment to study baryon acoustic oscillations (BAO) and the growth of structure with a wide-area galaxy and quasar redshift survey over 14,000 square degrees. It has been conditionally accepted by NOAO in response to a call for major new instrumentation and a high-impact science program for the 4-m Mayall telescope at Kitt Peak. The BigBOSS instrument is a robotically-actuated, fiber-fed spectrograph capable of taking 5000 simultaneous spectra over a wavelength range from 340 nm to 1060 nm, with a resolution R = 3000-4800. Using data from imaging surveys that are already underway, spectroscopic targets are selected that trace the underlying dark matter distribution. In particular, targets include luminous red galaxies (LRGs) up to z = 1.0, extending the BOSS LRG survey in both redshift and survey area. To probe the universe out to even higher redshift, BigBOSS will target bright [OII] emission line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are obtained to measure the BAO feature, trace the matter power spectrum at smaller scales, and detect redshift space distortions. BigBOSS will provide additional constraints on early dark energy and on the curvature of the universe by measuring the Ly-alpha forest in the spectra of over 600,000 2.2 < z < 3.5 quasars. BigBOSS galaxy BAO measurements combined with an analysis of the broadband power, including the Ly-alpha forest in BigBOSS quasar spectra, achieves a FOM of 395 with Planck plus Stage III priors. This FOM is based on conservative assumptions for the analysis of broad band power (kmax = 0.15), and could grow to over 600 if current work allows us to push the analysis to higher wave numbers (kmax = 0.3). BigBOSS will also place constraints on theories of modified gravity and inflation, and will measure the sum of neutrino masses to 0.024 eV accuracy.

The matter density field exhibits a nearly lognormal probability density distribution (PDF) after entering into the nonlinear regime. Recently, it has been shown that the shape of the power spectrum of a logarithmically transformed density field is very close to the linear density power spectrum, motivating an analytic study of it. In this paper, we develop cosmological perturbation theory for the power spectrum of this field. Our formalism is developed in the context of renormalized perturbation theory, which helps to regulate the convergence behavior of the perturbation series, and of the Taylor- series expansion we use of the logarithmic mapping. This approach allows us to handle the critical issue of density smoothing in a straightforward way. We also compare our perturbative results with simulation measurements.

We show that the "cold" Hubble flow observed for galaxies around the Milky Way does not represent a problem in cosmology but is due to the particular geometry and dynamics of our local wall. The behavior of the perturbed Hubble flow around the Milky Way is the result of two main factors: at small scales (R < 1 Mpc) the inflow is dominated by the gravitational influence of the Milky Way. At large scales (R > 1 Mpc) the out flow reflects the expansion of our local wall which "cools down" the peculiar velocities. This is an intrinsic property of walls and is independent of cosmology. We find the dispersion of the local Hubble flow (1 < R < 3 Mpc) around simulated "Milky Way" haloes located at the centre of low-density cosmological walls to be {\sigma}_H ~ 30 km/s, in excellent agreement with observations. The expansion of our local wall is also reflected in the value of the measured local Hubble constant. For "Milky Way" haloes inside walls, we find super-Hubble flows with h_local \simeq 0.77 - 1.13. The radius of equilibrium (R_0) depends not only on the mass of the central halo and the Hubble expansion but also on the dynamics given by the local LSS geometry. The super-Hubble flow inside our local wall has the effect of reducing the radius at which the local expansion balances the gravitational influence of the Milky Way. By ignoring the dynamical effect of the local wall, the mass of the Milky Way estimated from R_0 can be underestimated by as much as ~ 30%.

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

We find that, even in the presence of discreteness noise, a Gaussianizing transform (producing a more-Gaussian one-point distribution) reduces nonlinearities in the power spectra of cosmological matter and galaxy density fields, in many cases drastically. Although Gaussianization does increase the effective shot noise, it also increases the power spectrum's fidelity to the linear power spectrum on scales where the shot noise is negligible. Gaussianizing also increases the Fisher information in the power spectrum in all cases and resolutions, although the gains are smaller in redshift space than in real space. We also find that the gain in cumulative Fisher information from Gaussianizing peaks at a particular grid resolution that depends on the sampling level.

The inclination of galaxies induces both reddening and extinction to their observed spectral energy distribution, which in turn impact the derived properties of the galaxies. Here we report a significant dependence of the error in photometric redshift (photo-z) on the inclination of disk galaxies from the Sloan Digital Sky Survey. The bias in the photo-z based on the template-fitting approach increases from -0.015 in face-on to 0.021 in edge-on galaxies. A Principal Component Analysis on the full sample of photometry reveals the inclination of the galaxies to be represented by the 2nd mode. The corresponding eigenspectrum resembles an extinction curve. The isolation of the inclination effect in a low-order mode demonstrates the significant reddening induced on the observed colors, leading to the over-estimated photo-z in galaxies of high inclinations. We present approaches to correct the photo-z and the other properties of the disk galaxies against the inclination effect.

Morphology is a powerful indicator of a galaxy's dynamical and merger history. It is strongly correlated with many physical parameters, including mass, star formation history and the distribution of mass. The Galaxy Zoo project collected simple morphological classifications of nearly 900,000 galaxies drawn from the Sloan Digital Sky Survey, contributed by hundreds of thousands of volunteers. This large number of classifications allows us to exclude classifier error, and measure the influence of subtle biases inherent in morphological classification. This paper presents the data collected by the project, alongside measures of classification accuracy and bias. The data are now publicly available and full catalogues can be downloaded in electronic format from http://data.galaxyzoo.org.

We show that redshift-space distortions of galaxy correlations have a strong effect on correlation functions with distinct, localized features, like the signature of the baryon acoustic oscillations (BAO). Near the line of sight, the features become sharper as a result of redshift-space distortions. We demonstrate this effect by measuring the correlation function in Gaussian simulations and the Millennium Simulation. We also analyze the SDSS DR7 main-galaxy sample (MGS), splitting the sample into slices 2.5 degrees on the sky in various rotations. Measuring 2D correlation functions in each slice, we do see a sharp bump along the line of sight. Using Mexican-hat wavelets, we localize it to (110 +/- 10) Mpc/h. Averaging only along the line of sight, we estimate its significance at a particular wavelet scale and location at 2.2 sigma. In a flat angular weighting in the (pi,r_p) coordinate system, the noise level is suppressed, pushing the bump's significance to 4 sigma. We estimate that there is about a 0.2% chance of getting such a signal anywhere in the vicinity of the BAO scale from a power spectrum lacking a BAO feature. However, these estimates of the significances make some use of idealized Gaussian simulations, and thus are likely a bit optimistic.

We use GALEX ultraviolet (UV) and optical integrated photometry of the hosts of seventeen luminous supernovae (LSNe, having peak M_V < -21) and compare them to a sample of 26,000 galaxies from a cross-match between the SDSS DR4 spectral catalog and GALEX interim release 1.1. We place the LSNe hosts on the galaxy NUV-r versus M_r color magnitude diagram (CMD) with the larger sample to illustrate how extreme they are. The LSN hosts appear to favor low-density regions of the galaxy CMD falling on the blue edge of the blue cloud toward the low luminosity end. From the UV-optical photometry, we estimate the star formation history of the LSN hosts. The hosts have moderately low star formation rates (SFRs) and low stellar masses (M_*) resulting in high specific star formation rates (sSFR). Compared with the larger sample, the LSN hosts occupy low-density regions of a diagram plotting sSFR versus M_* in the area having higher sSFR and lower M_*. This preference for low M_*, high sSFR hosts implies the LSNe are produced by an effect having to do with their local environment. The correlation of mass with metallicity suggests that perhaps wind-driven mass loss is the factor that prevents LSNe from arising in higher-mass, higher-metallicity hosts. The massive progenitors of the LSNe (>100 M_sun), by appearing in low-SFR hosts, are potential tests for theories of the initial mass function that limit the maximum mass of a star based on the SFR.

We present the SpineWeb framework for the topological analysis of the Cosmic Web and the identification of its walls, filaments and cluster nodes. Based on the watershed segmentation of the cosmic density field, the SpineWeb method invokes the local adjacency properties of the boundaries between the watershed basins to trace the critical points in the density field and the separatrices defined by them. The separatrices are classified into walls and the spine, the network of filaments and nodes in the matter distribution. Testing the method with a heuristic Voronoi model yields outstanding results. Following the discussion of the test results, we apply the SpineWeb method to a set of cosmological N-body simulations. The latter illustrates the potential for studying the structure and dynamics of the Cosmic Web.

We present a novel method for computing the Minkowski Functionals from isodensity surfaces extracted directly from the Delaunay tessellation of a point distribution. This is an important step forward compared to the previous cosmological studies when the isodensity surface was built in the field on a uniform cubic grid and therefore having a uniform spatial resolution. The density field representing a particular interest in cosmology is the density of galaxies which is obtained from the highly nonuniform distribution of the galaxy positions. Therefore, the constraints caused by the spatially uniform grid put severe limitations on the studies of the geometry and shapes of the large-scale objects: superclusters and voids of galaxies. Our technique potentially is able to eliminate most of these limitations. The method is tested with some simple geometric models and an application to the density field from an N-body simulation is shown.

The cross-identification of sources in separate catalogs is one of the most basic tasks in observational astronomy. It is, however, surprisingly difficult and generally ill-defined. Recently Budav\'ari & Szalay (2008) formulated the problem in the realm of probability theory, and laid down the statistical foundations of an extensible methodology. In this paper, we apply their Bayesian approach to stars that, we know, can move measurably on the sky, with detectable proper motion, and show how to associate their observations. We study models on a sample of stars in the Sloan Digital Sky Survey, which allow for an unknown proper motion per object, and demonstrate the improvements over the analytic static model. Our models and conclusions are directly applicable to upcoming surveys such as PanSTARRS, the Dark Energy Survey, Sky Mapper, and the LSST, whose data sets will contain hundreds of millions of stars observed multiple times over several years.

Sky coverage is one of the most important pieces of information about astronomical observations. We discuss possible representations, and present algorithms to create and manipulate shapes consisting of generalized spherical polygons with arbitrary complexity and size on the celestial sphere. This shape specification integrates well with our Hierarchical Triangular Mesh indexing toolbox, whose performance and capabilities are enhanced by the advanced features presented here. Our portable implementation of the relevant spherical geometry routines comes with wrapper functions for database queries, which are currently being used within several scientific catalog archives including the Sloan Digital Sky Survey, the Galaxy Evolution Explorer and the Hubble Legacy Archive projects as well as the Footprint Service of the Virtual Observatory.

We measure the UV-optical color dependence of galaxy clustering in the local universe. Using the clean separation of the red and blue sequences made possible by the NUV - r color-magnitude diagram, we segregate the galaxies into red, blue and intermediate "green" classes. We explore the clustering as a function of this segregation by removing the dependence on luminosity and by excluding edge-on galaxies as a means of a non-model dependent veto of highly extincted galaxies. We find that \xi (r_p, \pi) for both red and green galaxies shows strong redshift space distortion on small scales -- the "finger-of-God" effect, with green galaxies having a lower amplitude than is seen for the red sequence, and the blue sequence showing almost no distortion. On large scales, \xi (r_p, \pi) for all three samples show the effect of large-scale streaming from coherent infall. On scales 1 Mpc/h < r_p < 10 Mpc/h, the projected auto-correlation function w_p(r_p) for red and green galaxies fits a power-law with slope \gamma ~ 1.93 and amplitude r_0 ~ 7.5 and 5.3, compared with \gamma ~ 1.75 and r_0 ~ 3.9 Mpc/h for blue sequence galaxies. Compared to the clustering of a fiducial L* galaxy, the red, green, and blue have a relative bias of 1.5, 1.1, and 0.9 respectively. The w_p(r_p) for blue galaxies display an increase in convexity at ~ 1 Mpc/h, with an excess of large scale clustering. Our results suggest that the majority of blue galaxies are likely central galaxies in less massive halos, while red and green galaxies have larger satellite fractions, and preferentially reside in virialized structures. If blue sequence galaxies migrate to the red sequence via processes like mergers or quenching that take them through the green valley, such a transformation may be accompanied by a change in environment in addition to any change in luminosity and color.

We present the fifth edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog, which is based upon the SDSS Seventh Data Release. The catalog, which contains 105,783 spectroscopically confirmed quasars, represents the conclusion of the SDSS-I and SDSS-II quasar survey. The catalog consists of the SDSS objects that have luminosities larger than M_i = -22.0 (in a cosmology with H_0 = 70 km/s/Mpc Omega_M = 0.3, and Omega_Lambda = 0.7) have at least one emission line with FWHM larger than 1000 km/s or have interesting/complex absorption features, are fainter than i > 15.0 and have highly reliable redshifts. The catalog covers an area of 9380 deg^2. The quasar redshifts range from 0.065 to 5.46, with a median value of 1.49; the catalog includes 1248 quasars at redshifts greater than four, of which 56 are at redshifts greater than five. The catalog contains 9210 quasars with i < 18; slightly over half of the entries have i< 19. For each object the catalog presents positions accurate to better than 0.1" rms per coordinate, five-band (ugriz) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains radio, near-infrared, and X-ray emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3800-9200 Ang. at a spectral resolution R = 2000 the spectra can be retrieved from the SDSS public database using the information provided in the catalog. Over 96% of the objects in the catalog were discovered by the SDSS. We also include a supplemental list of an additional 207 quasars with SDSS spectra whose archive photometric information is incomplete.

We present morphological classifications obtained using machine learning for objects in SDSS DR6 that have been classified by Galaxy Zoo into three classes, namely early types, spirals and point sources/artifacts. An artificial neural network is trained on a subset of objects classified by the human eye and we test whether the machine learning algorithm can reproduce the human classifications for the rest of the sample. We find that the success of the neural network in matching the human classifications depends crucially on the set of input parameters chosen for the machine-learning algorithm. The colours and parameters associated with profile-fitting are reasonable in separating the objects into three classes. However, these results are considerably improved when adding adaptive shape parameters as well as concentration and texture. The adaptive moments, concentration and texture parameters alone cannot distinguish between early type galaxies and the point sources/artifacts. Using a set of twelve parameters, the neural network is able to reproduce the human classifications to better than 90% for all three morphological classes. We find that using a training set that is incomplete in magnitude does not degrade our results given our particular choice of the input parameters to the network. We conclude that it is promising to use machine- learning algorithms to perform morphological classification for the next generation of wide-field imaging surveys and that the Galaxy Zoo catalogue provides an invaluable training set for such purposes.

We study the spectroscopic properties and environments of red spiral galaxies found by the Galaxy Zoo project. By carefully selecting face-on, disk dominated spirals we construct a sample of truly passive disks (not dust reddened, nor dominated by old stellar populations in a bulge). As such, our red spirals represent an interesting set of possible transition objects between normal blue spirals and red early types. We use SDSS data to investigate the physical processes which could have turned these objects red without disturbing their morphology. Red spirals prefer intermediate density regimes, however there are no obvious correlations between red spiral properties and environment - environment alone is not sufficient to determine if a spiral will become red. Red spirals are a small fraction of spirals at low masses, but are a significant fraction at large stellar masses - massive galaxies are red independent of morphology. We confirm that red spirals have older stellar popns and less recent star formation than the main spiral population. While the presence of spiral arms suggests that major star formation cannot have ceased long ago, we show that these are not recent post-starbursts, so star formation must have ceased gradually. Intriguingly, red spirals are ~4 times more likely than normal spirals to host optically identified Seyfert or LINER, with most of the difference coming from LINERs. We find a curiously large bar fraction in the red spirals suggesting that the cessation of star formation and bar instabilities are strongly correlated. We conclude by discussing the possible origins. We suggest they may represent the very oldest spiral galaxies which have already used up their reserves of gas - probably aided by strangulation, and perhaps bar instabilities moving material around in the disk.

We investigate the effect of dust on spiral galaxies by measuring the inclination-dependence of optical colours for 24,276 well-resolved SDSS galaxies visually classified in Galaxy Zoo. We find clear trends of reddening with inclination which imply a total extinction from face-on to edge-on of 0.7, 0.6, 0.5 and 0.4 magnitudes for the ugri passbands. We split the sample into "bulgy" (early-type) and "disky" (late-type) spirals using the SDSS fracdeV (or f_DeV) parameter and show that the average face-on colour of "bulgy" spirals is redder than the average edge-on colour of "disky" spirals. This shows that the observed optical colour of a spiral galaxy is determined almost equally by the spiral type (via the bulge-disk ratio and stellar populations), and reddening due to dust. We find that both luminosity and spiral type affect the total amount of extinction, with "disky" spirals at M_r ~ -21.5 mags having the most reddening. This decrease of reddening for the most luminous spirals has not been observed before and may be related to their lower levels of recent star formation. We compare our results with the latest dust attenuation models of Tuffs et al. We find that the model reproduces the observed trends reasonably well but overpredicts the amount of u-band attenuation in edge-on galaxies. We end by discussing the effects of dust on large galaxy surveys and emphasize that these effects will become important as we push to higher precision measurements of galaxy properties and their clustering.

We present a framework for the hierarchical identification and characterization of voids based on the Watershed Void Finder. The Hierarchical Void Finder is based on a generalization of the scale space of a density field invoked in order to trace the hierarchical nature and structure of cosmological voids. At each level of the hierarchy, the watershed transform is used to identify the voids at that particular scale. By identifying the overlapping regions between watershed basins in adjacent levels, the hierarchical void tree is constructed. Applications on a hierarchical Voronoi model and on a set of cosmological simulations illustrate its potential.

We use data from the Sloan Digital Sky Survey and visual classifications of morphology from the Galaxy Zoo project to study black hole growth in the nearby Universe (z < 0.05) and to break down the AGN host galaxy population by color, stellar mass and morphology. We find that black hole growth at luminosities L_OIII >1E40 erg/s in early- and late-type galaxies is fundamentally different. AGN host galaxies as a population have a broad range of stellar masses (1E10-1E11 Msun), reside in the green valley of the color-mass diagram and their central black holes have median masses around 1E6.5 Msun. However, by comparing early- and late-type AGN host galaxies to their non-active counterparts, we find several key differences: in early-type galaxies, it is preferentially the galaxies with the least massive black holes that are growing, while late-type galaxies, it is preferentially the most massive}black holes that are growing. The duty cycle of AGN in early-type galaxies is strongly peaked in the green valley below the low-mass end (1E10 Msun) of the red sequence at stellar masses where there is a steady supply of blue cloud progenitors. The duty cycle of AGN in late-type galaxies on the other hand peaks in massive (1E11 Msun) green and red late-types which generally do not have a corresponding blue cloud population of similar mass. At high Eddington ratios (L/L_Edd > 0.1), the only population with a substantial fraction of AGN are the low-mass green valley early-type galaxies. Finally, the Milky Way likely resides in the "sweet spot" on the color-mass diagram where the AGN duty cycle of late-type galaxies is highest. We discuss the implications of these results for our understanding of the role of AGN in the evolution of galaxies

We report on the finding of a correlation between galaxies' past star formation activity and the degree to which neighbouring galaxies rotation axes are aligned. This is obtained by cross-correlating star formation histories, derived with MOPED, and spin direction (chirality), as determined by the Galaxy Zoo project, for a sample of SDSS galaxies. Our findings suggest that spiral galaxies which formed the majority of their stars early (z > 2) tend to display coherent rotation over scales of ~10 Mpc/h. The correlation is weaker for galaxies with significant recent star formation. We find evidence for this alignment at more than the 5-sigma level, but no correlation with other galaxy stellar properties. This finding can be explained within the context of hierarchical tidal-torque theory if the SDSS galaxies harboring the majority of the old stellar population where formed in the past, in the same filament and at about the same time. Galaxies with significant recent star formation instead are in the field, thus influenced by the general tidal field that will align them in random directions or had a recent merger which would promote star formation, but deviate the spin direction.

Extinction in galaxies affects their observed properties. In scenarios describing the distribution of dust and stars in individual disk galaxies, the amplitude of the extinction can be modulated by the inclination of the galaxies. In this work we investigate the inclination dependency in composite spectra of star-forming disk galaxies from the Sloan Digital Sky Survey Data Release 5. In a volume-limited sample within a redshift range 0.065-0.075 and a r-band Petrosian absolute magnitude range -19.5 to -22 which exhibits a flat distribution of inclination, the inclined relative to face-on extinction in the stellar continuum is found empirically to increase with inclination in the g, r, and i bands. Within the central 0.5 intrinsic half-light radius of the galaxies, the g-band relative extinction in the stellar continuum for the highly-inclined objects (axis ratio b/a = 0.1) is 1.2 mag, agreeing with previous studies. The extinction curve of the disk galaxies is given in the restframe wavelengths 3700-8000 angstrom, identified with major optical emission and absorption lines in diagnostics. The Balmer decrement remains constant with inclination, suggesting a different kind of dust configuration and/or reddening mechanism in the HII region from that in the stellar continuum. One factor is shown to be the presence of spatially non-uniform interstellar extinction, presumably caused by clumped dust in the vicinity of the HII region.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We use multi-wavelength, matched aperture, integrated photometry from GALEX, SDSS and the RC3 to estimate the physical properties of 166 nearby galaxies hosting 168 well-observed Type Ia supernovae (SNe Ia). Our data corroborate well-known features that have been seen in other SN Ia samples. Specifically, hosts with active star formation produce brighter and slower SNe Ia on average, and hosts with luminosity-weighted ages older than 1 Gyr produce on average more faint, fast and fewer bright, slow SNe Ia than younger hosts. New results include that in our sample, the faintest and fastest SNe Ia occur only in galaxies exceeding a stellar mass threshhold of ~10^10 M_sun, indicating that their progenitors must arise in populations that are older and/or more metal rich than the general SN Ia population. A low host extinction sub-sample hints at a residual trend in peak luminosity with host age, after correcting for light-curve shape, giving the appearance that older hosts produce less-extincted SNe Ia on average. This has implications for cosmological fitting of SNe Ia and suggests that host age could be useful as a parameter in the fitting. Converting host mass to metallicity and computing 56Ni mass from the supernova light curves, we find that our local sample is consistent with a model that predicts a shallow trend between stellar metallicity and the 56Ni mass that powers the explosion, but we cannot rule out the absence of a trend. We measure a correlation between 56Ni mass and host age in the local universe that is shallower and not as significant as that seen at higher redshifts. The details of the age -- 56Ni mass correlations at low and higher redshift imply a luminosity-weighted age threshhold of ~3 Gyr for SN Ia hosts, above which they are less likely to produce SNe Ia with 56Ni masses above ~0.5 M_sun. (Abridged)

The spectroscopic Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) galaxy sample represents the final set of galaxies observed using the original SDSS target selection criteria. We analyse the clustering of galaxies within this sample, including both the Luminous Red Galaxy (LRG) and Main samples, and also include the 2-degree Field Galaxy Redshift Survey (2dFGRS) data. Baryon Acoustic Oscillations are observed in power spectra measured for different slices in redshift; this allows us to constrain the distance--redshift relation at multiple epochs. We achieve a distance measure at redshift z=0.275, of r_s(z_d)/D_V(0.275)=0.1390+/-0.0037 (2.7% accuracy), where r_s(z_d) is the comoving sound horizon at the baryon drag epoch, D_V(z)=[(1+z)^2D_A^2cz/H(z)]^(1/3), D_A(z) is the angular diameter distance and H(z) is the Hubble parameter. We find an almost independent constraint on the ratio of distances D_V(0.35)/D_V(0.2)=1.736+/-0.065, which is consistent at the 1.1sigma level with the best fit Lambda-CDM model obtained when combining our z=0.275 distance constraint with the WMAP 5-year data. The offset is similar to that found in previous analyses of the SDSS DR5 sample, but the discrepancy is now of lower significance, a change caused by a revised error analysis and a change in the methodology adopted, as well as the addition of more data. Using WMAP5 constraints on Omega_bh^2 and Omega_ch^2, and combining our BAO distance measurements with those from the Union Supernova sample, places a tight constraint on Omega_m=0.286+/-0.018 and H_0 = 68.2+/-2.2km/s/Mpc that is robust to allowing curvature and non-Lambda dark energy. This result is independent of the behaviour of dark energy at redshifts greater than those probed by the BAO and supernova measurements. (abridged)

Following the study of Darg et al. (2009; hereafter D09a) we explore the environments, optical colours, stellar masses, star formation and AGN activity in a sample of 3003 pairs of merging galaxies drawn from the SDSS using visual classifications from the Galaxy Zoo project. While D09a found that the spiral-to-elliptical ratio in (major) mergers appeared higher than that of the global galaxy population, no significant differences are found between the environmental distributions of mergers and a randomly selected control sample. This makes the high occurrence of spirals in mergers unlikely to be an environmental effect and must, therefore, arise from differing time-scales of detectability for spirals and ellipticals. We find that merging galaxies have a wider spread in colour than the global galaxy population, with a significant blue tail resulting from intense star formation in spiral mergers. Galaxies classed as star-forming using their emission-line properties have average star-formation rates approximately doubled by the merger process though star formation is negligibly enhanced in merging elliptical galaxies. We conclude that the internal properties of galaxies significantly affect the time-scales over which merging systems can be detected (as suggested by recent theoretical studies) which leads to spirals being `over-observed' in mergers. We also suggest that the transition mass $3\times10^{10}{M}_{\astrosun}$, noted by \citet{kauffmann1}, below which ellipticals are rare could be linked to disc survival/destruction in mergers.

We investigate the quality of associations of astronomical sources from multi-wavelength observations using simulated detections that are realistic in terms of their astrometric accuracy, small-scale clustering properties and selection functions. We present a general method to build such mock catalogs for studying associations, and compare the statistics of cross-identifications based on angular separation and Bayesian probability criteria. In particular, we focus on the highly relevant problem of cross-correlating the ultraviolet Galaxy Evolution Explorer (GALEX) and optical Sloan Digital Sky Survey (SDSS) surveys. Using refined simulations of the relevant catalogs, we find that the probability thresholds yield lower contamination of false associations, and are more efficient than angular separation. Our study presents a set of recommended criteria to construct reliable cross-match catalogs between SDSS and GALEX with minimal artifacts.

The Galaxy Zoo citizen science website invites anyone with an Internet connection to participate in research by classifying galaxies from the Sloan Digital Sky Survey. As of April 2009, more than 200,000 volunteers had made more than 100 million galaxy classifications. In this paper, we present results of a pilot study into the motivations and demographics of Galaxy Zoo volunteers, and define a technique to determine motivations from free responses that can be used in larger multiple-choice surveys with similar populations. Our categories form the basis for a future survey, with the goal of determining the prevalence of each motivation.

We present the largest, most homogeneous catalogue of merging galaxies in the nearby universe obtained through the Galaxy Zoo project - an interface on the world-wide web enabling large-scale morphological classification of galaxies through visual inspection of images from the Sloan Digital Sky Survey (SDSS). The method converts a set of visually-inspected classifications for each galaxy into a single parameter (the `weighted-merger-vote fraction,' $f_m$) which describes our confidence that the system is part of an ongoing merger. We describe how $f_m$ is used to create a catalogue of 3003 visually-selected pairs of merging galaxies from the SDSS in the redshift range $0.005 < z <0.1$. We use our merger sample and values of $f_m$ applied to the SDSS Main Galaxy Spectral sample (MGS) to estimate that the fraction of volume-limited ($M_r < -20.55$) major mergers ($1/3 < {M}^*_1/{M}^*_2 < 3$) in the nearby universe is $1 - 3 \times C%$ where $C \sim 1.5$ is a correction factor for spectroscopic incompleteness. Having visually classified the morphologies of the constituent galaxies in our mergers, we find that the spiral-to-elliptical ratio of galaxies in mergers is higher by a factor $\sim 2$ relative to the global population. In a companion paper, we examine the internal properties of these merging galaxies and conclude that this high spiral-to-elliptical ratio in mergers is due to a longer time-scale over which mergers with spirals are detectable compared to mergers with ellipticals.

We examine the correlation function \xi of the Sloan Digital Sky Survey (SDSS) Luminous Red Galaxy sample (LRG) at large scales (60<s<400 Mpc/h) using the final data release (DR7; 105,831 LRGs between 0.16<z<0.47). Using mock catalogs, we demonstrate that the observed baryonic acoustic peak and larger scale signal are consistent with LCDM at the 1.5\sigma level. The signal at 155<s<200 Mpc/h tends to be high relative to theoretical expectations; this slight deviation can be attributed to a bright subsample of the LRGs. Fitting data to a non-linear, redshift-space, template based-model, we constrain the peak position at s_p=103.6+3.6-2.4 Mpc/h when fitting the range 60<s<150 Mpc/h (1\sigma uncertainties measured from the mocks. This redshift-space distance s_p is related to the comoving sound horizon scale r_s after taking into account matter clustering non-linearities, redshift distortions and galaxy clustering bias. Mock catalogs show that the probability that a DR7-sized sample would not have an identifiable peak is at least 10%. As a consistency check of a fiducial cosmology, we use the observed s_p to obtain the distance D_V=[(1+z)^2D_A^2cz/H(z)]^(1/3) relative to the acoustic scale. We find r_s/D_V(z=0.278)=0.1394+-0.0049. This result is in excellent agreement with Percival et. al (2009), who examine roughly the same data set, but using the power spectrum. Comparison with other determinations in the literature are also in very good agreement. We have tested our results against a battery of possible systematic effects, finding all effects are smaller than our estimated sample variance.

We investigate a class of rapidly growing emission line galaxies, known as "Green Peas", first noted by volunteers in the Galaxy Zoo project because of their peculiar bright green colour and small size, unresolved in SDSS imaging. Their appearance is due to very strong optical emission lines, namely [O III] 5007 A, with an unusually large equivalent width of up to ~1000 A. We discuss a well-defined sample of 251 colour-selected objects, most of which are strongly star forming, although there are some AGN interlopers including 8 newly discovered narrow Line Seyfert 1 galaxies. The star-forming Peas are low mass galaxies (M~10^8.5 - 10^10 M_sun) with high star formation rates (~10 M_sun/yr), low metallicities (log[O/H] + 12 ~ 8.7) and low reddening (E(B-V) < 0.25) and they reside in low density environments. They have some of the highest specific star formation rates (up to ~10^{-8} yr^{-1}) seen in the local Universe, yielding doubling times for their stellar mass of hundreds of Myrs. The few star-forming Peas with HST imaging appear to have several clumps of bright star-forming regions and low surface density features that may indicate recent or ongoing mergers. The Peas are similar in size, mass, luminosity and metallicity to Luminous Blue Compact Galaxies. They are also similar to high redshift UV-luminous galaxies, e.g., Lyman-break galaxies and Lyman-alpha emitters, and therefore provide a local laboratory with which to study the extreme star formation processes that occur in high-redshift galaxies. Studying starbursting galaxies as a function of redshift is essential to understanding the build up of stellar mass in the Universe.

We present the power spectrum of the reconstructed halo density field derived from a sample of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey Seventh Data Release (DR7). The halo power spectrum has a direct connection to the underlying dark matter power for k <= 0.2 h/Mpc, well into the quasi-linear regime. This enables us to use a factor of ~8 more modes in the cosmological analysis than an analysis with kmax = 0.1 h/Mpc, as was adopted in the SDSS team analysis of the DR4 LRG sample (Tegmark et al. 2006). The observed halo power spectrum for 0.02 < k < 0.2 h/Mpc is well-fit by our model: chi^2 = 39.6 for 40 degrees of freedom for the best fit LCDM model. We find \Omega_m h^2 * (n_s/0.96)^0.13 = 0.141^{+0.009}_{-0.012} for a power law primordial power spectrum with spectral index n_s and \Omega_b h^2 = 0.02265 fixed, consistent with CMB measurements. The halo power spectrum also constrains the ratio of the comoving sound horizon at the baryon-drag epoch to an effective distance to z=0.35: r_s/D_V(0.35) = 0.1097^{+0.0039}_{-0.0042}. Combining the halo power spectrum measurement with the WMAP 5 year results, for the flat LCDM model we find \Omega_m = 0.289 +/- 0.019 and H_0 = 69.4 +/- 1.6 km/s/Mpc. Allowing for massive neutrinos in LCDM, we find \sum m_{\nu} < 0.62 eV at the 95% confidence level. If we instead consider the effective number of relativistic species Neff as a free parameter, we find Neff = 4.8^{+1.8}_{-1.7}. Combining also with the Kowalski et al. (2008) supernova sample, we find \Omega_{tot} = 1.011 +/- 0.009 and w = -0.99 +/- 0.11 for an open cosmology with constant dark energy equation of state w.

We analyze the environmental dependence of galaxy morphology and colour with two-point clustering statistics, using data from the Galaxy Zoo, the largest sample of visually classified morphologies yet compiled, extracted from the Sloan Digital Sky Survey. We present two-point correlation functions of spiral and early-type galaxies, and we quantify the correlation between morphology and environment with marked correlation functions. These yield clear and precise environmental trends across a wide range of scales, analogous to similar measurements with galaxy colours, indicating that the Galaxy Zoo classifications themselves are very precise. We measure morphology marked correlation functions at fixed colour and find that they are relatively weak, with the only residual correlation being that of red galaxies at small scales, indicating a morphology gradient within haloes for red galaxies. At fixed morphology, we find that the environmental dependence of colour remains strong, and these correlations remain for fixed morphology \textit{and} luminosity. An implication of this is that much of the morphology--density relation is due to the relation between colour and density. Our results also have implications for galaxy evolution: the morphological transformation of galaxies is usually accompanied by a colour transformation, but not necessarily vice versa. A spiral galaxy may move onto the red sequence of the colour-magnitude diagram without quickly becoming an early-type. We analyze the significant population of red spiral galaxies, and present evidence that they tend to be located in moderately dense environments and are often satellite galaxies in the outskirts of haloes. Finally, we combine our results to argue that central and satellite galaxies tend to follow different evolutionary paths.

We report the discovery of an unusual object near the spiral galaxy IC 2497, discovered by visual inspection of the Sloan Digital Sky Survey (SDSS) as part of the Galaxy Zoo project. The object, known as Hanny's Voorwerp, is bright in the SDSS g band due to unusually strong OIII 4959-5007 emission lines. We present the results of the first targeted observations of the object in the optical, UV and X-ray, which show that the object contains highly ionized gas. Although the line ratios are similar to extended emission-line regions near luminous AGN, the source of this ionization is not apparent. The emission-line properties, and lack of x-ray emission from IC 2497, suggest either a highly obscured AGN with a novel geometry arranged to allow photoionization of the object but not the galaxy's own circumnuclear gas, or, as we argue, the first detection of a quasar light echo. In this case, either the luminosity of the central source has decreased dramatically or else the obscuration in the system has increased within 10^5 years. This object may thus represent the first direct probe of quasar history on these timescales.

We explore the age distribution of the globular cluster (GC) system of the nearby elliptical galaxy NGC 5128 using ultraviolet (UV) photometry from Galaxy Evolution Explorer (GALEX) observations, with UV - optical colors used as the age indicator. Most GCs in NGC 5128 follow the general trends of GCs in M31 and Milky Way in UV - optical color-color diagram, which indicates that the majority of GCs in NGC 5128 are old similar to the age range of old GCs in M31 and Milky Way. A large fraction of spectroscopically identified intermediate-age GC (IAGC) candidates with ~ 3-8 Gyr are not detected in the FUV passband. Considering the nature of intermediate-age populations being faint in the far-UV (FUV) passband, we suggest that many of the spectroscopically identified IAGCs may be truly intermediate in age. This is in contrast to the case of M31 where a large fraction of spectroscopically suggested IAGCs are detected in FUV and therefore may not be genuine IAGCs but rather older GCs with developed blue horizontal branch stars. Our UV photometry strengthens the results previously suggesting the presence of GC and stellar subpopulation with intermediate age in NGC 5128. The existence of IAGCs strongly indicates the occurrence of at least one more major star formation episode after a starburst at high redshift.

Studying the environments of 0.4<z<1.2 UV-selected galaxies, as examples of extreme star-forming galaxies (with star formation rates in the range of 3-30 M_sol/yr), we explore the relationship between high rates of star-formation, host halo mass and pair fractions. We study the large-scale and small-scale environments of local Ultraviolet Luminous Galaxies (UVLGs) by measuring angular correlation functions. We cross-correlate these systems with other galaxy samples: a volume-limited sample (ALL), a Blue Luminous Galaxy sample (BLG) and a Luminous Red Galaxy sample (LRG). We determine the UVLG comoving correlation length to be r_0=4.8(+11.6/-2.4) h^-1 Mpc at <z> =1.0, which is unable to constrain the halo mass for this sample. However, we find that UVLGs form close (separation < 30 kpc) pairs with the ALL sample, but do not frequently form pairs with LRGs. A rare subset of UVLGs, those with the highest FUV surface brightnesses, are believed to be local analogs of high redshift Lyman Break Galaxies (LBGs) and are called Lyman Break Analogs (LBAs). LBGs and LBAs share similar characteristics (i.e., color, size, surface brightness, specific star formation rates, metallicities, and dust content). Recent HST images of z~0.2 LBAs show disturbed morphologies, signs of mergers and interactions. UVLGs may be influenced by interactions with other galaxies and we discuss this result in terms of other high star-forming, merging systems.

We find that nonlinearities in the dark-matter power spectrum are dramatically smaller if the density field first undergoes a logarithmic mapping. In the Millennium simulation, this procedure produces a power spectrum with a shape hardly departing from the linear power spectrum for k <~ 1 h/Mpc at all redshifts. Also, this procedure unveils pristine Fisher information on a range of scales reaching a factor of 2-3 smaller than in the standard power spectrum, yielding 10 times more cumulative signal-to-noise at z=0.

We measure the projected spatial correlation function w_p(r_p) from a large sample combining GALEX ultraviolet imaging with the SDSS spectroscopic sample. We study the dependence of the clustering strength for samples selected on (NUV - r)_abs color, specific star formation rate (SSFR), and stellar mass. We find that there is a smooth transition in the clustering of galaxies as a function of this color from weak clustering among blue galaxies to stronger clustering for red galaxies. The clustering of galaxies within the "green valley" has an intermediate strength, and is consistent with that expected from galaxy groups. The results are robust to the correction for dust extinction. The comparison with simple analytical modeling suggests that the halo occupation number increases with older star formation epochs. When splitting according to SSFR, we find that the SSFR is a more sensitive tracer of environment than stellar mass.

We present a detailed study of the Galaxy Evolution Explorer's photometric catalogs with special focus on the statistical properties of the All-sky and Medium Imaging Surveys. We introduce the concept of primaries to resolve the issue of multiple detections and follow a geometric approach to define clean catalogs with well-understood selection functions. We cross-identify the GALEX sources (GR2+3) with Sloan Digital Sky Survey (DR6) observations, which indirectly provides an invaluable insight about the astrometric model of the UV sources and allows us to revise the band merging strategy. We derive the formal description of the GALEX footprints as well as their intersections with the SDSS coverage along with analytic calculations of their areal coverage. The crossmatch catalogs are made available for the public. We conclude by illustrating the implementation of typical selection criteria in SQL for catalog subsets geared toward statistical analyses, e.g., correlation and luminosity function studies.

We present a novel technique to overcome the limitations of the applicability of Principal Component Analysis to typical real-life data sets, especially astronomical spectra. Our new approach addresses the issues of outliers, missing information, large number of dimensions and the vast amount of data by combining elements of robust statistics and recursive algorithms that provide improved eigensystem estimates step-by-step. We develop a generic mechanism for deriving reliable eigenspectra without manual data censoring, while utilising all the information contained in the observations. We demonstrate the power of the methodology on the attractive collection of the VIMOS VLT Deep Survey spectra that manifest most of the challenges today, and highlight the improvements over previous workarounds, as well as the scalability of our approach to collections with sizes of the Sloan Digital Sky Survey and beyond.

Under the unified model for active galactic nuclei (AGNs), narrow-line (Type 2) AGNs are, in fact, broad-line (Type 1) AGNs but each with a heavily obscured accretion disk. We would therefore expect the optical continuum emission from Type 2 AGN to be composed mainly of stellar light and non-variable on the time-scales of months to years. In this work we probe the spectroscopic variability of galaxies and narrow-line AGNs using the multi-epoch data in the Sloan Digital Sky Survey (SDSS) Data Release 6. The sample contains 18,435 sources for which there exist pairs of spectroscopic observations (with a maximum separation in time of ~700 days) covering a wavelength range of 3900-8900 angstrom. To obtain a reliable repeatability measurement between each spectral pair, we consider a number of techniques for spectrophotometric calibration resulting in an improved spectrophotometric calibration of a factor of two. From these data we find no obvious continuum and emission-line variability in the narrow-line AGNs on average -- the spectroscopic variability of the continuum is 0.07+/-0.26 mag in the g band and, for the emission-line ratios log10([NII]/Halpha) and log10([OIII]/Hbeta), the variability is 0.02+/-0.03 dex and 0.06+/-0.08 dex, respectively. From the continuum variability measurement we set an upper limit on the ratio between the flux of varying spectral component, presumably related to AGN activities, and that of host galaxy to be ~30%. We provide the corresponding upper limits for other spectral classes, including those from the BPT diagram, eClass galaxy classification, stars and quasars.

Over the last decade optical spectroscopic surveys have characterized the low redshift galaxy population and uncovered populations of star-forming galaxies back to z ~ 7. This work has shown that the primary epoch of galaxy building and black hole growth occurs at redshifts of 2 to 3. The establishment of the concordance LCDM cosmology shifted the focus of galaxy population studies from constraining cosmological parameters to characterizing the processes which regulate the formation and evolution of galaxies.In the next decade, high redshift observers will attempt to formulate a coherent evolutionary picture connecting galaxies in the high redshift Universe to galaxies today. In order to link galaxy populations at different redshifts, we must not only characterize their evolution in a systematic way, we must establish which physical processes are responsible for it. Considerable progress has already been made in understanding how galaxies evolved from z ~ 1 to the present day. Large spectroscopic surveys in the near infrared are required to push these studies back towards the main epoch of galaxy building. Only then will we understand the full story of the formation of L* galaxies like our own Milky Way. A large near-IR spectroscopic survey will also provide the calibration needed to avoid systematics in the large photometric programs proposed to study the nature of dark matter and dark energy. We provide an outline design for a multi-object 0.4 to 1.8 micron spectrograph, which could be placed on an existing telescope, and which would allow a full characterization of the galaxy population out to z ~ 2. We strongly recommend a serious further study to design a real instrument, which will be required for galaxy formation studies to advance to the next frontier.

We report the discovery of a population of nearby, blue early-type galaxies with high star formation rates (0.5 < SFR < 50 Msun/yr). They are identified by their visual morphology as provided by Galaxy Zoo for SDSS DR6 and their u-r colour. We select a volume-limited sample in the redshift range 0.02 < z < 0.05, corresponding to luminosities of approximately L* and above, and with u-r colours significantly bluer than the red sequence. We confirm the early-type morphology of the objects in this sample and investigate their environmental dependence and star formation properties. Blue early-type galaxies tend to live in lower-density environments than `normal' red sequence early-types and make up 5.7 +/-0.4% of the low-redshift early-type galaxy population. We find that such blue early-type galaxies are virtually absent at high velocity dispersions above 200 km/s. Our analysis uses emission line diganostic diagrams and we find that ~25% of them are actively starforming, while another ~25% host both star formation and an AGN. Another ~12% are AGN. The remaining 38% show no strong emission lines. When present and uncontaminated by an AGN contribution, the star formation is generally intense. We consider star formation rates derived from Halpha, u-band and infrared luminosities, and radial colour profiles, and conclude that the star formation is spatially extended. Of those objects that are not currently undergoing star formation must have ceased doing so recently in order to account for their blue optical colours. The gas phase metallicity of the actively starforming blue early-types galaxies is supersolar in all cases. We discuss the place of these objects in the context of galaxy formation. A catalogue of all 204 blue early-type galaxies in our sample, including star formation rates and emission line classification, is provided.

We investigate the nature of the star formation law at low gas surface densities using a sample of 19 low surface brightness (LSB) galaxies with existing HI maps in the literature, UV imaging from the Galaxy Evolution Explorer satellite, and optical images from the Sloan Digital Sky Survey. All of the LSB galaxies have (NUV-r) colors similar to those for higher surface brightness star-forming galaxies of similar luminosity indicating that their average star formation histories are not very different. Based upon four LSB galaxies with both UV and FIR data, we find FIR/UV ratios significantly less than one, implying low amounts of internal UV extinction in LSB galaxies. We use the UV images and HI maps to measure the star formation rate and hydrogen gas surface densities within the same region for all of the galaxies. The LSB galaxy star formation rate surface densities lie below the extrapolation of the power law fit to the star formation rate surface density as a function of the total gas density for higher surface brightness galaxies. Although there is more scatter, the LSB galaxies also lie below a second version of the star formation law in which the star formation rate surface density is correlated with the gas density divided by the orbital time in the disk. The downturn seen in both star formation laws is consistent with theoretical models that predict lower star formation efficiencies in LSB galaxies due to the declining molecular fraction with decreasing density.

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 analyse the relationships between galaxy morphology, colour, environment and stellar mass using data for over 100,000 objects from Galaxy Zoo, the largest sample of visually classified morphologies yet compiled. We conclusively show that colour and morphology fractions are very different functions of environment. Both are sensitive to stellar mass; however, at fixed stellar mass, while colour is also highly sensitive to environment, morphology displays much weaker environmental trends. Only a small part of both relations can be attributed to variation in the stellar mass function with environment. Galaxies with high stellar masses are mostly red, in all environments and irrespective of their morphology. Low stellar-mass galaxies are mostly blue in low-density environments, but mostly red in high-density environments, again irrespective of their morphology. The colour-density relation is primarily driven by variations in colour fractions at fixed morphology, in particular the fraction of spiral galaxies that have red colours, and especially at low stellar masses. We demonstrate that our red spirals primarily include galaxies with true spiral morphology. We clearly show there is an environmental dependence for colour beyond that for morphology. Before using the Galaxy Zoo morphologies to produce the above results, we first quantify a luminosity-, size- and redshift-dependent classification bias that affects this dataset, and probably most other studies of galaxy population morphology. A correction for this bias is derived and applied to produce a sample of galaxies with reliable morphological type likelihoods, on which we base our analysis.

We describe features of the LSST science database that are amenable to scientific data mining, object classification, outlier identification, anomaly detection, image quality assurance, and survey science validation. The data mining research agenda includes: scalability (at petabytes scales) of existing machine learning and data mining algorithms; development of grid-enabled parallel data mining algorithms; designing a robust system for brokering classifications from the LSST event pipeline (which may produce 10,000 or more event alerts per night); multi-resolution methods for exploration of petascale databases; indexing of multi-attribute multi-dimensional astronomical databases (beyond spatial indexing) for rapid querying of petabyte databases; and more.

Galaxy Zoo is the first study of nearby galaxies that contains reliable information about the spiral sense of rotation of galaxy arms for a sizeable number of galaxies. We measure the correlation function of spin chirality (the sense in which galaxies appear to be spinning) of face-on spiral galaxies in angular, real and projected spaces. Our results indicate a hint of positive correlation at separations less than ~0.5 Mpc at a statistical significance of 2-3 sigma. This is the first experimental evidence for chiral correlation of spins. Within tidal torque theory it indicates that the inertia tensors of nearby galaxies are correlated. This is complementary to the studies of nearby spin axis correlations that probe the correlations of the tidal field. Theoretical interpretation is made difficult by the small distances at which the correlations are detected, implying that substructure might play a significant role, and our necessary selection of face-on spiral galaxies, rather than a general volume-limited sample.

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.

We re-examine the evidence for a violation of large-scale statistical isotropy in the distribution of projected spin vectors of spiral galaxies. We have a sample of $\sim 37,000$ spiral galaxies from the Sloan Digital Sky Survey, with their line of sight spin direction confidently classified by members of the public through the online project Galaxy Zoo. After establishing and correcting for a certain level of bias in our handedness results we find the winding sense of the galaxies to be consistent with statistical isotropy. In particular we find no significant dipole signal, and thus no evidence for overall preferred handedness of the Universe. We compare this result to those of other authors and conclude that these may also be affected and explained by a bias effect.

In order to understand the formation and subsequent evolution of galaxies one must first distinguish between the two main morphological classes of massive systems: spirals and early-type systems. This paper introduces a project, Galaxy Zoo, which provides visual morphological classifications for nearly one million galaxies, extracted from the Sloan Digital Sky Survey (SDSS). This achievement was made possible by inviting the general public to visually inspect and classify these galaxies via the internet. The project has obtained more than 40,000,000 individual classifications made by ~100,000 participants. We discuss the motivation and strategy for this project, and detail how the classifications were performed and processed. We find that Galaxy Zoo results are consistent with those for subsets of SDSS galaxies classified by professional astronomers, thus demonstrating that our data provides a robust morphological catalogue. Obtaining morphologies by direct visual inspection avoids introducing biases associated with proxies for morphology such as colour, concentration or structual parameters. In addition, this catalogue can be used to directly compare SDSS morphologies with older data sets. The colour--magnitude diagrams for each morphological class are shown, and we illustrate how these distributions differ from those inferred using colour alone as a proxy for morphology.

We study the large-scale anisotropic two-point correlation function using 46,760 luminous red galaxies at redshifts 0.16 -- 0.47 from the Sloan Digital Sky Survey. We measure the correlation function as a function of separations parallel and perpendicular to the line-of-sight in order to take account of anisotropy of the large-scale structure in redshift space. We find a slight signal of baryonic features in the anisotropic correlation function, i.e., a ``baryon ridge'' which corresponds to a baryon acoustic peak in the spherically averaged correlation function which has already been reported using the same sample. The baryon ridge has primarily a spherical structure with a known radius in comoving coordinates. It enables us to divide the redshift distortion effects into dynamical and geometrical components and provides further constraints on cosmological parameters, including the dark energy equation-of-state. With an assumption of a flat $\Lambda$ cosmology, we find the best-fit values of $\Omega_{\rm m} = 0.218^{+0.047}_{-0.037}$ and $\Omega_{\rm b} = 0.047^{+0.016}_{-0.016}$ (68% C.L.) when we use the overall shape of the anisotropic correlation function of $40<s<200\himpc$ including a scale of baryon acoustic oscillations. When an additional assumption $\Omega_{\rm b}h^2=0.024$ is adopted, we obtain $\Omega_{\rm DE}=0.770^{+0.051}_{-0.040}$ and $w=-0.93^{+0.45}_{-0.35}$. These constraints are estimated only from our data of the anisotropic correlation function, and they agree quite well with values both from the cosmic microwave background (CMB) anisotropies and from other complementary statistics using the LRG sample. With the CMB prior from the 3 year WMAP results, we give stronger constraints on those parameters.

The GALEX (Galaxy Evolution Explorer) spectroscopic survey mode, with a resolution of about 8 A in the FUV (1350 - 1750 A) and about 20 A in the NUV (1950 - 2750 A) is used for a systematic search of Ly-a emitting galaxies at low redshift. This aims at filling a gap between high-redshift surveys and a small set of objects studied in detail in the nearby universe. A blind search of 7018 spectra extracted in 5 deep exposures (5.65 sq.deg) has resulted in 96 Ly-a emitting galaxy candidates in the FUV domain, after accounting for broad-line AGNs. The Ly-a EWs (equivalent width) are consistent with stellar population model predictions and show no trends as a function of UV color or UV luminosity, except a possible decrease in the most luminous that may be due to small-number statistics. Their distribution in EW is similar to that at z about 3 but their fraction among star-forming galaxies is smaller. Avoiding uncertain candidates, a sub-sample of 66 objects in the range 0.2 < z < 0.35 has been used to build a Ly-a LF (luminosity function). The incompleteness due to objects with significant Ly-a emission but a UV continuum too low for spectral extraction has been evaluated. A comparison with H-a LF in the same redshift domain is consistent with an average Ly-a/H-a of about 1 in about 15 % of the star-forming galaxies. A comparison with high-redshift Ly-a LFs implies an increase of the Ly-a luminosity density by a factor of about 16 from z about 0.3 to z about 3. By comparison with the factor 5 increase of the UV luminosity density in the same redshift range, this suggests an increase of the average Ly-a escape fraction with redshift.

We present a general probabilistic formalism for cross-identifying astronomical point sources in multiple observations. Our Bayesian approach, symmetric in all observations, is the foundation of a unified framework for object matching, where not only spatial information, but physical properties, such as colors, redshift and luminosity, can also be considered in a natural way. We provide a practical recipe to implement an efficient recursive algorithm to evaluate the Bayes factor over a set of catalogs with known circular errors in positions. This new methodology is crucial for studies leveraging the synergy of today's multi-wavelength observations and to enter the time-domain science of the upcoming survey telescopes.

The meaningful comparison of models of galaxy evolution to observations is critically dependent on the accurate treatment of dust attenuation. To investigate dust absorption and emission in galaxies we have assembled a sample of ~1000 galaxies with ultraviolet (UV) through infrared (IR) photometry from GALEX, SDSS, and Spitzer and optical spectroscopy from SDSS. The ratio of IR to UV emission (IRX) is used to constrain the dust attenuation in galaxies. We use the 4000A break as a robust and useful, although coarse, indicator of star formation history (SFH). We examine the relationship between IRX and the UV spectral slope (a common attenuation indicator at high-redshift) and find little dependence of the scatter on 4000A break strength. We construct average UV through far-IR spectral energy distributions (SEDs) for different ranges of IRX, 4000A break strength, and stellar mass (M_*) to show the variation of the entire SED with these parameters. When binned simultaneously by IRX, 4000A break strength, and M_* these SEDs allow us to determine a low resolution average attenuation curve for different ranges of M_*. The attenuation curves thus derived are consistent with a lambda^{-0.7} attenuation law, and we find no significant variations with M_*. Finally, we show the relationship between IRX and the global stellar mass surface density and gas-phase-metallicity. Among star forming galaxies we find a strong correlation between IRX and stellar mass surface density, even at constant metallicity, a result that is closely linked to the well-known correlation between IRX and star-formation rate.

We have initiated a search for extended ultraviolet disk (XUV-disk) galaxies in the local universe. Herein, we compare GALEX UV and visible--NIR images of 189 nearby (D$<$40 Mpc) S0--Sm galaxies included in the GALEX Atlas of Nearby Galaxies and present the first catalogue of XUV-disk galaxies. We find that XUV-disk galaxies are surprisingly common but have varied relative (UV/optical) extent and morphology. Type~1 objects ($\ga$20% incidence) have structured, UV-bright/optically-faint emission features in the outer disk, beyond the traditional star formation threshold. Type~2 XUV-disk galaxies ($\sim$10% incidence) exhibit an exceptionally large, UV-bright/optically-low-surface-brightness (LSB) zone having blue $UV-K_s$ outside the effective extent of the inner, older stellar population, but not reaching extreme galactocentric distance. If the activity occuring in XUV-disks is episodic, a higher fraction of present-day spirals could be influenced by such outer disk star formation. Type~1 disks are associated with spirals of all types, whereas Type~2 XUV-disks are predominantly found in late-type spirals. Type~2 XUV-disks are forming stars quickly enough to double their [presently low] stellar mass in the next Gyr (assuming a constant SF rate). XUV-disk galaxies of both types are systematically more gas-rich than the general galaxy population. Minor external perturbation may stimulate XUV-disk incidence, at least for Type~1 objects. XUV-disks are the most actively evolving galaxies growing via inside-out disk formation in the current epoch, and may constitute a segment of the galaxy population experiencing significant, continued gas accretion from the intergalactic medium or neighboring objects.

The color of galaxies is a fundamental property, easily measured, that constrains models of galaxies and their evolution. Dust attenuation and star formation history (SFH) are the dominant factors affecting the color of galaxies. Here we explore the empirical relation between SFH, attenuation, and color for a wide range of galaxies, including early types. These galaxies have been observed by GALEX, SDSS, and Spitzer, allowing the construction of measures of dust attenuation from the ratio of infrared (IR) to ultraviolet (UV) flux and measures of SFH from the strength of the 4000A break. The empirical relation between these three quantities is compared to models that separately predict the effects of dust and SFH on color. This comparison demonstrates the quantitative consistency of these simple models with the data and hints at the power of multiwavelength data for constraining these models. The UV color is a strong constraint; we find that a Milky Way extinction curve is disfavored, and that the UV emission of galaxies with large 4000A break strengths is likely to arise from evolved populations. We perform fits to the relation between SFH, attenuation, and color. This relation links the production of starlight and its absorption by dust to the subsequent reemission of the absorbed light in the IR. Galaxy models that self-consistently treat dust absorption and emission as well as stellar populations will need to reproduce these fitted relations in the low-redshift universe.

We use the UV-optical color magnitude diagram in combination with spectroscopic and photometric measurements derived from the SDSS spectroscopic sample to measure the distribution of galaxies in the local universe (z<0.25) and their physical properties as a function of specific star formation rate (SSFR) and stellar mass. Throughout this study our emphasis is on the properties of galaxies on and off of a local "star-forming sequence." We discuss how the physical characteristics of galaxies along this sequence are related to scaling relations typically derived for galaxies of different morphological types. We find, among other trends that our measure of the star formation rate surface density is nearly constant along this sequence. We discuss this result and implications for galaxies at higher redshift. For the first time, we report on measurements of the local UV luminosity function versus galaxy structural parameters as well as inclination. We also split our sample into disk-dominated and bulge-dominated subsamples using the i-band Sersic index and find that disk-dominated galaxies occupy a very tight locus in SSFR vs. stellar mass space while bulge-dominated galaxies display a much larger spread of SSFR at fixed stellar mass. A significant fraction of galaxies with SSFR and SF surface density above those on the "star-forming sequence" are bulge-dominated. We can use our derived distribution functions to ask whether a significant fraction of these galaxies may be experiencing a final episode of star formation (possibly induced by a merger or other burst), soon to be quenched, by determining whether this population can explain the growth rate of the non-star-forming galaxies on the "red sequence." (Abridged)

Given multiband photometric data from the SDSS DR6, we estimate galaxy redshifts. We employ a Random Forest trained on color features and spectroscopic redshifts from 80,000 randomly chosen primary galaxies yielding a mapping from color to redshift such that the difference between the estimate and the spectroscopic redshift is small. Our methodology results in tight RMS scatter in the estimates limited by photometric errors. Additionally, this approach yields an error distribution that is nearly Gaussian with parameter estimates giving reliable confidence intervals unique to each galaxy photometric redshift.

We have studied ~2100 early-type galaxies in the SDSS DR3 which have been detected by the GALEX Medium Imaging Survey (MIS), in the redshift range 0 < z < 0.11. Combining GALEX UV photometry with corollary optical data from the SDSS, we find that, at a 95 percent confidence level, at least ~30 percent of galaxies in this sample have UV to optical colours consistent with some recent star formation within the last Gyr. In particular, galaxies with a NUV - r colour less than 5.5 are very likely to have experienced such recent star formation, taking into account the possibility of a contribution to NUV flux from the UV upturn phenomenon. We find quantitative agreement between the observations and the predictions of a semi-analytical LCDM hierarchical merger model and deduce that early-type galaxies in the redshift range 0 < z < 0.11 have ~1 to 3 percent of their stellar mass in stars less than 1 Gyr old. The average age of this recently formed population is ~300 to 500 Myrs. We also find that monolithically evolving galaxies, where recent star formation can be driven solely by recycled gas from stellar mass loss, cannot exhibit the blue colours (NUV - r < 5.5) seen in a significant fraction (~30 percent) of our observed sample.

We report on the properties of a sample of ultraviolet luminous galaxies (UVLGs) selected by matching the Galaxy Evolution Explorer (GALEX) Surveys with the Sloan Digital Sky Survey Third Data Release. Out of 25362 galaxies between 0.0<z<0.3 detected by GALEX, there are 215 galaxies with L>2x10^10 L_solar at 1530 Angstroms (observed wavelength). The properties of this population are well correlated with ultraviolet surface brightness. We find that the galaxies with low UV surface brightness are primarily large spiral systems with a mixture of old and young stellar populations, while the high surface brightness galaxies consist primarily of compact starburst systems. In terms of the behavior of surface brightness with luminosity, size with luminosity, the mass-metallicity relation, and other parameters, the compact UVLGs clearly depart from the trends established by the full sample of galaxies. The subset of compact UVLGs with the highest surface brightness (``supercompact UVLGs'') have characteristics that are remarkably similar to Lyman Break Galaxies at higher redshift. They are much more luminous than typical local ultraviolet-bright starburst galaxies and blue compact dwarf galaxies. They have metallicities that are systematically lower than normal galaxies of the same stellar mass, indicating that they are less chemically evolved. In all these respects, they are the best local analogs for Lyman Break Galaxies.

We present the first measurements of the angular correlation function of galaxies selected in the far (1530 A) and near (2310 A) Ultraviolet from the GALEX survey fields overlapping SDSS DR5 in low galactic extinction regions. The area used covers 120 sqdeg (GALEX - MIS) down to magnitude AB = 22, yielding a total of 100,000 galaxies. The mean correlation length is ~ 3.7 \pm 0.6 Mpc and no significant trend is seen for this value as a function of the limiting apparent magnitude or between the GALEX bands. This estimate is close to that found from samples of blue galaxies in the local universe selected in the visible, and similar to that derived at z ~ 3 for LBGs with similar rest frame selection criteria. This result supports models that predict anti-biasing of star forming galaxies at low redshift, and brings an additional clue to the downsizing of star formation at z<1.

We introduce a method to constrain general cosmological models using Baryon Acoustic Oscillation (BAO) distance measurements from galaxy samples covering different redshift ranges, and apply this method to analyse samples drawn from the SDSS and 2dFGRS. BAO are detected in the clustering of the combined 2dFGRS and SDSS main galaxy samples, and measure the distance--redshift relation at z=0.2. BAO in the clustering of the SDSS luminous red galaxies measure the distance--redshift relation at z=0.35. The observed scale of the BAO calculated from these samples and from the combined sample are jointly analysed using estimates of the correlated errors, to constrain the form of the distance measure D_V(z)=[(1+z)^2D_A^2cz/H(z)]^(1/3). Here D_A is the angular diameter distance, and H(z) is the Hubble parameter. This gives r_s/D_V(0.2)=0.1980+/-0.0058 and r_s/D_V(0.35)=0.1094+/-0.0033 (1sigma errors), with correlation coefficient of 0.39, where r_s is the comoving sound horizon scale at recombination. Matching the BAO to have the same measured scale at all redshifts then gives D_V(0.35)/D_V(0.2)=1.812+/-0.060. The recovered ratio is roughly consistent with that predicted by the higher redshift SNLS supernovae data for Lambda cosmologies, but does require slightly stronger cosmological acceleration at low redshift. If we force the cosmological model to be flat with constant w, then we find Om_m=0.249+/-0.018 and w=-1.004+/-0.089 after combining with the SNLS data, and including the WMAP measurement of the apparent acoustic horizon angle in the CMB.

We use a new stacking technique to obtain mean mid IR and far IR to far UV flux ratios over the rest near-UV/near-IR color-magnitude diagram. We employ COMBO-17 redshifts and COMBO-17 optical, GALEX far and near UV, Spitzer IRAC and MIPS Mid IR photometry. This technique permits us to probe infrared excess (IRX), the ratio of far IR to far UV luminosity, and specific star formation rate (SSFR) and their co-evolution over two orders of magnitude of stellar mass and redshift 0.1<z<1.2. We find that the SSFR and the characteristic mass (M_0) above which the SSFR drops increase with redshift (downsizing). At any given epoch, IRX is an increasing function of mass up to M_0. Above this mass IRX falls, suggesting gas exhaustion. In a given mass bin below M_0 IRX increases with time in a fashion consistent with enrichment. We interpret these trends using a simple model with a Schmidt-Kennicutt law and extinction that tracks gas density and enrichment. We find that the average IRX and SSFR follows a galaxy age parameter which is determined mainly by the galaxy mass and time since formation. We conclude that blue sequence galaxies have properties which show simple, systematic trends with mass and time such as the steady build-up of heavy elements in the interstellar media of evolving galaxies and the exhaustion of gas in galaxies that are evolving off the blue sequence. The IRX represents a tool for selecting galaxies at various stages of evolution.

We study the ultraviolet to far-infrared (hereafter UV-to-IR) SEDs of a sample of intermediate redshift (0.2 < z < 0.7) UV-selected galaxies from the ELAIS-N1 and ELAIS-N2 fields by fitting a multi-wavelength dataset to a library of GRASIL templates. Star formation related properties of the galaxies are derived from the library of models by using the Bayesian statistics. We find a decreasing presence of galaxies with low attenuation and low total luminosity as redshift decreases, which does not hold for high total luminosity galaxies. In addition the dust attenuation of low mass galaxies increases as redshift decreases, and this trend seems to disappear for galaxies with M* > 10^11 M_sun. This result is consistent with a mass dependent evolution of the dust to gas ratio, which could be driven by a mass dependent efficiency of star formation in star forming galaxies. The specific star formation rates (SSFR) decrease with increasing stellar mass at all redshifts, and for a given stellar mass the SSFR decreases with decreasing redshift. The differences in the slope of the M*--SSFR relation found between this work and others at similar redshift could be explained by the adopted selection criteria of the samples which, for a UV selected sample, favours blue, star forming galaxies.

Using a sample of galaxies from the Sloan Digital Sky Survey spectroscopic catalog with measured star-formation rates (SFRs) and ultraviolet (UV) photometry from the GALEX Medium Imaging Survey, we derived empirical linear correlations between the SFR to UV luminosity ratio and the UV-optical colors of blue sequence galaxies. The relations provide a simple prescription to correct UV data for dust attenuation that best reconciles the SFRs derived from UV and emission line data. The method breaks down for the red sequence population as well as for very blue galaxies such as the local ``supercompact'' UV luminous galaxies and the majority of high redshift Lyman Break Galaxies which form a low attenuation sequence of their own.

Ultraviolet Luminous Galaxies (UVLGs) have been identified as intensely star-forming, nearby galaxies. A subset of these, the supercompact UVLGs, are believed to be local analogs of high redshift Lyman Break Galaxies. Here we investigate the radio continuum properties of this important population for the first time. We have observed 42 supercompact UVLGs with the VLA, all of which have extensive coverage in the UV/optical by GALEX and SDSS. Our analysis includes comparison samples of multiwavelength data from the Spitzer First Look Survey and from the SDSS-Galex matched catalogs. In addition we have Spitzer MIPS data for 24 of our galaxies and find that they fall on the radio-FIR correlation of normal star-forming galaxies. We find that our galaxies have lower radio-to-UV ratios and lower Balmer decrements than other local galaxies with similar (high) star formation rates. Optical spectra show they have lower Dn(4000) and HdeltaA indices, higher Hbeta emission-line equivalents widths, and higher [OIII]5007/Hbeta emission-line ratios than normal star forming galaxies. Comparing these results to galaxy spectral evolution models we conclude that supercompact UVLGs are distinguished from normal star forming galaxies firstly by their high specific star formation rates. Moreover, compared to other types of galaxies with similar star formation rates, they have significantly less dust attenuation. In both regards they are similar to Lyman Break Galaxies. This suggests that the process that causes star formation in the supercompact UVLGs differs from other local star forming galaxies, but may be similar to Lyman Break Galaxies.

We investigate the production of nitrogen in star forming galaxies with ultraviolet (UV) radiation detected by the Galaxy Evolution Explorer Satellite (GALEX). We use a sample of 8,745 GALEX emission line galaxies matched to the Sloan Digital Sky Survey (SDSS) spectroscopic sample. We derive both gas-phase oxygen and nitrogen abundances for the sample, and apply stellar population synthesis models to derive stellar masses and star formation histories of the galaxies. We compare oxygen abundances derived using three different diagnostics. We derive the specific star formation rates of the galaxies by modeling the 7-band GALEX+SDSS photometry. We find that galaxies that have log SFR/M$_*$ > -10.0 typically have values of log N/O ~0.05 dex less than galaxies with log SFR/M$_*$ < -10.0 and similar oxygen abundances.

We have analyzed the bivariate distribution of galaxies as a function of ultraviolet-optical colors and absolute magnitudes in the local universe. The sample consists of galaxies with redshifts and optical photometry from the Sloan Digital Sky Survey (SDSS) main galaxy sample matched with detections in the near-ultraviolet (NUV) and far-ultraviolet (FUV) bands in the Medium Imaging Survey being carried out by the Galaxy Evolution Explorer (GALEX) satellite. In the (NUV-r)_{0.1} vs. M_{r,0.1} galaxy color-magnitude diagram, the galaxies separate into two well-defined blue and red sequences. The (NUV-r)_{0.1} color distribution at each M_{r,0.1} is not well fit by the sum of two Gaussians due to an excess of galaxies in between the two sequences. The peaks of both sequences become redder with increasing luminosity with a distinct blue peak visible up to M_{r,0.1}\sim-23. The r_{0.1}-band luminosity functions vary systematically with color, with the faint end slope and characteristic luminosity gradually increasing with color. After correcting for attenuation due to dust, we find that approximately one quarter of the color variation along the blue sequence is due to dust with the remainder due to star formation history and metallicity. Finally, we present the distribution of galaxies as a function of specific star formation rate and stellar mass. The specific star formation rates imply that galaxies along the blue sequence progress from low mass galaxies with star formation rates that increase somewhat with time to more massive galaxies with a more or less constant star formation rate. Above a stellar mass of ~10^10.5 M_{sun}, galaxies with low ratios of current to past averaged star formation rate begin to dominate.

We report results from a pilot program to obtain spectroscopy for objects detected in the Galaxy Evolution Explorer (GALEX) Medium Imaging Survey (MIS). Our study examines the properties of galaxies detected by GALEX fainter than the Sloan Digital Sky Survey (SDSS) spectroscopic survey. This is the first study to extend the techinques of Salim et al. 2005 to estimate stellar masses, star formation rates (SFR) and the b (star formation history) parameter for star-forming galaxies out to z~0.7. We obtain redshifts for 50 GALEX MIS sources reaching NUV=23.9 (AB mag), having counterparts in the SDSS Data Release 4 (DR4). Of our sample, 43 are starforming galaxies with z<0.7, 3 have emission line ratios indicative of AGN with z<0.7, and 4 objects with z~1 are QSOs, 3 of which are not previously cataloged. We compare our sample to a much larger sample of ~50,000 matched GALEX/SDSS galaxies with SDSS spectroscopy; while our survey is shallow, the optical counterparts to our sources reach ~3 magnitudes fainter in SDSS r magnitude than the SDSS spectroscopic sample. We use emission line diagnostics for the galaxies to determine that the sample contains mostly star-forming galaxies. The galaxies in the sample populate the blue sequence in the NUV-r vs M_r color-magnitude diagram. Our sample has SFRs, luminosities, and velocity dispersions that are similar to the samples of faint compact blue galaxies studied previously in the same redshift range by Koo et. al 1995, Guzman et. al 1996 & Phillips et. al 1997. However, our sample is ~2 mag fainter in surface brightness than the compact blue galaxies. We find that the star-formation histories for a majority of the galaxies are consistent with a recent starburst within the last 100 Myr.

We analyze the clustering properties of ultraviolet selected galaxies by using GALEX-SDSS data at z<0.6 and CFHTLS deep u' imaging at z=1. These datasets provide a unique basis at z< 1 which can be directly compared with high redshift samples built with similar selection criteria. We discuss the dependence of the correlation function parameters (r0, delta) on the ultraviolet luminosity as well as the linear bias evolution. We find that the bias parameter shows a gradual decline from high (b > 2) to low redshift (b ~ 0.79^{+0.1}_{-0.08}). When accounting for the fraction of the star formation activity enclosed in the different samples, our results suggest that the bulk of star formation migrated from high mass dark matter halos at z>2 (10^12 < M_min < 10^13 M_sun, located in high density regions), to less massive halos at low redshift (M_min < 10^12 M_sun, located in low density regions). This result extends the ``downsizing'' picture (shift of the star formation activity from high stellar mass systems at high z to low stellar mass at low z) to the dark matter distribution.

We describe the calibration status and data products pertaining to the GR2 and GR3 data releases of the Galaxy Evolution Explorer (GALEX). These releases have identical pipeline calibrations that are significantly improved over the GR1 data release. GALEX continues to survey the sky in the Far Ultraviolet (FUV, ~154 nm) and Near Ultraviolet (NUV, ~232 nm) bands, providing simultaneous imaging with a pair of photon counting, microchannel plate, delay line readout detectors. These 1.25 degree field-of-view detectors are well-suited to ultraviolet observations because of their excellent red rejection and negligible background. A dithered mode of observing and photon list output pose complex requirements on the data processing pipeline, entangling detector calibrations and aspect reconstruction algorithms. Recent improvements have achieved photometric repeatability of 0.05 and 0.03 mAB in the FUV and NUV, respectively. We have detected a long term drift of order 1% FUV and 6% NUV over the mission. Astrometric precision is of order 0.5" RMS in both bands. In this paper we provide the GALEX user with a broad overview of the calibration issues likely to be confronted in the current release. Improvements are likely as the GALEX mission continues into an extended phase with a healthy instrument, no consumables, and increased opportunities for guest investigations.

We measure star formation rates of ~50,000 optically-selected galaxies in the local universe (z~0.1), spanning a range from gas-rich dwarfs to massive ellipticals. We obtain dust-corrected SFRs by fitting the GALEX (UV) and SDSS (optical) photometry to a library of population synthesis models that include dust attenuation. For star-forming galaxies, our UV-based SFRs compare remarkably well with those derived from SDSS H alpha. Deviations from perfect agreement between these two methods are due to differences in the dust attenuation estimates. In contrast to H alpha, UV provides reliable SFRs for galaxies with weak or no H alpha emission, and where H alpha is contaminated with an emission from an AGN. We use full-SED SFRs to calibrate a simple prescription that uses GALEX UV magnitudes to produce good SFRs for normal star-forming galaxies. The specific SFR is considered as a function of stellar mass for (1) star-forming galaxies with no AGN, (2) those hosting an AGN, and for (3) galaxies without H alpha emission. We find that the three have distinct star formation histories, with AGN lying intermediate between the star-forming and the quiescent galaxies. Normal star forming galaxies (without an AGN) lie on a relatively narrow linear sequence. Remarkably, galaxies hosting a strong AGN appear to represent the massive continuation of this sequence. Weak AGN, while also massive, have lower SFR, sometimes extending to the realm of quiescent galaxies. We propose an evolutionary sequence for massive galaxies that smoothly connects normal star-forming galaxies to quiescent (red sequence) galaxies via strong and weak AGN. We confirm that some galaxies with no H alpha emission show signs of SF in the UV. We derive a UV-based cosmic SFR density at z=0.1 with smaller total error than previous measurements (abridged).

We present the fourth edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog. The catalog contains 77,429 objects; this is an increase of over 30,000 entries since the previous edition. The catalog consists of the objects in the SDSS Fifth Data Release that have luminosities larger than M_i = -22.0 (in a cosmology with H_0 = 70 km/s/Mpc, Omega_M = 0.3, and Omega_Lambda = 0.7) have at least one emission line with FWHM larger than 1000 km/s, or have interesting/complex absorption features, are fainter than i=15.0, and have highly reliable redshifts. The area covered by the catalog is 5740 sq. deg. The quasar redshifts range from 0.08 to 5.41, with a median value of 1.48; the catalog includes 891 quasars at redshifts greater than four, of which 36 are at redshifts greater than five. Approximately half of the catalog quasars have i < 19; nearly all have i < 21. For each object the catalog presents positions accurate to better than 0.2 arcsec. rms per coordinate, five-band (ugriz) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains basic radio, near-infrared, and X-ray emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3800--9200A at a spectral resolution of ~2000. The spectra can be retrieved from the public database using the information provided in the catalog. The average SDSS colors of quasars as a function of redshift, derived from the catalog entries, are presented in tabular form. Approximately 96% of the objects in the catalog were discovered by the SDSS.

We present the GALEX UV photometry of the elliptical galaxies in Abell clusters at moderate redshifts (z < 0.2) for the study of the look-back time evolution of the UV upturn phenomenon. The brightest elliptical galaxies (M_r < -22) in 12 remote clusters are compared with the nearby giant elliptical galaxies of comparable optical luminosity in the Fornax and Virgo clusters. The sample galaxies presented here appear to be quiescent without signs of massive star formation or strong nuclear activity, and show smooth, extended profiles in their UV images indicating that the far-UV (FUV) light is mostly produced by hot stars in the underlying old stellar population. Compared to their counterparts in nearby clusters, the FUV flux of cluster giant elliptical galaxies at moderate redshifts fades rapidly with ~ 2 Gyrs of look-back time, and the observed pace in FUV - V color evolution agrees reasonably well with the prediction from the population synthesis models where the dominant FUV source is hot horizontal-branch stars and their progeny. A similar amount of color spread (~ 1 mag) in FUV - V exists among the brightest cluster elliptical galaxies at z ~ 0.1, as observed among the nearby giant elliptical galaxies of comparable optical luminosity.

We introduce a new quantity, the mass flux density of galaxies evolving from the blue sequence to the red sequence. We propose a simple technique for constraining this mass flux using the volume corrected number density in the extinction-corrected UV-optical color magnitude distribution, the stellar age indexes H-delta-a and D_n(4000), and a simple prescription for spectral evolution using a quenched star formation history. We exploit the excellent separation of red and blue sequences in the NUV-r band color-magnitude diagram. The final value we measure, 0.033 M_sun yr^-1 Mpc^-3, is strictly speaking an upper limit due to the possible contributions of bursting, composite, and extincted galaxies. However, it compares favorably with estimates of the average mass flux that we make based on the red luminosity function evolution derived from the DEEPII and COMBO-17 surveys (Bell et al 2004; Faber et al. 2005), 0.034 M_sun yr^-1 Mpc^-3. We find that the blue sequence mass has remained roughly constant since z=1 (mass flux 0.01 M_sun yr^-1 Mpc^-3) but the average on-going star formation of 0.037$ M_sun yr^-1 Mpc^-3 over 0<z<1 is balanced by mass flux off the blue sequence. We explore the nature of the galaxies in the transition zone with particular attention to the frequency and impact of AGNs. The AGN fraction peaks in the transition zone. We find circumstantial, albeit weak evidence that the quench rates are higher in higher luminosity AGNs.

We present the results from the analysis of optical spectra of 31 Halpha-selected regions in the extended UV (XUV) disks of M83 (NGC5236) and NGC4625 recently discovered by GALEX. The spectra were obtained using IMACS at Las Campanas Observatory 6.5m Magellan I telescope and COSMIC at the Palomar 200-inch telescope, respectively for M83 and NGC4625. The line ratios measured indicate nebular oxygen abundances (derived from the R23 parameter) of the order of Zsun/5-Zsun/10. For most emission-line regions analyzed the line fluxes and ratios measured are best reproduced by models of photoionization by single stars with masses in the range 20-40 Msun and oxygen abundances comparable to those derived from the R23 parameter. We find indications for a relatively high N/O abundance ratio in the XUV disk of M83. Although the metallicities derived imply that these are not the first stars formed in the XUV disks, such a level of enrichment could be reached in young spiral disks only 1 Gyr after these first stars would have formed. The amount of gas in the XUV disks allow maintaining the current level of star formation for at least a few Gyr.

(Abridged) We study the two-point correlation function of a uniformly selected sample of 4,426 luminous optical quasars with redshift $2.9 \le z\le 5.4$ selected over 4041 deg$^2$ from the Fifth Data Release of the Sloan Digital Sky Survey. For a real-space correlation function of the form $\xi(r)=(r/r_0)^{-\gamma}$, the fitted parameters in comoving coordinates are $r_0 = 15.2 \pm 2.7 h^{-1}$ Mpc and $\gamma = 2.0 \pm 0.3$, over a scale range $4\le r_p\le 150 h^{-1}$ Mpc. Thus high-redshift quasars are appreciably more strongly clustered than their $z \approx 1.5$ counterparts, which have a comoving clustering length $r_0 \approx 6.5 h^{-1}$ Mpc. Dividing our sample into two redshift bins: $2.9\le z\le 3.5$ and $z\ge 3.5$, and assuming a power-law index $\gamma=2.0$, we find a correlation length of $r_0 = 16.9 \pm 1.7 h^{-1}$ Mpc for the former, and $r_0 = 24.3 \pm 2.4 h^{-1}$ Mpc for the latter. Following Martini & Weinberg, we relate the clustering strength and quasar number density to the quasar lifetimes and duty cycle. Using the Sheth & Tormen halo mass function, the quasar lifetime is estimated to lie in the range $4\sim 50$ Myr for quasars with $2.9\le z\le 3.5$; and $30\sim 600$ Myr for quasars with $z\ge 3.5$. The corresponding duty cycles are $0.004\sim 0.05$ for the lower redshift bin and $0.03\sim 0.6$ for the higher redshift bin. The minimum mass of halos in which these quasars reside is $2-3\times 10^{12}\ h^{-1}M_\odot$ for quasars with $2.9\le z\le 3.5$ and $4-6\times 10^{12}\ h^{-1}M_\odot$ for quasars with $z\ge 3.5$.

Using subspace methods, we study the distribution of physical components of galaxies in wavelength space. We find that it is valid to assume that the stellar and the gaseous components of galaxies span complementary subspaces. To first order, stellar and gaseous spectral features can be extracted from galaxy spectra through a simple matrix multiplication. By comparing the stellar continua obtained respectively using the model-based and the empirical approach through a commonality measure, we conclude that the latter may lose higher-order spectral features.

We present a study of the morphological nature of redshift z~0.7 star-forming galaxies using a combination of HST/ACS, GALEX and ground-based images of the COSMOS field. Our sample consists of 8,146 galaxies, 5,777 of which are detected in the GALEX near-ultraviolet band down to a limiting magnitude of 25.5 (AB). We make use of the UV to estimate star formation rates, correcting for the effect of dust using the UV-slope, and compute, from the ACS F814W images, the C,A,S,G,M20 morphological parameters for all objects in our sample. We observe a morphological bimodality in the galaxy population and show that it has a strong correspondence with the FUV - g color bimodality. We conclude that UV-optical color predominantly evolves concurrently with morphology. We observe many of the most star-forming galaxies to have morphologies approaching that of early-type galaxies, and interpret this as evidence that strong starburst events are linked to bulge growth and constitute a process through which galaxies can be brought from the blue to the red sequence while simultaneously modifying their morphology accordingly. We conclude that the red sequence has continued growing at z~<0.7. We also observe z~0.7 galaxies to have physical properties similar to that of local galaxies, except for higher star formation rates. Whence we infer that the dimming of star-forming galaxies is responsible for most of the evolution in the star formation rate density of the Universe since that redshift, although our data are also consistent with a mild number evolution. [abridged]

We present ultraviolet (UV) photometry of M31 globular clusters (GCs) found in 23 Galaxy Evolution Explorer (GALEX) images covering the entirety of M31. We detect 485 and 273 GCs (and GC candidates) in the near-ultraviolet (NUV; 2267 A) and far-ultraviolet (FUV; 1516 A), respectively. Comparing M31 data with those of Galactic GCs in the UV with the aid of population models, we find that the age ranges of old GCs in M31 and the Galactic halo are similar. Three metal-rich ([Fe/H]>-1) GCs in M31 produce significant FUV flux making their FUV-V colors unusually blue for their metallicities. These are thought to be analogs of the two peculiar Galactic GCs NGC 6388 and NGC 6441 with extended blue HB stars. Based on the models incorporating helium enriched subpopulations in addition to the majority of the population that have a normal helium abundance, we suggest that even small fraction of super-helium-rich subpopulations in GCs can reproduce the observed UV bright metal-rich GCs. Young clusters in M31 show distinct UV and optical properties from GCs in Milky Way. Population models indicate that their typical age is less than ~ 2 Gyrs. A large fraction of young GCs have the kinematics of the thin, rapidly rotating disk component. However, a subset of the old GCs also shares the thin-disk kinematics of the younger clusters. The existence of young GCs on the outskirts of M31 disk suggests the occurrence of a significant recent star formation in the thin-disk of M31. Old thin-disk GCs may set constraints on the epoch of early formation of the M31 thin-disk. We detect 12 (10) intermediate-age GC candidates in NUV (FUV). We suggest that some of spectroscopically identified intermediate-age GCs may not be truly intermediate in age, but rather older GCs that possess developed HB.

We use the Galaxy Evolution Explorer (GALEX) Medium and All-Sky-Imaging Survey (MIS & AIS) data from the first public data release (GR1), matched to the Sloan Digital Sky Survey (SDSS) DR3 catalog, to perform source classification. The GALEX surveys provide photometry in far- and near-UV bands and the SDSS in five optical bands (u,g,r,i,z). The GR1/DR3 overlapping areas are 363[83]deg^2 for the GALEX AIS[MIS], for sources within the 0.5deg central area of the GALEX fields. Our sample covers mostly |b|>30deg galactic latitudes. We present statistical properties of the GALEX/SDSS matched sources catalog, containing >2x10^6 objects detected in at least one UV band. We classify the matched sources by comparing the seven-band photometry to model colors constructed for different classes of astrophysical objects. For sources with photometric errors <0.3 mag, the corresponding typical AB-magnitude limits are m_FUV~21.5, m_NUV~22.5 for AIS, and m_FUV~24, m_NUV~24.5 for MIS. At AIS depth, the number of Galactic and extragalactic objects are comparable, but the latter predominate in the MIS. Based on our stellar models, we estimate the GALEX surveys detect hot White Dwarfs throughout the Milky Way halo (down to a radius of 0.04 R_sun at MIS depth), providing an unprecedented improvement in the Galactic WD census. Their observed surface density is consistent with Milky Way model predictions. We also select low-redshift QSO candidates, extending the known QSO samples to lower magnitudes, and providing candidates for detailed z~1 follow-up investigations. SDSS optical spectra available for a large subsample confirm the classification for the photometrically selected candidates with 97% purity for single hot stars, ~45%(AIS)/31%(MIS) for binaries containing a hot star and a cooler companion, and about 85% for QSOs.

We measure the cosmological matter density by observing the positions of baryon acoustic oscillations in the clustering of galaxies in the Sloan Digital Sky Survey (SDSS). We jointly analyse the main galaxies and LRGs in the SDSS DR5 sample, using over half a million galaxies in total. The oscillations are detected with 99.74% confidence (3.0sigma assuming Gaussianity) compared to a smooth power spectrum. When combined with the observed scale of the peaks within the CMB, we find a best-fit value of Omega_m=0.256+0.029-0.024 (68% confidence interval), for a flat Lambda cosmology when marginalising over the Hubble parameter and the baryon density. This value of the matter density is derived from the locations of the baryon oscillations in the galaxy power spectrum and in the CMB, and does not include any information from the overall shape of the power spectra. This is an extremely clean cosmological measurement as the physics of the baryon acoustic oscillation production is well understood, and the positions of the oscillations are expected to be independent of systematics such as galaxy bias.

We measure the large-scale real-space power spectrum P(k) using luminous red galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS) and use this measurement to sharpen constraints on cosmological parameters from the Wilkinson Microwave Anisotropy Probe (WMAP). We employ a matrix-based power spectrum estimation method using Pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 20 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.01h/Mpc < k < 0.2h/Mpc. Results from the LRG and main galaxy samples are consistent, with the former providing higher signal-to-noise. Our results are robust to omitting angular and radial density fluctuations and are consistent between different parts of the sky. They provide a striking confirmation of the predicted large-scale LCDM power spectrum. Combining only SDSS LRG and WMAP data places robust constraints on many cosmological parameters that complement prior analyses of multiple data sets. The LRGs provide independent cross-checks on Om and the baryon fraction in good agreement with WMAP. Within the context of flat LCDM models, our LRG measurements complement WMAP by sharpening the constraints on the matter density, the neutrino density and the tensor amplitude by about a factor of two, giving Omega_m=0.24+-0.02 (1 sigma), sum m_nu < 0.9 eV (95%) and r<0.3 (95%). Baryon oscillations are clearly detected and provide a robust measurement of the comoving distance to the median survey redshift z=0.35 independent of curvature and dark energy properties. Within the LCDM framework, our power spectrum measurement improves the evidence for spatial flatness, sharpening the curvature constraint Omega_tot=1.05+-0.05 from WMAP alone to Omega_tot=1.003+-0.010. Assuming Omega_tot=1, the equation of state parameter is constrained to w=-0.94+-0.09, indicating the potential for more ambitious future LRG measurements to provide precision tests of the nature of dark energy. All these constraints are essentially independent of scales k>0.1h/Mpc and associated nonlinear complications, yet agree well with more aggressive published analyses where nonlinear modeling is crucial.

We present a Fourier analysis of the clustering of galaxies in the combined Main galaxy and Luminous Red Galaxy (LRG) Sloan Digital Sky Survey (SDSS) Data Release 5 (DR5) sample. The aim of our analysis is to consider how well we can measure the cosmological matter density using the signature of the horizon at matter-radiation equality embedded in the large-scale power spectrum. The new data constrains the power spectrum on scales 100--600h^-1Mpc with significantly higher precision than previous analyses of just the SDSS Main galaxies, due to our larger sample and the inclusion of the LRGs. This improvement means that we can now reveal a discrepancy between the shape of the measured power and linear CDM models on scales 0.01<k<0.15hMpc^-1, with linear model fits favouring a lower matter density (Omega_m=0.22+/-0.04) on scales 0.01<k<0.06hMpc^-1 and a higher matter density (Omega_m=0.32+/-0.01) when smaller scales are included, assuming a flat LCDM model with h=0.73 and n_s=0.96. This discrepancy could be explained by scale-dependent bias and, by analysing subsamples of galaxies, we find that the ratio of small-scale to large-scale power increases with galaxy luminosity, so all of the SDSS galaxies cannot trace the same power spectrum shape over 0.01<k<0.2hMpc^-1. However, the data are insufficient to clearly show a luminosity-dependent change in the largest scale at which a significant increase in clustering is observed, although they do not rule out such an effect. Significant scale-dependent galaxy bias on large-scales, which changes with the r-band luminosity of the galaxies, could potentially explain differences in our Omega_m estimates and differences previously observed between 2dFGRS and SDSS power spectra and the resulting parameter constraints.

We select far-infrared (FIR-60 microns) and far-ultraviolet (FUV-1530 A) samples of nearby galaxies in order to discuss the biases encountered by monochromatic surveys (FIR or FUV). Very different volumes are sampled by each selection and much care is taken to apply volume corrections to all the analyses. The distributions of the bolometric luminosity of young stars are compared for both samples: they are found to be consistent with each other for galaxies of intermediate luminosities but some differences are found for high (>5 10^{10} L_sun) luminosities. The shallowness of the IRAS survey prevents us from securing comparison at low luminosities (<2 10^9 L_sun). The ratio of the total infrared (TIR) luminosity to the FUV luminosity is found to increase with the bolometric luminosity in a similar way for both samples up to 5 10^{10} L_sun. Brighter galaxies are found to have a different behavior according to their selection: the L_TIR/L_FUV ratio of the FUV-selected galaxies brighter than 5 10^{10} L_sun reaches a plateau whereas L_TIR/L_FUV continues to increase with the luminosity of bright galaxies selected in FIR. The volume-averaged specific star formation rate (SFR per unit galaxy stellar mass, SSFR) is found to decrease toward massive galaxies within each selection. The SSFR is found to be larger than that measured for optical and NIR-selected sample over the whole mass range for the FIR selection, and for masses larger than 10^{10} M_sun for the FUV selection. Luminous and massive galaxies selected in FIR appear as active as galaxies with similar characteristics detected at z ~ 0.7.

We analyze a volume-limited sample of massive bulge-dominated galaxies with data from both the Sloan Digital Sky Survey and the Galaxy Evolution Explorer (GALEX) satellite. The galaxies have central velocity dispersions greater than 100 km/s and stellar surface mass densities that lie above the value where galaxies transition from actively star forming to passive systems. The sample is limited to redshifts 0.03<z<0.07. At these distances, the SDSS spectra sample the light from the bulge-dominated central regions of the galaxies. The GALEX NUV data provide high sensitivity to low rates of global star formation in these systems. Our sample of bulge-dominated galaxies exhibits a much larger dispersion in NUV-r colour than in optical g-r colour. Nearly all of the galaxies with bluer NUV-r colours are AGN. Both GALEX images and SDSS colour profiles demonstrate that the excess UV light is associated with an extended disk. We find that galaxies with red outer regions almost never have a young bulge or a strong AGN. Galaxies with blue outer regions have bulges and black holes that span a wide range in age and accretion rate. Galaxies with young bulges and strongly accreting black holes almost always have blue outer disks. Our suggested scenario is one in which the source of gas that builds the bulge and black hole is a low mass reservoir of cold gas in the disk.The presence of this gas is a necessary, but not sufficient condition for bulge and black hole growth. Some mechanism must transport this gas inwards in a time variable way. As the gas in the disk is converted into stars, the galaxies will turn red, but further gas infall can bring them back into the blue NUV-r sequence.(Abridged)

For a sample of 43 nearby, late-type galaxies, we have investigated the radial variation of both the current star formation rate and the dust-induced UV light attenuation. To do this we have cross-correlated IRAS images and GALEX observations for each of these galaxies, and compiled observations of the gas (CO and HI) and metal-abundance gradients found in the literature. We find that attenuation correlates with metallicity. We then use the UV profiles, corrected for attenuation, to study several variants of the Schmidt law and conclude that our results are compatible with a simple law similar to the one of Kennicutt but extending smoothly to lower surface densities, but with considerable scatter. We do not detect an abrupt break in the UV light at the threshold radius derived from H-alpha data (at which the H-alpha profile shows a break and beyond which only a few HII regions are usually found). We interpret the H-alpha sudden break not as a change in the star formation regime (as often suggested) but as the vanishingly small number of ionizing stars corresponding to low levels of star formation.

We use GALEX/optical photometry to construct color-color relationships for early-type galaxies sorted by morphological type. We have matched objects in the GALEX GR1 public release and the first IR1.1 internal release, with the RC3 early-type galaxies having a morphological type -5.5<T<-1.5 with mean error in T<1.5, and mean error on (B-V)T<0.05. After visual inspection of each match, we are left with 130 galaxies with a reliable GALEX pipeline photometry in the far-UV and near-UV bands. This sample is divided into Ellipticals (-5.5<T<-3.5) and Lenticulars (-3.5<T<-1.5). After correction for the Galactic extinction, the color-color diagrams FUV-NUV vs. (B-V)_{Tc} are plotted for the two subsamples. We find a tight anti-correlation between the FUV-NUV and (B-V)_{Tc} colors for Ellipticals, the UV color getting bluer when the (B-V)_{Tc} get redder. This relationship very likely is an extension of the color-metallicity relationship into the GALEX NUV band. We suspect that the main source of the correlation is metal line blanketing in the NUV band. The FUV-NUV vs B-V correlation has larger scatter for lenticular galaxies; we speculate this reflects the presence of low level star formation. If the latter objects (i.e. those that are blue both in FUV-NUV and B-V) are interpreted as harboring recent star formation activity, this would be the case for a few percent (~4%) of Ellipticals and ~15% of Lenticulars; this would make about 10% of early-type galaxies with residual star formation in our full sample of 130 early-type galaxies. We also plot FUV-NUV vs. the Mg_2 index and central velocity dispersion. We find a tight anti-correlation between FUV-NUV and the Mg_2 index(...).

Detailed high-resolution observations of the innermost regions of nearby galaxies have revealed the presence of supermassive black holes1. These black holes may interact with their host galaxies by means of 'feedback' in the form of energy and material jets; this feedback affects the evolution of the host and gives rise to observed relations between the black hole and the host. Here we report observations of the ultraviolet emissions of massive early-type galaxies. We derive an empirical relation for a critical black-hole mass (as a function of velocity dispersion) above which the outflows from these black holes suppress star formation in their hosts by heating and expelling all available cold gas. Supermassive black holes are negligible in mass compared to their hosts but nevertheless seem to play a critical role in the star formation history of galaxies.

We present an analysis of the mid-infrared (MIR) and optical properties of type 1 (broad-line) quasars detected by the Spitzer Space Telescope. The MIR color-redshift relation is characterized to z=3, with predictions to z=7. We demonstrate how combining MIR and optical colors can yield even more efficient selection of active galactic nuclei (AGN) than MIR or optical colors alone. Composite spectral energy distributions (SEDs) are constructed for 259 quasars with both Sloan Digital Sky Survey and Spitzer photometry, supplemented by near-IR, GALEX, VLA and ROSAT data where available. We discuss how the spectral diversity of quasars influences the determination of bolometric luminosities and accretion rates; assuming the mean SED can lead to errors as large as a factor of 2 for individual quasars. Finally, we show that careful consideration of the shape of the mean quasar SED and its redshift dependence leads to a lower estimate of the fraction of reddened/obscured AGNs missed by optical surveys as compared to estimates derived from a single mean MIR to optical flux ratio.

We present images, integrated photometry, surface-brightness and color profiles for a total of 1034 nearby galaxies recently observed by the GALEX satellite in its far-ultraviolet (FUV; 1516A) and near-ultraviolet (NUV; 2267A) bands. (...) This data set has been complemented with archival optical, near-infrared, and far-infrared fluxes and colors. We find that the integrated (FUV-K) color provides robust discrimination between elliptical and spiral/irregular galaxies and also among spiral galaxies of different sub-types. Elliptical galaxies with brighter K-band luminosities (i.e. more massive) are redder in (NUV-K) color but bluer in (FUV-NUV) than less massive ellipticals. In the case of the spiral/irregular galaxies our analysis shows the presence of a relatively tight correlation between the (FUV-NUV) color and the total infrared-to-UV ratio. The correlation found between (FUV-NUV) color and K-band luminosity (with lower luminosity objects being bluer than more luminous ones) can be explained as due to an increase in the dust content with galaxy luminosity. The images in this Atlas along with the profiles and integrated properties are publicly available through a dedicated web page at http://nedwww.ipac.caltech.edu/level5/GALEX_Atlas/

We combine data from SDSS and the GALEX and Spitzer observatories to create a sample of galaxies observed homogeneously from the UV to the Far-IR. This sample, consisting of ~460 galaxies observed spectroscopically by SDSS provides a multiwavelength (0.15-24 micron) view of obscured and unobscured star formation in nearby (z<0.3) galaxies with SFRs ranging from 0.01 to 100 M_solar/yr. We calculate a robust dust measure from the infrared to UV ratio (IRX) and explore the influence of star formation history (SFH) on the dust-UV color relation (i.e. the IRX-beta relation). We find that the UV colors of galaxies are only weakly dependent on their SFH as measured by the 4000A break. However, we find that the contributions of dust and SFH are distinguishable when colors at widely separated wavelengths (e.g. 0.23-3.6 micron) are introduced. We show this explicitly by recasting the IRX-beta relation as a more general IRX-SFH-color relation, which we examine in different projections. We also determine simple fits to this relation.

We study two samples of local galaxies, one is UV (GALEX) selected and the other FIR (IRAS) selected, to address the question whether UV and FIR surveys see the two sides ('bright' and 'dark') of the star formation of the same population of galaxies or two different populations of star forming galaxies. No significant difference between the L$_{tot}$ ($=L_{60}+L_{FUV}$) luminosity functions of the UV and FIR samples is found. Also, after the correction for the `Malmquist bias' (bias for flux limited samples), the FIR-to-UV ratio v.s. L$_{tot}$ relations of the two samples are consistent with each other. In the range of $9 \la \log(L_{tot}/L_\sun) \la 12$, both can be approximated by a simple linear relation of $\log (L_{60}/L_{FUV})=\log(L_{tot}/L_\sun)-9.66$. These are consistent with the hypothesis that the two samples represent the same population of star forming galaxies, and their well documented differences in L$_{tot}$ and in FIR-to-UV ratio are due only to the selection effect. A comparison between the UV luminosity functions shows marginal evidence for a population of faint UV galaxies missing in the FIR selected sample. The contribution from these 'FIR-quiet' galaxies to the overall UV population is insignificant, given that the K-band luminosity functions (i.e. the stellar mass functions) of the two samples do not show any significant difference.

[Abridged] We use mock catalogues based on the GALICS model (Hatton et al. 03) to explore the nature of galaxy clustering observed in the SDSS. We measure low and high order angular clustering statistic from these mock catalogues, after selecting galaxies the same way as for observations, and compare them directly to estimates from SDSS data. Note that we also present measurements of S3-S5 on the SDSS DR1. We find that our model is in general good agreement with observations in the scale/luminosity range where we can trust the predictions. This range is found to be limited (i) by the size of the dark matter simulation used -- which introduces finite volume effects at large scales -- and by the mass resolution of this simulation -- which introduces incompleteness at apparent magnitudes fainter than $r\sim 20$. We then focus on the small scale clustering properties of galaxies and investigate the behaviour of three different prescriptions for positioning galaxies within haloes of dark matter. We show that galaxies are poor tracers both of DM particles or DM sub-structures, within groups and clusters. Instead, SDSS data tells us that the distribution of galaxies lies somewhat in between these two populations. This confirms the general theoretical expectation from numerical simulations and semi-analytic modelling.

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 determine the number counts and z=0-5 luminosity function for a well-defined, homogeneous sample of quasars from the Sloan Digital Sky Survey (SDSS). We conservatively define the most uniform statistical sample possible, consisting of 15,343 quasars within an effective area of 1622 deg^2 that was derived from a parent sample of 46,420 spectroscopically confirmed broad-line quasars in the 5282 deg^2 of imaging data from SDSS Data Release Three. The sample extends from i=15 to i=19.1 at z<3 and to i=20.2 for z>3. The number counts and luminosity function agree well with the results of the 2dF QSO Survey, but the SDSS data probe to much higher redshifts than does the 2dF sample. The number density of luminous quasars peaks between redshifts 2 and 3, although uncertainties in the selection function in this range do not allow us to determine the peak redshift more precisely. Our best fit model has a flatter bright end slope at high redshift than at low redshift. For z<2.4 the data are best fit by a redshift-independent slope of beta = -3.1 (Phi(L) propto L^beta). Above z=2.4 the slope flattens with redshift to beta=-2.37 at z=5. This slope change, which is significant at a >5-sigma level, must be accounted for in models of the evolution of accretion onto supermassive black holes.

We use \textit{GALEX} (Galaxy Evolution Explorer) near-UV (NUV) photometry of a sample of early-type galaxies selected in \textit{SDSS} (Sloan Digital Sky Survey) to study the UV color-magnitude relation (CMR). $NUV-r$ color is an excellent tracer of even small amounts ($\sim 1$% mass fraction) of recent ($\la 1$ Gyr) star formation and so the $NUV-r$ CMR allows us to study the effect of environment on the recent star formation history. We analyze a volume-limited sample of 839 visually-inspected early-type galaxies in the redshift range $0.05 < z < 0.10$ brighter than $M_{r}$ of -21.5 with any possible emission-line or radio-selected AGN removed to avoid contamination. We find that contamination by AGN candidates and late-type interlopers highly bias any study of recent star formation in early-type galaxies and that, after removing those, our lower limit to the fraction of massive early-type galaxies showing signs of recent star formation is roughly $30 \pm 3%$ This suggests that residual star formation is common even amongst the present day early-type galaxy population. We find that the fraction of UV-bright early-type galaxies is 25% higher in low-density environments. However, the density effect is clear only in the lowest density bin. The blue galaxy fraction for the subsample of the brightest early-type galaxies however shows a very strong density dependence, in the sense that the blue galaxy fraction is lower in a higher density region.

This work presents the main ultraviolet (UV) and far-infrared (FIR) properties of two samples of nearby galaxies selected from the GALEX ($\lambda = 2315$\AA, hereafter NUV) and IRAS ($\lambda = 60\mu$m) surveys respectively. They are built in order to get detection at both wavelengths for most of the galaxies. Star formation rate (SFR) estimators based on the UV and FIR emissions are compared. Systematic differences are found between the SFR estimators for individual galaxies based on the NUV fluxes corrected for dust attenuation and on the total IR luminosity. A combined estimator based on NUV and IR luminosities seems to be the best proxy over the whole range of values of SFR. Although both samples present similar average values of the birthrate parameter b, their star-formation-related properties are substantially different: NUV-selected galaxies tend to show larger values of $b$ for lower masses, SFRs and dust attenuations, supporting previous scenarios for the star formation history (SFH). Conversely, about 20% of the FIR-selected galaxies show high values of $b$, SFR and NUV attenuation. These galaxies, most of them being LIRGs and ULIRGs, break down the downsizing picture for the SFH, however their relative contribution per unit volume is small in the local Universe. Finally, the cosmic SFR density of the local Universe is estimated in a consistent way from the NUV and IR luminosities.

We present the first large, unbiased sample of Lyman Break Galaxies (LBGs) at z ~ 1. Far ultraviolet-dropout (1530 A) galaxies in the Chandra Deep Field South have been selected using GALEX data. This first large sample in the z ~ 1 universe provides us with a high quality reference sample of LBGs. We analyzed the sample from the UV to the IR using GALEX, SPITZER, ESO and HST data. The morphology (obtained from GOODS data) of 75 % of our LBGs is consistent with a disk. The vast majority of LBGs with an IR detection are also Luminous Infrared Galaxies (LIRGs). As a class, the galaxies not detected at 24 microns are an order of magnitude fainter relative to the UV compared with those detected individually, suggesting that there may be two types of behavior within the sample. For the IR-bright galaxies, there is an apparent upper limit for the UV dust attenuation and this upper limit is anti-correlated with the observed UV luminosity. Previous estimates of dust attenuations based on the ultraviolet slope are compared to new ones based on the FIR/UV ratio (for LBGs detected at 24 microns), which is usually a more reliable estimator. Depending on the calibration we use to estimate the total IR luminosity, beta-based attenuations A_{FUV} are larger by 0.2 to 0.6 mag. than the ones estimated from FIR/UV ratio. Finally, for IR-bright LBGs, median estimated beta-based SFRs are 2-3 times larger than the total SFRs estimated as SFR_{TOT} = SFR_{UV} + SFR_{IR} while IR-based SFRs provide values below SFR_{TOT} by 15 - 20 %. We use a stacking method to statistically constrain the 24 microns flux of LBGs non individually detected. The results suggest that these LBGs do not contain large amounts of dust.

In this paper we describe a detailed analysis of the photometric uncertainties present within the Sloan Digital Sky Survey (SDSS) imaging survey based on repeat observations of approximately 200 square degrees of the sky. We show that, for the standard SDSS aperture systems (petrocounts, counts_model, psfcounts and cmodel_counts), the errors generated by the SDSS photometric pipeline under-estimate the observed scatter in the individual bands. The degree of disagreement is a strong function of aperture and magnitude (ranging from 20% to more than a factor of 2). We also find that the photometry in the five optical bands can be highly correlated for both point sources and galaxies, although the correlation for point sources is almost entirely due to variable objects. Without correcting for this covariance the SDSS color errors could be in over-estimated by a factor of two to three. Combining these opposing effects, the SDSS errors on the colors differ from the observed color variation by approximately 10-20% for most apertures and magnitudes. We provide a prescription to correct the errors derived from the SDSS photometric pipeline as a function of magnitude and a semi-analytic method for generating the appropriate covariance between the different photometric passbands. Given the intrinsic nature of these correlations, we expect that all current and future multi-band surveys will also observe strongly covariant magnitudes. The ability of these surveys to complete their science goals is largely dependent on color-based target selection and photometric redshifts; these results show the importance of spending a significant fraction of early survey operations on re-imaging to empirically determine the photometric covariance of any observing/reduction pipeline.

Using an eigenspectrum decomposition technique, we separate the host galaxy from the broad line active galactic nucleus (AGN) in a set of 4666 spectra from the Sloan Digital Sky Survey (SDSS), from redshifts near zero up to about 0.75. The decomposition technique uses separate sets of galaxy and quasar eigenspectra to efficiently and reliably separate the AGN and host spectroscopic components. The technique accurately reproduces the host galaxy spectrum, its contributing fraction, and its classification. We show how the accuracy of the decomposition depends upon S/N, host galaxy fraction, and the galaxy class. Based on the eigencoefficients, the sample of SDSS broad-line AGN host galaxies spans a wide range of spectral types, but the distribution differs significantly from inactive galaxies. In particular, post-starburst activity appears to be much more common among AGN host galaxies. The luminosities of the hosts are much higher than expected for normal early-type galaxies, and their colors become increasingly bluer than early-type galaxies with increasing host luminosity. Most of the AGNs with detected hosts are emitting at between 1% and 10% of their estimated Eddington luminosities, but the sensitivity of the technique usually does not extend to the Eddington limit. There are mild correlations among the AGN and host galaxy eigencoefficients, possibly indicating a link between recent star formation and the onset of AGN activity. The catalog of spectral reconstruction parameters is available as an electronic table.

We present a first analysis of the clustering of SDSS galaxies using the distribution function of the sum of Fourier phases. This statistic was recently proposed by one of authors as a new method to probe phase correlations of cosmological density fields. Since the Fourier phases are statistically independent of Fourier amplitudes, the phase statistic plays a complementary role to the conventional two-point statistics of galaxy clustering. In particular, we focus on the distribution functions of phase sum over three closed wavevectors as a function of the configurations of triangle wavevectors. We find that the observed distribution functions of phase sum are in good agreement with the lowest-order approximation from perturbation theory. For direct comparison with observations, we construct mock catalogs from N-body simulations taking account of several observational effects such as the survey geometry, the redshift distortion, and the discreteness due to the limited number of data. Indeed the observed phase correlations for the galaxies in the range of absolute magnitude -22<Mr<-18 agree well with those for Lambda CDM predictions with sigma8=0.9 evolved from Gaussian initial condition. This agreement implies that the galaxy biasing is approximately linear in redshift space. Instead, assuming that the galaxy biasing is described by a quadratic deterministic function at scale larger than 30Mpc/h, we can constrain the ratio of the quadratic biasing parameter b2 to the linear biasing parameter b1. We find that resulting b2/b1 is mainly dependent on sigma8 and that it is well fitted by b2/b1=0.54(+-0.06)-0.62(+-0.08)sigma8 in CDM models. Indeed, b2/b1 is nearly zero when sigma8=0.9.

The Galaxy Evolution Explorer (GALEX) satellite has obtained high time resolution ultraviolet photometry during a large flare on the M4 dwarf star GJ 3685A. Simultaneous NUV (1750 - 2800A) and FUV (1350 - 1750A) time-tagged photometry with time resolution better than 0.1 s shows that the overall brightness in the FUV band increased by a factor of 1000 in 200 s. Under the assumption that the NUV emission is mostly due to a stellar continuum, and that the FUV flux is shared equally between emission lines and continuum, then there is evidence for two distinct flare components for this event. The first flare type is characterized by an exponential increase in flux with little or no increase in temperature. The other involves rapid increases in both temperature and flux. While the decay time for the first flare component may be several hours, the second flare event decayed over less than 1 minute, suggesting that there was little or no confinement of the heated plasma.

We report the discovery of a new two-image gravitationally lensed quasar, SDSS J024634.11-082536.2 (SDSS J0246-0825). This object was selected as a lensed quasar candidate from the Sloan Digital Sky Survey (SDSS) by the same algorithm that was used to discover other SDSS lensed quasars (e.g., SDSS J0924+0219). Multicolor imaging with the Magellan Consortium's Walter Baade 6.5-m telescope and the spectroscopic observations using the W. M. Keck Observatory's Keck II telescope confirm that SDSS J0246-0825 consists of two lensed images ($\Delta{\theta}=$1\farcs04) of a source quasar at z=1.68. Imaging observations with the Keck telescope and the Hubble Space Telescope reveal an extended object between the two quasar components, which is likely to be a lensing galaxy of this system. From the absorption lines in the spectra of quasar components and the apparent magnitude of the galaxy, combined with the expected absolute magnitude from the Faber-Jackson relation, we estimate the redshift of the lensing galaxy to be z=0.724. A highly distorted ring is visible in the Hubble Space Telescope images, which is likely to be the lensed host galaxy of the source quasar. Simple mass modeling predicts the possibility that there is a small (faint) lensing object near the primary lensing galaxy.

The DEEP2 and COMBO-17 surveys are used to study the evolution of the luminosity function of red and blue galaxies to $z \sim 1$. Schechter function fits show that, since $z = 1$, $M^*_B$ dims by $\sim$ 1.3 mag per unit redshift for both color classes, $\phi^*$ of blue galaxies shows little change, while $\phi^*$ for red galaxies has formally nearly quadrupled. At face value, the number density of blue galaxies has remained roughly constant since $ z = 1$, whereas that of red galaxies has been rising. Luminosity densities support both conclusions, but we note that most red-galaxy evolution occurs between our data and local surveys and in our highest redshift bin, where the data are weakest. We discuss the implications of having most red galaxies emerge after $z = 1$ from precursors among the blue population, taking into account the properties of local and distant E/S0s. We suggest a ``mixed'' scenario in which some blue galaxies have their star-formation quenched in gas-rich mergers, migrate to the red sequence with a variety of masses, and merge further on the red sequence in one or more purely stellar mergers. E/S0s of a given mass today will have formed via different routes, in a manner that may help to explain the fundamental plane and other local scaling laws.

We discuss the construction of a photometric redshift catalogue of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS), emphasizing the principal steps necessary for constructing such a catalogue -- (i) photometrically selecting the sample, (ii) measuring photometric redshifts and their error distributions, (iii) and estimating the true redshift distribution. We compare two photometric redshift algorithms for these data and find that they give comparable results. Calibrating against the SDSS and SDSS-2dF spectroscopic surveys, we find that the photometric redshift accuracy is $\sigma \sim 0.03$ for redshifts less than 0.55 and worsens at higher redshift ($\sim 0.06$). These errors are caused by photometric scatter, as well as systematic errors in the templates, filter curves, and photometric zeropoints. We also parametrize the photometric redshift error distribution with a sum of Gaussians, and use this model to deconvolve the errors from the measured photometric redshift distribution to estimate the true redshift distribution. We pay special attention to the stability of this deconvolution, regularizing the method with a prior on the smoothness of the true redshift distribution. The methods we develop are applicable to general photometric redshift surveys.

We report the first result of the clustering analysis of Sloan Digital Sky Survey (SDSS) quasars. We compute the two-point correlation function (2PCF) of SDSS quasars in redshift space at $8h^{-1}{\rm Mpc} < s < 500h^{-1}{\rm Mpc}$, with particular attention to its baryonic signature. Our sample consists of 19986 quasars extracted from the SDSS Data Release 4 (DR4). The redshift range of the sample is $0.72 \le z \le 2.24$ (the mean redshift is $\bar z = 1.46$) and the reddening-corrected $i$-band apparent magnitude range is $15.0 \le m_{i,{\rm rc}} \le 19.1$. Due to the relatively low number density of the quasar sample, the bump in the power spectrum due to the baryon density, $\Omega_{\rm b}$, is not clearly visible. The effect of the baryon density is, however, to distort the overall shape of the 2PCF.The degree of distortion makes it an interesting alternate measure of the baryonic signature. Assuming a scale-independent linear bias and the spatially flat universe, i.e., $\Omega_{\rm b} + \Omega_{\rm d} + \Omega_\Lambda =1$, where $\Omega_{\rm d}$ and $\Omega_\Lambda$ denote the density parameters of dark matter and the cosmological constant, we combine the observed quasar 2PCF and the predicted matter 2PCF to put constraints on $\Omega_{\rm b}$ and $\Omega_\Lambda$. Our result is fitted as $0.80- 2.8\Omega_{\rm b} < \Omega_\Lambda < 0.90 - 1.4\Omega_{\rm b}$ at the 2$\sigma$ confidence level, which is consistent with results from other cosmological observations such as WMAP. (abridged)

We present Version 1.0 of the NASA Galaxy Evolution Explorer (GALEX) ultraviolet variability catalog (GUVV) that contains information on 84 time-variable and transient sources gained with simultaneous near and far ultraviolet photometric observations. These time-variable sources were serendipitously revealed in the various 1.2 degree star fields currently being surveyed by the GALEX satellite in two ultraviolet bands (NUV 1750-2750A, FUV 1350-1750A) with limiting AB magnitudes of 23-25. The largest-amplitude variable objects presently detected by GALEX are M-dwarf flare stars, which can brighten by 5-10 mag in both the NUV and FUV bands during short duration (< 500s) outbursts. Other types of large-amplitude ultraviolet variable objects include ab-type RR Lyrae stars, which can vary periodically by 2-5mag in the GALEX FUV band. This first GUVV catalog lists galactic positions and possible source identifications in order to provide the astronomical community with a list of time-variable objects that can now be repeatedly observed at other wavelengths. We expect the total number of time-variable source detections to increase as the GALEX mission progresses, such that later version numbers of the GUVV catalog will contain substantially more variable sources.

The ELAIS S1 field was observed by GALEX in both its Wide Spectroscopic and Deep Imaging Survey modes. This field was previously observed by the Infrared Space Observatory and we made use of the catalogue of multi-wavelength data published by the ELAIS consortium to select galaxies common to the two samples. Among the 959 objects with GALEX spectroscopy, 88 are present in the ELAIS catalog and 19 are galaxies with an optical spectroscopic redshift. The distribution of redshifts covers the range $0<z<1.6$. The selected galaxies have bolometric IR luminosities $10<Log(L_{IR})<13$ (deduced from the $15 \mu m$ flux using ISOCAM) which means that we cover a wide range of galaxies from normal to Ultra Luminous IR Galaxies. The mean ($\sigma$) UV luminosity (not corrected for extinction) amounts to $Log(\lambda.L_{1530}) = 9.8 (0.6)$ $L_\sun$ for the low-z ($z \le 0.35$) sample. The UV slope $\beta$ (assuming $f_\lambda \propto \lambda^\beta$) correlates with the GALEX FUV-NUV color if the sample is restricted to galaxies below $z < 0.1$. Taking advantage of the UV and IR data, we estimate the dust attenuation from the IR/UV ratio and compare it to the UV slope $\beta$. We find that it is not possible to uniquely estimate the dust attenuation from $\beta$ for our sample of galaxies. These galaxies are highly extinguished with a median value $A_{FUV} = 2.7 \pm 0.8$. Once the dust correction applied, the UV- and IR-based SFRs correlate. For the closest galaxy with the best quality spectrum, we see a feature consistent with being produced by a bump near 220nm in the attenuation curve.

We have used the 2dF instrument on the AAT to obtain redshifts of a sample of z<3, 18.0<g<21.85 quasars selected from SDSS imaging. These data are part of a larger joint programme: the 2dF-SDSS LRG and QSO Survey (2SLAQ). We describe the quasar selection algorithm and present the resulting luminosity function of 5645 quasars in 105.7 deg^2. The bright end number counts and luminosity function agree well with determinations from the 2dF QSO Redshift Survey (2QZ) data to g\sim20.2. However, at the faint end the 2SLAQ number counts and luminosity function are steeper than the final 2QZ results from Croom et al. (2004), but are consistent with the preliminary 2QZ results from Boyle et al. (2000). Using the functional form adopted for the 2QZ analysis, we find a faint end slope of beta=-1.78+/-0.03 if we allow all of the parameters to vary and beta=-1.45+/-0.03 if we allow only the faint end slope and normalization to vary. Our maximum likelihood fit to the data yields 32% more quasars than the final 2QZ parameterization, but is not inconsistent with other g>21 deep surveys. The 2SLAQ data exhibit no well defined ``break'' but do clearly flatten with increasing magnitude. The shape of the quasar luminosity function derived from 2SLAQ is in good agreement with that derived from type I quasars found in hard X-ray surveys. [Abridged]

We present the third edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog. The catalog consists of the 46,420 objects in the SDSS Third Data Release that have luminosities larger than M_i = -22 (in a cosmology with H_0 = 70 km/s/Mpc, Omega_M = 0.3, and Omega_Lambda = 0.7), have at least one emission line with FWHM larger than 1000 km/s or are unambiguously broad absorption line quasars, are fainter than i = 15.0, and have highly reliable redshifts. The area covered by the catalog is 4188 sq. deg. The quasar redshifts range from 0.08 to 5.41, with a median value of 1.47; the high-redshift sample includes 520 quasars at redshifts greater than four, of which 17 are at redshifts greater than five. For each object the catalog presents positions accurate to better than 0.2 arcsec. rms per coordinate, five-band (ugriz) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains radio, near-infrared, and X-ray emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3800--9200A at a spectral resolution about 2000; the spectra can be retrieved from the public database using the information provided in the catalog. A total of 44,221 objects in the catalog were discovered by the SDSS; 28,400 of the SDSS discoveries are reported here for the first time.

The Deep Extragalactic Exploratory Probe (DEEP) is a multi-phase research program dedicated to the study of the formation and evolution of galaxies and of large scale structure in the distant Universe. This paper describes the first five-year phase, denoted DEEP1. A series of ten DEEP1 papers will discuss a range of scientific topics (e.g., the study of photometric and spectral properties of a general distant galaxy survey, the evolution observed in galaxy populations of varied morphologies). The observational basis for these studies is the Groth Survey Strip field, a 127 square arcminute region which has been observed with the Hubble Space Telescope in both broad I-band and V-band optical filters and with the Low Resolution Imaging Spectrograph on the Keck Telescopes. Catalogs of photometric and structural parameters have been constructed for 11,547 galaxies and stars at magnitudes brighter than 29, and spectroscopy has been conducted for a magnitude-color weighted subsample of 818 objects. We evaluate three independent techniques for constructing an imaging catalog for the field from the HST data, and discuss the depth and sampling of the resultant catalogs. The selection of the spectroscopic subsample is discussed, and we describe the multifaceted approach taken to prioritizing objects of interest for a variety of scientific subprograms. A series of Monte Carlo simulations then demonstrates that the spectroscopic subsample can be adequately modeled as a simple function of magnitude and color cuts in the imaging catalog.

We present the GALEX NUV (2310 A) and FUV (1530 A) galaxy luminosity functions of the nearby cluster of galaxies A1367 in the magnitude range -20.3< M_AB < -13.3. The luminosity functions are consistent with previous (~ 2 mag shallower) estimates based on the FOCA and FAUST experiments, but display a steeper faint-end slope than the GALEX luminosity function for local field galaxies. Using spectro-photometric optical data we select out star-forming systems from quiescent galaxies and study their separate contributions to the cluster luminosity function. We find that the UV luminosity function of cluster star-forming galaxies is consistent with the field. The difference between the cluster and field LF is entirely due to the contribution at low luminosities (M_AB >-16 mag) of non star-forming, early-type galaxies that are significantly over dense in clusters.

We present GALEX NUV (2310 A) and FUV (1530 A) images of the interacting galaxy NGC 4438 (Arp 120) in the center of the Virgo cluster. These images show an extended (20 kpc) tidal tail at the north-west edge of the galaxy previously undetected at other wavelengths, at 15-25 kpc from its nucleus. Except in the nucleus, the UV morphology of NGC 4438 is totally different from the Halpha+[NII] one, more similar to the X-ray emission, confirming its gas cooling origin. We study the star formation history of NGC 4438 combining spectro-photometric data in the UV-visible-near-IR wavelength range with population synthesis and galaxy evolution models. The data are consistent with a recent (~ 10 Myr), instantaneous burst of star formation in the newly discovered UV north-western tail which is significantly younger than the age of the tidal interaction with NGC 4435, dated by dynamical models at ~ 100 Myr ago. Recent star formation events are also present at the edge of the northern arm and in the southern tail, while totally lacking in the other regions, which are dominated by the old stellar population perturbed during the dynamical interaction with NGC 4435. The contribution of this recent starburst to the total galaxy stellar mass is lower than 0.1%, an extremely low value for such a violent interaction. High-velocity, off-center tidal encounters such as that observed in Arp 120 are thus not sufficient to significantly increase the star formation activity of cluster galaxies.

We present the large-scale correlation function measured from a spectroscopic sample of 46,748 luminous red galaxies from the Sloan Digital Sky Survey. The survey region covers 0.72 h^{-3} Gpc^3 over 3816 square degrees and 0.16<z<0.47, making it the best sample yet for the study of large-scale structure. We find a well-detected peak in the correlation function at 100h^{-1} Mpc separation that is an excellent match to the predicted shape and location of the imprint of the recombination-epoch acoustic oscillations on the low-redshift clustering of matter. This detection demonstrates the linear growth of structure by gravitational instability between z=1000 and the present and confirms a firm prediction of the standard cosmological theory. The acoustic peak provides a standard ruler by which we can measure the ratio of the distances to z=0.35 and z=1089 to 4% fractional accuracy and the absolute distance to z=0.35 to 5% accuracy. From the overall shape of the correlation function, we measure the matter density Omega_mh^2 to 8% and find agreement with the value from cosmic microwave background (CMB) anisotropies. Independent of the constraints provided by the CMB acoustic scale, we find Omega_m = 0.273 +- 0.025 + 0.123 (1+w_0) + 0.137 Omega_K. Including the CMB acoustic scale, we find that the spatial curvature is Omega_K=-0.010+-0.009 if the dark energy is a cosmological constant. More generally, our results provide a measurement of cosmological distance, and hence an argument for dark energy, based on a geometric method with the same simple physics as the microwave background anisotropies. The standard cosmological model convincingly passes these new and robust tests of its fundamental properties.

We have used the first matched set of GALEX and SDSS data to investigate the properties of a sample of 74 nearby galaxies with far-ultraviolet luminosities chosen to overlap the luminosity range of typical high-z Lyman Break Galaxies (LBGs). GALEX deep surveys have shown that ultraviolet-luminous galaxies (UVLGs) similar to these are the fastest evolving component of the UV galaxy population. Model fits to the combined GALEX and SDSS photometry yield typical FUV extinctions similar to LBGs. The implied star formation rates are SFR ~ 3 to 30 solar mass per year. This overlaps the range of SFRs for LBGs. We find a strong inverse correlation between galaxy mass and far-ultraviolet surface brightness, and on this basis divide the sample into ``large'' and ``compact'' UVLGs. The compact UVLGs have half-light radii of a few kpc or less (similar to LBGs). They are relatively low mass galaxies (~10 billion solar masses) with typical velocity dispersions of 60 to 150 km/s. They span a range in metallicity from 0.3 to 1 times solar, have blue optical-UV colors, and are forming stars at a rate sufficient to build the present galaxy in ~a Gigayear. In all these respects they appear similar to the LBG population. These ``living fossils'' may therefore provide an opportunity for detailed investigation of the physical processes occurring in typical star forming galaxies in the early universe.

We use the GALEX (Galaxy Evolution Explorer) Medium Imaging Survey (MIS) and All-Sky Imaging Survey (AIS) data available in the first internal release, matched to the SDSS catalogs in the overlapping regions, to classify objects by comparing the multi-band photometry to model colors. We show an example of the advantage of such broad wavelength coverage (GALEX far-UV and near-UV, SDSS ugriz) in classifying objects and augmenting the existing samples and catalogs. From the MIS [AIS] sample over an area of 75 [92] square degrees, we select a total of 1736 [222] QSO candidates at redshift less than 2, significantly extending the number of fainter candidates, and moderately increasing the number of bright objects in the SDSS list of spectroscopically confirmed QSO. Numerous hot stellar objects are also revealed by the UV colors, as expected.

The GALEX (Galaxy Evolution Explorer) Nearby Galaxies Survey is providing deep far-UV and near-UV imaging for a representative sample of galaxies in the local universe. We present early results for M51 and M101, from GALEX UV imaging and SDSS optical data in five bands. The multi-band photometry of compact stellar complexes in M101 is compared to population synthesis models, to derive ages, reddening, reddening-corrected luminosities and current/initial masses. The GALEX UV photometry provides a complete census of young compact complexes on a approximately 160pc scale. A galactocentric gradient of the far-UV - near-UV color indicates younger stellar populations towards the outer parts of the galaxy disks, the effect being more pronounced in M101 than in M51.

New GALEX observations have detected significant FUV (1530 Angstroms) and NUV (2310 Angstroms) emission from stellar substructures within the tidal tails of four ongoing galaxy mergers. The UV-bright regions are optically faint and are coincident with HI density enhancements. FUV emission is detected at any location where the HI surface density exceeds ~2 M_sun pc^{-2}, and is often detected in the absence of visible wavelength emission. UV luminosities of the brighter regions of the tidal tails imply masses of 10^6 M_sun to ~10^9 M_sun in young stars in the tails, and HI luminosities imply similar HI masses. UV-optical colors of the tidal tails indicate stellar populations as young as a few Myr, and in all cases ages < 400Myr. Most of the young stars in the tails formed in single bursts rather than resulting from continuous star formation, and they formed *in situ* as the tails evolved. Star formation appears to be older near the parent galaxies and younger at increasing distances from the parent galaxy. This could be because the star formation occurs progressively along the tails, or because the star formation has been inhibited near the galaxy/tail interface. The youngest stellar concentrations, usually near the ends of long tidal tails, have masses comparable to confirmed tidal dwarf galaxies and may be newly forming galaxies undergoing their first burst of star formation.

We use the far-UV (FUV) and near-UV (NUV) images of M83 obtained by GALEX to compute the radial profile of the UV spectral slope in the star forming disk. We briefly present a model of its chemical evolution which allows us to obtain realistic intrinsic properties of the stellar populations. Using corollary data, we also compute the profiles of H-alpha/H-beta and of the total IR (TIR) to FUV ratio. Both data and model are used to estimate and compare the extinction gradients at the FUV wavelength obtained from these various indicators. We discuss the implications for the determination of the star formation rate.

We present ultraviolet photometry for globular clusters (GCs) in M31 from 15 square deg of imaging using the Galaxy Evolution Explorer (GALEX). We detect 200 and 94 GCs with certainty in the near-ultraviolet (NUV; 1750 - 2750 Angstroms) and far-ultraviolet (FUV; 1350 - 1750 Angstroms) bandpasses, respectively. Our rate of detection is about 50% in the NUV and 23% in the FUV, to an approximate limiting V magnitude of 19. Out of six clusters with [Fe/H]>-1 seen in the NUV, none is detected in the FUV bandpass. Furthermore, we find no candidate metal-rich clusters with significant FUV flux, because of the contribution of blue horizontal-branch (HB) stars, such as NGC 6388 and NGC 6441, which are metal-rich Galactic GCs with hot HB stars. We show that our GALEX photometry follows the general color trends established in previous UV studies of GCs in M31 and the Galaxy. Comparing our data with Galactic GCs in the UV and with population synthesis models, we suggest that the age range of M31 and Galactic halo GCs are similar.

We have used the GALEX UV photometric data to construct a first near-ultraviolet (NUV) color-magnitude relation (CMR) for the galaxies pre-classified as early-type by SDSS studies. The NUV CMR is a powerful tool for tracking the recent star formation history in early-type galaxies, owing to its high sensitivity to the presence of young stellar populations. Our NUV CMR for UV-weak galaxies shows a well-defined slope and thus will be useful for interpreting the restframe NUV data of distant galaxies and studying their star formation history. Compared to optical CMRs, the NUV CMR shows a substantially larger scatter, which we interpret as evidence of recent star formation activities. Roughly 15% of the recent epoch (z < 0.13) bright (M[r] < -22) early-type galaxies show a sign of recent (< 1Gyr) star formation at the 1-2% level (lower limit) in mass compared to the total stellar mass. This implies that low level residual star formation was common during the last few billion years even in bright early-type galaxies.

We derive a variety of physical parameters including star formation rates (SFRs), dust attenuation and burst mass fractions for 6472 galaxies observed by the Galaxy Evolution Explorer (GALEX) and present in the SDSS DR1 main spectroscopic sample. Parameters are estimated in a statistical way by comparing each observed broad-band SED (two GALEX and five SDSS bands) with an extensive library of model galaxy SEDs, which cover a wide range of star formation histories and include stochastic starbursts. We compare the constraints derived using SDSS bands only with those derived using the combination of SDSS and GALEX photometry. We find that the addition of the GALEX bands leads to significant improvement in the estimation of both the dust optical depth and the star formation rate over timescales of 100 Myr to 1 Gyr in a galaxy. We are sensitive to SFRs as low as 10^{-3} M_sun/yr, and we find that low levels of star formation (SF) are mostly associated with early-type, red galaxies. The least massive galaxies have ratios of current to past-averaged SF rates (b-parameter) consistent with constant SF over a Hubble time. For late-type galaxies, this ratio on average decreases with mass. We find that b correlates tightly with NUV-r color, implying that the SF history of a galaxy can be constrained on the basis of the NUV-r color alone. The fraction of galaxies that have undergone a significant starburst episode within the last 1 Gyr steeply declines with mass-from ~20% for galaxies with ~10^8 M_sun to ~5% for ~10^11 M_sun galaxies.

We compare the dust attenuation properties of two samples of galaxies purely selected in the near-ultraviolet (NUV) band (1750-2750 A, lambda_m = 2310 A) and in the far-infrared (FIR) at 60micron. These samples are built using the GALEX and IRAS sky surveys over ~600 square degrees. The NUV selected sample contains 95 galaxies detected down to NUV=16mag (AB system). 83 galaxies in this sample are spirals or irregulars and only two of them are not detected at 60micron. The FIR selected sample is built from the IRAS PSCz catalog complete down to 0.6Jy. Among the 163 sources, we select 118 star forming galaxies well measured by IRAS, all but 1 are detected in NUV and 14 galaxies are not detected in the far-ultraviolet (FUV) band (1350-1750 A, lambda_m = 1530 A). The dust to ultraviolet (NUV and FUV) flux ratio is calibrated to estimate the dust attenuation at both wavelengths. The median value of the attenuation in NUV is found to be ~1 mag for the NUV selected sample versus ~2 mag for the FIR selected one. Within both samples, the dust attenuation is found to correlate with the luminosity of the galaxies. Almost all the NUV selected galaxies and 2/3 of the FIR selected sample exhibit a lower dust attenuation than that expected from the tight relation found previously for starburst galaxies between the dust attenuation and the slope of the ultraviolet continuum. The situation is inverse for one third of the FIR selected galaxies: their extinction is higher than that deduced from their FUV-NUV color and the relation valid for starbursts.

We present ultraviolet photometry for a sample of morphologically early-type galaxies selected by matching the Sloan Digital Sky Survey Data Release 1 with the GALEX Medium and All-sky Imaging Surveys. We obtain a working sample of 1032 early-type galaxies with GALEX FUV detections, SDSS spectroscopy, and $z<0.2$. Using the SDSS spectra to identify galaxies with even weak star formation or evidence of AGN, and further removing galaxies with any evidence of non early-type morphology, we derive a final sample of 172 red quiescent early-type galaxies. We find that the $FUV-r$ color has a full range of 5 mag. Plotting against the $FUV-r$ color the metallicity sensitive Lick $\rm Mg_2$ and D4000 indices, and the stellar velocity dispersion, we find no correlation between our measurement of UV rising flux, and any parameter sensitive to metallicity.

We present the results of a determination of the galaxy luminosity function at ultraviolet wavelengths at redshifts of $z=0.0-0.1$ from GALEX data. We determined the luminosity function in the GALEX FUV and NUV bands from a sample of galaxies with UV magnitudes between 17 and 20 that are drawn from a total of 56.73 deg^2 of GALEX fields overlapping the b_j-selected 2dF Galaxy Redshift Survey. The resulting luminosity functions are fainter than previous UV estimates and result in total UV luminosity densities of 10^(25.55+/-0.12) ergs s^-1 Hz^-1 Mpc^-3 and 10^(25.72+/-0.12) ergs s^-1 Hz^-1 Mpc^-3 at 1530 Ang. and 2310 Ang., respectively. This corresponds to a local star formation rate density in agreement with previous estimates made with H-alpha-selected data for reasonable assumptions about the UV extinction.

The total ultraviolet (UV) flux (from 1412 to 2718 AA) of M 101 is compared on a pixel-to-pixel basis with the total far-infrared (FIR) flux (from 60 to 170 micron), using the maps of the galaxy taken by GALEX in the near-UV and far-UV and by ISOPHOT at 60, 100 and 170 micron. The main result of this investigation is the discovery of a tight dependence of the FIR/UV ratio on radius, with values monotonically decreasing from around 4 in the nuclear region to nearly zero towards the edge of the optical disk. Although the tightness of this dependence is in part attributable to resolution effects, the result is consistent with the presence of a large-scale distribution of diffuse dust having a face-on optical depth which decreases with radius and which dominates over the more localized variations in opacity between the arm and interarm regions. We also find a trend for the FIR/UV ratio to take on higher values in the regions of diffuse interarm emission than in the spiral-arm regions, at a given radius. This is interpreted quantitatively in terms of the escape probability of UV photons from spiral arms and their subsequent scattering in the interarm regions, and in terms of the larger relative contribution of optical photons to the heating of the dust in the interarm regions.

We present measurements of the UV galaxy luminosity function and the evolution of luminosity density from GALEX observations matched to the Sloan Digital Sky Survey (SDSS). We analyze galaxies in the Medium Imaging Survey overlapping the SDSS DR1 with a total coverage of 44 deg^2. Using the combined GALEX+SDSS photometry, we compute photometric redshifts and study the LF in three redshift shells between z=0.07 and 0.25. The Schechter function fits indicate that the faint-end slope alpha is consistent with -1.1 at all redshifts but the characteristic UV luminosity M* brightens by 0.2 mag from z=0.07 to 0.25. In the lowest redshift bin, early and late type galaxies are studied separately and we confirm that red galaxies tend to be brighter and have a shallower slope alpha than blue ones. The derived luminosity densities are consistent with other GALEX results based on a local spectroscopic sample from 2dF and the evolution follows the trend reported by deeper studies.

We present Galaxy Evolution Explorer (GALEX) far-UV and near-UV mosaic observations covering the entirety of M31 and M33. For both targets, we measure the decline of surface brightness (in FUV and NUV) and changes in FUV--NUV color as a function of galactocentric radius. These UV radial profiles are compared to the distribution of ionized gas traced by H-alpha emission. We find that the extent of the UV emission, in both targets, is greater than the extent of the observed HII regions and diffuse ionized gas. We determine the ultraviolet diffuse fraction in M33 using our FUV observations and compare it to the H-alpha diffuse fraction obtained from wide-field narrow-band imaging. The FUV diffuse fraction appears to be remarkably constant near 0.65 over a large range in galactocentric radius, with departures to higher values in circumnuclear regions and, most notably, at the limit of the H-alpha disk. We suggest that the increase in FUV diffuse fraction at large galactocentric radii could indicate that a substantial portion of the diffuse emission beyond this point is not generated in situ but rather scattered from dust, after originating in the vicinity of the disk's outermost HII regions. Radial variation of the H-alpha diffuse fraction was also measured. We found the H-alpha diffuse fraction generally near 0.4 but rising toward the galaxy center, up to 0.6. We made no attempt to correct our diffuse fraction measurements for position-dependent extinction, so the quoted values are best interpreted as upper limits given the plausibly higher extinction for stellar clusters relative to their surroundings.

we present the first GALEX UV images of the well known interacting group of galaxies, Stephan's Quintet (SQ). We detect widespread UV emission throughout the group. However, there is no consistent coincidence between UV structure and emission in the optical, H\alpha, or HI. Excluding the foreground galaxy NGC7320 (Sd), most of the UV emission is found in regions associated with the two spiral members of the group, NGC7319 and NGC7318b, and the intragroup medium starburst SQ-A. The extinction corrected UV data are analyzed to investigate the overall star formation activity in SQ. It is found that the total star formation rate (SFR) of SQ is 6.69+-0.65 M_\sun/yr. Among this, 1.34+-0.16 M_sun/yr is due to SQ-A. This is in excellent agreement with that derived from extinction corrected H\alpha luminosity of SQ-A. The SFR in regions related to NGC 7319 is 1.98+-0.58 M_\sun/yr, most of which(68%) is contributed by the disk. The contribution from the 'young tail' is only 15%. In the UV, the 'young tail' is more extended (~100 kpc) and shows a loop-like structure, including the optical tail, the extragalactic HII regions recently discovered in H\alpha, and other UV emission regions discovered for the first time. The UV and optical colors of the 'old tail' are consistent with a single stellar population of age t ~10^{8.5+-0.4} yrs. The UV emission associated with NGC 7318b is found in a very large (~80 kpc) disk, with a net SFR of 3.37+-0.25 M_sun/yr. Several large UV emission regions are 30 -- 40 kpc away from the nucleus of NGC7318b. Although both NGC7319 and NGC7318b show peculiar UV morphology, their SFR is consistent with that of normal Sbc galaxies, indicating that the strength of star formation activity is not enhenced by interactions.

We present Galaxy Evolution Explorer (GALEX) images of the prototypical edge-on starburst galaxies M82 and NGC253. Our initial analysis is restricted to the complex of ultraviolet (UV) filaments in the starburst-driven outflows in the galaxy halos. The UV luminosities in the halo are too high to be provided by shock-heated or photoionized gas except perhaps in the brightest filaments in M82, suggesting that most of the UV light is the stellar continuum of the starburst scattered into our line of sight by dust in the outflow. This interpretation agrees with previous results from optical imaging polarimetry in M82. The morphology of the UV filaments in both galaxies shows a high degree of spatial correlation with H-alpha and X-ray emission. This indicates that these outflows contain cold gas and dust, some of which may be vented into the intergalactic medium (IGM). UV light is seen in the ``H-alpha cap'' 11 kpc North of M82. If this cap is a result of the wind fluid running into a pre-existing gas cloud, the gas cloud contains dust and is not primordial in nature but was probably stripped from M82 or M81. If starburst winds efficiently expel dust into the IGM, this could have significant consequences for the observation of cosmologically distant objects.

We present measurements of the FUV (1530A) and NUV (2310A) galaxy luminosity functions (LF) at low redshift (z<0.2) from GALEX observations matched to the 2dF Galaxy Redshift Survey. We split our FUV and NUV samples into two UV-bj color bins and two redshift bins. As observed at optical wavelengths, the local LF of the bluest galaxies tend to have steeper faint end slopes and fainter characteristic magnitudes M* than the reddest subsamples. We find evidence for color dependent evolution at very low redshift in both bands, with bright blue galaxies becoming dominant in the highest redshift bin. The evolution of the total LF is consistent with an 0.3 magnitude brightening between z=0 and 0.13, in agreement with the first analysis of deeper GALEX fields probing adjacent and higher redshifts.

We have derived the bivariate luminosity function for the far ultraviolet (1530Angstroms) and far infrared (60 microns). We used matched GALEX and IRAS data, and redshifts from NED and PSC-z. We have derived a total star formation luminosity function phi(L_{tot}), with L_{tot} = L_{FUV}+L_{FIR}. Using these, we determined the cosmic ``star formation rate'' function and density for the local universe. The total SFR function is fit very well by a log-normal distribution over five decades of luminosity. We find that the bivariate luminosity function phi(L_{FUV},L_{FIR}) shows a bimodal behavior, with L_{FIR} tracking L_{FUV} for L_{TOT}< 10^10 L_sun, and L_{FUV} saturating at 10^10 L_sun, while L_{TOT} L_{FIR} for higher luminosities. We also calculate the SFR density and compare it to other measurements.

Ultraviolet imaging with the Galaxy Evolution Explorer (GALEX) has revealed an extensive sample of UV-bright stellar complexes in the extreme outer disk of M83, extending to about four times the radius where the majority of HII regions are detected (R_HII = 5.1' or 6.6 kpc). These sources are typically associated with large-scale filamentary HI structures in the warped outer disk of M83, and are distributed beyond the galactocentric radii at which molecular ISM has yet been detected. We present measured properties of these stellar complexes, including FUV and NUV magnitudes and local gas surface density. Only a subset of the outer disk UV sources have corresponding HII regions detected in H-alpha imaging, consistent with a sample of mixed age in which some sources are a few Myr old and others are much more evolved (~ 10^8 yr).

In order to investigate the origin of the far-UV (FUV) flux from the early-type galaxies, Galaxy Evolution Explorer (GALEX) is collecting the UV data for the elliptical-rich clusters at moderate redshifts (z < 0.2) where the dominant FUV source is predicted to be hot horizontal-branch (HB) stars and their post-HB progeny. Here we present our first result for the early-type galaxies in Abell 2670 at z = 0.076. Compared to NGC 1399, a nearby giant elliptical galaxy in the Fornax cluster, it appears that the rest-frame FUV - V color of the giant ellipticals gets redder by ~ 0.7 mag at the distance of Abell 2670 (z = 0.076; look-back time ~ 1.0 Gyr). Although a detailed comparison with the models is postponed until more cluster data are accumulated, it is interesting to note that this value is consistent with the variation predicted by the population synthesis models where the mean temperature of HB stars declines rapidly with increasing look-back time.

M32, the compact elliptical-galaxy companion to the Andromeda spiral galaxy has been imaged by the Galaxy Evolution Explorer (GALEX) in two ultraviolet bands, centered at ~1500 (FUV) and 2300 Angstroms (NUV). The imaging data have been carefully decomposed so as to properly account for the complicated background contamination from the disk of M31. We have derived the surface brightness and color profiles finding a slightly positive color gradient of Delta(FUV-B)/Delta log(r)=+0.15+/-0.03 within one effective radius. Earlier data from the Ultraviolet Imaging Telescope suggested that M32 had an extremely large (negative) FUV-optical color gradient (Delta(FUV-B)/Delta log(r)<-2), inverted with respect to the majority of gradients seen in giant elliptical galaxies. Our new results show that, despite of its very low UV-upturn, M32 has similar UV properties to those observed in luminous elliptical galaxies.

We give an overview of the Galaxy Evolution Explorer (GALEX), a NASA Explorer Mission launched on April 28, 2003. GALEX is performing the first space UV sky-survey, including imaging and grism surveys in two bands (1350-1750 Angstroms and 1750-2750 Angstroms). The surveys include an all-sky imaging survey (m[AB] ~ 20.5), a medium imaging survey of 1000 square degrees (m[AB] ~ 23), a deep imaging survey of 100 square degrees (m[AB] ~ 25), and a nearby galaxy survey. Spectroscopic grism surveys (R=100-200) are underway with various depths and sky coverage. Many targets overlap existing or planned surveys. We will use the measured UV properties of local galaxies, along with corollary observations, to calibrate the UV-global star formation rate relationship in local galaxies. We will apply this calibration to distant galaxies discovered in the deep imaging and spectroscopic surveys to map the history of star formation in the universe over the redshift range 0 < z < 1.5, and probe the physical drivers of star formation in galaxies. The GALEX mission includes a Guest Investigator program supporting the wide variety of programs made possible by the first UV sky survey.

We report the first year on-orbit performance results for the Galaxy Evolution Explorer (GALEX), a NASA Small Explorer that is performing a survey of the sky in two ultraviolet bands. The instrument comprises a 50 cm diameter modified Ritchey-Chretien telescope with a 1.25 degree field of view, selectable imaging and objective grism spectroscopic modes, and an innovative optical system with a thin-film multilayer dichroic beam splitter that enables simultaneous imaging by a pair of photon counting, microchannel plate, delay line readout detectors. Initial measurements demonstrate that GALEX is performing well, meeting its requirements for resolution, efficiency, astrometry, bandpass definition and survey sensitivity.

The NASA Galaxy Evolution Explorer (GALEX) satellite has obtained simultaneous near and far ultraviolet light curves of the ROTSE-I Catalog RR Lyrae ab-type variable star J143753.84+345924.8. A series of 38 GALEX Deep Imaging Survey observations well distributed in phase within the star's 0.56432d period shows an AB=4.9mag variation in the far UV (1350-1750A) band and an AB=1.8mag variation in the near UV (1750-2750A) band, compared with only a 0.8mag variation in the broad, unfiltered ROTSE-I (4500-10000A) band. These GALEX UV observations are the first to reveal a large RR Lyrae amplitude variation at wavelengths below 1800A. We compare the GALEX and ROTSE-I observations to predictions made by recent Kurucz stellar atmosphere models. We use published physical parameters for the comparable period (0.57433d), well-observed RR Lyrae star WY Antliae to compute predicted FUV, NUV, and ROTSE-I light curves for J143753.84+345924.8. The observed light curves agree with the Kurucz predictions for [Fe/H]=-1.25 to within AB=0.2mag in the GALEX NUV and ROTSE-I bands, and within 0.5mag in the FUV. At all metallicities between solar and one hundredth solar, the Kurucz models predict 6-8mag of variation at wavelengths between 1000-1700A. Other variable stars with similar temperature variations, such as Cepheids, should also have large-amplitude FUV light curves, observable during the ongoing GALEX imaging surveys.

The Deep Extragalactic Evolutionary Probe (DEEP) is a series of spectroscopic surveys of faint galaxies, targeted at the properties and clustering of galaxies at redshifts z ~ 1. We present the redshift catalog of the DEEP 1 GSS pilot phase of this project, a Keck/LRIS survey in the HST/WFPC2 Groth Survey Strip. The redshift catalog and data, including reduced spectra, are publicly available through a Web-accessible database. The catalog contains 658 secure galaxy redshifts with a median z=0.65, and shows large-scale structure walls to z = 1. We find a bimodal distribution in the galaxy color-magnitude diagram which persists to z = 1. A similar color division has been seen locally by the SDSS and to z ~ 1 by COMBO-17. For red galaxies, we find a reddening of only 0.11 mag from z ~ 0.8 to now, about half the color evolution measured by COMBO-17. We measure structural properties of the galaxies from the HST imaging, and find that the color division corresponds generally to a structural division. Most red galaxies, ~ 75%, are centrally concentrated, with a red bulge or spheroid, while blue galaxies usually have exponential profiles. However, there are two subclasses of red galaxies that are not bulge-dominated: edge-on disks and a second category which we term diffuse red galaxies (DIFRGs). The distant edge-on disks are similar in appearance and frequency to those at low redshift, but analogs of DIFRGs are rare among local red galaxies. DIFRGs have significant emission lines, indicating that they are reddened mainly by dust rather than age. The DIFRGs in our sample are all at z>0.64, suggesting that DIFRGs are more prevalent at high redshifts; they may be related to the dusty or irregular extremely red objects (EROs) beyond z>1.2 that have been found in deep K-selected surveys. (abridged)

We summarize the Sloan Digital Sky Survey data acquisition and processing steps, and describe runQA, a pipeline designed for automated data quality assessment. In particular, we show how the position of the stellar locus in color-color diagrams can be used to estimate the accuracy of photometric zeropoint calibration to better than 0.01 mag in 0.03 deg2 patches. Using this method, we estimate that typical photometric zeropoint calibration errors for SDSS imaging data are not larger than ~0.01 mag in the g, r, and i bands, 0.02 mag in the z band, and 0.03 mag in the u band (root-mean-scatter for zeropoint offsets).

We study 16,707 quasar spectra from the SDSS using the Karhunen-Lo\`eve (KL) transform (or Principal Component Analysis, PCA). The quasar eigenspectra of the full catalog reveal the following: 1st order - the mean spectrum; 2nd order - a host-galaxy component; 3rd order - the UV-optical continuum slope; 4th order - the correlations of Balmer emission lines. We find that the spectral classification of quasars is redshift and luminosity dependent, as such there does not exist a compact set (< 10 modes) of eigenspectra (covering 900 to 8000 angstrom) which can describe most variations (> 95%) of the entire catalog. We therefore construct several sets of eigenspectra in different redshift and luminosity bins. From these eigenspectra we find that quasar spectra can be classified (by the first two eigenspectra) into a sequence that is defined by a simple progression in the steepness of the slope of the continuum. We also find a dependence on redshift and luminosity in the eigencoefficients. The dominant redshift effect is a result of the evolution of the blended FeII emission (optical) and the Balmer continuum (the ``small bump''). A luminosity dependence is also present in the eigencoefficients and is related to the Baldwin effect. The eigenspectra constructed from the wavelength-selected SDSS spectra are found to agree with the principal components by Francis et al. (1992) and the well-known ``Eigenvector-1'' (Boroson and Green 1992). (Abridged)

We present measurements of the redshift--space three-point correlation function of galaxies in the Sloan Digital Sky Survey (SDSS). For the first time, we analyze the dependence of this statistic on galaxy morphology, color and luminosity. In order to control systematics due to selection effects, we used $r$--band, volume-limited samples of galaxies, constructed from the magnitude-limited SDSS data ($14.5<r<17.5$), and further divided the samples into two morphological types (early and late) or two color populations (red and blue). The three-point correlation function of SDSS galaxies follow the hierarchical relation well and the reduced three-point amplitudes in redshift--space are almost scale-independent ($Q_z=0.5\sim1.0$). In addition, their dependence on the morphology, color and luminosity is not statistically significant. Given the robust morphological, color and luminosity dependences of the two-point correlation function, this implies that galaxy biasing is complex on weakly non-linear to non-linear scales. We show that simple deterministic linear relation with the underlying mass could not explain our measurements on these scales.

We perform an objective classification of 170,000 galaxy spectra in the Sloan Digital Sky Survey (SDSS) using the Karhunen-Lo\`eve (KL) transform. With about one-sixth of the total set of galaxy spectra which will be obtained by the survey, we are able to carry out the most extensive analysis of its kind to date. The formalism proposed by Connolly and Szalay is adopted to correct for gappy regions in the spectra, and to derive eigenspectra and eigencoefficients. From this analysis, we show that this gap-correction formalism leads to a converging set of eigenspectra and KL-repaired spectra. Furthermore, KL eigenspectra of galaxies are found to be convergent not only as a function of iteration, but also as a function of the number of randomly selected galaxy spectra used in the analysis. From these data a set of ten eigenspectra of galaxy spectra are constructed, with rest-wavelength coverage 3450-8350 angstrom. The eigencoefficients describing these galaxies naturally place the spectra into several classes defined by the plane formed by the first three eigencoefficients of each spectrum. Spectral types, corresponding to different Hubble-types and galaxies with extreme emission lines, are identified for the 170,000 spectra and are shown to be complementary to existing spectral classifications. From a non-parametric classification technique, we find that the population of galaxies can be divided into three classes which correspond to early late- through to intermediate late-types galaxies. (Abridged)

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}.

Galaxy-galaxy lensing has emerged as a powerful probe of the dark matter halos of galaxies, but is subject to contamination if intrinsically aligned satellites of the lens galaxy are used as part of the source sample. We present a measurement of this intrinsic shear using 200,747 lens galaxies from the Sloan Digital Sky Survey (SDSS) spectroscopic sample and a sample of satellites selected using photometric redshifts. The mean intrinsic shear at transverse separations of 30--446$h^{-1}$ kpc is constrained to be $-0.0062<\Delta\gamma<+0.0066$ (99.9 per cent confidence, including identified systematics), which limits contamination of the galaxy-galaxy lensing signal to at most $\sim 15$ per cent on these scales. We present these limits as a function of transverse separation and lens luminosity. We furthermore investigate shear calibration biases in the SDSS, which can also affect galaxy-galaxy lensing, and conclude that the shear amplitude is calibrated to better than 18 per cent. This includes noise-induced calibration biases in the ellipticity, which are small for the sample considered here, but which can be more important if low signal-to-noise or poorly resolved source galaxies are used.

We construct eigenspectra from the DR1 quasars in the SDSS using the Karhunen-Lo\`eve (KL) transform (or Principal Component Analysis, PCA) in different redshift and luminosity bins. We find that the quasar spectra can be classified, by the first two eigenspectra, into a continuous sequence in the variation of the spectral slope. We also find a dependence on redshift and luminosity in the eigencoefficients. The dominant redshift effect is the evolution of the blended Fe~II emission (optical) and the Balmer continuum (the ``small-bump'', $\lambda_{rest} \approx 2000-4000$\AA), while the luminosity effect is related to the Baldwin effect. Correlations among several major broad emission lines are found, including the well-known ``Eigenvector-1''.

The advent of wide-area multicolour synoptic sky surveys is leading to data sets unprecedented in size, complexity and data throughput. VO technology offers a way to exploit these to the full but requires changes in design philosophy. The Palomar-QUEST survey is a major new survey being undertaken by Caltech, Yale, JPL and Indiana University to repeatedly observe 1/3 of the sky (~15000 sq. deg. between -27 < Dec <27 in seven passbands. Utilising the 48-inch Oschin Schmidt Telescope at the Palomar Observatory with the 112-CCD QUEST camera covering the full 4 x 4 sq. deg. field of view, it will generate \~1TB of data per month. In this paper, we review the design of QUEST as a VO resource, a federated data set and an exemplar of VO standards.

We present estimates of cosmological parameters from the application of the Karhunen-Loeve transform to the analysis of the 3D power spectrum of density fluctuations using Sloan Digital Sky Survey galaxy redshifts. We use Omega_m*h and f_b = Omega_b/Omega_m to describe the shape of the power spectrum, sigma8 for the (linearly extrapolated) normalization, and beta to parametrize linear theory redshift space distortions. On scales k < 0.16 h/Mpc, our maximum likelihood values are Omega_m*h = 0.264 +/-0.043, f_b = 0.286 +/- 0.065, sigma8 = 0.966 +/- 0.048, and beta = 0.45 +/- 0.12. When we take a prior on Omega_b from WMAP, we find Omega_m*h = 0.207 +/- 0.030, which is in excellent agreement with WMAP and 2dF. This indicates that we have reasonably measured the gross shape of the power spectrum but we have difficulty breaking the degeneracy between Omega_m*h and f_b because the baryon oscillations are not resolved in the current spectroscopic survey window function.

We measure cosmological parameters using the three-dimensional power spectrum P(k) from over 200,000 galaxies in the Sloan Digital Sky Survey (SDSS) in combination with WMAP and other data. Our results are consistent with a ``vanilla'' flat adiabatic Lambda-CDM model without tilt (n=1), running tilt, tensor modes or massive neutrinos. Adding SDSS information more than halves the WMAP-only error bars on some parameters, tightening 1 sigma constraints on the Hubble parameter from h~0.74+0.18-0.07 to h~0.70+0.04-0.03, on the matter density from Omega_m~0.25+/-0.10 to Omega_m~0.30+/-0.04 (1 sigma) and on neutrino masses from <11 eV to <0.6 eV (95%). SDSS helps even more when dropping prior assumptions about curvature, neutrinos, tensor modes and the equation of state. Our results are in substantial agreement with the joint analysis of WMAP and the 2dF Galaxy Redshift Survey, which is an impressive consistency check with independent redshift survey data and analysis techniques. In this paper, we place particular emphasis on clarifying the physical origin of the constraints, i.e., what we do and do not know when using different data sets and prior assumptions. For instance, dropping the assumption that space is perfectly flat, the WMAP-only constraint on the measured age of the Universe tightens from t0~16.3+2.3-1.8 Gyr to t0~14.1+1.0-0.9 Gyr by adding SDSS and SN Ia data. Including tensors, running tilt, neutrino mass and equation of state in the list of free parameters, many constraints are still quite weak, but future cosmological measurements from SDSS and other sources should allow these to be substantially tightened.

Eigenmode analysis is one of the most promising methods of analyzing large data sets in ongoing and near-future galaxy surveys. In such analyses, a fast evaluation of the correlation matrix in arbitrary cosmological models is crucial. The observational effects, including peculiar velocity distortions in redshift space, light-cone effects, selection effects, and effects of the complex shape of the survey geometry, should be taken into account in the analysis. In the framework of the linear theory of gravitational instability, we provide the methodology to quickly compute the correlation matrix. Our methods are not restricted to shallow redshift surveys, arbitrarily deep samples can be dealt with as well. Therefore, our methods are useful in constraining the geometry of the universe and the dark energy component, as well as the power spectrum of galaxies, since ongoing and near-future galaxy surveys probe the universe at intermediate to deep redshifts, z ~ 0.2--5. In addition to the detailed methods to compute the correlation matrix in 3-dimensional redshift surveys, methods to calculate the matrix in 2-dimensional projected samples are also provided. Prospects of applying our methods to likelihood estimation of the cosmological parameters are discussed.

We report the discovery of the two-image gravitationally lensed quasar SDSS J1335+0118. The object was selected as a lens candidate from the Sloan Digital Sky Survey. The imaging and spectroscopic follow-up observations confirm that the system exhibits two gravitationally lensed images of a quasar at z=1.57. The image separation is 1.56''. We also detect an extended component between the two quasar images, likely the lensing galaxy. Preliminary mass modeling predicts the differential time delay \Delta t ~ 30 h^{-1} day assuming the redshift of the lens galaxy is 0.5.

We study the recently discovered gravitational lens SDSS J1004+4112, the first quasar lensed by a cluster of galaxies. It consists of four images with a maximum separation of 14.62''. The system has been confirmed as a lensed quasar at z=1.734 on the basis of deep imaging and spectroscopic follow-up observations. We present color-magnitude relations for galaxies near the lens plus spectroscopy of three central cluster members, which unambiguously confirm that a cluster at z=0.68 is responsible for the large image separation. We find a wide range of lens models consistent with the data, but they suggest four general conclusions: (1) the brightest cluster galaxy and the center of the cluster potential well appear to be offset by several kpc; (2) the cluster mass distribution must be elongated in the North--South direction, which is consistent with the observed distribution of cluster galaxies; (3) the inference of a large tidal shear (~0.2) suggests significant substructure in the cluster; and (4) enormous uncertainty in the predicted time delays between the images means that measuring the delays would greatly improve constraints on the models. We also compute the probability of such large separation lensing in the SDSS quasar sample, on the basis of the CDM model. The lack of large separation lenses in previous surveys and the discovery of one in SDSS together imply a mass fluctuation normalization \sigma_8=1.0^{+0.4}_{-0.2} (95% CL), if cluster dark matter halos have an inner slope -1.5. Shallower profiles would require higher values of \sigma_8. Although the statistical conclusion might be somewhat dependent on the degree of the complexity of the lens potential, the discovery is consistent with the predictions of the abundance of cluster-scale halos in the CDM scenario. (Abridged)

Compact groups (CGs) of galaxies -- relatively poor groups of galaxies in which the typical separations between members is of the order of a galaxy diameter -- offer an exceptional laboratory for the study of dense galaxian environments with short (<1Gyr) dynamical time-scales. In this paper, we present an objectively defined catalog of CGs in 153 sq deg of the Sloan Digital Sky Survey Early Data Release (SDSS EDR). To identify CGs, we applied a modified version of Hickson's (1982) criteria aimed at finding the highest density CGs and thus reducing the number of chance alignments. Our catalog contains 175 CGs down to a limiting galaxy magnitude of r* = 21. The resulting catalog has a median depth of approximately z = 0.13, substantially deeper than previous CG catalogs. Since the SDSS will eventually image up to one quarter of the celestial sphere, we expect our final catalog, based upon the completed SDSS, will contain on the order of 5,000 - 10,000 CGs. This catalog will be useful for conducting studies of the general characteristics of CGs, their environments, and their component galaxies.

We present a new, and we believe arguably correct, algorithm for producing Red-Green-Blue (RBG) composites from 3-band astronomical images. Our method ensures that an object with a specified astronomical color (e.g. g-r and r-i) has a unique color in the RGB image, as opposed to the burnt-out white stars to which we are accustomed. A natural consequence of this is that we can use the same colors to code color-magnitude diagrams, providing a natural `index' to our images. We also introduce the use of an asinh stretch, which allows us to show faint objects while simultaneously preserving the structure of brighter objects in the field, such as the spiral arms of large galaxies. We believe that, in addition to their aesthetic value, our images convey far more information than do the traditional ones, and provide examples from Sloan Digital Sky Survey (SDSS) imaging, the Hubble Deep Field (HDF), and Chandra to support our claims. More examples are available at http://www.astro.princeton.edu/~rhl/PrettyPictures

We measure the large-scale real-space power spectrum P(k) using a sample of 205,443 galaxies from the Sloan Digital Sky Survey, covering 2417 square degrees with mean redshift z~0.1. We employ a matrix-based method using pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 22 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.02 h/Mpc < k < 0.3h/Mpc. We pay particular attention to modeling, quantifying and correcting for potential systematic errors, nonlinear redshift distortions and the artificial red-tilt caused by luminosity-dependent bias. Our final result is a measurement of the real-space matter power spectrum P(k) up to an unknown overall multiplicative bias factor. Our calculations suggest that this bias factor is independent of scale to better than a few percent for k<0.1h/Mpc, thereby making our results useful for precision measurements of cosmological parameters in conjunction with data from other experiments such as the WMAP satellite. As a simple characterization of the data, our measurements are well fit by a flat scale-invariant adiabatic cosmological model with h Omega_m =0.201+/- 0.017 and L* galaxy sigma_8=0.89 +/- 0.02 when fixing the baryon fraction Omega_b/Omega_m=0.17 and the Hubble parameter h=0.72; cosmological interpretation is given in a companion paper.

We analyse the bivariate distribution, in color versus absolute magnitude (u-r vs. M_r), of a low redshift sample of galaxies from the Sloan Digital Sky Survey (SDSS; 2400 deg^2, 0.004<z<0.08, -23.5<M_r<-15.5). We trace the bimodality of the distribution from luminous to faint galaxies by fitting double-Gaussians to the color functions separated in absolute magnitude bins. Color-magnitude (CM) relations are obtained for red and blue distributions (early- and late-type, predominantly field, galaxies) without using any cut in morphology. Instead, the analysis is based on the assumption of normal Gaussian distributions in color. We find that the CM relations are well fit by a straight line plus a tanh function. Both relations can be described by a shallow CM trend (slopes of about -0.04, -0.05) plus a steeper transition in the average galaxy properties over about two magnitudes. The midpoints of the transitions (M_r=-19.8 and -20.8 for the red and blue distributions, respectively) occur around 2x10^10 M_solar after converting luminosities to stellar mass. Separate luminosity functions are obtained for the two distributions. The red distribution has a more luminous characteristic magnitude and a shallower faint-end slope (M^*=-21.5, alpha=-0.8) compared to the blue distribution (alpha=-1.3 depending on the parameterization). These are approximately converted to galaxy stellar mass functions. The red distribution galaxies have a higher number density per magnitude for masses greater than about 3x10^10 M_solar. Using a simple merger model, we show that the differences between the two functions are consistent with the red distribution being formed from major galaxy mergers.

We analyze the photometric information contained in individual pixels of galaxies in the Hubble Deep Field North (HDFN) using a new technique, _pixel-z_, that combines predictions of evolutionary synthesis models with photometric redshift template fitting. Each spectral energy distribution template is a result of modeling of the detailed physical processes affecting gas properties and star formation efficiency. The criteria chosen to generate the SED templates is that of sampling a wide range of physical characteristics such as age, star formation rate, obscuration and metallicity. A key feature of our method is the sophisticated use of error analysis to generate error maps that define the reliability of the template fitting on pixel scales and allow for the separation of the interplay among dust, metallicity and star formation histories. This technique offers a number of advantages over traditional integrated color studies. As a first application, we derive the star formation and metallicity histories of galaxies in the HDFN. Our results show that the comoving density of star formation rate, determined from the UV luminosity density of sources in the HDFN, increases monotonically with redshift out to at least redshift of 5. This behavior can plausibly be explained by a smooth increase of the UV luminosity density with redshift coupled with an increase in the number of star forming regions as a function of redshift. We also find that the information contained in individual pixels in a galaxy can be linked to its morphological history. Finally, we derive the metal enrichment rate history of the universe and find it in good agreement with predictions based on the evolving HI content of Lyman-alpha QSO absorption line systems.

We report the discovery of a new gravitationally lensed quasar from the Sloan Digital Sky Survey, SDSS J090334.92+502819.2. This object was targeted for SDSS spectroscopy as a Luminous Red Galaxy (LRG), but manual examination of the spectrum showed the presence of a quasar at z= 3.6 in addition to a red galaxy at z=0.388, and the SDSS image showed a second possible quasar image nearby. Follow-up imaging and spectroscopy confirmed the lensing hypothesis. In images taken at the ARC 3.5-meter telescope, two quasars are separated by 2.8 arc-seconds; the lensing galaxy is clearly seen and is blended with one of the quasar images. Spectroscopy taken at the Keck II telescope shows that the quasars have identical redshifts of z=3.6 and both show the presence of the same broad absorption line-like troughs. We present simple lens models which account for the geometry and magnifications. The lens galaxy lies near two groups of galaxies and may be a part of them. The models suggest that the groups may contribute considerable shear and may have a strong effect on the lens configuration.

We present measurements of the angular cross-correlation between luminous red galaxies from the Sloan Digital Sky Survey and the cosmic microwave background temperature maps from the Wilkinson Microwave Anisotropy Probe. We find a statistically significant achromatic positive correlation between these two data sets, which is consistent with the expected signal from the late Integrated Sachs-Wolfe (ISW) effect. We do not detect any anti-correlation on small angular scales as would be produced from a large Sunyaev-Zel'dovich (SZ) effect, although we do see evidence for some SZ effect for our highest redshift samples. Assuming a flat universe, our preliminary detection of the ISW effect provides independent physical evidence for the existence of dark energy.

Understanding the clustering of galaxies has long been a goal of modern observational cosmology. Utilizing our photometric redshift technique a volume limited sample containing more than 2 million galaxies is constructed from the SDSS galaxy catalog. In the largest such analysis to date, we study the angular clustering as a function of luminosity and spectral type. Using Limber's equation we calculate the clustering length for the full data set as r0=5.77+/-0.10 Mpc/h. We find that r0 increases with luminosity by a factor of 1.6 over the sampled luminosity range, in agreement with previous redshift surveys. We also find that both the clustering length and the slope of the correlation function depend on the galaxy type. In particular, by splitting the galaxies in four groups by their rest-frame type we find a bimodal behavior in their clustering properties. Galaxies with spectral types similar to elliptical galaxies have a correlation length of 6.59 +/- 0.17 Mpc/h and a slope of the angular correlation function of 0.96 +/- 0.05 while blue galaxies have a clustering length of 4.51 +/- 0.19 Mpc/h and a slope of 0.68 +/- 0.09. The two intermediate color groups behave like their more extreme 'siblings', rather than showing a gradual transition in slope. We discuss these correlations in the context of current cosmological models for structure formation.

We show the apparent redshift-space clustering of galaxies in redshift range of 0.2--0.4 provides surprisingly useful constraints on dark energy component in the universe, because of the right balance between the density of objects and the survey depth. We apply Fisher matrix analysis to the the Luminous Red Galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS), as a concrete example. Possible degeneracies in the evolution of the equation of state (EOS) and the other cosmological parameters are clarified.

We combine Sloan Digital Sky Survey spectra of 22,000 luminous, red, bulge-dominated galaxies to get high S/N average spectra in the rest-frame optical and ultraviolet (2600A to 7000A). The average spectra of these massive, quiescent galaxies are early-type with weak emission lines and with absorption lines indicating an apparent excess of alpha elements over solar abundance ratios. We make average spectra of subsamples selected by luminosity, environment and redshift. The average spectra are remarkable in their similarity. What variations do exist in the average spectra as a function of luminosity and environment are found to form a nearly one-parameter family in spectrum space. We present a high signal-to-noise ratio spectrum of the variation. We measure the properties of the variation with a modified version of the Lick index system and compare to model spectra from stellar population syntheses. The variation may be a combination of age and chemical abundance differences, but the conservative conclusion is that the quality of the data considerably exceeds the current state of the models.

The Early Data Release from the Sloan Digital Sky survey provides one of the largest multicolor photometric catalogs currently available to the astronomical community. In this paper we present the first application of photometric redshifts to the $\sim 6$ million extended sources within these data (with 1.8 million sources having $r' < 21$). Utilizing a range of photometric redshift techniques, from empirical to template and hybrid techniques, we investigate the statistical and systematic uncertainties present within the redshift estimates for the EDR data. For $r'<21$ we find that the redshift estimates provide realistic redshift histograms with an rms uncertainty in the photometric redshift relation of 0.035 at $r'<18$ and rising to 0.1 at $r'<21$. We conclude by describing how these photometric redshifts and derived quantities, such as spectral type, restframe colors and absolute magnitudes, are stored within the SDSS database. We provide sample queries for searching on photometric redshifts and list the current caveats and issues that should be understood before using these photometric redshifts in statistical analyses of the SDSS galaxies.

The angular distribution of galaxies encodes a wealth of information about large scale structure. Ultimately, the Sloan Digital Sky Survey (SDSS) will record the angular positions of order 10^8 galaxies in five bands, adding significantly to the cosmological constraints. This is the first in a series of papers analyzing a rectangular stripe 2.5x90 degrees from early SDSS data. We present the angular correlation function for galaxies in four separate magnitude bins on angular scales ranging from 0.003 degrees to 15 degrees. Much of the focus of this paper is on potential systematic effects. We show that the final galaxy catalog -- with the mask accounting for regions of poor seeing, reddening, bright stars, etc. -- is free from external and internal systematic effects for galaxies brighter than r* = 22. Our estimator of the angular correlation function includes the effects of the integral constraint and the mask. The full covariance matrix of errors in these estimates is derived using mock catalogs with further estimates using a number of other methods.

A precision method for determining the value of Beta= Omega_m^{0.6}/b, where b is the galaxy bias parameter, is presented. In contrast to other existing techniques that focus on estimating this quantity by measuring distortions in the redshift space galaxy-galaxy correlation function or power spectrum, this method removes the distortions by reconstructing the real space density field and determining the value of Beta that results in a symmetric signal. To remove the distortions, the method modifies the amplitudes of a Fourier plane-wave expansion of the survey data parameterized by Beta. This technique is not dependent on the small-angle/plane-parallel approximation and can make full use of large redshift survey data. It has been tested using simulations with four different cosmologies and returns the value of Beta to +/- 0.031, over a factor of two improvement over existing techniques.

We estimate how clustering in large-scale redshift surveys can constrain various cosmological parameters. Depth and sky coverage of modern redshift surveys are greater than ever, opening new possibilities for statistical analysis. We have constructed a novel maximum likelihood technique applicable to deep redshift surveys of wide sky coverage by taking into account the effects of both curvature and linear velocity distortions. The Fisher information matrix is evaluated numerically to show the bounds derived from a given redshift sample. We find that intermediate-redshift galaxies, such as the Luminous Red Galaxies (LRGs) in the Sloan Digital Sky Survey, can constrain cosmological parameters, including the cosmological constant, unexpectedly well. The importance of the dense as well as deep sampling in designing redshift surveys is emphasized.

We compute the angular power spectrum C_l from 1.5 million galaxies in early SDSS data on large angular scales, l<600. The data set covers about 160 square degrees, with a characteristic depth of order 1 Gpc/h in the faintest (21<r<22) of our four magnitude bins. Cosmological interpretations of these results are presented in a companion paper by Dodelson et al (2001). The data in all four magnitude bins are consistent with a simple flat ``concordance'' model with nonlinear evolution and linear bias factors of order unity. Nonlinear evolution is particularly evident for the brightest galaxies. A series of tests suggest that systematic errors related to seeing, reddening, etc., are negligible, which bodes well for the sixtyfold larger sample that the SDSS is currently collecting. Uncorrelated error bars and well-behaved window functions make our measurements a convenient starting point for cosmological model fitting.

We present initial results for counts in cells statistics of the angular distribution of galaxies in early data from the Sloan Digital Sky Survey (SDSS). We analyze a rectangular stripe $2.5^\circ$ wide, covering approximately 160 sq. degrees, containing over $10^6$ galaxies in the apparent magnitude range $18 < r^\prime < 22$, with areas of bad seeing, contamination from bright stars, ghosts, and high galactic extinction masked out. This survey region, which forms part of the SDSS Early Data Release, is the same as that for which two-point angular clustering statistics have recently been computed. The third and fourth moments of the cell counts, $s_3$ (skewness) and $s_4$ (kurtosis), constitute the most accurate measurements to date of these quantities (for $r^\prime < 21$) over angular scales $0.015^\circ-0.3^\circ$. They display the approximate hierarchical scaling expected from non-linear structure formation models and are in reasonable agreement with the predictions of $\Lambda$-dominated cold dark matter models with galaxy biasing that suppresses higher order correlations at small scales. The results are in general consistent with previous measurements in the APM, EDSGC, and Deeprange surveys. These results suggest that the SDSS imaging data are free of systematics to a high degree and will therefore enable determination of the skewness and kurtosis to 1% and less then 10%, as predicted by Colombi, Szapudi, & Szalay (1998).

We present moderate resolution Keck spectroscopy of quasars at z=5.82, 5.99 and 6.28, discovered by the Sloan Digital Sky Survey (SDSS). We find that the Ly Alpha absorption in the spectra of these quasars evolves strongly with redshift. To z~5.7, the Ly Alpha absorption evolves as expected from an extrapolation from lower redshifts. However, in the highest redshift object, SDSSp J103027.10+052455.0 (z=6.28), the average transmitted flux is 0.0038+-0.0026 times that of the continuum level over 8450 A < lambda < 8710 A (5.95<z(abs)<6.16), consistent with zero flux. Thus the flux level drops by a factor of >150, and is consistent with zero flux in the Ly Alpha forest region immediately blueward of the Ly Alpha emission line, compared with a drop by a factor of ~10 at z(abs)~5.3. A similar break is seen at Ly Beta; because of the decreased oscillator strength of this transition, this allows us to put a considerably stronger limit, tau(eff) > 20, on the optical depth to Ly Alpha absorption at z=6. This is a clear detection of a complete Gunn-Peterson trough, caused by neutral hydrogen in the intergalactic medium. Even a small neutral hydrogen fraction in the intergalactic medium would result in an undetectable flux in the Ly Alpha forest region. Therefore, the existence of the Gunn-Peterson trough by itself does not indicate that the quasar is observed prior to the reionization epoch. However, the fast evolution of the mean absorption in these high-redshift quasars suggests that the mean ionizing background along the line of sight to this quasar has declined significantly from z~5 to 6, and the universe is approaching the reionization epoch at z~6.

As a scientific discipline, Astronomy is rather unique. We only have one laboratory, the Universe, and we cannot, of course, change the initial conditions and study the resulting effects. On top of this, acquiring Astronomical data has historically been a very labor-intensive effort. As a result, data has traditionally been preserved for posterity. With recent technological advances, however, the rate at which we acquire new data has grown exponentially, which has generated a Data Tsunami, whose wave train threatens to overwhelm the field. In this conference proceedings, we present and define the concept of virtual observatories, which we feel is the only logical answer to this dilemma.

We describe the target selection and resulting properties of a spectroscopic sample of luminous, red galaxies (LRG) from the imaging data of the Sloan Digital Sky Survey (SDSS). These galaxies are selected on the basis of color and magnitude to yield a sample of luminous, intrinsically red galaxies that extends fainter and further than the main flux-limited portion of the SDSS galaxy spectroscopic sample. The sample is designed to impose a passively-evolving luminosity and rest-frame color cut to a redshift of 0.38. Additional, yet more luminous, red galaxies are included to a redshift of 0.5. Approximately 12 of these galaxies per square degree are targeted for spectroscopy, so the sample will number over 100,000 with the full survey. SDSS commissioning data indicate that the algorithm efficiently selects luminous (M_g=-21.4), red galaxies, that the spectroscopic success rate is very high, and that the resulting set of galaxies is approximately volume-limited out to z=0.38. When the SDSS is complete, the LRG spectroscopic sample will fill over 1h^-3 Gpc^3 with an approximately homogeneous population of galaxies and will therefore be well suited to studies of large-scale structure and clusters out to z=0.5.

We present measurements of parameters of the 3-dimensional power spectrum of galaxy clustering from 222 square degrees of early imaging data in the Sloan Digital Sky Survey. The projected galaxy distribution on the sky is expanded over a set of Karhunen-Loeve eigenfunctions, which optimize the signal-to-noise ratio in our analysis. A maximum likelihood analysis is used to estimate parameters that set the shape and amplitude of the 3-dimensional power spectrum. Our best estimates are Gamma=0.188 +/- 0.04 and sigma_8L = 0.915 +/- 0.06 (statistical errors only), for a flat Universe with a cosmological constant. We demonstrate that our measurements contain signal from scales at or beyond the peak of the 3D power spectrum. We discuss how the results scale with systematic uncertainties, like the radial selection function. We find that the central values satisfy the analytically estimated scaling relation. We have also explored the effects of evolutionary corrections, various truncations of the KL basis, seeing, sample size and limiting magnitude. We find that the impact of most of these uncertainties stay within the 2-sigma uncertainties of our fiducial result.

The Sloan Digital Sky Survey is one of the first multicolor photometric and spectroscopic surveys designed to measure the statistical properties of galaxies within the local Universe. In this Letter we present some of the initial results on the angular 2-point correlation function measured from the early SDSS galaxy data. The form of the correlation function, over the magnitude interval 18<r*<22, is shown to be consistent with results from existing wide-field, photographic-based surveys and narrower CCD galaxy surveys. On scales between 1 arcminute and 1 degree the correlation function is well described by a power-law with an exponent of ~ -0.7. The amplitude of the correlation function, within this angular interval, decreases with fainter magnitudes in good agreement with analyses from existing galaxy surveys. There is a characteristic break in the correlation function on scales of approximately 1-2 degrees. On small scales, < 1', the SDSS correlation function does not appear to be consistent with the power-law form fitted to the 1'< theta <0.5 deg data. With a data set that is less than 2% of the full SDSS survey area, we have obtained high precision measurements of the power-law angular correlation function on angular scales 1' < theta < 1 deg, which are robust to systematic uncertainties. Because of the limited area and the highly correlated nature of the error covariance matrix, these initial results do not yet provide a definitive characterization of departures from the power-law form at smaller and larger angles. In the near future, however, the area of the SDSS imaging survey will be sufficient to allow detailed analysis of the small and large scale regimes, measurements of higher-order correlations, and studies of angular clustering as a function of redshift and galaxy type.

Early photometric data from the Sloan Digital Sky Survey (SDSS) contain angular positions for 1.5 million galaxies. In companion papers, the angular correlation function $w(\theta)$ and 2D power spectrum $C_l$ of these galaxies are presented. Here we invert Limber's equation to extract the 3D power spectrum from the angular results. We accomplish this using an estimate of $dn/dz$, the redshift distribution of galaxies in four different magnitude slices in the SDSS photometric catalog. The resulting 3D power spectrum estimates from $w(\theta)$ and $C_l$ agree with each other and with previous estimates over a range in wavenumbers $0.03 < k/{\rm h Mpc}^{-1} < 1$. The galaxies in the faintest magnitude bin ($21 < \rstar < 22$, which have median redshift $z_m=0.43$) are less clustered than the galaxies in the brightest magnitude bin ($18 < \rstar < 19$ with $z_m=0.17$), especially on scales where nonlinearities are important. The derived power spectrum agrees with that of Szalay et al. (2001) who go directly from the raw data to a parametric estimate of the power spectrum. The strongest constraints on the shape parameter $\Gamma$ come from the faintest galaxies (in the magnitude bin $21 < \rstar < 22$), from which we infer $\Gamma = 0.14^{+0.11}_{-0.06}$ (95% C.L.).

Astronomy has a long history of acquiring, systematizing, and interpreting large quantities of data. Starting from the earliest sky atlases through the first major photographic sky surveys of the 20th century, this tradition is continuing today, and at an ever increasing rate. Like many other fields, astronomy has become a very data-rich science, driven by the advances in telescope, detector, and computer technology. Numerous large digital sky surveys and archives already exist, with information content measured in multiple Terabytes, and even larger, multi-Petabyte data sets are on the horizon. Systematic observations of the sky, over a range of wavelengths, are becoming the primary source of astronomical data. Numerical simulations are also producing comparable volumes of information. Data mining promises to both make the scientific utilization of these data sets more effective and more complete, and to open completely new avenues of astronomical research. Technological problems range from the issues of database design and federation, to data mining and advanced visualization, leading to a new toolkit for astronomical research. This is similar to challenges encountered in other data-intensive fields today. These advances are now being organized through a concept of the Virtual Observatories, federations of data archives and services representing a new information infrastructure for astronomy of the 21st century. In this article, we provide an overview of some of the major datasets in astronomy, discuss different techniques used for archiving data, and conclude with a discussion of the future of massive datasets in astronomy.

We propose a novel type filter for multicolor imaging to improve on the photometric redshift estimation of galaxies. An extra filter - specific to a certain photometric system - may be utilized with high efficiency. We present a case study of the Hubble Space Telescope's Advanced Camera for Surveys and show that one extra exposure could cut down the mean square error on photometric redshifts by 34% over the z<1.3 redshift range.

From SDSS commissioning photometric and spectroscopic data, we investigate the utility of photometric redshift techniques to the task of estimating QSO redshifts. We consider empirical methods (e.g. nearest-neighbor searches and polynomial fitting), standard spectral template fitting and hybrid approaches (i.e. training spectral templates from spectroscopic and photometric observations of QSOs). We find that in all cases, due to the presence of strong emission-lines within the QSO spectra, the nearest-neighbor and template fitting methods are superior to the polynomial fitting approach. Applying a novel reconstruction technique, we can, from the SDSS multicolor photometry, reconstruct a statistical representation of the underlying SEDs of the SDSS QSOs. Although, the reconstructed templates are based on only broadband photometry the common emission lines present within the QSO spectra can be recovered in the resulting spectral energy distributions. The technique should be useful in searching for spectral differences among QSOs at a given redshift, in searching for spectral evolution of QSOs, in comparing photometric redshifts for objects beyond the SDSS spectroscopic sample with those in the well calibrated photometric redshifts for objects brighter than 20th magnitude and in searching for systematic and time variable effects in the SDSS broad band photometric and spectral photometric calibrations.

We show how the cosmological constant can be estimated from redshift surveys at different redshifts, using maximum-likelihood techniques. The apparent redshift-space clustering on large scales (\simgt 20 \himpc) are affected in the radial direction by infall, and curvature influences the apparent correlations in the transverse direction. The relative strengths of the two effects will strongly vary with redshift. Using a simple idealized survey geometry, we compute the smoothed correlation matrix of the redshift-space correlation function, and the Fisher matrix for $\omm$ and $\oml$. These represent the best possible measurement of these parameters given the geometry. We find that the likelihood contours are turning, according to the behavior of the angular-diameter distance relation. The clustering measures from redshift surveys at intermediate-to-high redshifts can provide a surprisingly tight constraint on $\oml$. We also estimate confidence contours for real survey geometries, using the SDSS LRG and QSO surveys as specific examples. We believe that this method will become a practical tool to constrain the nature of the dark energy.

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 propose and implement a fast, universally applicable method for extracting the angular power spectrum C_l from CMB temperature maps by first estimating the correlation function \xi(\theta). Our procedure recovers the C_l's using N^2 (but potentially N logN), operations, where N is the number of pixels. This is in contrast with standard maximum likelihood techniques which require N^3 operations. Our method makes no special assumptions about the map, unlike present fast techniques which rely on symmetries of the underlying noise matrix, sky coverage, scanning strategy, and geometry. This enables for the first time the analysis of megapixel maps without symmetries. The key element of our technique is the accurate multipole decomposition of \xi(\theta). The C_l error bars and cross-correlations are found by a Monte-Carlo approach. We applied our technique to a large number of simulated maps with Boomerang sky coverage in 81000 pixels. We used a diagonal noise matrix, with approximately the same amplitude as Boomerang. These studies demonstrate that our technique provides an unbiased estimator of the C_l's. Even though our method is approximate, the error bars obtained are nearly optimal, and converged only after few tens of Monte-Carlo realizations. Our method is directly applicable for the non-diagonal noise matrix. This, and other generalizations, such as minimum variance weighting schemes, polarization, and higher order statistics are also discussed.

We present an analysis of the star formation rate in galaxies between 0.7<z<1.8 using Near Infrared Camera and Multi-Object Spectrograph (NICMOS) grism spectral observations. We detect 163 galaxies in an area of ~4.4 square arcminutes, 37 of which show possible H-alpha emission. We extend the observed H-alpha luminosity function (LF) in this redshift range to luminosities a factor of two fainter than earlier work, and are consistent in the region of overlap. Using the H-alpha LF, we estimate a star formation rate (SFR) density in this redshift range of 0.166 M_sun/yr/Mpc^3 (H0=75 km/s/Mpc), consistent with other estimates based on emission lines, and supporting the order of magnitude increase in SFR density between z=0 and z=1. Our measurement of SFR density is a factor of ~2-3 greater than that estimated from UV data, comparable to the factor observed locally, implying little evolution in the relative extinctions between UV and H-alpha out to z~1.3.

Measuring the evolution in the clustering of galaxies over a large redshift range is a challenging problem. We have developed a new technique which uses photometric redshifts to measure the angular correlation function in redshift shells. This novel approach minimizes the galaxy projection effect inherent in standard angular correlation measurements, and allows for a measurement of the evolution in the galaxy correlation strength with redshift. In this paper, we present new results which utilize more accurate photometric redshifts, which are derived from a multi-band dataset (U,B,R, and I) covering almost two hundred square arcminutes to B_{AB} ~ 26.5, to quantify the evolution in the clustering of galaxies for z < 1. We also extend our technique to incorporate absolute magnitudes, which provides a simultaneous measurement of the evolution of clustering with both redshift and intrinsic luminosity. Specifically, we find a gradual decline in the strength of clustering with redshift out to z ~ 1, as predicted by semi-analytic models of structure formation. Furthermore, we find that r_0(z=0) ~ 4.0 h^{-1} Mpc for the predictions of linear theory in an Omega_0 = 0.1 universe.

We report the discovery of five quasars with redshifts of 4.67 - 5.27 and z'-band magnitudes of 19.5-20.7 M_B ~ -27. All were originally selected as distant quasar candidates in optical/near-infrared photometry from the Sloan Digital Sky Survey (SDSS), and most were confirmed as probable high-redshift quasars by supplementing the SDSS data with J and K measurements. The quasars possess strong, broad Lyman-alpha emission lines, with the characteristic sharp cutoff on the blue side produced by Lyman-alpha forest absorption. Three quasars contain strong, broad absorption features, and one of them exhibits very strong N V emission. The amount of absorption produced by the Lyman-alpha forest increases toward higher redshift, and that in the z=5.27 object (D_A ~ 0.7) is consistent with a smooth extrapolation of the absorption seen in lower redshift quasars. The high luminosity of these objects relative to most other known objects at z >~ 5 makes them potentially valuable as probes of early quasar properties and of the intervening intergalactic medium.