The classical gravitational force on a torus is anisotropic and always lower than Newton's $1/r^2$ law. We demonstrate the effects of periodicity in dark matter only $N$-body simulations of spherical collapse and standard $\Lambda$CDM initial conditions. Periodic boundary conditions cause an overall negative and anisotropic bias in cosmological simulations of cosmic structure formation. The effect is most significant when the largest mildly non-linear scales are comparable to the linear size of the simulation box, as often is the case for high-resolution hydrodynamical simulations. Spherical collapse morphs into a shape similar to an octahedron. The anisotropic growth distorts the large-scale $\Lambda$CDM dark matter structures. We introduce the direction-dependent power spectrum invariant under the octahedral group of the simulation volume and show that the results break spherical symmetry.

The standard cosmological paradigm narrates a reassuring story of a universe currently dominated by an enigmatic dark energy component. Disquietingly, its universal explaining power has recently been challenged by, above all, the $\sim4\sigma$ tension in the values of the Hubble constant. A more durable yet less studied other anomaly is the repeated observation of integrated Sachs-Wolfe imprints $\sim5\times$ stronger than expected in the $\Lambda$CDM model from $R_v > 100\,h^{-1}Mpc$ super-structures. Here we show that the inhomogeneous AvERA model of emerging curvature is capable of telling a plausible albeit radically different story that explains both observational anomalies without dark energy. We demonstrate that while stacked imprints of $R_v > 100\,h^{-1}Mpc$ supervoids in cosmic microwave background temperature maps can discriminate between the AvERA and $\Lambda$CDM models, their characteristic differences may remain hidden using alternative void definitions and stacking methodologies. Testing the extremes, we then also show that the CMB Cold Spot can plausibly be explained in the AvERA model as an ISW imprint. We hence conclude that the anomalous imprint of supervoids may well be the canary in the coal mine, and existing observational evidence for dark energy should be re-interpreted to further test alternative models.

We present our results from training and evaluating a convolutional neural network (CNN) to predict galaxy shapes from wide-field survey images of the first data release of the Dark Energy Survey (DES DR1). We use conventional shape measurements as ground truth from an overlapping, deeper survey with less sky coverage, the Canada-France Hawaii Telescope Lensing Survey (CFHTLenS). We demonstrate that CNN predictions from single band DES images reproduce the results of CFHTLenS at bright magnitudes and show higher correlation with CFHTLenS at fainter magnitudes than maximum likelihood model fitting estimates in the DES Y1 im3shape catalogue. Prediction of shape parameters with a CNN is also extremely fast, it takes only 0.2 milliseconds per galaxy, improving more than 4 orders of magnitudes over forward model fitting. The CNN can also accurately predict shapes when using multiple images of the same galaxy, even in different color bands, with no additional computational overhead. The CNN is again more precise for faint objects, and the advantage of the CNN is more pronounced for blue galaxies than red ones when compared to the DES Y1 metacalibration catalogue, which fits a single Gaussian profile using riz band images. We demonstrate that CNN shape predictions within the metacalibration self-calibrating framework yield shear estimates with negligible multiplicative bias, $ m < 10^{-3}$, and no significant PSF leakage. Our proposed setup is applicable to current and next generation weak lensing surveys where higher quality ground truth shapes can be measured in dedicated deep fields.

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.

Weak gravitational lensing is one of the most promising cosmological probes of the late universe. Several large ongoing (DES, KiDS, HSC) and planned (LSST, EUCLID, WFIRST) astronomical surveys attempt to collect even deeper and larger scale data on weak lensing. Due to gravitational collapse, the distribution of dark matter is non-Gaussian on small scales. However, observations are typically evaluated through the two-point correlation function of galaxy shear, which does not capture non-Gaussian features of the lensing maps. Previous studies attempted to extract non-Gaussian information from weak lensing observations through several higher-order statistics such as the three-point correlation function, peak counts or Minkowski-functionals. Deep convolutional neural networks (CNN) emerged in the field of computer vision with tremendous success, and they offer a new and very promising framework to extract information from 2 or 3-dimensional astronomical data sets, confirmed by recent studies on weak lensing. We show that a CNN is able to yield significantly stricter constraints of ($\sigma_8, \Omega_m$) cosmological parameters than the power spectrum using convergence maps generated by full N-body simulations and ray-tracing, at angular scales and shape noise levels relevant for future observations. In a scenario mimicking LSST or Euclid, the CNN yields 2.4-2.8 times smaller credible contours than the power spectrum, and 3.5-4.2 times smaller at noise levels corresponding to a deep space survey such as WFIRST. We also show that at shape noise levels achievable in future space surveys the CNN yields 1.4-2.1 times smaller contours than peak counts, a higher-order statistic capable of extracting non-Gaussian information from weak lensing maps.

Dark matter cannot be observed directly, but its weak gravitational lensing slightly distorts the apparent shapes of background galaxies, making weak lensing one of the most promising probes of cosmology. Several observational studies have measured the effect, and there are currently running, and planned efforts to provide even larger, and higher resolution weak lensing maps. Due to nonlinearities on small scales, the traditional analysis with two-point statistics does not fully capture all the underlying information. Multiple inference methods were proposed to extract more details based on higher order statistics, peak statistics, Minkowski functionals and recently convolutional neural networks (CNN). Here we present an improved convolutional neural network that gives significantly better estimates of $\Omega_m$ and $\sigma_8$ cosmological parameters from simulated convergence maps than the state of art methods and also is free of systematic bias. We show that the network exploits information in the gradients around peaks, and with this insight, we construct a new, easy-to-understand, and robust peak counting algorithm based on the 'steepness' of peaks, instead of their heights. The proposed scheme is even more accurate than the neural network on high-resolution noiseless maps. With shape noise and lower resolution its relative advantage deteriorates, but it remains more accurate than peak counting.

We present a novel $N$-body simulation method that compactifies the infinite spatial extent of the Universe into a finite sphere with isotropic boundary conditions to follow the evolution of the large-scale structure. Our approach eliminates the need for periodic boundary conditions, a mere numerical convenience which is not supported by observation and which modifies the law of force on large scales in an unrealistic fashion. We demonstrate that our method outclasses standard simulations executed on workstation-scale hardware in dynamic range, it is balanced in following a comparable number of high and low $k$ modes and, its fundamental geometry and topology match observations. Our approach is also capable of simulating an expanding, infinite universe in static coordinates with Newtonian dynamics. The price of these achievements is that most of the simulated volume has smoothly varying mass and spatial resolution, an approximation that carries different systematics than periodic simulations. Our initial implementation of the method is called StePS which stands for Stereographically Projected Cosmological Simulations. It uses stereographic projection for space compactification and naive $\mathcal{O}(N^2)$ force calculation which is nevertheless faster to arrive at a correlation function of the same quality than any standard (tree or P$^3$M) algorithm with similar spatial and mass resolution. The $N^2$ force calculation is easy to adapt to modern graphics cards, hence our code can function as a high-speed prediction tool for modern large-scale surveys. To learn about the limits of the respective methods, we compare StePS with GADGET-2 \citep{Gadget2_2005MNRAS.364.1105S} running matching initial conditions.

