The advent of new time domain surveys and the imminent increase in astronomical data expose the shortcomings in traditional time series analysis (such as power spectra analysis) in characterising the abundantly varied, complex and stochastic light curves of Active Galactic Nuclei (AGN). Recent applications of novel methods from non-linear dynamics have shown promise in characterising higher modes of variability and time-scales in AGN. Recurrence analysis in particular can provide complementary information about characteristic time-scales revealed by other methods, as well as probe the nature of the underlying physics in these objects. Recurrence analysis was developed to study the recurrences of dynamical trajectories in phase space, which can be constructed from one-dimensional time series such as light curves. We apply the methods of recurrence analysis to two optical light curves of Kepler-monitored AGN. We confirm the detection and period of an optical quasi-periodic oscillation in one AGN, and confirm multiple other time-scales recovered from other methods ranging from 5 days to 60 days in both objects. We detect regions in the light curves that deviate from regularity, provide evidence of determinism and non-linearity in the mechanisms underlying one light curve (KIC 9650712), and determine a linear stochastic process recovers the dominant variability in the other light curve (Zwicky 229--015). We discuss possible underlying processes driving the dynamics of the light curves and their diverse classes of variability.

This work develops application techniques for stochastic modelling of Active Galactic Nuclei (AGN) variability as a probe of accretion disk physics. Stochastic models, specifically Continuous Auto-Regressive Moving Average (CARMA) models, characterize lightcurves by estimating delay timescales that describe movements away from and toward equilibrium (mean flux) as well as an amplitude and frequency of intrinsic perturbations to the AGN flux. We begin this tutorial by reviewing discrete auto-regressive (AR) and moving-average (MA) processes, we bridge these components to their continuous analogs, and lastly we investigate the significance of timescales from direct stochastic modelling of a lightcurve projected in power spectrum (PSD) and structure function (SF) space. We determine that higher order CARMA models, for example the Damped Harmonic Oscillator (DHO or CARMA(2,1)) are more sensitive to deviations from a single-slope power-law description of AGN variability; unlike Damped Random Walks (DRW or CAR(1)) where the PSD slope is fixed, the DHO slope is not. Higher complexity stochastic models than the DRW capture additional covariance in data and output additional characteristic timescales that probe the driving mechanisms of variability.

We study how the void environment affects galactic chemical evolution by comparing the oxygen and nitrogen abundances of dwarf galaxies in voids with dwarf galaxies in denser regions. Using spectroscopic observations from SDSS DR7, we estimate oxygen, nitrogen, and neon abundances of 889 void dwarf galaxies and 672 dwarf galaxies in denser regions. A substitute for the [OII] 3727 doublet is developed, permitting oxygen abundance estimates of SDSS dwarf galaxies at all redshifts with the Direct Te method. We find that void dwarf galaxies have about the same oxygen abundance and Ne/O ratio, slightly higher neon abundances, and slightly lower nitrogen abundance and N/O ratio than dwarf galaxies in denser environments. We conclude that the void environment has a slight influence on dwarf galaxy chemical evolution. Our mass-N/O relationship shows that the secondary production of nitrogen commences at a lower stellar mass in void dwarf galaxies than in dwarf galaxies in denser environments. Our dwarf galaxy sample demonstrates a strong anti-correlation between the sSFR and N/O ratio, providing evidence that oxygen is produced in higher mass stars than those which synthesize nitrogen. The lower N/O ratios and smaller stellar mass for secondary nitrogen production seen in void dwarf galaxies may indicate both delayed star formation and a dependence of cosmic downsizing on the large-scale environment. A shift toward slightly higher oxygen abundances in void dwarf galaxies could be evidence of larger ratios of dark matter halo mass to stellar mass in voids than in denser regions.

AGN exhibit rapid, high amplitude stochastic flux variations across the entire electromagnetic spectrum on timescales ranging from hours to years. The cause of this variability is poorly understood. We present a Green's Function-based method for using variability to (1) measure the time-scales on which flux perturbations evolve and (2) characterize the driving flux perturbations. We model the observed light curve of an AGN as a linear differential equation driven by stochastic impulses. We analyze the light curve of the Kepler AGN Zw 229-15 and find that the observed variability behavior can be modeled as a damped harmonic oscillator perturbed by a colored noise process. The model powerspectrum turns over on time-scale $385$~d. On shorter time-scales, the log-powerspectrum slope varies between $2$ and $4$, explaining the behavior noted by previous studies. We recover and identify both the $5.6$~d and $67$~d timescales reported by previous work using the Green's Function of the C-ARMA equation rather than by directly fitting the powerspectrum of the light curve. These are the timescales on which flux perturbations grow, and on which flux perturbations decay back to the steady-state flux level respectively. We make the software package KALI used to study light curves using our method available to the community.

We study how the cosmic environment affects galaxy evolution in the Universe by comparing the metallicities of dwarf galaxies in voids with dwarf galaxies in more dense regions. Ratios of the fluxes of emission lines, particularly those of the forbidden [O III] and [S II] transitions, provide estimates of a region's electron temperature and number density. From these two quantities and the emission line fluxes [O II] 3727, [O III] 4363, and [O III] 4959,5007, we estimate the abundance of oxygen with the Direct Te method. We estimate the metallicity of 42 blue, star-forming void dwarf galaxies and 89 blue, star-forming dwarf galaxies in more dense regions using spectroscopic observations from the Sloan Digital Sky Survey Data Release 7, as re-processed in the MPA-JHU value-added catalog. We find very little difference between the two sets of galaxies, indicating little influence from the large-scale environment on their chemical evolution. Of particular interest are a number of extremely metal-poor dwarf galaxies that are less prevalent in voids than in the denser regions.

We examine how the cosmic environment affects the chemical evolution of galaxies in the Universe by comparing the N/O ratio of dwarf galaxies in voids with dwarf galaxies in more dense regions. Ratios of the forbidden [O III] and [S II] transitions provide estimates of a region's electron temperature and number density. We estimate the abundances of oxygen and nitrogen using these temperature and density estimates and the emission line fluxes [O II] 3727, [O III] 4959, 5007, and [N II] 6548, 6584 with the direct Te method. Using spectroscopic observations from the Sloan Digital Sky Survey Data Release 7, we are able to estimate the N/O ratio in 42 void dwarf galaxies and 89 dwarf galaxies in more dense regions. The N/O ratio for void dwarfs (Mr > -17) is slightly lower (12%) than for dwarf galaxies in denser regions. We also estimate the nitrogen and oxygen abundances of 2050 void galaxies and 3883 galaxies in more dense regions with Mr > -20. These somewhat brighter galaxies (but still fainter than L*) also display similar minor shifts in the N/O ratio. The shifts in the average and median element abundance values in all absolute magnitude bins studied are in the same direction, suggesting that the large-scale environment may influence the chemical evolution of galaxies. We discuss possible causes of such a large-scale environmental dependence of the chemical evolution of galaxies, including retarded star formation and a higher dark matter halo mass to stellar mass ratio in void galaxies.

We measure the star formation properties of two large samples of galaxies from the SDSS in large-scale cosmic voids on time scales of 10 Myr and 100 Myr, using H$\alpha$ emission line strengths and GALEX FUV fluxes, respectively. The first sample consists of 109,818 optically selected galaxies. We find that void galaxies in this sample have higher specific star formation rates (SSFRs; star formation rates per unit stellar mass) than similar stellar mass galaxies in denser regions. The second sample is a subset of the optically selected sample containing 8070 galaxies with reliable HI detections from ALFALFA. For the full HI detected sample, SSFRs do not vary systematically with large-scale environment. However, investigating only the HI detected dwarf galaxies reveals a trend towards higher SSFRs in voids. Furthermore, we estimate the star formation rate per unit HI mass (known as the star formation efficiency; SFE) of a galaxy, as a function of environment. For the overall HI detected population, we notice no environmental dependence. Limiting the sample to dwarf galaxies again reveals a trend towards higher SFEs in voids. These results suggest that void environments provide a nurturing environment for dwarf galaxy evolution allowing for higher specific star formation rates and efficiencies.