The recent AvERA cosmological simulation of R\'acz et al. (2017) has a $\Lambda \mathrm{CDM}$-like expansion history and removes the tension between local and Planck (cosmic microwave background) Hubble constants. We contrast the AvERA prediction of the integrated Sachs--Wolfe (ISW) effect with that of $\Lambda \mathrm{CDM}$. The linear ISW effect is proportional to the derivative of the growth function, thus it is sensitive to small differences in the expansion histories of the respective models. We create simulated ISW maps tracing the path of light-rays through the Millennium XXL cosmological simulation, and perform theoretical calculations of the ISW power spectrum. AvERA predicts a significantly higher ISW effect than $\Lambda \mathrm{CDM}$, $A=1.93-5.29$ times larger depending on the $l$ index of the spherical power spectrum, which could be utilized to definitively differentiate the models. We also show that AvERA predicts an opposite-sign ISW effect in the redshift range $z \approx 1.5 - 4.4$, in clear contrast with $\Lambda \mathrm{CDM}$. Finally, we compare our ISW predictions with previous observations. While at present these cannot distinguish between the two models due to large error bars, and lack of internal consistency suggesting systematics, ISW probes from future surveys will tightly constrain the models.

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.

According to the separate universe conjecture, spherically symmetric sub-regions in an isotropic universe behave like mini-universes with their own cosmological parameters. This is an excellent approximation in both Newtonian and general relativistic theories. We estimate local expansion rates for a large number of such regions, and use a scale parameter calculated from the volume-averaged increments of local scale parameters at each time step in an otherwise standard cosmological $N$-body simulation. The particle mass, corresponding to a coarse graining scale, is an adjustable parameter. This mean field approximation neglects tidal forces and boundary effects, but it is the first step towards a non-perturbative statistical estimation of the effect of non-linear evolution of structure on the expansion rate. Using our algorithm, a simulation with an initial $\Omega_m=1$ Einstein--de~Sitter setting closely tracks the expansion and structure growth history of the $\Lambda$CDM cosmology. Due to small but characteristic differences, our model can be distinguished from the $\Lambda$CDM model by future precision observations. Moreover, our model can resolve the emerging tension between local Hubble constant measurements and the Planck best-fitting cosmology. Further improvements to the simulation are necessary to investigate light propagation and confirm full consistency with cosmic microwave background observations.

The Fermi collaboration identified a possible electromagnetic counterpart of the gravitational wave event of September 14, 2015. Our goal is to provide an unsupervised data analysis algorithm to identify similar events in Fermi's Gamma-ray Burst Monitor CTTE data stream. We are looking for signals that are typically weak. Therefore, they can only be found by a careful analysis of count rates of all detectors and energy channels simultaneously. Our Automatized Detector Weight Optimization (ADWO) method consists of a search for the signal, and a test of its significance. We developed ADWO, a virtual detector analysis tool for multi-channel multi-detector signals, and performed successful searches for short transients in the data-streams. We have identified GRB150522B, as well as possible electromagnetic candidates of the transients GW150914 and LVT151012. ADWO is an independently developed, unsupervised data analysis tool that only relies on the raw data of the Fermi satellite. It can therefore provide a strong, independent test to any electromagnetic signal accompanying future gravitational wave observations.

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

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.

Position angle measurements of Sloan Digital Sky Survey (SDSS) galaxies, as measured by the surface brightness profile fitting code of the SDSS photometric pipeline (Lupton 2001), are known to be strongly biased, especially in the case of almost face-on and highly inclined galaxies. To address this issue we developed a reliable algorithm which determines position angles by means of isophote fitting. In this paper we present our algorithm and a catalogue of position angles for 26397 SDSS galaxies taken from the deep co-added Stripe 82 (equatorial stripe) images.

Getting spectra at good signal-to-noise ratios takes orders of magnitudes more time than photometric observations. Building on the technique developed for photometric redshift estimation of galaxies, we develop and demonstrate a non-parametric photometric method for estimating the chemical composition of galactic stars. We investigate the efficiency of our method using spectroscopically determined stellar metallicities from SDSS DR7. The technique is generic in the sense that it is not restricted to certain stellar types or stellar parameter ranges and makes it possible to obtain metallicities and error estimates for a much larger sample than spectroscopic surveys would allow. We find that our method performs well, especially for brighter stars and higher metallicities and, in contrast to many other techniques, we are able to reliably estimate the error of the predicted metallicities.

More than half of the sources identified by recent radio sky surveys have not been detected by wide-field optical surveys. We present a study based on our co-added image stacking technique, in which our aim is to detect the optical emission from unresolved, isolated radio sources of the Very Large Array (VLA) Faint Images of the Radio Sky at Twenty-cm (FIRST) survey that have no identified optical counterparts in the Sloan Digital Sky Survey (SDSS) Stripe 82 co-added data set. From the FIRST catalogue, 2116 such radio point sources were selected, and cut-out images, centred on the FIRST coordinates, were generated from the Stripe 82 images. The already co-added cut-outs were stacked once again to obtain images of high signal-to-noise ratio, in the hope that optical emission from the radio sources would become detectable. Multiple stacks were generated, based on the radio luminosity of the point sources. The resulting stacked images show central peaks similar to point sources. The peaks have very red colours with steep optical spectral energy distributions. We have found that the optical spectral index alpha_nu falls in the range -2.9 < alpha_nu < -2.2, depending only weakly on the radio flux. The total integration times of the stacks are between 270 and 300 h, and the corresponding 5 sigma detection limit is estimated to be about m_r = 26.6 mag. We argue that the detected light is mainly from the central regions of dust-reddened Type 1 active galactic nuclei. Dust-reddened quasars might represent an early phase of quasar evolution, and thus they can also give us an insight into the formation of massive galaxies. The data used in the paper are available on-line at http://www.vo.elte.hu/doublestacking.

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.

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.

We apply order statistics (OS) to the bright end ($M_r < -22$) of the luminosity distribution of early-type galaxies spectroscopically identified in the SDSS DR7 catalog. We calculate the typical OS quantities of this distribution numerically, measuring the expectation value and variance of the $k^{th}$ most luminous galaxy in a sample with cardinality $N$ over a large ensemble of such samples. From these statistical quantities we explain why and in what limit the $k^{th}$ most luminous galaxies can be used as standard candles for cosmological studies. Since our sample contains all bright galaxies including the brightest cluster galaxies (BCG), based on OS we argue that BCGs can be considered as statistical extremes of a well-established Schechter luminosity distribution when galaxies are binned by redshift and not cluster-by-cluster. We presume that the reason behind this might be that luminous red ellipticals in galaxy clusters are \em not random \em samples of an overall luminosity distribution but biased by the fact that they are in a cluster containing the BCG. We show that a simple statistical toy model can reproduce the well-known magnitude gap between the BCG and the second brightest galaxy of the clusters.

We determine the evolution of the co-moving density of the most massive ($M_* \geq 10^{12} M_\odot$) early-type galaxy population in the redshift range of $z = 0.15$ - 0.45 in different stellar mass ranges using data from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) catalog. We find that the co-moving number density of these galaxies grew exponentially, weakly depending on the stellar mass range, as a function of cosmic time with a time-scale of $\tau \simeq 1.16 \pm 0.16$ Gyr for at least 4 Gyr ending around $z \simeq 0.15$. This is about a factor of ten of growth between $z=0.5$ - 0.15. Since $z \simeq 0.15$ a constant co-moving number density can be measured. According to theoretical models the most massive early-type galaxies gain most of their stellar mass via dry merging but the major merger rate measured by others cannot account for the high growth in number density we measured thus, stellar mass gain from minor mergers and slow, smooth accretion seems to play an important role. We outline a simple analytic model that explains the observed evolution based on the exponential decline of the luminosity function and sets constraints on the time dependence of the close-pair fraction of merger candidate galaxies.

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.

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.

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.

Wide-angle surveys have been an engine for new discoveries throughout the modern history of astronomy, and have been among the most highly cited and scientifically productive observing facilities in recent years. This trend is likely to continue over the next decade, as many of the most important questions in astrophysics are best tackled with massive surveys, often in synergy with each other and in tandem with the more traditional observatories. We argue that these surveys are most productive and have the greatest impact when the data from the surveys are made public in a timely manner. The rise of the "survey astronomer" is a substantial change in the demographics of our field; one of the most important challenges of the next decade is to find ways to recognize the intellectual contributions of those who work on the infrastructure of surveys (hardware, software, survey planning and operations, and databases/data distribution), and to make career paths to allow them to thrive.

The low energy spectra of some gamma-ray bursts' show excess components beside the power-law dependence. The consequences of such a feature allows to estimate the gamma photometric redshift of the long gamma-ray bursts in the BATSE Catalog. There is good correlation between the measured optical and the estimated gamma photometric redshifts. The estimated redshift values for the long bright gamma-ray bursts are up to z=4, while for the the faint long bursts - which should be up to z=20 - the redshifts cannot be determined unambiguously with this method. The redshift distribution of all the gamma-ray bursts with known optical redshift agrees quite well with the BATSE based gamma photometric redshift distribution.

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 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 measure the mean halo-mass of z=0.5 MgII absorbers using the cross-correlation (over co-moving scales 0.05-13h^{-1}Mpc) between 1806 MgII quasar absorption systems and ~250,000 Luminous Red Galaxies (LRGs), both selected from the SDSS DR3. The MgII systems have rest-frame equivalent widths W_r(2796)>=0.3A. From the ratio of the MgII-LRG cross-correlation to the LRG-LRG auto-correlation, we find that the bias ratio between MgII absorbers and LRGs is 0.65+/-0.08, which implies that the absorber host-galaxies have a mean halo-mass 20-40 times smaller than that of the LRGs; the MgII absorbers have halos of mean mass <logM_h(M_sol)>=11.94+/-0.31(stat)+/-0.25(sys). We demonstrate that this statistical technique, which does not require any spectroscopic follow-up, does not suffer from contaminants such as stars or foreground and background galaxies. Finally, we find that the absorber halo-mass is anti-correlated with the equivalent width. If MgII absorbers were virialized in galaxy halos a positive M_h-W_r correlation would have been observed since W_r(2796) is a direct measure of the velocity spread of the MgII sub-components. Thus, our results demonstrate that the individual clouds of a MgII system are not virialized in the gaseous halos of the host-galaxies. We review past results in the literature on the statistics of MgII absorbers and find that they too require an M_h-W_r anti-correlation. When combined with measurements of the equivalent width distribution, the M_h-W_r anti-correlation naturally explains why absorbers with W_r(2796)>=2A are not seen at large impact parameters. We interpret the M_h-W_r anti-correlation within the starburst scenario where strong MgII absorbers are produced by supernovae-driven winds.