We present a study of spectral properties of galaxies in underdense large-scale structures, voids. Our void galaxy sample (75,939 galaxies) is selected from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) with $\rm z < 0.107$. We find that there are no significant differences in the luminosities, stellar masses, stellar populations, and specific star formation rates between void galaxies of specific spectral types and their wall counterparts. However, the fraction of star-forming galaxies in voids is significantly higher ($\ge 9\%$) than that in walls. Void galaxies, when considering all spectral types, are slightly fainter, less massive, have younger stellar populations and of higher specific star formation rates than wall galaxies. These minor differences are totally caused by the higher fraction of star-forming galaxies in voids. We confirm that AGNs exist in voids, already found by \cite{co08}, with similar abundance as in walls. Type I AGNs contribute $\sim$ 1\%-2\% of void galaxies, similar to their fraction in walls. The intrinsic [O III] luminosities , spanning over $10^6\ L_{\sun} \sim 10^9\ L_{\sun}$, and Eddington ratios are similar comparing our void AGNs versus wall AGNs. Small scale statistics show that all spectral types of void galaxies are less clustered than their counterparts in walls. Major merger may not be the dominant trigger of black hole accretion in the luminosity range we probe. Our study implies that the growth of black holes relies weakly on large scale structures.

We measure the r-band galaxy luminosity function (LF) across environments over the redshift range 0<$z$<0.107 using the SDSS. We divide our sample into galaxies residing in large scale voids (void galaxies) and those residing in denser regions (wall galaxies). The best fitting Schechter parameters for void galaxies are: log$\Phi^*$= -3.40$\pm$0.03 log(Mpc$^{-3}$), $M^*$= -19.88$\pm$0.05, and $\alpha$=-1.20$\pm$0.02. For wall galaxies, the best fitting parameters are: log$\Phi^*$=-2.86$\pm$0.02 log(Mpc$^{-3}$), $M^*$=-20.80$\pm$0.03, and $\alpha$=-1.16$\pm$0.01. We find a shift in the characteristic magnitude, $M^*$, towards fainter magnitudes for void galaxies and find no significant difference between the faint-end slopes of the void and wall galaxy LFs. We investigate how low surface brightness selections effects can affect the galaxy LF. To attempt to examine a sample of galaxies that is relatively free of surface brightness selection effects, we compute the optical galaxy LF of galaxies detected by the blind HI survey, ALFALFA. We find that the global LF of the ALFALFA sample is not well fit by a Schechter function, because of the presence of a wide dip in the LF around $M_r$=-18 and an upturn at fainter magnitudes ($\alpha$~-1.47). We compare the HI selected r-band LF to various LFs of optically selected populations to determine where the HI selected optical LF obtains its shape. We find that sample selection plays a large role in determining the shape of the LF.

We gauge the impact of spacecraft-induced effects on the inferred variability properties of the light curve of the Seyfert 1 AGN Zw 229-15 observed by \Kepler. We compare the light curve of Zw 229-15 obtained from the Kepler MAST database with a re-processed light curve constructed from raw pixel data (Williams & Carini, 2015). We use the first-order structure function, $SF(\delta t)$, to fit both light curves to the damped power-law PSD of Kasliwal, Vogeley & Richards, 2015. On short timescales, we find a steeper log-PSD slope ($\gamma = 2.90$ to within $10$ percent) for the re-processed light curve as compared to the light curve found on MAST ($\gamma = 2.65$ to within $10$ percent)---both inconsistent with a damped random walk which requires $\gamma = 2$. The log-PSD slope inferred for the re-processed light curve is consistent with previous results (Carini & Ryle, 2012, Williams & Carini, 2015) that study the same re-processed light curve. The turnover timescale is almost identical for both light curves ($27.1$ and $27.5$~d for the reprocessed and MAST database light curves). Based on the obvious visual difference between the two versions of the light curve and on the PSD model fits, we conclude that there remain significant levels of spacecraft-induced effects in the standard pipeline reduction of the Kepler data. Re-processing the light curves will change the model inferenced from the data but is unlikely to change the overall scientific conclusion reached by Kasliwal et al. 2015---not all AGN light curves are consistent with the DRW.

We test the consistency of active galactic nuclei (AGN) optical flux variability with the $\textit{damped random walk}$ (DRW) model. Our sample consists of 20 multi-quarter $\textit{Kepler}$ AGN light curves including both Type 1 and 2 Seyferts, radio-loud and -quiet AGN, quasars, and blazars. $\textit{Kepler}$ observations of AGN light curves offer a unique insight into the variability properties of AGN light curves because of the very rapid ($11.6-28.6$ min) and highly uniform rest-frame sampling combined with a photometric precision of $1$ part in $10^{5}$ over a period of 3.5 yr. We categorize the light curves of all 20 objects based on visual similarities and find that the light curves fall into 5 broad categories. We measure the first order structure function of these light curves and model the observed light curve with a general broken power-law PSD characterized by a short-timescale power-law index $\gamma$ and turnover timescale $\tau$. We find that less than half the objects are consistent with a DRW and observe variability on short timescales ($\sim 2$ h). The turnover timescale $\tau$ ranges from $\sim 10-135$ d. Interesting structure function features include pronounced dips on rest-frame timescales ranging from $10-100$ d and varying slopes on different timescales. The range of observed short-timescale PSD slopes and the presence of dip and varying slope features suggests that the DRW model may not be appropriate for all AGN. We conclude that AGN variability is a complex phenomenon that requires a more sophisticated statistical treatment.

We measure the HI mass function (HIMF) and velocity width function (WF) across environments over a range of masses $7.2<\log(M_{HI}/M_{\odot})<10.8$, and profile widths $1.3\log(km/s)<\log(W)<2.9\log(km/s)$, using a catalog of ~7,300 HI-selected galaxies from the ALFALFA Survey, located in the region of sky where ALFALFA and SDSS (Data Release 7) North overlap. We divide our galaxy sample into those that reside in large-scale voids (void galaxies) and those that live in denser regions (wall galaxies). We find the void HIMF to be well fit by a Schechter function with normalization $\Phi^*=(1.37\pm0.1)\times10^{-2} h^3Mpc^{-3}$, characteristic mass $\log(M^*/M_{\odot})+2\log h_{70}=9.86\pm0.02$, and low-mass-end slope $\alpha=-1.29\pm0.02$. Similarly, for wall galaxies, we find best-fitting parameters $\Phi^*=(1.82\pm0.03)\times10^{-2} h^3Mpc^{-3}$, $\log(M^*/M_{\odot})+2\log h_{70}=10.00\pm0.01$, and $\alpha=-1.35\pm0.01$. We conclude that void galaxies typically have slightly lower HI masses than their non-void counterparts, which is in agreement with the dark matter halo mass function shift in voids assuming a simple relationship between DM mass and HI mass. We also find that the low-mass slope of the void HIMF is similar to that of the wall HIMF suggesting that there is either no excess of low-mass galaxies in voids or there is an abundance of intermediate HI mass galaxies. We fit a modified Schechter function to the ALFALFA void WF and determine its best-fitting parameters to be $\Phi^*=0.21\pm0.1 h^3Mpc^{-3}$, $\log(W^*)=2.13\pm0.3$, $\alpha=0.52\pm0.5$ and high-width slope $\beta=1.3\pm0.4$. For wall galaxies, the WF parameters are: $\Phi^*=0.022\pm0.009 h^3Mpc^{-3}$, $\log(W^*)=2.62\pm0.5$, $\alpha=-0.64\pm0.2$ and $\beta=3.58\pm1.5$. Because of large uncertainties on the void and wall width functions, we cannot conclude whether the WF is dependent on the environment.