We synthesize the rest-frame Stroemgren colors using SDSS spectra for 99,088 galaxies selected from Data Release 1. This narrow-band ~200 AA photometric system (uz, vz, bz, yz), first designed for the determination of effective temperature, metallicity and gravity of stars, measures the continuum spectral slope of galaxies in the rest-frame 3200-5800 AA wavelength range. Galaxies form a remarkably narrow locus (~0.03 mag) in the resulting color-color diagram. The Bruzual & Charlot population synthesis models suggest that the position of a galaxy along this locus is controlled by a degenerate combination of metallicity and age of the dominant stellar population. Galaxy distribution along the locus is bimodal, with the local minimum corresponding to an ~1 Gyr old single stellar population. The position perpendicular to the locus is independent of metallicity and age, and reflects the galaxy's dust content, as implied by both the models and the statistics of IRAS detections. A comparison of this locus with the galaxy locus in the H_delta-D_n(4000) diagram, utilized by Kauffmann et al. (2003) to estimate stellar masses, reveals a tight correlation, although the two analyzed spectral ranges barely overlap. Overall, the galaxy spectral energy distribution in the entire UV to near-IR range can be described as a single-parameter family with an accuracy of 0.1 mag, or better. This nearly one-dimensional distribution of galaxies in the multi-dimensional space of measured parameters strongly supports the conclusion of Yip et al. (2004), based on a principal component analysis, that SDSS galaxy spectra can be described by a small number of eigenspectra. Apparently, the contributions of stellar populations that dominate the optical emission from galaxies are combined in a simple and well-defined way.

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 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 study the cross-correlation between 716 MgII quasar absorption systems and about 100,000 Luminous Red Galaxies (LRGs) selected from the Sloan Digital Sky Survey Data Release 3 in the redshift range 0.4<z<0.8. The MgII systems were selected to have 2796 & 2803 rest-frame equivalent widths greater than 1.0 \AA and identifications confirmed by the FeII 2600 or MgI 2852 lines. Over co-moving scales 0.2--13/h Mpc, the MgII--LRG cross-correlation has an amplitude 0.69+/-0.09 times that of the LRG--LRG auto-correlation. Since LRGs have halo-masses of 10^{13} \msun, this strong cross-correlation implies that the absorber host-galaxies have halo-masses 1--2 times 10^{12} \msun.

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.

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 construct and analyze a u-band selected galaxy sample from the SDSS Southern Survey, which covers 275 sq. deg. The sample includes 43223 galaxies with spectroscopic redshifts in the range 0.005<z<0.3 and with 14.5<u<20.5. The S/N in the u-band Petrosian aperture is improved by coadding multiple epochs of imaging data and by including sky-subtraction corrections. Luminosity functions for the near-UV ^{0.1}u band (322+-26 nm) are determined in redshift slices of width 0.02, which show a highly significant evolution in M^{star} of -0.8+-0.1 mag between z=0 and 0.3; with M^{star} = -18.84+-0.05 (AB mag), log phi^{star} = -2.06+-0.03 (Mpc^{-3}) and log rho_L = 19.11+-0.02 (W Hz^{-1} Mpc^{-3}) at z=0.1. The faint-end slope determined for z<0.06 is given by alpha=-1.05+-0.08. This is in agreement with recent determinations from GALEX at shorter wavelengths. Comparing our z<0.3 luminosity density measurements with 0.2<z<1.2 from COMBO-17, we find that the 280-nm density evolves as rho_L proportional to (1+z)^{beta} with beta=2.1+-0.2; and find no evidence for any change in slope over this redshift range. By comparing with other measurements of cosmic star formation history, we estimate that the effective dust attenuation at 280 nm has increased by 0.8+-0.3 mag between z=0 and 1.

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)

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 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 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 investigate the continuum and emission line properties of 4576 SDSS quasars as a function of their optical/UV SEDs. The optical/UV color distribution of our sample is roughly Gaussian, but with a red tail; we distinguish between 1) intrinsically blue (optically flat) quasars, 2) intrinsically red (optically steep) quasars, and 3) the 273 (6%) of our quasars whose continua are inconsistent with a single power-law and appear redder due to SMC-like dust reddening rather than synchrotron emission. The color distribution suggests that the population of moderately dust reddened broad-line quasars is smaller than that of unobscured quasars, but we estimate that a further 10% of the luminous quasar population is missing from the SDSS sample because of dust extinction with E(B-V)<0.5. We also investigate the emission and absorption line properties of these quasars as a function of color with regard to Boroson & Green type eigenvectors. Intrinsically red (optically steep) quasars tend to have narrower Balmer lines and weaker CIV, CIII], HeII and 3000A bump emission as compared with bluer (optically flatter) quasars. The change in strength of the 3000A bump appears to be dominated by the Balmer continuum and not by FeII emission. The dust reddened quasars have even narrower Balmer lines and weaker 3000A bumps, in addition to having considerably larger equivalent widths of [OII] and [OIII] emission. The fraction of broad absorption line quasars (BALQSOs) increases from ~3.4% for the bluest quasars to perhaps as large as 20% for the dust reddened quasars, but the intrinsic color distribution is affected by dust reddening. (abridged)

abridged: We use a complete sample of about 140,000 galaxies from the Sloan Digital Sky Survey (SDSS) to study the size distribution of galaxies and its dependence on their luminosity, stellar mass, and morphological type. The large SDSS database provides statistics of unprecedented accuracy. For each type of galaxy, the size distribution at given luminosity (or stellar mass) is well described by a log-normal function, characterized by its median $\bar{R}$ and dispersion $\sigma_{\ln R}$. For late-type galaxies, there is a characteristic luminosity at $M_{r,0}\sim -20.5$ (assuming $h=0.7$) corresponding to a stellar mass $M_0\sim 10^{10.6}\Msun$. Galaxies more massive than $M_0$ have $\bar{R}\propto M^{0.4}$ and $\sigma_{\ln R}\sim 0.3$, while less massive galaxies have $\bar{R}\propto M^{0.15}$ and $\sigma_{\ln R}\sim 0.5$. For early-type galaxies, the $\bar{R}$ - $M$ relation is significantly steeper, $\bar{R}\propto M^{0.55}$, but the $\sigma_{\ln R}$ - $M$ relation is similar to that of late-type galaxies. Faint red galaxies have sizes quite independent of their luminosities.

Using the Sloan Digital Sky Survey (SDSS), we probe the halo mass distribution by studying the velocities of satellites orbiting isolated galaxies. In a subsample that covers 2500 sq. degrees on the sky, we detect about 3000 satellites with absolute blue magnitudes going down to M_B = -14; most of the satellites have M_B=-16 to -18, comparable to the magnitudes of M32 and the Magellanic Clouds. After a careful, model-independent removal of interlopers, we find that the line-of-sight velocity dispersion of satellites declines with distance to the primary. For an L* galaxy the r.m.s. line-of-sight velocity changes from ~120 km/s at 20 kpc to ~60 km/s at 350 kpc. This decline agrees remarkably well with theoretical expectations, as all modern cosmological models predict that the density of dark matter in the peripheral parts of galaxies declines as rho_DM propto r^{-3}. Thus, for the first time we find direct observational evidence of the density decline predicted by cosmological models; we also note that this result contradicts alternative theories of gravity such as MOND. We also find that the velocity dispersion of satellites within 100 kpc scales with the absolute magnitude of the central galaxy as sigma propto L^{0.3}; this is very close to the Tully-Fisher relation for normal spiral galaxies.

We present a determination of the `Cosmic Optical Spectrum' of the Universe, i.e. the ensemble emission from galaxies, as determined from the red-selected Sloan Digital Sky Survey main galaxy sample and compare with previous results of the blue-selected 2dF Galaxy Redshift Survey. Broadly we find good agreement in both the spectrum and the derived star-formation histories. If we use a power-law star-formation history model where star-formation rate $\propto (1+z)^\beta$ out to z=1, then we find that $\beta$ of 2 to 3 is still the most likely model and there is no evidence for current surveys missing large amounts of star formation at high redshift. In particular `Fossil Cosmology' of the local universe gives measures of star-formation history which are consistent with direct observations at high redshift. Using the photometry of SDSS we are able to derive the cosmic spectrum in absolute units (i.e.$ W \AA$^{-1}$ Mpc$^{-3}$) at 2--5\AA resolution and find good agreement with published broad-band luminosity densities. For a Salpeter IMF the best fit stellar mass/light ratio is 3.7--7.5 $\Msun/\Lsun$ in the r-band (corresponding to $\omstars h = 0.0025$--0.0055) and from both the stellar emission history and the H$\alpha$ luminosity density independently we find a cosmological star-formation rate of 0.03--0.04 h $\Msun$ yr$^{-1}$ Mpc$^{-3}$ today.