We explore the mid-infrared (mid-IR) through ultraviolet (UV) spectral energy distributions (SEDs) of 119,652 luminous broad-lined quasars with 0.064<z<5.46 using mid-IR data from Spitzer and WISE, near-infrared data from Two Micron All Sky Survey and UKIDSS, optical data from Sloan Digital Sky Survey, and UV data from Galaxy Evolution Explorer. The mean SED requires a bolometric correction (relative to 2500A) of BC=2.75+-0.40 using the integrated light from 1um-2keV, and we further explore the range of bolometric corrections exhibited by individual objects. In addition, we investigate the dependence of the mean SED on various parameters, particularly the UV luminosity for quasars with 0.5<z<3 and the properties of the UV emission lines for quasars with z>1.6; the latter is a possible indicator of the strength of the accretion disk wind, which is expected to be SED dependent. Luminosity-dependent mean SEDs show that, relative to the high-luminosity SED, low-luminosity SEDs exhibit a harder (bluer) far-UV spectral slope, a redder optical continuum, and less hot dust. Mean SEDs constructed instead as a function of UV emission line properties reveal changes that are consistent with known Principal Component Analysis (PCA) trends. A potentially important contribution to the bolometric correction is the unseen extream-UV (EUV) continuum. Our work suggests that lower-luminosity quasars and/or quasars with disk-dominated broad emission lines may require an extra continuum component in the EUV that is not present (or much weaker) in high-luminosity quasars with strong accretion disk winds. As such, we consider four possible models and explore the resulting bolometric corrections. Understanding these various SED-dependent effects will be important for accurate determination of quasar accretion rates.

Using the sample presented in Pan:2011, we analyse the photometric properties of 88,794 void galaxies and compare them to galaxies in higher density environments with the same absolute magnitude distribution. In Pan et al. (2011), we found a total of 1054 dynamically distinct voids in the SDSS with radius larger than 10h^-1 Mpc. The voids are underdense, with delta rho/rho < -0.9 in their centers. Here we study the photometric properties of these void galaxies. We look at the u - r colours as an indication of star formation activity and the inverse concentration index as an indication of galaxy type. We find that void galaxies are statistically bluer than galaxies found in higher density environments with the same magnitude distribution. We examine the colours of the galaxies as a function of magnitude, and we fit each colour distribution with a double-Gaussian model for the red and blue subpopulations. As we move from bright to dwarf galaxies, the population of red galaxies steadily decreases and the fraction of blue galaxies increases in both voids and walls, however the fraction of blue galaxies in the voids is always higher and bluer than in the walls. We also split the void and wall galaxies into samples depending on galaxy type. We find that late type void galaxies are bluer than late type wall galaxies and the same holds for early galaxies. We also find that early type, dwarf void galaxies are blue in colour. We also study the properties of void galaxies as a function of their distance from the center of the void. We find very little variation in the properties, such as magnitude, colour and type, of void galaxies as a function of their location in the void. The only exception is that the dwarf void galaxies may live closer to the center. The centers of voids have very similar density contrast and hence all void galaxies live in very similar density environments (ABRIDGED)

We study the distribution of cosmic voids and void galaxies using Sloan Digital Sky Survey Data Release 7 (SDSS DR7). Using the VoidFinder algorithm as described by Hoyle 2002, we identify 1054 statistically significant voids in the northern galactic hemisphere with radii > 10 h^{-1} Mpc. The filling factor of voids in the sample volume is 62%. The largest void is just over 30 h^{-1} Mpc in effective radius. The median effective radius is 17 h^{-1} Mpc. The voids are found to be significantly underdense, with density contrast \delta < -0.85 at the edges of the voids. The radial density profiles of these voids are similar to predictions of dynamically distinct underdensities in gravitational theory. We find 8,046 galaxies brighter than M_r = -20.09 within the voids, accounting for 7% of the galaxies. We compare the results of VoidFinder on SDSS DR7 to mock catalogs generated from a SPH halo model simulation as well as other \Lambda -CDM simulations and find similar void fractions and void sizes in the data and simulations. This catalog is made publicly available at http://www.physics.drexel.edu/~pan/voidcatalog.html for download.

Nowadays, science has been coming into a new paradigm, called data-intensive science. While current studies of the new phenomenon focused on building up infrastructure for this new paradigm, yet a few studies concern users of scientific data, particularly their usage practices in the newly emerging paradigm, even though the importance of understanding users' work flow and practices has been summoned. This study endeavors to improve our understanding of users' data usage behavior through a content analysis of publications in a frequently cited new paradigm-related project, Sloan Digital Sky Survey (SDSS). We found that (1) nearly half studies used one data source only. A few studies exploited three or more data sources; (2) the number of objects that were analyzed in SDSS publications is in all scales from one digit to millions; (3) different paper types may affect the data usage patterns; (4) Users are not only consumers of scientific data. They are producers too; (5) studies that can use multiple large scale data sources are relative rare. Issues of data provenance, trust, and usability may prevent researchers from doing this kind of research.

We measure the topology of the main galaxy distribution using the Seventh Data Release of the Sloan Digital Sky Survey, examining the dependence of galaxy clustering topology on galaxy properties. The observational results are used to test galaxy formation models. A volume-limited sample defined by $M_r<-20.19$ enables us to measure the genus curve with amplitude of $G=378$ at $6h^{-1}$Mpc smoothing scale, with 4.8\% uncertainty including all systematics and cosmic variance. The clustering topology over the smoothing length interval from 6 to $10 h^{-1}$Mpc reveals a mild scale-dependence for the shift ($\Delta\nu$) and void abundance ($A_V$) parameters of the genus curve. We find substantial bias in the topology of galaxy clustering with respect to the predicted topology of the matter distribution, which varies with luminosity, morphology, color, and the smoothing scale of the density field. The distribution of relatively brighter galaxies shows a greater prevalence of isolated clusters and more percolated voids. Even though early (late)-type galaxies show topology similar to that of red (blue) galaxies, the morphology dependence of topology is not identical to the color dependence. In particular, the void abundance parameter $A_V$ depends on morphology more strongly than on color. We test five galaxy assignment schemes applied to cosmological N-body simulations of a $\Lambda$CDM universe to generate mock galaxies: the Halo-Galaxy one-to-one Correspondence model, the Halo Occupation Distribution model, and three implementations of Semi-Analytic Models (SAMs). None of the models reproduces all aspects of the observed clustering topology; the deviations vary from one model to another but include statistically significant discrepancies in the abundance of isolated voids or isolated clusters and the amplitude and overall shift of the genus curve. (Abridged)

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

We review the current standard model for the evolution of the Universe from an early inflationary epoch to the complex hierarchy of structure seen today. We summarize and provide key references for the following topics: observations of the expanding Universe; the hot early Universe and nucleosynthesis; theory and observations of the cosmic microwave background; Big Bang cosmology; inflation; dark matter and dark energy; theory of structure formation; the cold dark matter model; galaxy formation; cosmological simulations; observations of galaxies, clusters, and quasars; statistical measures of large-scale structure; and measurement of cosmological parameters. We conclude with discussion of some open questions in cosmology. This review is designed to provide a graduate student or other new worker in the field an introduction to the cosmological literature.