We study the mean environments of galaxies in the Sloan Digital Sky Survey as a function of rest-frame luminosity and color. Overdensities in galaxy number are estimated in $8 h^{-1} \mathrm{Mpc}$ and $1 h^{-1} \mathrm{Mpc}$ spheres centered on $125,000$ galaxies taken from the SDSS spectroscopic sample. We find that, at constant color, overdensity is independent of luminosity for galaxies with the blue colors of spirals. This suggests that, at fixed star-formation history, spiral-galaxy mass is a very weak function of environment. Overdensity does depend on luminosity for galaxies with the red colors of early types; both low-luminosity and high-luminosity red galaxies are found to be in highly overdense regions.

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.

It is known that the soft tail of the gamma-ray bursts' spectra show excesses from the exact power-law dependence. In this article we show that this departure can be detected in the peak flux ratios of different BATSE DISCSC energy channels. This effect allows to estimate the redshift of the bright long gamma-ray bursts in the BATSE Catalog. A verification of these redshifts is obtained for the 8 GRB which have both BATSE DISCSC data and measured optical spectroscopic redshifts. There is good correlation between the measured and esti redshifts, and the average error is $\Delta z \approx 0.33$. The method is similar to the photometric redshift estimation of galaxies in the optical range, hence it can be called as "gamma photometric redshift estimation". The estimated redshifts for the long bright gamma-ray bursts are up to $z \simeq 4$. For the the faint long bursts - which should be up to $z \simeq 20$ - the redshifts cannot be determined unambiguously with this method.

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.

We develop a new method to constrain the star formation histories, dust attenuation and stellar masses of galaxies. It is based on two stellar absorption line indices, the 4000 Angstrom break strength and the Balmer absorption line index Hdelta_A. Together, these indices allow us to constrain the mean stellar ages of galaxies and the fractional stellar mass formed in bursts over the past few Gyr. A comparison with broad band photometry then yields estimates of dust attenuation and of stellar mass. We generate a large library of Monte Carlo realizations of different star formation histories, including starbursts of varying strength and a range of metallicities. We use this library to generate median likelihood estimates of burst mass fractions, dust attenuation strengths, stellar masses and stellar mass-to-light ratios for a sample of 122,208 galaxies drawn from the Sloan Digital Sky Survey. The typical 95% confidence range in our estimated stellar masses is +-40%. We study how the stellar mass-to-light ratios of galaxies vary as a function of absolute magnitude, concentration index and photometric pass-band and how dust attenuation varies as a function of absolute magnitude and 4000 Angstrom break strength. We also calculate how the total stellar mass of the present Universe is distributed over galaxies as a function of their mass, size, concentration, colour, burst mass fraction and surface mass density. We find that most of the stellar mass in the local Universe resides in galaxies that have stellar masses \~5\times 10^10 M_sol, half light radii ~3 kpc, and half- light surface mass densities ~10^9 M_sol/kpc^2. The distribution of D(4000) is strongly bimodal, showing a clear division between galaxies dominated by old stellar populations and galaxies with more recent star formation.

Using a catalog of 147,986 galaxy redshifts and fluxes from the Sloan Digital Sky Survey (SDSS) we measure the galaxy luminosity density at z=0.1 in five optical bandpasses corresponding to the SDSS bandpasses shifted to match their restframe shape at z=0.1. We denote the bands {0.1}{u}, {0.1}{g}, {0.1}{r}, {0.1}{i}, {0.1}{z}, with \lambda_{eff} = [3216, 4240, 5595, 6792, 8111] Angstroms respectively. We use a maximum likelihood method which allows for a general form for the shape of the luminosity function, simple luminosity and number evolution, incorporates flux uncertainties, and accounts for the flux limits of the survey. We find luminosity densities at z=0.1 in absolute AB magnitudes in a Mpc^3 of [-14.10 \pm 0.15, -15.18 \pm 0.03, -15.90 \pm 0.03, -16.24 \pm 0.03, -16.56 \pm 0.02] in [{0.1}{u}, {0.1}{g}, {0.1}{r}, {0.1}{i}, {0.1}{z}], respectively, using \Omega_0 =0.3, \Omega_\Lambda=0.7, and h=1, and using Petrosian magnitudes. Similar results are obtained using Sersic model magnitudes, suggesting that flux from outside the Petrosian apertures is not a major correction. In the {0.1}{r} band, the best fit Schechter function to our results has \phi_\ast = (1.49 \pm 0.04) \times 10^{-2} h^3 Mpc^{-3}, M_\ast - 5\log_{10} h = -20.44 \pm 0.01, and \alpha = -1.05\pm 0.01. In solar luminosities, the luminosity density in {0.1}{r} is (1.84 \pm 0.04) h 10^8 L_{{0.1}{r},\odot} Mpc^{-3}. Our results are consistent with other estimates of the luminosity density, from the Two-degree Field Galaxy Redshift Survey and the Millenium Galaxy Catalog. They represent a substantial change (\sim 0.5 mag) from earlier SDSS luminosity density results based on commissioning data, almost entirely because of the inclusion of evolution in the luminosity function model.

Broad-band measurements of flux for galaxies at different redshifts measure different regions of the rest-frame galaxy spectrum. Certain astronomical questions, such as the evolution of the luminosity function of galaxies, require transforming these magnitudes into redshift-independent quantities. To prepare to address these astronomical questions, investigated in detail in subsequent papers, we fit spectral energy distributions (SEDs) to broad band photometric observations, in the context of the optical observations of the Sloan Digital Sky Survey (SDSS). Linear combinations of four spectral templates can reproduce the five SDSS magnitudes of all galaxies to the precision of the photometry. Expressed in the appropriate coordinate system, the locus of the coefficients multiplying the templates is planar, and in fact nearly linear. The resulting reconstructed SEDs can be used to recover fixed frame magnitudes over a range of redshifts. This process yields consistent results, in the sense that within each sample the intrinsic colors of similar type galaxies are nearly constant with redshift. We compare our results to simpler interpolation methods and galaxy spectrophotometry from the SDSS. The software that generates these results is publicly available and easily adapted to handle a wide range of galaxy observations.

Using photometry and spectroscopy of 144,609 galaxies from the Sloan Digital Sky Survey, we present bivariate distributions of pairs of seven galaxy properties: four optical colors, surface brightness, radial profile shape as measured by the Sersic index, and absolute magnitude. In addition, we present the dependence of local galaxy density (smoothed on 8 h^{-1} Mpc scales) on all of these properties. Several classic, well-known relations among galaxy properties are evident at extremely high signal-to-noise ratio: the color-color relations of galaxies, the color-magnitude relations, the magnitude-surface brightness relation, and the dependence of density on color and absolute magnitude. We show that most of the i-band luminosity density in the universe is in the absolute magnitude and surface brightness ranges used. Some of the relationships between parameters, in particular the color--magnitude relations, show stronger correlations for exponential galaxies and concentrated galaxies taken separately than for all galaxies taken together. We provide a simple set of fits of the dependence of galaxy properties on luminosity for these two sets of galaxies.

Galaxy masses can be estimated by a variety of methods; each applicable in different circumstances, and each suffering from different systematic uncertainties. Confirmation of results obtained by one technique with analysis by another is particularly important. Recent SDSS weak lensing measurements of the projected-mass correlation function reveal a linear relation between galaxy luminosities and the depth of their dark matter halos (measured on 260 \hinv kpc scales). In this work we use an entirely independent dynamical method to confirm these results. We begin by assembling a sample of 618 relatively isolated host galaxies, surrounded by a total of 1225 substantially fainter satellites. We observe the mean dynamical effect of these hosts on the motions of their satellites by assembling velocity difference histograms. Dividing the sample by host properties, we find significant variations in satellite velocity dispersion with host luminosity. We quantify these variations using a simple dynamical model, measuring \mtsd a dynamical mass within 260 \hinv kpc. The appropriateness of this mass reconstruction is checked by conducting a similar analysis within an N-body simulation. Comparison between the dynamical and lensing mass-to-light scalings shows reasonable agreement, providing some quantitative confirmation for the lensing results.

We present multi-instrument optical observations of the High Energy Transient Explorer (HETE-2)/Interplanetary Network (IPN) error box of GRB 010921. This event was the first gamma ray burst (GRB) localized by HETE-2 which has resulted in the detection of an optical afterglow. In this paper we report the earliest known observations of the GRB010921 field, taken with the 0.11-m Livermore Optical Transient Imaging System (LOTIS) telescope, and the earliest known detection of the GRB010921 optical afterglow, using the 0.5-m Sloan Digital Sky Survey Photometric Telescope (SDSS PT). Observations with the LOTIS telescope began during a routine sky patrol 52 minutes after the burst. Observations were made with the SDSS PT, the 0.6-m Super-LOTIS telescope, and the 1.34-m Tautenburg Schmidt telescope at 21.3, 21.8, and 37.5 hours after the GRB, respectively. In addition, the host galaxy was observed with the USNOFS 1.0-m telescope 56 days after the burst. We find that at later times (t > 1 day after the burst), the optical afterglow exhibited a power-law decline with a slope of $\alpha = 1.75 \pm 0.28$. However, our earliest observations show that this power-law decline can not have extended to early times (t < 0.035 day).

We identify new structures in the halo of the Milky Way Galaxy from positions, colors and magnitudes of five million stars detected in the Sloan Digital Sky Survey. Most of these stars are within 1.26 degrees of the celestial equator. We present color-magnitude diagrams (CMDs) for stars in two previously discovered, tidally disrupted structures. The CMDs and turnoff colors are consistent with those of the Sagittarius dwarf galaxy, as had been predicted. In one direction, we are even able to detect a clump of red stars, similar to that of the Sagittarius dwarf, from stars spread across 110 square degrees of sky. Focusing on stars with the colors of F turnoff objects, we identify at least five additional overdensities of stars. Four of these may be pieces of the same halo structure, which would cover a region of the sky at least 40 degrees in diameter, at a distance of 11 kpc from the Sun (18 kpc from the center of the Galaxy). The turnoff is significantly bluer than that of thick disk stars, and closer to the Galactic plane than a power-law spheroid. We suggest two models to explain this new structure. One possibility is that this new structure could be a new dwarf satellite of the Milky Way, hidden in the Galactic plane, and in the process of being tidally disrupted. The other possibility is that it could be part of a disk-like distribution of stars which is metal-poor, with a scale height of approximately 2 kpc and a scale length of approximately 10 kpc. The fifth overdensity, which is 20 kpc away, is some distance from the Sagittarius dwarf streamer orbit and is not associated with any known structure in the Galactic plane. It is likely that there are many smaller streams of stars in the Galactic halo.

We identify close pairs of galaxies from 278 deg^2 of Sloan Digital Sky Survey commissioning imaging data. The pairs are drawn from a sample of 330,041 galaxies with 18 < r^* < 20. We determine the angular correlation function of galaxy pairs, and find it to be stronger than the correlation function of single galaxies by a factor of 2.9 +/- 0.4. The two correlation functions have the same logarithmic slope of 0.77. We invert Limber's equation to estimate the three-dimensional correlation functions; we find clustering lengths of r_0= 4.2 +/- 0.4 h^{-1} Mpc for galaxies and 7.8 +/- 0.7 h^{-1} Mpc for galaxy pairs. These results agree well with the global richness dependence of the correlation functions of galaxy systems.

A sample of nearly 9000 early-type galaxies, in the redshift range 0.01<z<0.3, was selected from the Sloan Digital Sky Survey using morphological and spectral criteria. The sample was used to study how early-type galaxy observables, including luminosity L, effective radius R, surface brightness I, color, and velocity dispersion V, are correlated with one another. Measurement biases are understood with mock catalogs which reproduce all of the observed scaling relations. At any given redshift, the intrinsic distribution of luminosities, sizes and velocity dispersions in our sample are all approximately Gaussian. In the r* band L ~ V^3.91, L ~ R^1.58, R ~ I^(-0.75), and the Fundamental Plane relation is R ~ V^(1.49) I^(-0.75). These relations are approximately the same in the g*, i* and z* bands. At fixed luminosity, the mass-to-light ratio scales as M/L ~ L^0.14. The g*-r* color scales as V^0.25. Color also correlates with magnitude and size, but these correlations are entirely due to the L-V and R-V relations. Chemical evolution and star formation histories are investigated using co-added spectra of similar objects in our sample. Chemical abundances correlate primarily with velocity dispersion. At fixed V, the higher redshift population is bluer, is weaker in Mg2, and is stronger in Hbeta than the population nearby. In addition, the population at higher redshifts is slightly more luminous. These differences are consistent with that of a passively evolving population which formed the bulk of its stars about 9 Gyrs ago. The Fundamental Plane suggests that galaxies in dense regions are slightly different from those in less dense regions, but the co-added spectra and color--magnitude relations show no statistically significant dependence on environment.