Despite a history that dates back at least a quarter of a century studies of voids in the large--scale structure of the Universe are bedevilled by a major problem: there exist a large number of quite different void--finding algorithms, a fact that has so far got in the way of groups comparing their results without worrying about whether such a comparison in fact makes sense. Because of the recent increased interest in voids, both in very large galaxy surveys and in detailed simulations of cosmic structure formation, this situation is very unfortunate. We here present the first systematic comparison study of thirteen different void finders constructed using particles, haloes, and semi--analytical model galaxies extracted from a subvolume of the Millennium simulation. The study includes many groups that have studied voids over the past decade. We show their results and discuss their differences and agreements. As it turns out, the basic results of the various methods agree very well with each other in that they all locate a major void near the centre of our volume. Voids have very underdense centres, reaching below 10 percent of the mean cosmic density. In addition, those void finders that allow for void galaxies show that those galaxies follow similar trends. For example, the overdensity of void galaxies brighter than $m_B = -20 $ is found to be smaller than about -0.8 by all our void finding algorithms.

We present a comprehensive study of accretion activity in the most underdense environments in the universe, the voids, based on the SDSS DR2 data. Based on investigations of multiple void regions, we show that AGN's occurrence rate and properties differ from those in walls. AGN are more common in voids than in walls, but only among moderately luminous and massive galaxies (M_r < -20, log M_*/M_sun < 10.5), and this enhancement is more pronounced for the weakly accreting systems (i.e., L_[O III] < 10^39 erg/s). Void AGN hosted by moderately massive and luminous galaxies are accreting at equal or lower rates than their wall counterparts, show less obscuration than in walls, and similarly aged stellar populations. The very few void AGN in massive bright hosts accrete more strongly, are more obscured, and are associated with younger stellar emission than wall AGN. Thus, accretion strength is probably connected to the availability of fuel supply, and accretion and star-formation co-evolve and rely on the same source of fuel. Nearest neighbor statistics indicate that the weak accretion activity (LINER-like) is not influenced by the local environment. However, H IIs, Seyferts, and Transition objects prefer more grouped small scale structures, indicating that the rate at which galaxies interact with each other affects their activity. These trends support a potential H II -> Seyfert/Transition Object -> LINER evolutionary sequence that we show is apparent in many properties of actively line-emitting galaxies, in both voids and walls. The subtle differences between void and wall AGN might be explained by a longer, less disturbed duty cycle of these systems in voids.

The current view of galaxy formation holds that all massive galaxies harbor a massive black hole at their center, but that these black holes are not always in an actively accreting phase. X-ray emission is often used to identify accreting sources, but for galaxies that are not harboring quasars (low-luminosity active galaxies), the X-ray flux may be weak, or obscured by dust. To aid in the understanding of weakly accreting black holes in the local universe, a large sample of galaxies with X-ray detections is needed. We cross-match the ROSAT All Sky Survey (RASS) with galaxies from the Sloan Digital Sky Survey Data Release 4 (SDSS DR4) to create such a sample. Because of the high SDSS source density and large RASS positional errors, the cross-matched catalog is highly contaminated by random associations. We investigate the overlap of these surveys and provide a statistical test of the validity of RASS-SDSS galaxy cross-matches. SDSS quasars provide a test of our cross-match validation scheme, as they have a very high fraction of true RASS matches. We find that the number of true matches between the SDSS main galaxy sample and the RASS is highly dependent on the optical spectral classification of the galaxy; essentially no star-forming galaxies are detected, while more than 0.6% of narrow-line Seyferts are detected in the RASS. Also, galaxies with ambiguous optical classification have a surprisingly high RASS detection fraction. This allows us to further constrain the SEDs of low-luminosity active galaxies. Our technique is quite general, and can be applied to any cross-matching between surveys with well-understood positional errors.

We examine volume-limited samples from the SDSS galaxies to look for relations among internal and collective physical parameters of galaxies as faint as M_r=- 17.5, which include morphology, luminosity, color, color gradient, concentration, size, velocity dispersion, equivalent width (EW) of H_alpha line,axis ratio, the luminosity and velocity dispersion functions. At fixed morphology and luminosity, we find that bright (M_r<-20) early-types show very small dispersions in color, color gradient, concentration, size, and velocity dispersion. These dispersions increase at fainter magnitudes, where the fraction of blue early-types increases. Late-types show wider dispersions in all physical parameters compared to early types at the same luminosity. Concentration indices of early-types are well-correlated with velocity dispersion, but are insensitive to luminosity and color for bright galaxies. The slope of the Faber-Jackson relation continuously changes from 4.6 +- 0.4 to 2.7+- 0.2 when luminosity changes from M_r = -22 to -20. The size of early- types is well-correlated with stellar velocity dispersion (for >100 km/s). We find that passive spirals are well-separated from star-forming late-types at EW (H_alpha) of about 4. An interesting finding is that many physical parameters of galaxies manifest different behaviors across the absolute magnitude of about M_* +- 1. The morphology fraction as a function of luminosity depends less sensitively on large scale structure than the luminosity function (LF) does, and thus seems to be more universal. The effects of internal extinction in late-types on the completeness of volume limited samples and on the LF and morphology fraction are found to be very important. (abridged)

We investigate the dependence of physical properties of galaxies brighter than M_r=-18.0 in the SDSS on environment, as measured by local density using an adaptive smoothing kernel. We find that variations of galaxy properties with environment are almost entirely due to the dependence of morphology and luminosity on environment. Because galaxy properties depend not only on luminosity but also on morphology, it is clear that galaxy properties cannot be determined solely by dark halo mass. When morphology and luminosity are fixed, other physical properties, such as color, color-gradient, concentration, size, velocity dispersion, and star formation rate, are nearly independent of local density. The only feature is the sharp decrease of late type fraction above the critical luminosity of about M_r=-21.3 in the morphology versus luminosity relation. Weak residual dependences on environment include that of the color of late-types (bluer at lower density) and of the L-sigma relation of early-types (larger dispersion at higher density for bright galaxies). The early-type fraction is a monotonically increasing function of local density and luminosity. The morphology-density- luminosity relation should be a key constraint on galaxy formation models. We demonstrate that the dependence on environment of the morphology of galaxie originates from variations in density on effective Gaussian smoothing scales much smaller than 12 h^{-1}Mpc. We find that galaxy morphology varies both with density measured on an effective Gaussian smoothing scale of 4.7 h^{-1} Mpc and on distance to the nearest bright galaxy (particularly, for about 0.2 h^{-1}Mpc). We propose a mechanism that the morphology of galaxies in galaxy systems is transformed by the tidal force.

We discuss here the spatial clustering of Seyferts and LINERs and consequences for their central engines. We show that Seyferts are less clustered than LINERs, and that this difference is not driven by the morphology-density relation, but it is related to the difference in clustering as a function of level of activity in these systems and the amount of fuel available for accretion. LINERs, which are the most clustered among AGN, show the lowest luminosities and obscuration levels, and relatively low gas densities, suggesting that these objects harbor black holes that are relatively massive yet weakly active or inefficient in their accretion, probably due to the insufficiency of their fuel supply. Seyferts, which are weakly clustered, are very luminous, show generally high gas densities and large quantities of obscuring material, suggesting that in these systems the black holes are less massive but abundantly fueled and therefore accrete quickly and probably efficiently enough to clearly dominate the ionization.

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

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

We measure the three-dimensional topology of large-scale structure in the Sloan Digital Sky Survey (SDSS). This allows the genus statistic to be measured with unprecedented statistical accuracy. The sample size is now sufficiently large to allow the topology to be an important tool for testing galaxy formation models. For comparison, we make mock SDSS samples using several state-of-the-art N-body simulations: the Millennium run of Springel et al. (2005)(10 billion particles), Kim & Park (2006) CDM models (1.1 billion particles), and Cen & Ostriker (2006) hydrodynamic code models (8.6 billion cell hydro mesh). Each of these simulations uses a different method for modeling galaxy formation. The SDSS data show a genus curve that is broadly characteristic of that produced by Gaussian random phase initial conditions. Thus the data strongly support the standard model of inflation where Gaussian random phase initial conditions are produced by random quantum fluctuations in the early universe. But on top of this general shape there are measurable differences produced by non-linear gravitational effects (cf. Matsubara 1994), and biasing connected with galaxy formation. The N-body simulations have been tuned to reproduce the power spectrum and multiplicity function but not topology, so topology is an acid test for these models. The data show a ``meatball'' shift (only partly due to the Sloan Great Wall of Galaxies; this shift also appears in a sub-sample not containing the Wall) which differs at the 2.5\sigma level from the results of the Millennium run and the Kim & Park dark halo models, even including the effects of cosmic variance.

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