The commissioning year of the Sloan Digital Sky Survey has demonstrated that many cataclysmic variables have been missed in previous surveys with brighter limits. We report the identification of 22 cataclysmic variables, of which 19 are new discoveries and 3 are known systems (SW UMa, BH Lyn and Vir4). A compendium of positions, colors and characteristics of these systems obtained from the SDSS photometry and spectroscopy is presented along with data obtained during follow-up studies with the Apache Point Observatory (APO) and Manastash Ridge Observatory (MRO) telescopes. We have determined orbital periods for 3 of the new systems: two show dwarf nova outbursts, and the third is a likely magnetic system with eclipses of its region of line emission. Based on these results, we expect the completed survey to locate at least 400 new CVs. Most of these will be faint systems with low accretion rates that will provide new constraints on binary evolution models.

Low dispersion optical spectra have been obtained with the Hobby-Eberly Telescope of 22 very red objects found in early imaging data from the Sloan Digital Sky Survey. The objects are assigned spectral types on the 2MASS system (Kirkpatrick et al. 1999) and are found to range from late M to late L. The red- and near-infrared colors from SDSS and 2MASS correlate closely with each other, and most of the colors are closely related to spectral type in this range; the exception is the (i^* - z^*) color, which appears to be independent of spectral type between about M7 and L4. The spectra suggest that this independence is due to the disappearance of the TiO and VO absorption in the i-band for later spectral types; to the presence of strong Na I and K I absorption in the i-band; and to the gradual disappearance of the 8400 Angstrom absorption of TiO and FeH in the z-band.

We present a comparison of three cluster finding algorithms from imaging data using Monte Carlo simulations of clusters embedded in a 25 deg^2 region of Sloan Digital Sky Survey (SDSS) imaging data: the Matched Filter (MF; Postman et al. 1996), the Adaptive Matched Filter (AMF; Kepner et al. 1999) and a color-magnitude filtered Voronoi Tessellation Technique (VTT). Among the two matched filters, we find that the MF is more efficient in detecting faint clusters, whereas the AMF evaluates the redshifts and richnesses more accurately, therefore suggesting a hybrid method (HMF) that combines the two. The HMF outperforms the VTT when using a background that is uniform, but it is more sensitive to the presence of a non-uniform galaxy background than is the VTT; this is due to the assumption of a uniform background in the HMF model. We thus find that for the detection thresholds we determine to be appropriate for the SDSS data, the performance of both algorithms are similar; we present the selection function for each method evaluated with these thresholds as a function of redshift and richness. For simulated clusters generated with a Schechter luminosity function (M_r^* = -21.5 and alpha = -1.1) both algorithms are complete for Abell richness >= 1 clusters up to z ~ 0.4 for a sample magnitude limited to r = 21. While the cluster parameter evaluation shows a mild correlation with the local background density, the detection efficiency is not significantly affected by the background fluctuations, unlike previous shallower surveys.

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.

We present new measurements of scaling laws relating the luminosity of galaxies to the amplitude and shape of their dark matter halos. Early imaging and spectroscopic data from the Sloan Digital Sky Survey are used to make weak lensing measurements of the surface mass density contrast Delta Sigma_+ around classes of lens objects. This surface mass density contrast as a function of radius is a measure of the galaxy-mass correlation function (GMCF). Because spectroscopic redshifts are available for all lens objects, the mass and distance scales are well constrained. The GMCF measured around ~31,000 lenses is well fit by a power law of the form Delta Sigma_+ = (2.5+0.7-0.6) (R/1 Mpc)^{-0.8+-0.2} h M_sun pc^-2. We compare this GMCF to galaxy luminosity, type, and environment, and find that it varies strongly with all three. We quantify these variations by comparing the normalization of a fit to the inner 260 h^-1 kpc, M_260, to the galaxy luminosity. While M_260 is not strongly related to luminosity in bluest band u', there is a simple, linear relation between M_260 and luminosity in redder bands (g', r', i', and z'). We test the universality of these mass-to-light scalings by independently measuring them for spiral and elliptical galaxies,and for galaxies in a variety of environments. We find remarkable consistency in these determinations in the red bands, especially i' and z'. This consistency across a wide range of systems suggests that the measured scaling represents an excellent cosmic average, and that the integrated star formation history of galaxies is strongly related to the dark matter environments in which they form.

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

The discovery of an optical counterpart to GRB010222 (detected by BeppoSAX; Piro 2001) was announced 4.4 hrs after the burst by Henden (2001a). The Sloan Digital Sky Survey's 0.5m photometric telescope (PT) and 2.5m survey telescope were used to observe the afterglow of GRB010222 starting 4.8 hours after the GRB. The 0.5m PT observed the afterglow in five, 300 sec g' band exposures over the course of half an hour, measuring a temporal decay rate in this short period of F_nu \propto t^{-1.0+/-0.5}. The 2.5m camera imaged the counterpart nearly simultaneously in five filters (u' g' r' i' z'), with r' = 18.74+/-0.02 at 12:10 UT. These multicolor observations, corrected for reddening and the afterglow's temporal decay, are well fit by the power-law F_nu \propto nu^{-0.90+/-0.03} with the exception of the u' band UV flux which is 20% below this slope. We examine possible interpretations of this spectral shape, including source extinction in a star forming region.

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 present a lensing study of 42 galaxy clusters imaged in Sloan Digital Sky Survey (SDSS) commissioning data. Cluster candidates are selected optically from SDSS imaging data and confirmed for this study by matching to X-ray sources found independently in the ROSAT all sky survey (RASS). Five color SDSS photometry is used to make accurate photometric redshift estimates that are used to rescale and combine the lensing measurements. The mean shear from these clusters is detected to 2 h-1 Mpc at the 7-sigma level, corresponding to a mass within that radius of 4.2 +/- 0.6 x 10^14 h-1 M_sun. The shear profile is well fit by a power law with index -0.9 +/- 0.3, consistent with that of an isothermal density profile. This paper demonstrates our ability to measure ensemble cluster masses from SDSS imaging data.

Early data taken during commissioning of the SDSS have resulted in the discovery of a very cool white dwarf. It appears to have stronger collision induced absorption from molecular hydrogen than any other known white dwarf, suggesting it has a cooler temperature than any other. While its distance is presently unknown, it has a surprisingly small proper motion, making it unlikely to be a halo star. An analysis of white dwarf cooling times suggests that this object may be a low-mass star with a helium core. The SDSS imaging and spectroscopy also recovered LHS 3250, the coolest previously known white dwarf, indicating that the SDSS will be an effective tool for identifying these extreme objects.

We report the discovery of two well-defined tidal tails emerging from the sparse remote globular cluster Palomar 5. These tails stretch out symmetrically to both sides of the cluster in the direction of constant Galactic latitude and subtend an angle of 2.6 degrees on the sky. The tails have been detected in commissioning data of the Sloan Digital Sky Survey (SDSS), providing deep five-color photometry in a 2.5 degrees wide band along the equator. The stars in the tails make up a substantial part (~1/3) of the current total population of cluster stars in the magnitude interval 19.5 < i* < 22.0. This reveals that the cluster is subject to heavy mass loss. The orientation of the tails provides an important key for the determination of the cluster's Galactic orbit.

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.

Understanding how the physical properties of galaxies (e.g. their spectral type or age) evolve as a function of redshift relies on having an accurate representation of galaxy spectral energy distributions. While it has been known for some time that galaxy spectra can be reconstructed from a handful of orthogonal basis templates, the underlying basis is poorly constrained. The limiting factor has been the lack of large samples of galaxies (covering a wide range in spectral type) with high signal-to-noise spectrophotometric observations. To alleviate this problem we introduce here a new technique for reconstructing galaxy spectral energy distributions directly from samples of galaxies with broadband photometric data and spectroscopic redshifts. Exploiting the statistical approach of the Karhunen-Loeve expansion, our iterative training procedure increasingly improves the eigenbasis, so that it provides better agreement with the photometry. We demonstrate the utility of this approach by applying these improved spectral energy distributions to the estimation of photometric redshifts for the HDF sample of galaxies. We find that in a small number of iterations the dispersion in the photometric redshifts estimator (a comparison between predicted and measured redshifts) can decrease by up to a factor of 2.

A sample of 4208 objects with magnitude 15 < g* < 22 and colors of main sequence A stars has been selected from 370 square degrees of Sloan Digital Sky Survey (SDSS) commissioning observations. The data is from two long, narrow stripes, each with an opening angle of greater than 60 deg, at Galactic latitudes 36 < abs(b) < 63 on the celestial equator. An examination of the sample's distribution shows that these stars trace considerable substructure in the halo. Large overdensities of A-colored stars in the North at (l,b,R) = (350, 50, 46 kpc) and in the South at (157, -58, 33 kpc) and extending over tens of degrees are present in the halo of the Milky Way. Using photometry to separate the stars by surface gravity, both structures are shown to contain a sequence of low surface gravity stars consistent with identification as a blue horizontal branch (BHB). Both structures also contain a population of high surface gravity stars two magnitudes fainter than the BHB stars, consistent with their identification as blue stragglers (BSs). From the numbers of detected BHB stars, lower limits to the implied mass of the structures are 6x10^6 M_sun and 2x10^6 M_sun. The fact that two such large clumps have been detected in a survey of only 1% of the sky indicates that such structures are not uncommon in the halo. Simple spheroidal parameters are fit to a complete sample of the remaining unclumped BHB stars and yield (at r < 40 kpc) a fit to a halo distribution with flattening (c/a = 0.65+/-0.2) and a density falloff exponent of alpha = -3.2+/-0.3.

We present a novel approach to photometric redshifts, one that merges the advantages of both the template fitting and empirical fitting algorithms, without any of their disadvantages. This technique derives a set of templates, describing the spectral energy distributions of galaxies, from a catalog with both multicolor photometry and spectroscopic redshifts. The algorithm is essentially using the shapes of the templates as the fitting parameters. From simulated multicolor data we show that for a small training set of galaxies we can reconstruct robustly the underlying spectral energy distributions even in the presence of substantial errors in the photometric observations. We apply these techniques to the multicolor and spectroscopic observations of the Hubble Deep Field building a set of template spectra that reproduced the observed galaxy colors to better than 10%. Finally we demonstrate that these improved spectral energy distributions lead to a photometric-redshift relation for the Hubble Deep Field that is more accurate than standard template-based approaches.

We present a novel method capable of creating optimal eigenspectra from multicolor redshift surveys for photometric redshift estimation. Our iterative training algorithm modifies the templates to represent the photometric measurements better. We present a short description of our algorithm here. We show that the corrected templates give more precise photometric redshifts, essentially a ``free'' feature, since we were not fitting for the redshifts themselves.

Results of a blind test of photometric redshift predictions against spectroscopic galaxy redshifts obtained in the Hubble Deep Field with the Keck Telescope are presented. The best photometric redshift schemes predict spectroscopic redshifts with a redshift accuracy of |Delta-z|<0.1 for more than 68 percent of sources and with |Delta-z|<0.3 for 100 percent, when single-feature spectroscopic redshifts are removed from consideration. This test shows that photometric redshift schemes work well at least when the photometric data are of high quality and when the sources are at moderate redshifts.

As a means of better understanding the evolution of optically selected galaxies we consider the distribution of galaxies within the multicolor space $U$, $B_J$, $R_F$ and $I_N$. We find that they form an almost planar distribution out to $B_J =22.5$ and $z<0.3$. The position of a galaxy within this plane is dependent on its redshift, luminosity and spectral type. While in the original $U$, $B_J$, $R_F$ and $I_N$ space these properties are highly correlated we can define an optimal rotation of the photometric axes that makes much of this information orthogonal. Fitting the observed spectroscopic redshifts with a quadratic function of the four magnitudes we show that redshifts for galaxies can be estimated to an accuracy better than $\Delta z =0.05$. This dispersion is due to the photometric uncertainties within the photographic data. Assuming no galaxy evolution we derive a set of simulated galaxy fluxes in the U, J, F and N passbands. Using these data we investigate how the redshift is encoded within the broadband magnitudes and the intrinsic dispersion of the photometric-redshift relation. We find that the signal that defines a galaxy's photometric redshift is not related to specific absorption or emission lines but comes from the break in the overall shape of the galaxy continuum at around 4000 \AA. Using high signal-to-noise photometric data we estimate that it is possible to achieve an intrinsic dispersion of less than $\Delta z =0.02$.