We present the first multi-parameter analysis of the narrow line AGN clustering properties. Estimates of the two-point correlation function (CF) based on SDSS DR2 data reveal that Seyferts are clearly less clustered than normal galaxies, while the clustering amplitude (r_0) of LINERs is consistent with that of the parent galaxy population. The similarities in the host properties (color and concentration index) of Seyferts and LINERs suggest that the difference in their r_0 is not driven by the morphology-density relation. We find that the luminosity of [O I] emission shows the strongest influence on AGN clustering, with low L([O I]) sources having the highest r_0. This trend is much stronger than the previously detected dependence on L([O III]), which we confirm. There is a strong correspondence between the clustering patterns of objects of given spectral type and their physical properties. LINERs, which exhibit high r_0, show the lowest luminosities and obscuration levels, and relatively low gas densities (n_e), suggesting that these objects harbor black holes that are relatively massive yet weakly active or inefficient in their accretion, probably due to the insufficiency of their fuel supply. Seyferts, which have low r_0, are luminous and show large n_e, suggesting that their black holes are less massive but accrete quickly and efficiently enough to clearly dominate the ionization. The low r_0 of the H II galaxies can be understood as a consequence of both the morphology-density and star formation rate-density relations, however, their spectral properties suggest that their centers hide amidst large amounts of obscuring material black holes of generally low mass whose activity remains relatively feeble. Our own Milky Way may be a typical such case.[abridged]

We have produced a new conformal map of the universe illustrating recent discoveries, ranging from Kuiper belt objects in the Solar system, to the galaxies and quasars from the Sloan Digital Sky Survey. This map projection, based on the logarithm map of the complex plane, preserves shapes locally, and yet is able to display the entire range of astronomical scales from the Earth's neighborhood to the cosmic microwave background. The conformal nature of the projection, preserving shapes locally, may be of particular use for analyzing large scale structure. Prominent in the map is a Sloan Great Wall of galaxies 1.37 billion light years long, 80% longer than the Great Wall discovered by Geller and Huchra and therefore the largest observed structure in the universe.

We study the luminosity and color dependence of the galaxy 2-point correlation function in the Sloan Digital Sky Survey, starting from a sample of 200,000 galaxies over 2500 deg^2. We concentrate on the projected correlation function w(r_p), which is directly related to the real space \xi(r). The amplitude of w(r_p) grows continuously with luminosity, rising more steeply above the characteristic luminosity L_*. Redder galaxies exhibit a higher amplitude and steeper correlation function at all luminosities. The correlation amplitude of blue galaxies increases continuously with luminosity, but the luminosity dependence for red galaxies is less regular, with bright red galaxies more strongly clustered at large scales and faint red galaxies more strongly clustered at small scales. We interpret these results using halo occupation distribution (HOD) models assuming concordance cosmological parameters. For most samples, an HOD model with two adjustable parameters fits the w(r_p) data better than a power-law, explaining inflections at r_p ~ 1-3 Mpc/h as the transition between the 1-halo and 2-halo regimes of \xi(r). The implied minimum mass for a halo hosting a central galaxy above a luminosity threshold L grows as M_min ~ L at low luminosities and more steeply above L_*. The mass at which an average halo has one satellite galaxy brighter than L is M_1 ~ 23 M_min(L). These results imply a conditional luminosity function (at fixed halo mass) in which central galaxies lie far above a Schechter function extrapolation of the satellite population. HOD models nicely explain the joint luminosity-color dependence of w(r_p) in terms of the color fractions of central and satellite populations as a function of halo mass. The inferred HOD properties are in good qualitative agreement with theoretical predictions.

We measure the topology of volume-limited galaxy samples selected from a parent sample of 314,050 galaxies in the Sloan Digital Sky Survey(SDSS), which is now complete enough to describe the fully three-dimensional topology and its dependence on galaxy properties. We compare the observed genus statistic G\nu_f) to predictions for a Gaussian random field and to the genus measured for mock surveys constructed from new large-volume simulations of the LambdaCDM cosmology. In this analysis we carefully examine the dependence of the observed genus statistic on the Gaussian smoothing scale R_G from 3.5 to 11 h^-1Mpc and on the luminosity of galaxies over the range -22.50<M_r<-18.5. The void multiplicity A_V is less than unity at all smoothing scales. Because A_V cannot become less than 1 through gravitational evolution, this result provides strong evidence for biased galaxy formation in low density environments. We also find clear evidence of luminosity bias of topology within the volume-limited sub-samples. The shift parameter Delta nu indicates that the genus of brighter galaxies shows a negative shift toward a ``meatball'' (i.e. cluster-dominated) topology, while faint galaxies show a positive shift toward a ``bubble'' (i.e. void-dominated) topology. The transition from negative to positive shift occurs approximately at the characteristic absolute magnitude M_{r*}=-20.4. Even in this analysis of the largest galaxy sample to date, we detect the influence of individual large-scale structures, as the shift parameter Delta nu and cluster multiplicity A_C reflect (at ~3 sigma) the presence of the Sloan Great Wall and a X-shaped structure which runs for several hundred Mpc across the survey volume.

We present a first morphometric investigation of a preliminary sample from the SDSS of 154287 galaxies with apparent magnitude 14.5<m_r<17.5 and redshift 0.001<z<0.4. We measure the Minkowski Functionals, which are a complete set of morphological descriptors. To account for the complicated wedge--like geometry of the present survey data, we construct isodensity contour surfaces from the galaxy positions in redshift space and employ two complementary methods of computing the Minkowski Functionals. We find that the observed Minkowski Functionals for SDSS galaxies are consistent with the prediction of a Lambda--dominated spatially--flat Cold Dark Matter model with random--Gaussian initial conditions, within the cosmic variance estimated from the corresponding mock catalogue. We expect that future releases of the SDSS survey will allow us to distinguish morphological differences in the galaxy distribution with regard to different morphological type and luminosity ranges.

We estimate the Mass Function of void galaxies in the second public data release of the Sloan Digital Sky Survey from a sample of 1000 galaxies with local density contrasts of delta_v < -0.6. The galaxy sample is split into ellipticals and spirals using a color-Sersic index criteria. We estimate the virial masses of ellipticals using the measured spectral line-widths along with the observed size. Projection effects and uncertainties in halo properties make mass estimates of spirals more difficult. We use an inversion of the Tully-Fisher relation to estimate the isothermal rotational velocity, and introduce a scaling factor to estimate the halo extent. We then fit the measured mass function against a theoretical Press-Schechter model, and find that the distribution of galaxies in voids appears to be nearly unbiased compared to the mass.

Using a nearest neighbor analysis, we construct a sample of void galaxies from the Sloan Digital Sky Survey (SDSS) and compare the photometric properties of these galaxies to the population of non-void (wall) galaxies. We trace the density field of galaxies using a volume-limited sample with z_{max}=0.089. Galaxies from the flux-limited SDSS with z\leq z_{max} and fewer than three volume-limited neighbors within 7h^{-1}Mpc are classified as void galaxies. This criterion implies a density contrast \delta \rho/ \rho < -0.6 around void galaxies. From 155,000 galaxies, we obtain a sub-sample of 13,742 galaxies with z\leq z_{max}, from which we identify 1,010 galaxies as void galaxies. To identify an additional 194 faint void galaxies from the SDSS in the nearby universe, r~ 72 h^{-1}Mpc, we employ volume-limited samples extracted from the Updated Zwicky Catalog and the Southern Sky Redshift Survey with z_{max}=0.025 to trace the galaxy distribution. Our void galaxies span a range of absolute magnitude from M_r=-13.5 to M_r=-22.5. Using SDSS photometry, we compare the colors, concentration indices, and Sersic indices of the void and wall samples. Void galaxies are significantly bluer than galaxies lying at higher density. The population of void galaxies with M_r ~ M* +1 and brighter is on average bluer and more concentrated (later type) than galaxies outside of voids. The latter behavior is only partly explained by the paucity of luminous red galaxies in voids. These results generally agree with the predictions of semi-analytic models for galaxy formation in cold dark matter models, which indicate that void galaxies should be relatively bluer, more disklike, and have higher specific star formation rates.

We study the spectroscopic properties of a sample of 10^3 void galaxies from the Sloan Digital Sky Survey (SDSS) and compare these with the properties of galaxies in higher density regions (wall galaxies). This sample of void galaxies covers the range of absolute magnitude from M_r=-13.5 to M_r=-22.5 in regions with density contrast delta < -0.6. In this paper we compare the equivalent widths of Halpha, [OII], [NII], Hbeta, and [OIII] of void and wall galaxies of similar luminosities and find that void galaxies have larger values, indicating that they are still forming stars at a high rate. A comparison of the Balmer break, as measured by the parameter Dn(4000), reveals that void galaxies have younger stellar populations than wall galaxies. Using standard techniques, we estimate Halpha and [OII] star formation rates of the void and wall galaxies and along with estimates of the stellar masses, we compute specific star formation rates. In most cases, we find that void galaxies have similar SFRs to wall galaxies but they are fainter and smaller mass. This means that, consistent with the EWs, void galaxies have higher specific star formation rates than wall galaxies.

Wide-angle, moderately deep redshift surveys such as that conducted as part of the Sloan Digital Sky Survey (SDSS) allow study of the relationship between the structural elements of the large-scale distribution of galaxies -- including groups, cluster, superclusters, and voids -- and the dependence of galaxy formation and evolution on these enviroments. We present a progress report on mapping efforts with the SDSS and discuss recently constructed catalogs of clusters, voids, and void galaxies, and evidence for a 420Mpc/h supercluster or ``Great Wall.'' Analysis of multi-band photometry and moderate-resolution spectroscopy from the SDSS reveals environmental dependence of the star formation history of galaxies that extends over more than a factor of 100 in density, from clusters all the way to the deep interiors of voids. On average, galaxies in the rarified environments of voids exhibit bluer colors, higher specific star formation rates, lower dust content, and more disk-like morphology than objects in denser regions. This trend persists in comparisons of samples in low vs. high-density regions with similar luminosity and morphology, thus this dependence is not simply an extension of the morphology-density relation. Large-scale modulation of the halo mass function and the temperature of the intergalactic medium might explain this dependence of galaxy evolution on the large-scale environment.

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

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

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

We present an analysis of voids in the 2dF Galaxy Redshift Survey (2dFGRS). This analysis includes identification of void regions and measurement of void statistics. The 2dFGRS is the largest completed redshift survey to date, including a total of 245,591 galaxies covering 1500 deg$^2$ to a median depth of z~0.11. We use the voidfinder algorithm to identify a total of 289 voids in the 2dFGRS with radius larger than 10h^{-1}Mpc. These voids have an average effective radius of 14.89+- 2.65 h^{-1}Mpc in the North Galactic Pole region (NGP) and 15.61+- 2.84h^{-1}Mpc in the South Galactic Pole region (SGP). These voids are extremely underdense, with average density contrast of delta rho/rho= -0.94+-0.02. The centers of voids are even emptier, because the few galaxies within the voids typically lie close to the edges. The total volume of the universe filled by these void regions is approximately 40%. We measure the Void Probability Function (VPF) of the 2dFGRS for volume-limited samples with limiting absolute magnitudes from -16 to -21 in b_J. We measure the Underdensity Probability Function for samples with limiting absolute magnitudes from -18 to -21. We find that the SGP is more underdense than the NGP for all but the brightest sample under consideration. There is good agreement between the VPF's of the Center for Astrophysics Survey and the 2dFGRS. Comparison of VPF's measured for the 2dFGRS with the distribution of simulated dark-matter halos of similar number density indicates that voids in the matter distribution in LCDM simulations are not empty enough. However, semi-analytic models of galaxy formation that include feedback effects yield VPF's that show excellent agreement with the data.

We present a novel method for simulation of the interior of large cosmic voids, suitable for study of the formation and evolution of objects lying within such regions. Following Birkhoff's theorem, void regions dynamically evolve as universes with cosmological parameters that depend on the underdensity of the void. We derive the values of $\Omega_M$, $\Omega_{\Lambda}$, and $H_0$ that describe this evolution. We examine how the growth rate of structure and scale factor in a void differ from the background universe. Together with a prescription for the power spectrum of fluctuations, these equations provide the initial conditions for running specialized void simulations. The increased efficiency of such simulations, in comparison with general-purpose simulations, allows an improvement of upwards of twenty in the mass resolution. As a sanity check, we run a moderate resolution simulation ($N=128^3$ particles) and confirm that the resulting mass function of void halos is consistent with other theoretical and numerical models.

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

(Abridged) The Sloan Digital Sky Survey now extends over a large enough region that galaxies in low density environments can be detected. Rojas et al. have developed a technique to extract more than 1000 void galaxies with delta rho/rho < -0.6 on scales r>7 h^{-1} Mpc from the SDSS. In this paper, we present the luminosity function of these galaxies and compare it to that of galaxies at higher densities. We obtain a measurement of the void galaxy luminosity function over the range -21.5<M_r*<-14.5 in the SDSS r-band. The void galaxy luminosity function (LF) is well fit by a Schechter function with values of Phi*=(0.19 +- 0.04)*0.01 h^3 Mpc^{-3}, M_r* - 5logh=-19.74 +- 0.11 and alpha=-1.18 +- 0.13, where as the LF of the wall galaxy has Schechter function parameters Phi*=(1.42 +- 0.03)*0.01 h^3 Mpc^{-3}, M_r* - 5logh=-20.62 +- 0.08 and alpha=-1.19 +- 0.07. The fainter value of M_r* is consistent with the work of Rojas et al.who found that void galaxies were fainter than galaxies at higher densities but the value of alpha is very similar to that of galaxies at higher densities. To investigate this further, the LF of wall and void galaxies as a function of density is considered. Only the wall galaxies in the highest density regions have a shallower faint end slope. The LF of the void galaxies is also compared to the LF of the wall galaxies with similar properties to the void galaxies. The LF of the bluest, Sersic index n<2 and EW(Halpha)>5A wall galaxies is very similar in shape to that of the void galaxies, with only slightly brighter M_r* values and similar values of alpha. The value of alpha ~ -1.2 found from the void galaxy LF is shallower that that found from CDM models, suggesting that voids are not filled with a large population of dwarf galaxies.

The mass function of clusters of galaxies is determined from 400 deg^2 of early commissioning imaging data of the Sloan Digital Sky Survey; ~300 clusters in the redshift range z = 0.1 - 0.2 are used. Clusters are selected using two independent selection methods: a Matched Filter and a red-sequence color magnitude technique. The two methods yield consistent results. The cluster mass function is compared with large-scale cosmological simulations. We find a best-fit cluster normalization relation of sigma_8*omega_m^0.6 = 0.33 +- 0.03 (for 0.1 ~< omega_m ~< 0.4), or equivalently sigma_8 = (0.16/omega_m)^0.6. The amplitude of this relation is significantly lower than the previous canonical value, implying that either omega_m is lower than previously expected (omega_m = 0.16 if sigma_8 = 1) or sigma_8 is lower than expected (sigma_8 = 0.7 if omega_m = 0.3). The best-fit mass function parameters are omega_m = 0.19 (+0.08,-0.07) and sigma_8 = 0.9 (+0.3,-0.2). High values of omega_m (>= 0.4) and low sigma_8 (=< 0.6) are excluded at >~ 2 sigma.

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.

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.

We present the first analysis of three-dimensional genus statistics for the SDSS EDR galaxy sample. Due to the complicated survey volume and the selection function, analytic predictions of the genus statistics for this sample are not feasible, therefore we construct extensive mock catalogs from N-body simulations in order to compare the observed data with model predictions. This comparison allows us to evaluate the effects of a variety of observational systematics on the estimated genus for the SDSS sample, including the shape of the survey volume, the redshift distortion effect, and the radial selection function due to the magnitude limit. The observed genus for the SDSS EDR galaxy sample is consistent with that predicted by simulations of a $\Lambda$-dominated spatially-flat cold dark matter model. Standard ($\Omega_0=1$) cold dark matter model predictions do not match the observations. We discuss how future SDSS galaxy samples will yield improved estimates of the genus.

We present predictions for numerous statistics related to the presence of voids in the distribution of galaxies in a cold dark matter model of structure formation using a semi-analytic model of galaxy formation. Our study is able to probe galaxies with masses as low as 10^9Msun/h corresponding to absolute magnitudes of M_bJ-5log(h)=-18.1 and M_r-5log(h)=-18.7. We quantify the void and underdense probability functions, distributions of nearest neighbour distances and void sizes and compute the density profiles of voids. These results are contrasted with the expectations for dark matter (and the difference examined in terms of the galaxy/dark matter biasing relation) and are compared to analytic predictions and observational data where available. The predicted void probability functions are consistent with those measured from the Center for Astrophysics redshift surveys given the rather large uncertainties in this relatively small (for studies of voids) observational sample. The size of the observational sample is too small to probe the bias between galaxies and dark matter that we predict. We also examine the predicted properties of galaxies living within voids and contrast these with the general galaxy population. Our predictions are aimed at forthcoming large galaxy redshift surveys which should for the first time provide statistically accurate measures of the void population.

We present the topology of a volume-limited sample of 11,884 galaxies, selected from an apparent-magnitude limited sample of over 100,000 galaxies observed as part of the Sloan Digital Sky Survey (SDSS). The data currently cover three main regions on the sky: one in the Galactic north and one in the south, both at zero degrees declination, and one area in the north at higher declination. Each of these areas covers a wide range of survey longitude but a narrow range of survey latitude, allowing the two dimensional genus to be measured. The genus curves of the SDSS sub-samples are similar, after appropriately normalizing these measurements for the different areas. We sum the genus curves from the three areas to obtain the total genus curve of the SDSS. The total curve has a shape similar to the genus curve derived from mock catalogs drawn from the Hubble Volume Lambda CDM simulation and is similar to that of a Gaussian random field. Likewise, comparison with the genus of the 2dFGRS, after normalization for the difference in area, reveals remarkable similarity in the topology of these samples. We test for the effects of galaxy type segregation by splitting the SDSS data into thirds, based on the u^*- r^* colors of the galaxies, and measure the genus of the reddest and bluest sub-samples. This red/blue split in u^*- r^* is essentially a split by morphology (Strateva et al. 2001). We find that the genus curve for the reddest galaxies exhibits a ``meatball'' shift of the topology -- reflecting the concentration of red galaxies in high density regions -- compared to the bluest galaxies and the full sample, in agreement with predictions from simulations.

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

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

We describe an algorithm to detect voids in galaxy redshift surveys. The method is based on the void finder algorithm of El-Ad & Piran. We apply a series of tests to determine how accurately we are able to recover the volumes of voids using our detection method. We simulate voids of different ellipticity and find that if voids are approximately spherical, our algorithm will recover 100% of the volume of the void. The more elliptical the void, the smaller the fraction of the volume we can recover. We insist that voids lie completely within the survey. Voids close to the edge of the survey will therefore be underestimated in volume. By considering a deeper sample, we estimate the maximal sphere diameters are correct to within 30%. We apply the algorithm to the Point Source Catalogue Survey (PSCz) and the Updated Zwicky Catalog (UZC). The PSCz survey is an almost all-sky survey with objects selected from the IRAS catalog. The UZC covers a smaller area of sky but is optically selected and samples the structures more densely. We detect 35 voids in the PSCz and 19 voids in the UZC with diameter larger than 20h^{-1}Mpc. Using this minimum size threshold, voids have an average effective diameter of 29.8+-3.5 h^{-1}Mpc (PSCz) and 29.2+-2.7 h^{-1}Mpc (UZC) and that they are underdense regions with delta rho /rho values of -0.92+-0.03 (PSCz) and -0.96+-0.01(UZC) respectively. Using this quite stringent threshold for void definition, voids fill up to 40% of the volume of the universe.

We study the topology of the publicly available data released by the 2dFGRS. The 2dFGRS data contains over 100,000 galaxy redshifts with a magnitude limit of b_J=19.45 and is the largest such survey to date. The data lie over a wide range of right ascension (75 degree strips) but only within a narrow range of declination (10 degree and 15 degree strips). This allows measurements of the two-dimensional genus to be made. The NGP displays a slight meatball shift topology, whereas the SGP displays a bubble like topology. The current SGP data also have a slightly higher genus amplitude. In both cases, a slight excess of overdense regions are found over underdense regions. We assess the significance of these features using mock catalogs drawn from the Virgo Consortium's Hubble Volume LCDM z=0 simulation. We find that differences between the NGP and SGP genus curves are only significant at the 1 sigma level. The average genus curve of the 2dFGRS agrees well with that extracted from the LCDM mock catalogs. We compare the amplitude of the 2dFGRS genus curve to the amplitude of a Gaussian random field with the same power spectrum as the 2dFGRS and find, contradictory to results for the 3D genus of other samples, that the amplitude of the GRF genus curve is slightly lower than that of the 2dFGRS. This could be due to a a feature in the current data set or the 2D genus may not be as sensitive as the 3D genus to non-linear clustering due to the averaging over the thickness of the slice in 2D. (Abridged)

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

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

One of the most striking features in galaxy redshift surveys is the ubiquitous presence of voids. However, voids have not been extensively studied due to observational limitations. Until recently, galaxy redshift surveys included only a few voids of diameter > 30 h-1Mpc. The first step in studying voids is the identification process. We outline a method for detecting voids (based on the method of El-Ad and Piran 1997, EP97). We apply it to the PSCz survey, and the Updated Zwicky Catalog. We find that voids have typical diameters of ~30 h-1Mpc and are very underdense regions with (delta rho)/rho ~ -0.95. Up to 40% of volume of the Universe is occupied by voids. We discuss the results from these surveys and our detection algorithm's usefulness for future surveys. The next generation of surveys (such as the Sloan Digital Sky Survey and 2dFGRS) will improve this situation. The surveys will extend to greater depths, allowing a more statistically complete sample of voids to be obtained and the accompanying digital imaging will provide accurate photometry of fainter objects. This will allow us to gain greater insight into the nature of voids, search for possible void galaxies and maybe place constraints on cosmological parameters as voids play a critical role in the evolution of large scale structure.

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

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

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

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

We use weighted mean and median statistics techniques to combine individual cosmic microwave background (CMB) anisotropy detections and determine binned, multipole-space, CMB anisotropy power spectra. The resultant power spectra are peaked. The derived weighted-mean CMB anisotropy power spectrum is not a good representation of the individual measurements in a number of multipole-space bins, if the CMB anisotropy is Gaussian and correlations between individual measurements are small. This could mean that some observational error bars are underestimated, possibly as a consequence of undetected systematic effects. Discarding the most discrepant 5% of the measurements alleviates but does not completely resolve this problem. The median-statistics power spectrum of this culled data set is not as constraining as the weighted-mean power spectrum. Nevertheless it indicates that there is more power at multipoles $\ell \sim 150-250$ than is expected in an open cold dark matter (CDM) model, and it is more consistent with a flat CDM model. Unlike the weighted-mean power spectrum, the median-statistics power spectrum at $\ell \sim 400-500$ does not exclude a second peak in the flat CDM model.

(Abridged) We develop median statistics that provide powerful alternatives to chi-squared likelihood methods and require fewer assumptions about the data. Applying median statistics to Huchra's compilation of nearly all estimates of the Hubble constant, we find a median value H_0=67 km/s/Mpc. Median statistics assume only that the measurements are independent and free of systematic errors. This estimate is arguably the best summary of current knowledge because it uses all available data and, unlike other estimates, makes no assumption about the distribution of measurement errors. The 95% range of purely statistical errors is +/- 2 km/s/Mpc. The statistical precision of this result leads us to analyze the range of possible systematic errors in the median, which we estimate to be roughly +/- 5 km/s/Mpc (95% limits), dominating over the statistical errors. A Bayesian median statistics treatment of high-redshift Type Ia supernovae (SNe Ia) apparent magnitude versus redshift data from Riess et al. yields a posterior probability that the cosmological constant Lambda > 0 of 70 or 89%, depending on the prior information used. The posterior probability of an open universe is about 47%. Analysis of the Perlmutter et al. high-redshift SNe Ia data show the best-fit flat-Lambda model favored over the best-fit Lambda = 0 open model by odds of 366:1; corresponding Riess et al. odds are 3:1 (assuming prior odds of 1:1).Median statistics analyses of the SNe Ia data do not rule out a time-variable Lambda model, and may even favor it over a time-independent Lambda and a Lambda = 0 open model.

The Sloan Digital Sky Survey has obtained 450 square degrees of photometric scan data, in five bands (u',g',r',i',z'), which we use to identify clusters of galaxies. We illustrate how we do star-galaxy separation, and present a simple and elegant method of detecting overdensities in the galaxy distribution, using the Voronoi Tessellation.

I review some results of estimation of the power spectrum of density fluctuations from galaxy redshift surveys and discuss advances that may be possible with the Sloan Digital Sky Survey. I then examine the realities of power spectrum estimation in the presence of Galactic extinction, photometric errors, galaxy evolution, clustering evolution, and uncertainty about the background cosmology.

Precision measurements of the galaxy power spectrum P(k) require a data analysis pipeline that is both fast enough to be computationally feasible and accurate enough to take full advantage of high-quality data. We present a rigorous discussion of different methods of power spectrum estimation, with emphasis on the traditional Fourier method, the linear (Karhunen-Loeve; KL), and quadratic data compression schemes, showing in what approximations they give the same result. To improve speed, we show how many of the advantages of KL data compression and power spectrum estimation may be achieved with a computationally faster quadratic method. To improve accuracy, we derive analytic expressions for handling the integral constraint, since it is crucial that finite volume effects are accurately corrected for on scales comparable to the depth of the survey. We also show that for the KL and quadratic techniques, multiple constraints can be included via simple matrix operations, thereby rendering the results less sensitive to galactic extinction and mis-estimates of the radial selection function. We present a data analysis pipeline that we argue does justice to the increases in both quality and quantity of data that upcoming redshift surveys will provide. It uses three analysis techniques in conjunction: a traditional Fourier approach on small scales, a pixelized quadratic matrix method on large scales and a pixelized KL eigenmode analysis to probe anisotropic effects such as redshift-space distortions.

(Abridged) Statistical analysis of the unresolved light in the Hubble Deep Field (HDF) strongly constrains possible sources of the optical Extragalactic Background Light (EBL). To test for the statistical signature of previously undetected sources, we estimate the auto, cross, and color correlations of the ``sky'' in the HDF that remains after masking objects brighter than I=30 mag. This measurement yields the most stringent limits to date on small-scale structure in the night sky; analysis of shallower imaging would be dominated by galaxies now detected by the HDF. Unless there is a truly uniform optical background, the mean EBL is likely to be within a small fraction of the surface brightness from detected galaxies. No currently plausible sources of additional EBL satisfy the constraints that they (1) would not have already been detected, (2) contribute EBL comparable to that from detected galaxies, and (3) do not produce EBL fluctuations in excess of the upper limits set by correlations in the HDF.

We describe a method for estimating the power spectrum of density fluctuations from galaxy redshift surveys that yields improvement in both accuracy and resolution over direct Fourier analysis. The key feature of this analysis is expansion of the observed density field in the unique set of statistically orthogonal spatial functions which obtains for a given survey's geometry and selection function and the known properties of galaxy clustering (the Karhunen-Loeve transform). Each of these eigenmodes of the observed density field optimally weights the data to yield the cleanest (highest signal/noise) possible measure of clustering power as a function of wavelength scale for any survey. Using Bayesian methods, we simultaneously estimate the mean density, power spectrum of density fluctuations, and redshift distortion parameters that best fit the observed data. This method is particularly important for analysis of surveys with small sky coverage, that are comprised of disjoint regions (e.g., an ensemble of pencil beams or slices), or that have large fluctuations in sampling density. We present algorithms for practical application of this technique to galaxy survey data.

The power spectrum of density fluctuations measured from galaxy redshift surveys provides important constraints on models for the formation of large-scale structure. I review current results for the 3-D power spectrum and examine the limitations of current measurements and estimation techniques. To span the decade of wavelength between the scales probed by galaxy surveys and COBE, measure the detailed shape of the power spectrum, and accurately examine the dependence of clustering on galaxy species, we require deeper samples with carefully controlled selection criteria and improved techniques for power spectrum estimation. I describe a new method for estimating the power spectrum that optimally treats survey data with arbitrary geometry and sampling.

The power spectrum of density fluctuations measured from galaxy redshift surveys provides important constraints on models for the formation of large-scale structure. I present new results for a redshift sample of 15,000 galaxies, and review the limitations of current measurements. To span the decade of wavelength between the scales probed by galaxy surveys and COBE, measure the detailed shape of the power spectrum, and accurately examine the dependence of clustering on galaxy species, we require deeper samples with carefully controlled selection criteria and improved techniques for power spectrum estimation. I review plans for the Sloan Digital Sky Survey and describe a new method for estimating the power spectrum that optimally treats survey data with arbitrary geometry and sampling.

We compute the power spectrum of galaxy density fluctuations in a recently completed redshift survey of optically-selected galaxies in the southern hemisphere (SSRS2). The amplitude and shape of the SSRS2 power spectrum are consistent with results of the Center for Astrophysics redshift survey of the northern hemisphere (CfA2), including the abrupt change of slope on a scale of 30-50Mpc/h; these results are reproducible for independent volumes of space and variations are consistent with the errors estimated from mock surveys. Taken together, the SSRS2 and CfA2 form a complete sample of 14,383 galaxies which covers one-third of the sky. The power spectrum of this larger sample continues to rise on scales up to ~ 200Mpc/h, with weak evidence for flattening on the largest scales. The SSRS2+CfA2 power spectrum and the power spectrum constraints implied by COBE are well-matched by an Omega*h ~ 0.2, Omega+lambda_0=1 CDM model with minimal biasing of optically-selected galaxies.