Halo models provide a simple and computationally inexpensive way to investigate the connection between galaxies and their dark matter haloes. However, these models rely on the assumption that the role of baryons can be easily parametrized in the modelling procedure. We aim to examine the ability of halo occupation distribution (HOD) modelling to reproduce the galaxy clustering found in two different hydrodynamic simulations, Illustris and EAGLE. For each simulation, we measure several galaxy clustering statistics on two different luminosity threshold samples. We then apply a simple five parameter HOD, which was fit to each simulation separately, to the corresponding dark matter only simulations, and measure the same clustering statistics. We find that the halo mass function is shifted to lower masses in the hydrodynamic simulations, resulting in a galaxy number density that is too high when an HOD is applied to the dark matter only simulation. However, the exact way in which baryons alter the mass function is remarkably different in the two simulations. After applying a correction to the halo mass function in each simulation, the HOD is able to accurately reproduce all clustering statistics for the high luminosity sample of galaxies. For the low luminosity sample, we find evidence that in addition to correcting the halo mass function, including spatial, velocity, and assembly bias parameters in the HOD is necessary to accurately reproduce clustering statistics.

We study the present-day connection between galaxy morphology and angular momentum using the {\sc Dark Sage} semi-analytic model of galaxy formation. For galaxies between $ 10^{11}-10^{12} \mathrm{M}_{\odot}$ in stellar mass, the model successfully predicts the observed trend whereby galaxies with more prominent disks exhibit higher {\em stellar} disk specific angular momentum ($j_{\rm stellar, disk}$) at fixed stellar mass. However, when we include the gas in the disk, bulge-dominated galaxies have the highest {\em total} disk specific angular momentum ($j_{\rm total, disk}$). We attribute this to a large contribution from an extended disk of cold gas in typical bulge-dominated galaxies. We find the relationship between $j_{\rm dark matter}$ and morphology to be quite complex. Surprisingly, in this stellar mass range, not only do bulge-dominated galaxies tend to live in halos with higher $j_{\rm dark matter}$ than disk-dominated galaxies, but intermediate galaxies (those with roughly equal fractions of bulge and disk mass) have the lowest $j_{\rm dark matter}$ of all. Yet, when controlling for halo mass, rather than stellar mass, the relationship between $j_{\rm dark matter}$ and morphology vanishes. Based on these results, halo mass rather than angular momentum is the main driver of the predicted morphology sequence at high masses. In fact, in our stellar mass range, disk-dominated galaxies live in dark matter halos that are roughly 1/10th the mass of their bulge-dominated counterparts.

We present a machine learning (ML) approach for the prediction of galaxies' dark matter halo masses that achieves an improved performance over conventional methods. We train three ML algorithms (\texttt{XGBoost}, Random Forests, and neural network) to predict halo masses using a set of synthetic galaxy catalogues that are built by populating dark matter haloes in N-body simulations with galaxies, and that match both the clustering and the joint-distributions of properties of galaxies in the Sloan Digital Sky Survey (SDSS). We explore the correlation of different galaxy- and group-related properties with halo mass, and extract the set of nine features that contribute the most to the prediction of halo mass. We find that mass predictions from the ML algorithms are more accurate than those from halo abundance matching (\texttt{HAM}) or dynamical mass (\texttt{DYN}) estimates. Since the danger of this approach is that our training data might not accurately represent the real Universe, we explore the effect of testing the model on synthetic catalogues built with different assumptions than the ones used in the training phase. We test a variety of models with different ways of populating dark matter haloes, such as adding velocity bias for satellite galaxies. We determine that, though training and testing on different data can lead to systematic errors in predicted masses, the ML approach still yields substantially better masses than either \texttt{HAM}or \texttt{DYN}. Finally, we apply the trained model to a galaxy and group catalogue from the SDSS DR7 and present the resulting halo masses.

We investigate the role of angular momentum in the clustering of dark matter haloes. We make use of data from two high-resolution N-body simulations spanning over four orders of magnitude in halo mass, from $10^{9.8}$ to $10^{14}\ h^{-1}\ \text{M}_\odot$. We explore the hypothesis that mass accretion in filamentary environments alters the angular momentum of a halo, thereby driving a correlation between the spin parameter $\lambda$ and the strength of clustering. However, we do not find evidence that the distribution of matter on large scales is related to the spin of haloes. We find that a halo's spin is correlated with its age, concentration, sphericity, and mass accretion rate. Removing these correlations strongly affects the strength of secondary spin bias at low halo masses. We also find that high spin haloes are slightly more likely to be found near another halo of comparable mass. These haloes that are found near a comparable mass neighbour - a \textit{twin} - are strongly spatially biased. We demonstrate that this \textit{twin bias}, along with the relationship between spin and mass accretion rates, statistically accounts for halo spin secondary bias.

A small survey of the UV-absorbing gas in 12 low-$z$ galaxy groups has been conducted using the Cosmic Origins Spectrograph (COS) on-board the Hubble Space Telescope (HST). Targets were selected from a large, homogeneously-selected sample of groups found in the Sloan Digital Sky Survey (SDSS). A critical selection criterion excluded sight lines that pass close ($<1.5$ virial radii) to a group galaxy, to ensure absorber association with the group as a whole. Deeper galaxy redshift observations are used both to search for closer galaxies and also to characterize these $10^{13.5}$ to $10^{14.5} M_{\odot}$ groups, the most massive of which are highly-virialized with numerous early-type galaxies (ETGs). This sample also includes two spiral-rich groups, not yet fully-virialized. At group-centric impact parameters of 0.3-2 Mpc, these $\mathrm{S/N}=15$-30 spectra detected HI absorption in 7 of 12 groups; high (OVI) and low (SiIII) ion metal lines are present in 2/3 of the absorption components. None of the three most highly-virialized, ETG-dominated groups are detected in absorption. Covering fractions $\gtrsim50$% are seen at all impact parameters probed, but do not require large filling factors despite an enormous extent. Unlike halo clouds in individual galaxies, group absorbers have radial velocities which are too low to escape the group potential well without doubt. This suggests that these groups are "closed boxes" for galactic evolution in the current epoch. Evidence is presented that the cool and warm group absorbers are not a pervasive intra-group medium (IGrM), requiring a hotter ($T\sim10^6$ to $10^7$ K) IGrM to be present to close the baryon accounting.

Galactic conformity is the phenomenon whereby galaxy properties exhibit excess correlations across distance than that expected if these properties only depended on halo mass. We perform a comprehensive study of conformity at low redshift using a galaxy group catalogue from the SDSS DR7 spectroscopic sample. We study correlations both between central galaxies and their satellites (1-halo), and between central galaxies in separate haloes (2-halo). We use the quenched fractions and the marked correlation function (MCF), to probe for conformity in three galaxy properties, $(g-r)$ colour, specific star formation rate (sSFR), and morphology. We assess the statistical significance of conformity signals with a suite of mock galaxy catalogues that have no built-in conformity, but contain the same group-finding and mass assignment errors as the real data. In the case of 1-halo conformity, quenched fractions show strong signals at all group masses. However, these signals are equally strong in mock catalogues, indicating that the conformity signal is spurious and likely entirely caused by group-finding systematics, calling into question previous claims of 1-halo conformity detection. The MCF reveals a significant detection of radial segregation within massive groups, but no evidence of conformity. In the case of 2-halo conformity, quenched fractions show no significant evidence of conformity in colour or sSFR once compared with mock catalogues, but a clear signal using morphology. In contrast, the MCF reveals a small, yet highly significant signal for all three properties in low mass groups and scales of $0.8-4\ h^{-1}\textrm{Mpc}$, possibly representing the first robust detection of 2-halo conformity.

Interpreting the small-scale clustering of galaxies with halo models can elucidate the connection between galaxies and dark matter halos. Unfortunately, the modelling is typically not sufficiently accurate for ruling out models statistically. It is thus difficult to use the information encoded in small scales to test cosmological models or probe subtle features of the galaxy-halo connection. In this paper, we attempt to push halo modelling into the "accurate" regime with a fully numerical mock-based methodology and careful treatment of statistical and systematic errors. With our forward-modelling approach, we can incorporate clustering statistics beyond the traditional two-point statistics. We use this modelling methodology to test the standard $\Lambda\mathrm{CDM}$ + halo model against the clustering of SDSS DR7 galaxies. Specifically, we use the projected correlation function, group multiplicity function and galaxy number density as constraints. We find that while the model fits each statistic separately, it struggles to fit them simultaneously. Adding group statistics leads to a more stringent test of the model and significantly tighter constraints on model parameters. We explore the impact of varying the adopted halo definition and cosmological model and find that changing the cosmology makes a significant difference. The most successful model we tried (Planck cosmology with Mvir halos) matches the clustering of low luminosity galaxies, but exhibits a 2.3$\sigma$ tension with the clustering of luminous galaxies, thus providing evidence that the "standard" halo model needs to be extended. This work opens the door to adding interesting freedom to the halo model and including additional clustering statistics as constraints.

Standard present day large-scale structure (LSS) analyses make a major assumption in their Bayesian parameter inference --- that the likelihood has a Gaussian form. For summary statistics currently used in LSS, this assumption, even if the underlying density field is Gaussian, cannot be correct in detail. We investigate the impact of this assumption on two recent LSS analyses: the Beutler et al. (2017) power spectrum multipole ($P_\ell$) analysis and the Sinha et al. (2017) group multiplicity function ($\zeta$) analysis. Using non-parametric divergence estimators on mock catalogs originally constructed for covariance matrix estimation, we identify significant non-Gaussianity in both the $P_\ell$ and $\zeta$ likelihoods. We then use Gaussian mixture density estimation and Independent Component Analysis on the same mocks to construct likelihood estimates that approximate the true likelihood better than the Gaussian $pseudo$-likelihood. Using these likelihood estimates, we accurately estimate the true posterior probability distribution of the Beutler et al. (2017) and Sinha et al. (2017) parameters. Likelihood non-Gaussianity shifts the $f\sigma_8$ constraint by $-0.44\sigma$, but otherwise, does not significantly impact the overall parameter constraints of Beutler et al. (2017). For the $\zeta$ analysis, using the pseudo-likelihood significantly underestimates the uncertainties and biases the constraints of Sinha et al. (2017) halo occupation parameters. For $\log M_1$ and $\alpha$, the posteriors are shifted by $+0.43\sigma$ and $-0.51\sigma$ and broadened by $42\%$ and $66\%$, respectively. The divergence and likelihood estimation methods we present provide a straightforward framework for quantifying the impact of likelihood non-Gaussianity and deriving more accurate parameter constraints.

We explore the phenomenon commonly known as halo assembly bias, whereby dark matter halos of the same mass are found to be more or less clustered when a second halo property is considered, for halos in the mass range $3.7 \times 10^{11} \; h^{-1} \mathrm{M_{\odot}} - 5.0 \times 10^{13} \; h^{-1} \mathrm{M_{\odot}}$. Using the Large Suite of Dark Matter Simulations (LasDamas) we consider nine commonly used halo properties and find that a clustering bias exists if halos are binned by mass or by any other halo property. This secondary bias implies that no single halo property encompasses all the spatial clustering information of the halo population. The mean values of some halo properties depend on their halo's distance to a more massive neighbor. Halo samples selected by having high values of one of these properties therefore inherit a neighbor bias such that they are much more likely to be close to a much more massive neighbor. This neighbor bias largely accounts for the secondary bias seen in halos binned by mass and split by concentration or age. However, halos binned by other mass-like properties still show a secondary bias even when the neighbor bias is removed. The secondary bias of halos selected by their spin behaves differently than that for other halo properties, suggesting that the origin of the spin bias is different than of other secondary biases.

We examine the z = 0 group-integrated stellar and cold baryonic (stars + cold atomic gas) mass functions (group SMF and CBMF) and the baryonic collapse efficiency (group cold baryonic to dark matter halo mass ratio) using the RESOLVE and ECO survey galaxy group catalogs and a galform semi-analytic model (SAM) mock catalog. The group SMF and CBMF fall off more steeply at high masses and rise with a shallower low-mass slope than the theoretical halo mass function (HMF). The transition occurs at group-integrated cold baryonic mass M_coldbary ~ 10^11 Msun. The SAM, however, has significantly fewer groups at the transition mass ~ 10^11 Msun and a steeper low-mass slope than the data, suggesting that feedback is too weak in low-mass halos and conversely too strong near the transition mass. Using literature prescriptions to include hot halo gas and potential unobservable galaxy gas produces a group BMF with slope similar to the HMF even below the transition mass. Its normalization is lower by a factor of ~2, in agreement with estimates of warm-hot gas making up the remaining difference. We compute baryonic collapse efficiency with the halo mass calculated two ways, via halo abundance matching (HAM) and via dynamics (extended all the way to three-galaxy groups using stacking). Using HAM, we find that baryonic collapse efficiencies reach a flat maximum for groups across the halo mass range of M_halo ~ 10^11.4-12 Msun, which we label "nascent groups." Using dynamics, however, we find greater scatter in baryonic collapse efficiencies, likely indicating variation in group hot-to-cold baryon ratios. Similarly, we see higher scatter in baryonic collapse efficiencies in the SAM when using its true groups and their group halo masses as opposed to friends-of-friends groups and HAM masses.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

In preparation for a Hubble Space Telescope (HST) observing project using the Cosmic Origins Spectrograph (COS), the positions of all AGN targets having high-S/N far-UV G130M spectra were cross-correlated with a large catalog of low-redshift galaxy groups homogeneously selected from the spectroscopic sample of the Sloan Digital Sky Survey (SDSS). Searching for targets behind only those groups at z = 0.1-0.2 (which places the OVI doublet in the wavelength region of peak COS sensitivity) we identified only one potential S/N = 15-20 target, FBQS 1010+3003. An OVI-only absorber was found in its G130M spectrum at z = 0.11326, close to the redshift of a foreground small group of luminous galaxies at z = 0.11685. Because there is no associated Lyalpha absorption, any characterization of this absorber is necessarily minimal; however, the OVI detection likely traces "warm" gas in collisional ionization equilibrium at T ~ 300,000 K. While this discovery is consistent with being interface gas between cooler, photoionized clouds and a hotter intra-group medium, it could also be warm, interface gas associated with the circum-galactic medium (CGM) of the single closest galaxy. In this case a detailed analysis of the galaxy distribution (complete to 0.2 L*) strongly favors the individual galaxy association. This analysis highlights the necessity of both high-S/N > 20 COS data and a deep galaxy redshift survey of the region in order to test more rigorously the association of OVI-absorbing gas with a galaxy group. A Cycle 23 HST/COS program currently is targeting 10 UV-bright AGN behind 12 low-redshift galaxy groups to test the warm, group gas hypothesis.

We have searched for the signature of cosmic voids in the CMB, in both the Planck temperature and lensing-convergence maps; voids should give decrements in both. We use zobov voids from the DR12 SDSS CMASS galaxy sample. We base our analysis on N-body simulations, to avoid a posteriori bias. For the first time, we detect the signature of voids in CMB lensing: the significance is $3.2\sigma$, close to $\Lambda$CDM in both amplitude and projected density-profile shape. A temperature dip is also seen, at modest significance ($2.3\sigma$), with amplitude about 6 times the prediction. This temperature signal is induced mostly by voids with radius between 100 and 150 Mpc/h, while the lensing signal is mostly contributed by smaller voids -- as expected; lensing relates directly to density, while ISW depends on gravitational potential. The void abundance in observations and simulations agree, as well. We also repeated the analysis excluding lower-significance voids: no lensing signal is detected, with an upper limit of about twice the $\Lambda$CDM prediction. But the mean temperature decrement now becomes non-zero at the $3.7\sigma$ level (similar to that found by Granett et al.), with amplitude about 20 times the prediction. However, the observed dependence of temperature on void size is in poor agreement with simulations, whereas the lensing results are consistent with $\Lambda$CDM theory. Thus, the overall tension between theory and observations does not favour non-standard theories of gravity, despite the hints of an enhanced amplitude for the ISW effect from voids.

We present the HI mass inventory for the RESOLVE survey, a volume-limited, multi-wavelength census of >1500 z=0 galaxies spanning diverse environments and complete in baryonic mass down to dwarfs of 10^9 Msun. This first 21cm data release provides robust detections or strong upper limits (1.4M_HI < 5 to 10% of stellar mass M_stars) for 94% of RESOLVE. We examine global atomic gas-to-stellar mass ratios (G/S) in relation to galaxy environment using several metrics: group dark matter halo mass M_h , central/satellite designation, relative mass density of the cosmic web, and distance to nearest massive group. We find that at fixed M_stars, satellites have decreasing G/S with increasing M_h starting clearly at M_h = 10^12 Msun, suggesting the presence of starvation and/or stripping mechanisms associated with halo gas heating in intermediate-mass groups. The analogous relationship for centrals is uncertain because halo abundance matching builds in relationships between central G/S, stellar mass, and halo mass, which depend on the integrated group property used as a proxy for halo mass (stellar or baryonic mass). On larger scales G/S trends are less sensitive to the abundance matching method. At fixed M_h < 10^12 Msun, the fraction of gas-poor centrals increases with large-scale structure density. In overdense regions, we identify a rare population of gas-poor centrals in low-mass (M_h < 10^11.4 Msun) halos primarily located within 1.5 times the virial radius of more massive (M_h > 10^12 Msun) halos, suggesting that gas stripping and/or starvation may be induced by interactions with larger halos or the surrounding cosmic web. We find that the detailed relationship between G/S and environment varies when we examine different subvolumes of RESOLVE independently, which we suggest may be a signature of assembly bias.

We apply the Alcock-Paczynski (AP) test to the stacked voids identified using the large-scale structure galaxy catalog from the Baryon Oscillation Spectroscopic Survey (BOSS). This galaxy catalog is part of the Sloan Digital Sky Survey (SDSS) Data Release 12 and is the final catalog of SDSS-III. We also use 1000 mock galaxy catalogs that match the geometry, density, and clustering properties of the BOSS sample in order to characterize the statistical uncertainties of our measurements and take into account systematic errors such as redshift space distortions. For both BOSS data and mock catalogs, we use the ZOBOV algorithm to identify voids, we stack together all voids with effective radii of 30-100Mpc/h in the redshift range 0.43-0.7, and we accurately measure the shape of the stacked voids. Our tests with the mock catalogs show that we measure the stacked void ellipticity with a statistical precision of 2.6%. We find that the stacked voids in redshift space are slightly squashed along the line of sight, which is consistent with previous studies. We repeat this measurement of stacked void shape in the BOSS data assuming several values of Omega_m within the flat LCDM model, and we compare to the mock catalogs in redshift space in order to perform the AP test. We obtain a constraint of $\Omega_m = 0.38^{+0.18}_{-0.15}$ at the 68% confidence level from the AP test. We discuss the various sources of statistical and systematic noise that affect the constraining power of this method. In particular, we find that the measured ellipticity of stacked voids scales more weakly with cosmology than the standard AP prediction, leading to significantly weaker constraints. We discuss how AP constraints will improve in future surveys with larger volumes and densities.

In this work, we present galaxy stellar and baryonic (stars plus cold gas) mass functions (SMF and BMF) and their halo mass dependence for two volume-limited data sets. The first, RESOLVE-B, coincides with the Stripe 82 footprint and is extremely complete down to baryonic mass Mbary ~ 10^9.1 Msun, probing the gas-rich dwarf regime below Mbary ~ 10^10 Msun. The second, ECO, covers a ~40 times larger volume (containing RESOLVE-A) and is complete to Mbary ~ 10^9.4 Msun. To construct the SMF and BMF we implement a new "cross-bin sampling" technique with Monte Carlo sampling from the full likelihood distributions of stellar or baryonic mass. Our SMFs exhibit the "plateau" feature starting below Mstar ~ 10^10 Msun that has been described in prior work. However, the BMF fills in this feature and rises as a straight power law below ~10^10 Msun, as gas-dominated galaxies become the majority of the population. Nonetheless, the low-mass slope of the BMF is not as steep as that of the theoretical dark matter halo MF. Moreover, we assign group halo masses by abundance matching, finding that the SMF and BMF separated into four physically motivated halo mass regimes reveal complex structure underlying the simple shape of the overall MFs. In particular, the satellite MFs are depressed below the central galaxy MF "humps" in groups with mass <10^13.5 Msun yet rise steeply in clusters. Our results suggest that satellite destruction and/or stripping are active from the point of nascent group formation. We show that the key role of groups in shaping MFs enables reconstruction of a given survey's SMF or BMF based on its group halo mass distribution.

We present a cosmic void catalog using the large-scale structure galaxy catalog from the Baryon Oscillation Spectroscopic Survey (BOSS). This galaxy catalog is part of the Sloan Digital Sky Survey (SDSS) Data Release 12 and is the final catalog of SDSS-III. We take into account the survey boundaries, masks, and angular and radial selection functions, and apply the ZOBOV void finding algorithm to the galaxy catalog. We identify a total of 10,643 voids. After making quality cuts to ensure that the voids represent real underdense regions, we obtain 1,228 voids with effective radii spanning the range 20-100Mpc/h and with central densities that are, on average, 30% of the mean sample density. We release versions of the catalogs both with and without quality cuts. We discuss the basic statistics of voids, such as their size and redshift distributions, and measure the radial density profile of the voids via a stacking technique. In addition, we construct mock void catalogs from 1000 mock galaxy catalogs, and find that the properties of BOSS voids are in good agreement with those in the mock catalogs. We compare the stellar mass distribution of galaxies living inside and outside of the voids, and find no significant difference. These BOSS and mock void catalogs are useful for a number of cosmological and galaxy environment studies.

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

We study the relationships between galaxy environments and galaxy properties related to disk (re)growth, considering two highly complete samples that are approximately baryonic mass limited into the high-mass dwarf galaxy regime, the Environmental COntext (ECO) catalog (data release herein) and the B-semester region of the REsolved Spectroscopy Of a Local VolumE (RESOLVE) survey. We quantify galaxy environments using both group identification and smoothed galaxy density field methods. We use by-eye and quantitative morphological classifications plus atomic gas content measurements and estimates. We find that blue early-type (E/S0) galaxies, gas-dominated galaxies, and UV-bright disk host galaxies all become distinctly more common below group halo mass ~10^11.5 Msun, implying that this low group halo mass regime may be a preferred regime for significant disk growth activity. We also find that blue early-type and blue late-type galaxies inhabit environments of similar group halo mass at fixed baryonic mass, consistent with a scenario in which blue early types can regrow late-type disks. In fact, we find that the only significant difference in the typical group halo mass inhabited by different galaxy classes is for satellite galaxies with different colors, where at fixed baryonic mass red early and late types have higher typical group halo masses than blue early and late types. More generally, we argue that the traditional morphology-environment relation (i.e., that denser environments tend to have more early types) can be largely attributed to the morphology-galaxy mass relation for centrals and the color-environment relation for satellites.

We measure the two-point correlation function of G-dwarf stars within 1-3 kpc of the Sun in multiple lines-of-sight using the Schlesinger et al. G-dwarf sample from the SDSS SEGUE survey. The shapes of the correlation functions along individual SEGUE lines-of-sight depend sensitively on both the stellar-density gradients and the survey geometry. We fit smooth disk galaxy models to our SEGUE clustering measurements, and obtain strong constraints on the thin- and thick-disk components of the Milky Way. Specifically, we constrain the values of the thin- and thick-disk scale heights with 3% and 2% precision, respectively, and the values of the thin- and thick-disk scale lengths with 20% and 8% precision, respectively. Moreover, we find that a two-disk model is unable to fully explain our clustering measurements, which exhibit an excess of clustering at small scales (< 50 pc). This suggests the presence of small-scale substructure in the disk system of the Milky Way.

We measure the angular clustering of galaxies from the Sloan Digital Sky Survey Data Release 7 in order to probe the spatial distribution of satellite galaxies within their dark matter halos. Specifically, we measure the angular correlation function on very small scales (7-320") in a range of luminosity threshold samples (absolute r-band magnitudes of -18 up to -21) that are constructed from the subset of SDSS that has been spectroscopically observed more than once (the so-called plate overlap region). We choose to measure angular clustering in this reduced survey footprint in order to minimize the effects of fiber collision incompleteness, which are otherwise substantial on these small scales. We model our clustering measurements using a fully numerical halo model that populates dark matter halos in N-body simulations to create realistic mock galaxy catalogs. The model has free parameters that specify both the number and spatial distribution of galaxies within their host halos. We adopt a flexible density profile for the spatial distribution of satellite galaxies that is similar to the dark matter Navarro-Frenk-White (NFW) profile, except that the inner slope is allowed to vary. We find that the angular clustering of our most luminous samples (Mr< -20 and -21) suggests that luminous satellite galaxies have substantially steeper inner density profiles than NFW. Lower luminosity samples are less constraining, however, and are consistent with satellite galaxies having shallow density profiles. Our results confirm the findings of Watson et al. 2012 while using different clustering measurements and modeling methodology.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

We investigate the ability of current implementations of galaxy group finders to recover colour-dependent halo occupation statistics. To test the fidelity of group catalogue inferred statistics, we run three different group finders used in the literature over a mock that includes galaxy colours in a realistic manner. Overall, the resulting mock group catalogues are remarkably similar, and most colour-dependent statistics are recovered with reasonable accuracy. However, it is also clear that certain systematic errors arise as a consequence of correlated errors in group membership determination, central/satellite designation, and halo mass assignment. We introduce a new statistic, the halo transition probability (HTP), which captures the combined impact of all these errors. As a rule of thumb, errors tend to equalize the properties of distinct galaxy populations (i.e. red vs. blue galaxies or centrals vs. satellites), and to result in inferred occupation statistics that are more accurate for red galaxies than for blue galaxies. A statistic that is particularly poorly recovered from the group catalogues is the red fraction of central galaxies as function of halo mass. Group finders do a good job in recovering galactic conformity, but also have a tendency to introduce weak conformity when none is present. We conclude that proper inference of colour-dependent statistics from group catalogues is best achieved using forward modelling (i.e., running group finders over mock data), or by implementing a correction scheme based on the HTP, as long as the latter is not too strongly model-dependent.

In this paper, we test the age matching hypothesis that the star formation rate (SFR) of a galaxy of fixed stellar mass is determined by its dark matter halo formation history, and as such, that more quiescent galaxies reside in older halos. This simple model has been remarkably successful at predicting color-based galaxy statistics at low redshift as measured in the Sloan Digital Sky Survey (SDSS). To further test this method with observations, we present new SDSS measurements of the galaxy two-point correlation function and galaxy-galaxy lensing as a function of stellar mass and SFR, separated into quenched and star-forming galaxy samples. We find that our age matching model is in excellent agreement with these new measurements. We also employ a galaxy group finder and show that our model is able to predict: (1) the relative SFRs of central and satellite galaxies, (2) the SFR-dependence of the radial distribution of satellite galaxy populations within galaxy groups, rich groups, and clusters and their surrounding larger scale environments, and (3) the interesting feature that the satellite quenched fraction as a function of projected radial distance from the central galaxy exhibits an ~ r^-.15 slope, independent of environment. The accurate prediction for the spatial distribution of satellites is intriguing given the fact that we do not explicitly model satellite-specific processes after infall, and that in our model the virial radius does not mark a special transition region in the evolution of a satellite, contrary to most galaxy evolution models. The success of the model suggests that present-day galaxy SFR is strongly correlated with halo mass assembly history.

We investigate the luminosity and colour dependence of clustering of CMASS galaxies in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Tenth Data Release. The halo occupation distribution framework is adopted to model the projected two-point correlation function measurements on small and intermediate scales (from $0.02$ to $60\,h^{-1}{\rm {Mpc}}$) and to interpret the observed trends and infer the connection of galaxies to dark matter halos. We find that luminous red galaxies reside in massive halos of mass $M{\sim}10^{13}$--$10^{14}\,h^{-1}{\rm M_\odot}$ and more luminous galaxies are more clustered and hosted by more massive halos. The strong small-scale clustering requires a fraction of these galaxies to be satellites in massive halos, with the fraction at the level of 5--8 per cent and decreasing with luminosity. The characteristic mass of a halo hosting on average one satellite galaxy above a luminosity threshold is about a factor $8.7$ larger than that of a halo hosting a central galaxy above the same threshold. At a fixed luminosity, progressively redder galaxies are more strongly clustered on small scales, which can be explained by having a larger fraction of these galaxies in the form of satellites in massive halos. Our clustering measurements on scales below $0.4\,h^{-1}{\rm {Mpc}}$ allow us to study the small-scale spatial distribution of satellites inside halos. While the clustering of luminosity-threshold samples can be well described by a Navarro-Frenk-White (NFW) profile, that of the reddest galaxies prefers a steeper or more concentrated profile. Finally, we also use galaxy samples of constant number density at different redshifts to study the evolution of luminous galaxies, and find the clustering to be consistent with passive evolution in the redshift range of $0.5 \lesssim z \lesssim 0.6$.

We use cosmological N-body simulations to investigate whether measurements of the moments of large-scale structure can yield constraints on primordial non-Gaussianity. We measure the variance, skewness, and kurtosis of the evolved density field from simulations with Gaussian and three different non-Gaussian initial conditions: a local model with f_NL=100, an equilateral model with f_NL=-400, and an orthogonal model with f_NL=-400. We show that the moments of the dark matter density field differ significantly between Gaussian and non-Gaussian models. We also make the measurements on mock galaxy catalogs that contain galaxies with clustering properties similar to those of luminous red galaxies (LRGs). We find that, in the case of skewness and kurtosis, galaxy bias reduces the detectability of non-Gaussianity, though we can still clearly discriminate between different models in our simulation volume. However, in the case of the variance, galaxy bias greatly amplifies the detectability of non-Gaussianity. In all cases we find that redshift distortions do not significantly affect the detectability. When we restrict our measurements to volumes equivalent to the Sloan Digital Sky Survey II (SDSS-II) or Baryon Oscillation Spectroscopic Survey (BOSS) samples, the probability of detecting a departure from the Gaussian model is high by using measurements of the variance, but very low by using only skewness and kurtosis measurements. We find that skewness and kurtosis measurements are never likely to yield useful constraints on primordial non-Gaussianity, but future surveys should be large enough to place meaningful constraints using measurements of the galaxy variance.(Abridged)

The two main advantages of space-based observation of extreme-energy ($\gtrsim 10^{19}$~eV) cosmic-rays (EECRs) over ground-based observatories are the increased field of view, and the all-sky coverage with nearly uniform systematics of an orbiting observatory. The former guarantees increased statistics, whereas the latter enables a partitioning of the sky into spherical harmonics. We have begun an investigation, using the spherical harmonic technique, of the reach of \J\ into potential anisotropies in the extreme-energy cosmic-ray sky-map. The technique is explained here, and simulations are presented. The discovery of anisotropies would help to identify the long-sought origin of EECRs.

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

The age matching model has recently been shown to predict correctly the luminosity L and g-r color of galaxies residing within dark matter halos. The central tenet of the model is intuitive: older halos tend to host galaxies with older stellar populations. In this paper, we demonstrate that age matching also correctly predicts the g-r color trends exhibited in a wide variety of statistics of the galaxy distribution for stellar mass M* threshold samples. In particular, we present new measurements of the galaxy two-point correlation function and the galaxy-galaxy lensing signal as a function of M* and g-r color from the Sloan Digital Sky Survey, and show that age matching exhibits remarkable agreement with these and other statistics of low-redshift galaxies. In so doing, we also demonstrate good agreement between the galaxy-galaxy lensing observed by SDSS and the signal predicted by abundance matching, a new success of this model. We describe how age matching is a specific example of a larger class of Conditional Abundance Matching models (CAM), a theoretical framework we introduce here for the first time. CAM provides a general formalism to study correlations at fixed mass between any galaxy property and any halo property. The striking success of our simple implementation of CAM provides compelling evidence that this technique has the potential to describe the same set of data as alternative models, but with a dramatic reduction in the required number of parameters. CAM achieves this reduction by exploiting the capability of contemporary N-body simulations to determine dark matter halo properties other than mass alone, which distinguishes our model from conventional approaches to the galaxy-halo connection.

This work investigates the preferred environment of local Luminous IR Galaxies (LIRGs) in terms of the host halos that they inhabit, and in comparison to a control galaxy sample. The LIRGs are drawn from the IRAS Point Source Catalogue redshift survey (PSCz), while the control sample is drawn from the 2MASS redshift survey (2MRS). A friends-of-friends algorithm was run on the 2MRS sample to identify galaxies living in the same dark matter halos and the PSCz galaxies were then associated with these identified halos. We show that the relative probability of finding local LIRGs with respect to 2MASS galaxies is largest in approximately group size halos (Mhalo~10^13M_sun), and declines both in the cluster regime and in smaller halos. This confirms, using a different technique than in previous work, that local LIRGs are indeed more abundant in group environments than elsewhere. We also find a trend between the LIR values of LIRGs and their location within their host dark matter halos, such that the average location of LIRGs with high IR luminosity is closer to the halo centre than for low IR luminosity galaxies. Moreover, this trend does not seem to depend on halo mass.

We construct mock catalogs of galaxy groups using subhalo abundance matching (SHAM) and undertake several new tests of the SHAM prescription for the galaxy-dark matter connection. All SHAM models we studied exhibit significant tension with galaxy groups observed in the Sloan Digital Sky Survey (SDSS). The SHAM prediction for the field galaxy luminosity function (LF) is systematically too dim, and the group galaxy LF systematically too bright, regardless of the details of the SHAM prescription. SHAM models connecting r-band luminosity, Mr, to Vacc, the maximum circular velocity of a subhalo at the time of accretion onto the host, faithfully reproduce galaxy group abundance as a function of richness, g(N). However, SHAM models connecting Mr with Vpeak, the peak value of Vmax over the entire merger history of the halo, over-predict galaxy group abundance. Our results suggest that no SHAM model can simultaneously reproduce the observed g(N) and two-point projected galaxy clustering. Nevertheless, we also report a new success of SHAM: an accurate prediction for Phi(m12), the abundance of galaxy groups as a function of magnitude gap m12, defined as the difference between the r-band absolute magnitude of the two brightest group members. We show that it may be possible to use joint measurements of g(N) and Phi(m12) to tightly constrain the details of the SHAM implementation. Additionally, we show that the hypothesis that the luminosity gap is constructed via random draws from a universal LF provides a poor description of the data, contradicting recent claims in the literature. Finally, we test a common assumption of the Conditional Luminosity Function (CLF) formalism, that the satellite LF need only be conditioned by the brightness of the central galaxy. We find this assumption to be well-supported by the observed Phi(m12).

We introduce a new technique that uses galaxy clustering to constrain how satellite galaxies lose stellar mass and contribute to the diffuse "intrahalo light" (IHL). We implement two models that relate satellite galaxy stellar mass loss to the detailed knowledge of subhalo dark matter mass loss. Model 1 assumes that the fractional stellar mass loss of a galaxy is proportional to the fractional amount of dark matter mass loss of its subhalo. Model 2 accounts for a delay in the time that stellar mass is lost since the galaxy resides deep in the potential well of the subhalo which may experience dark matter mass loss for some time before the galaxy is affected. We use these models to predict the stellar masses of a population of galaxies and use abundance matching to predict the clustering of several r-band luminosity threshold samples from the Sloan Digital Sky Survey. Abundance matching assuming no stellar mass loss (akin to abundance matching at the time of subhalo infall) over-estimates the correlation function on small scales (<~ 1 Mpc), while allowing too much stellar mass loss leads to an under-estimate. For each sample, we are thus able to constrain the amount of stellar mass loss required to match the observed clustering. We find that less luminous satellite galaxies experience more efficient stellar mass loss than luminous satellites. From these models, we can infer the amount of stellar mass that is deposited into the IHL. We find that both of our model predictions for the mean amount of IHL as a function of halo mass are consistent with current observational measurements. However, our two models predict a different amount of scatter in the IHL from halo to halo, with Model 2 being favored by observations. This demonstrates that a comparison to IHL measurements provides independent verification of our stellar mass loss models. (Abridged)

We evaluate the exposure during nadir observations with JEM-EUSO, the Extreme Universe Space Observatory, on-board the Japanese Experiment Module of the International Space Station. Designed as a mission to explore the extreme energy Universe from space, JEM-EUSO will monitor the Earth's nighttime atmosphere to record the ultraviolet light from tracks generated by extensive air showers initiated by ultra-high energy cosmic rays. In the present work, we discuss the particularities of space-based observation and we compute the annual exposure in nadir observation. The results are based on studies of the expected trigger aperture and observational duty cycle, as well as, on the investigations of the effects of clouds and different types of background light. We show that the annual exposure is about one order of magnitude higher than those of the presently operating ground-based observatories.

We investigate the potential to improve optical tracers of cluster mass by exploiting measurements of the magnitude gap, m12, defined as the difference between the r-band absolute magnitude of the two brightest cluster members. We find that in a mock sample of galaxy groups and clusters constructed from the Bolshoi simulation, the scatter about the mass-richness relation decreases by 15-20% when magnitude gap information is included. A similar trend is evident in a volume-limited, spectroscopic sample of galaxy groups observed in the Sloan Digital Sky Survey (SDSS). We find that SDSS groups with small magnitude gaps are richer than large-gap groups at fixed values of the one-dimensional velocity dispersion among group members sigma_v, which we use as a mass proxy. We demonstrate explicitly that m12 contains information about cluster mass that supplements the information provided by group richness and the luminosity of the brightest cluster galaxy, L_bcg. In so doing, we show that the luminosities of the members of a group with richness N are inconsistent with the distribution of luminosities that results from N random draws from the global galaxy luminosity function. As the cosmological constraining power of galaxy clusters is limited by the precision in cluster mass determination, our findings suggest a new way to improve the cosmological constraints derived from galaxy clusters.

We measure the luminosity and color dependence and the redshift evolution of galaxy clustering in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Ninth Data Release. We focus on the projected two-point correlation function (2PCF) of subsets of its CMASS sample, which includes about 260,000 galaxies over ~3,300 sq. deg in the redshift range 0.43<z<0.7. To minimize the selection effect on galaxy clustering, we construct well-defined luminosity and color subsamples by carefully accounting for the CMASS galaxy selection cuts. The 2PCF of the whole CMASS sample, if approximated by a power-law, has a correlation length of r_0=7.93\pm0.06Mpc/h and an index of \gamma=1.85\pm0.01. Clear dependences on galaxy luminosity and color are found for the projected 2PCF in all redshift bins, with more luminous and redder galaxies generally exhibiting stronger clustering and steeper 2PCF. The color dependence is also clearly seen for galaxies within the red sequence, consistent with the behavior of SDSS-II main sample galaxies at lower redshifts. At a given luminosity (k+e corrected), no significant evolution of the projected 2PCFs with redshift is detected for red sequence galaxies. We also construct galaxy samples of fixed number density at different redshifts, using redshift-dependent magnitude thresholds. The clustering of these galaxies in the CMASS redshift range is found to be consistent with that predicted by passive evolution. Our measurements of the luminosity and color dependence and redshift evolution of galaxy clustering will allow for detailed modeling of the relation between galaxies and dark matter halos and new constraints on galaxy formation and evolution.

We report on the small scale (0.5<r<40h^-1 Mpc) clustering of 78895 massive (M*~10^11.3M_sun) galaxies at 0.2<z<0.4 from the first two years of data from the Baryon Oscillation Spectroscopic Survey (BOSS), to be released as part of SDSS Data Release 9 (DR9). We describe the sample selection, basic properties of the galaxies, and caveats for working with the data. We calculate the real- and redshift-space two-point correlation functions of these galaxies, fit these measurements using Halo Occupation Distribution (HOD) modeling within dark matter cosmological simulations, and estimate the errors using mock catalogs. These galaxies lie in massive halos, with a mean halo mass of 5.2x10^13 h^-1 M_sun, a large scale bias of ~2.0, and a satellite fraction of 12+/-2%. Thus, these galaxies occupy halos with average masses in between those of the higher redshift BOSS CMASS sample and the original SDSS I/II LRG sample.

The Baryon Oscillation Spectroscopic Survey (BOSS) is designed to measure the scale of baryon acoustic oscillations (BAO) in the clustering of matter over a larger volume than the combined efforts of all previous spectroscopic surveys of large scale structure. BOSS uses 1.5 million luminous galaxies as faint as i=19.9 over 10,000 square degrees to measure BAO to redshifts z<0.7. Observations of neutral hydrogen in the Lyman alpha forest in more than 150,000 quasar spectra (g<22) will constrain BAO over the redshift range 2.15<z<3.5. Early results from BOSS include the first detection of the large-scale three-dimensional clustering of the Lyman alpha forest and a strong detection from the Data Release 9 data set of the BAO in the clustering of massive galaxies at an effective redshift z = 0.57. We project that BOSS will yield measurements of the angular diameter distance D_A to an accuracy of 1.0% at redshifts z=0.3 and z=0.57 and measurements of H(z) to 1.8% and 1.7% at the same redshifts. Forecasts for Lyman alpha forest constraints predict a measurement of an overall dilation factor that scales the highly degenerate D_A(z) and H^{-1}(z) parameters to an accuracy of 1.9% at z~2.5 when the survey is complete. Here, we provide an overview of the selection of spectroscopic targets, planning of observations, and analysis of data and data quality of BOSS.

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

Contributions of the JEM-EUSO Collaboration to the 32nd International Cosmic Ray Conference, Beijing, August, 2011.

We present a fast method of producing mock galaxy catalogues that can be used to compute covariance matrices of large-scale clustering measurements and test the methods of analysis. Our method populates a 2nd-order Lagrangian Perturbation Theory (2LPT) matter field, where we calibrate masses of dark matter halos by detailed comparisons with N-body simulations. We demonstrate the clustering of halos is recovered at ~10 per cent accuracy. We populate halos with mock galaxies using a Halo Occupation Distribution (HOD) prescription, which has been calibrated to reproduce the clustering measurements on scales between 30 and 80 Mpc/h. We compare the sample covariance matrix from our mocks with analytic estimates, and discuss differences. We have used this method to make catalogues corresponding to Data Release 9 of the Baryon Oscillation Spectroscopic Survey (BOSS),producing 600 mock catalogues of the "CMASS" galaxy sample. These mocks enabled detailed tests of methods and errors that formed an integral part of companion analyses of these galaxy data.

We investigate the radial distribution of galaxies within their host dark matter halos by modeling their small-scale clustering, as measured in the Sloan Digital Sky Survey. Specifically, we model the Jiang et al. (2011) measurements of the galaxy two-point correlation function down to very small projected separations (10 < r < 400 kpc/h), in a wide range of luminosity threshold samples (absolute r-band magnitudes of -18 up to -23). We use a halo occupation distribution (HOD) framework with free parameters that specify both the number and spatial distribution of galaxies within their host dark matter halos. We assume that the first galaxy in each halo lives at the halo center and that additional satellite galaxies follow a radial density profile similar to the dark matter Navarro-Frenk-White (NFW) profile, except that the concentration and inner slope are allowed to vary. We find that in low luminosity samples, satellite galaxies have radial profiles that are consistent with NFW. M_r < -20 and brighter satellite galaxies have radial profiles with significantly steeper inner slopes than NFW (we find inner logarithmic slopes ranging from -1.6 to -2.1, as opposed to -1 for NFW). We define a useful metric of concentration, M_(1/10), which is the fraction of satellite galaxies (or mass) that are enclosed within one tenth of the virial radius of a halo. We find that M_(1/10) for low luminosity satellite galaxies agrees with NFW, whereas for luminous galaxies it is 2.5-4 times higher, demonstrating that these galaxies are substantially more centrally concentrated within their dark matter halos than the dark matter itself. Our results therefore suggest that the processes that govern the spatial distribution of galaxies, once they have merged into larger halos, must be luminosity dependent, such that luminous galaxies become poor tracers of the underlying dark matter.

We provide realistic mock-catalogs of cosmic rays above 40 EeV, for a pure proton composition, assuming their sources are a random subset of ordinary galaxies in a simulated, volume-limited survey, for various choices of source density: 10^-3.5 Mpc^-3, 10^-4.0 Mpc^-3 and 10^-4.5 Mpc^-3. The spectrum at the source is taken to be E^-2.3 and the effects of cosmological redshifting as well as photo-pion and e^+ e^- energy losses are included.

This work studies the environments and star formation relationships of local luminous infrared galaxies (LIRG) in comparison to other types of local and distant (z~1) galaxies. The infrared (IR) galaxies are drawn from the IRAS sample. The density of the environment is quantified using 6dF and Point Source Catalogue redshift survey (PSCz) galaxies in a cylinder of 2h^-1 Mpc radius and 10h^-1 Mpc length. Our most important result shows the existence of a dramatic density difference between local LIRGs and local non-LIRG IR galaxies. LIRGs live in denser environments than non-LIRG IR galaxies implying that L_IR=10^11 h^-2 L_sun marks an important transition point among IR-selected local galaxies. We also find that there is a strong correlation between the densities around LIRGs and their L_IR luminosity, while the IR-activity of non-LIRG IR galaxies does not show any dependence on environment. This trend is independent of mass-bin selection. The SF-density trend in local LIRGs is similar to that found in some studies of blue cloud galaxies at z~1 which show a correlation between star formation and local density (the reversal of the relation seen for local galaxies). This, together with the rapid decline of the number count of LIRGs since z~1, could mean that local LIRGs are survivors of whatever process transformed blue cloud galaxies at z~1 to the present day or local LIRGs came into existence by similar process than high redshift LIRGs but at later stage.

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

We model the evolution of galaxy clustering through cosmic time to investigate the nature of the power-law shape of xi(r), the galaxy two-point correlation function. Departures from a power law are mainly governed by galaxy pair counts on small scales, subject to non-linear dynamics. We assume that galaxies reside within dark matter halos and subhalos and use a semi-analytic substructure evolution model to study subhalo populations within host halos. We find that tidal mass loss and, to a lesser degree, dynamical friction deplete the number of subhalos within larger host halos over time by ~90%, just the right amount for achieving a power-law xi(r) at z = 0. We find that xi(r) breaks from a power law at high masses, implying that only galaxies of luminosities <= Lstar should exhibit power-law clustering. We also demonstrate that xi(r) evolves from being far from a power law at high redshift, toward a near power law at z = 0 and deviates in the future. This is mainly caused by the evolving competition between the accretion and destruction rates of subhalos, which happen to strike just the right balance at z~0. We show that key ingredients determining the shape of xi(r) are the fraction of galaxies that are satellites, the relative difference in mass between the halos of isolated galaxies and halos that contain a single satellite on average, and the rareness of halos that host galaxies. These pieces are intertwined and we find no simple, universal rule for which a power-law xi(r) will occur. However, we do show that the physics responsible for setting the galaxy content of halos do not care about the conditions needed to achieve a power law xi(r) and these conditions are met only in a narrow mass and redshift range. We conclude that the power-law nature of xi(r) for Lstar and fainter galaxy samples at low redshift is a cosmic coincidence. (Abridged)

We measure the luminosity and color dependence of galaxy clustering in the SDSS DR7 main galaxy sample, focusing on the projected correlation function w_p(r_p) of volume-limited samples. We interpret our measurements using halo occupation distribution (HOD) modeling assuming a Lambda-CDM cosmology. The amplitude of w_p(r_p) grows slowly with luminosity for L < L_* and increases sharply at higher luminosities, with bias factor b(>L)=1.06+0.23(L/L_*)^{1.12}. At fixed luminosity, redder galaxies have a stronger and steeper w_p(r_p), a trend that runs steadily from the bluest galaxies to the reddest galaxies. The individual luminosity trends for the red and blue galaxy populations are strikingly different. Blue galaxies show a slow but steady increase of w_p(r_p) with luminosity, at all scales. The large-scale clustering of red galaxies shows little luminosity dependence until a sharp increase at L > 4L_*, but the lowest luminosity red galaxies (0.04-0.25 L_*) show very strong clustering on scales r_p < 2 Mpc/h. Most of the observed trends can be naturally understood within the LCDM+HOD framework. The growth of w_p(r_p) with luminosity reflects an overall shift in the halo mass scale, in particular an increase in the minimum host halo mass Mmin. The mass at which a halo has, on average, one satellite galaxy brighter than L is M_1 ~ 17 Mmin(L) over most of the luminosity range. The growth and steepening of w_p(r_p) for redder galaxies reflects the increasing fraction of galaxies that are satellite systems in high mass halos instead of central systems in low mass halos, a trend that is especially marked at low luminosities. Our extensive measurements, provided in tabular form, will allow detailed tests of theoretical models of galaxy formation, a firm grounding of semi-empirical models of the galaxy population, and new cosmological tests.

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

We calculate the real- and redshift-space clustering of massive galaxies at z~0.5 using the first semester of data by the Baryon Oscillation Spectroscopic Survey (BOSS). We study the correlation functions of a sample of 44,000 massive galaxies in the redshift range 0.4<z<0.7. We present a halo-occupation distribution modeling of the clustering results and discuss the implications for the manner in which massive galaxies at z~0.5 occupy dark matter halos. The majority of our galaxies are central galaxies living in halos of mass 10^{13}Msun/h, but 10% are satellites living in halos 10 times more massive. These results are broadly in agreement with earlier investigations of massive galaxies at z~0.5. The inferred large-scale bias (b~2) and relatively high number density (nbar=3e-4 h^3 Mpc^{-3}) imply that BOSS galaxies are excellent tracers of large-scale structure, suggesting BOSS will enable a wide range of investigations on the distance scale, the growth of large-scale structure, massive galaxy evolution and other topics.

We constrain the linear and quadratic bias parameters from the configuration dependence of the three-point correlation function (3PCF) in both redshift and projected space, utilizing measurements of spectroscopic galaxies in the Sloan Digital Sky Survey (SDSS) Main Galaxy Sample. We show that bright galaxies (M_r < -21.5) are biased tracers of mass, measured at a significance of 4.5 sigma in redshift space and 2.5 sigma in projected space by using a thorough error analysis in the quasi-linear regime (9-27 Mpc/h). Measurements on a fainter galaxy sample are consistent with an unbiased model. We demonstrate that a linear bias model appears sufficient to explain the galaxy-mass bias of our samples, although a model using both linear and quadratic terms results in a better fit. In contrast, the bias values obtained from the linear model appear in better agreement with the data by inspection of the relative bias, and yield implied values of sigma_8 that are more consistent with current constraints. We investigate the covariance of the 3PCF, which itself is a measurement of galaxy clustering. We assess the accuracy of our error estimates by comparing results from mock galaxy catalogs to jackknife re-sampling methods. We identify significant differences in the structure of the covariance. However, the impact of these discrepancies appears to be mitigated by an eigenmode analysis that can account for the noisy, unresolved modes. Our results demonstrate that using this technique is sufficient to remove potential systematics even when using less-than-ideal methods to estimate errors.

We analyse the Fundamental Plane (FP) relation of $39,993$ early-type galaxies (ETGs) in the optical (griz) and $5,080$ ETGs in the Near-Infrared (YJHK) wavebands, forming an optical$+$NIR sample of $4,589$ galaxies. We focus on the analysis of the FP as a function of the environment where galaxies reside. We characterise the environment using the largest group catalogue, based on 3D data, generated from SDSS at low redshift ($z < 0.1$). We find that the intercept $``c''$ of the FP decreases smoothly from high to low density regions, implying that galaxies at low density have on average lower mass-to-light ratios than their high-density counterparts. The $``c''$ also decreases as a function of the mean characteristic mass of the parent galaxy group. However, this trend is weak and completely accounted for by the variation of $``c''$ with local density. The variation of the FP offset is the same in all wavebands, implying that ETGs at low density have younger luminosity-weighted ages than cluster galaxies, consistent with the expectations of semi-analytical models of galaxy formation. We measure an age variation of $\sim 0.048$~dex ($\sim 11\%$) per decade of local galaxy density. This implies an age difference of about $32 \%$ ($\sim 3 \, Gyr$) between galaxies in the regions of highest density and the field. We find the metallicity decreasing, at $\sim 2$~$\sigma$, from low to high density. We also find $2.5 \, \sigma$ evidence that the variation in age per decade of local density augments, up to a factor of two, for galaxies residing in massive relative to poor groups. (abridged)

We analyze the line-of-sight baryonic acoustic feature in the two-point correlation function {\xi} of the Sloan Digital Sky Survey (SDSS) luminous red galaxy (LRG) sample (0.16 < z < 0.47). By defining a narrow line-of-sight region, rp < 5.5 Mpc/h, where rp is the transverse separation component, we measure a strong excess of clustering at ~ 110 Mpc/h, as previously reported in the literature. We also test these results in an alternative coordinate system, by defining the line-of-sight as {\theta} < 3{\deg}, where {\theta} is the opening angle. This clustering excess appears much stronger than the feature in the better-measured monopole. A fiducial {\Lambda}CDM non-linear model in redshift-space predicts a much weaker signature. We use realistic mock catalogs to model the expected signal and noise. We find that the line-of-sight measurements can be explained well by our mocks as well as by a featureless {\xi} = 0. We conclude that there is no convincing evidence that the strong clustering measurement is the line-of-sight baryonic acoustic feature. We also evaluate how detectable such a signal would be in the upcoming Baryon Oscillation Spectroscopic Survey LRG volume (BOSS). Mock LRG catalogs (z < 0.6) suggest that: (i) the narrow line- of-sight cylinder and cone defined above probably will not reveal a detectable acoustic feature in BOSS; (ii) a clustering measurement as high as that in the current sample can be ruled out (or confirmed) at a high confidence level using a BOSS-sized data set; and (iii) an analysis with wider angular cuts, which provide better signal-to-noise ratios, can nevertheless be used to compare line-of-sight and transverse distances, and thereby constrain the expansion rate H(z) and diameter distance DA(z).

We confirm the UHECR horizon established by the Pierre Auger Observatory using the heterogeneous Veron-Cetty Veron (VCV) catalog of AGNs, by performing a redshift-angle-IR luminosity scan using PSCz galaxies having infrared luminosity greater than 10^{10}L_sun. The strongest correlation -- for z < 0.016, psi = 2.1 deg, and L_ir > 10^{10.5}L_sun -- arises in fewer than 0.3% of scans with isotropic source directions. When we apply a penalty for using the UHECR energy threshold that was tuned to maximize the correlation with VCV, the significance degrades to 1.1%. Since the PSCz catalog is complete and volume-limited for these parameters, this suggests that the UHECR horizon discovered by the Pierre Auger Observatory is not an artifact of the incompleteness and other idiosyncrasies of the VCV catalog. The strength of the correlation between UHECRs and the nearby highest-IR-luminosity PSCz galaxies is stronger than in about 90% percent of trials with scrambled luminosity assignments for the PSCz galaxies. If confirmed by future data, this result would indicate that the sources of UHECRs are more strongly associated with luminous IR galaxies than with ordinary, lower IR luminosity galaxies.

Any multivariate distribution can be uniquely decomposed into marginal (1-point) distributions, and a function called the copula, which contains all of the information on correlations between the distributions. The copula provides an important new methodology for analyzing the density field in large-scale structure. We derive the empirical 2-point copula for the evolved dark matter density field. We find that this empirical copula is well-approximated by a Gaussian copula. We consider the possibility that the full n-point copula is also Gaussian and describe some of the consequences of this hypothesis. Future directions for investigation are discussed.

We model the small-scale clustering of luminous red galaxies (LRGs; Masjedi et al. 2006) in the Sloan Digital Sky Survey (SDSS). Specifically, we use the halo occupation distribution (HOD) formalism to model the projected two-point correlation function of LRGs on scales well within the sizes of their host halos (0.016 Mpc/h < r < 0.42 Mpc/h). We start by varying P(N|M), the probability distribution that a dark matter halo of mass M contains N LRGs, and assuming that the radial distribution of satellite LRGs within halos traces the NFW dark matter density profile. We find that varying P(N|M) alone is not sufficient to match the small-scale data. We next allow the concentration of satellite LRG galaxies to differ from that of dark matter and find that this is also not sufficient. Finally, we relax the assumption of an NFW profile and allow the inner slope of the density profile to vary. We find that this model provides a good fit to the data and the resulting value of the slope is -2.17 +/- 0.12. The radial density profile of satellite LRGs within halos is thus not compatible with that of the underlying dark matter, but rather is closer to an isothermal distribution.

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

We investigate several aspects of the correlations reported by the Pierre Auger Observatory between the highest energy cosmic rays (UHECRs) and galaxies in the Veron-Cetty Veron (VCV) catalog of AGNs. First, we quantify the extent of the inhomogeneity and impurity of the VCV catalog. Second, we determine how the correlation between the highest energy Auger UHECRs and VCV galaxies is modified when only optically-identified AGNs are used. Finally, we measure the correlation between the published Auger UHECRs and the distribution of matter. Our most important finding is that the correlation between UHECRs and AGNs is too strong to be explained purely by AGNs tracing the large scale distribution of matter, indicating that (barring the correlation being a statistical fluke) some substantial fraction of UHECRs are produced by AGNs. We also find that once we take into account the heavy oversampling of the VCV catalog in the Virgo region, the lack of UHECR events from that region is not incompatible with UHECR having AGN sources.

Galaxies are seen from different viewing angles and their observed properties change as a function of viewing angle. In many circumstances we would rather know the intrinsic properties of galaxies - those properties that do not depend on viewing angle. For a large sample of galaxies it is possible to recover the intrinsic properties of galaxies, statistically, by looking for correlations with galaxy inclination, and then applying a correction to remove those dependencies. Studying the intrinsic properties of galaxies can give a different impression of the galaxy population and help avoid the mistake of connecting observed properties to quantities that don't depend on inclination like halo mass.

The observed properties of galaxies vary with inclination; for most applications we would rather have properties that are independent of inclination, intrinsic properties. One way to determine inclination corrections is to consider a large sample of galaxies, study how the observed properties of these galaxies depend on inclination and then remove this dependence to recover the intrinsic properties. We perform such an analysis for galaxies selected from the Sloan Digital Sky Survey which have been matched to galaxies from the Two-Micron All Sky Survey. We determine inclination corrections for these galaxies as a function of galaxy luminosity and Sersic index. In the g-band these corrections reach as as high as 1.2 mag and have a median value of 0.3 mag for all galaxies in our sample. We find that the corrections show little dependence on galaxy luminosity, except in the $u$ band, but are strongly dependent on galaxy Sersic index. We find that the ratio of red-to-blue galaxies changes from 1:1 to 1:2 when going from observed to intrinsic colors for galaxies in the range -22.75 < M_K < -17.75. We also discuss how survey completeness and photometric redshifts should be determined when taking into account that observed and intrinsic properties differ. Finally, we examine whether previous determinations of stellar mass give an intrinsic quantity or one that depends on galaxy inclination.

We show that future Ultra-High Energy Cosmic Ray samples should be able to distinguish whether the sources of UHECRs are hosted by galaxy clusters or ordinary galaxies, or whether the sources are uncorrelated with the large-scale structure of the universe. Moreover, this is true independently of arrival direction uncertainty due to magnetic deflection or measurement error. The reason for this is the simple property that the strength of large-scale clustering for extragalactic sources depends on their mass, with more massive objects, such as galaxy clusters, clustering more strongly than lower mass objects, such as ordinary galaxies.

(Abridged) We present a semi-analytic model for Cold Dark Matter halo substructure that can be used as a framework for studying galaxy formation and an ingredient in halo models of galaxy clustering. We perform a comprehensive comparison of the model to the results of a suite of high-resolution cosmological simulations. The comparisons reveal that subhalo statistics, such as velocity and mass functions, radial distributions, and the halo occupation distributions, agree well over three orders of magnitude in host halo mass at various redshifts. Both in the simulations and in our model, the radial distributions of subhalos are significantly shallower than that of the dark matter density. The abundance of subhalos in a host is set by competition between destruction and new accretion. High mass halos and halos at high redshift tend to host more subhalos because the subhalos have, on average, been accreted more recently. Similarly, at fixed mass, halos that formed more recently host more subhalos. Observed "fossil groups" may represent an extreme tail of this correlation. We find a related correlation between host halo concentration and subhalo abundance at fixed host mass, Nsat ~ c^-a, where a changes with redshift and host-to-subhalo mass ratio. Lastly, we use our model to populate host halos in one of our high-resolution cosmological simulations, replacing the subhalos resolved in the simulation with subhalos computed according to our model. We show that the resulting 2pt correlation function of such a "hybrid" halo ensemble is indistinguishable from that measured directly in the simulation. This supports one of the key tenets of the standard halo model -- the assumption that the halo occupation distribution is statistically independent of host halo environment.

We investigate the dependences of galaxy star-formation history and galaxy morphology on environment, using color and H-alpha equivalent width as star-formation history indicators, using concentration and central surface brightness as morphological indicators, and using clustocentric distance as an environment indicator. Clustocentric distance has the virtue that it can be measured with very high precision over a large dynamic range. We find the following asymmetry between morphological and star-formation history parameters: star-formation history parameters relate directly to the clustocentric distance while morphological parameters relate to the clustocentric distance only indirectly through their relationships with star-formation history. This asymmetry has important implications for the role that environment plays in shaping galaxy properties and it places strong constraints on theoretical models of galaxy formation. Current semi-analytic models do not reproduce this effect.

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 investigate the clustering of galaxy groups and clusters in the SDSS using the Berlind et al. (2006) group sample, which is designed to identify galaxy systems that each occupy a single dark matter halo. We estimate group masses from their abundances, and measure their relative large-scale bias as a function of mass. Our measurements are in agreement with the theoretical halo bias function, given a standard LCDM cosmological model, and they tend to favor a low value of the power spectrum amplitude sigma_8. We search for a residual dependence of clustering on other group properties at fixed mass, and find the strongest signal for central galaxy color in high mass groups. Massive groups with less red central galaxies are more biased on large scales than similar mass groups with redder central galaxies. We show that this effect is unlikely to be caused by errors in our mass estimates, and is most likely observational evidence of recent theoretical findings that halo bias depends on a ``second parameter'' other than mass, such as age or concentration. To compare with the data, we study the bias of massive halos in N-body simulations and quantify the strength of the relation between halo bias and concentration at fixed mass. In addition to confirming a non-trivial prediction of the LCDM cosmological model, these results have important implications for the role that environment plays in shaping galaxy properties.

We determine what aspects of the density field surrounding galaxies most affect their properties. For Sloan Digital Sky Survey galaxies, we measure the group environment, meaning the host group luminosity and the distance from the group center (hereafter, ``groupocentric distance''). For comparison, we measure the surrounding density field on scales ranging from 100 kpc/h to 10 Mpc/h. We use the relationship between color and group environment to test the null hypothesis that only the group environment matters, searching for a residual dependence of properties on the surrounding density. Generally, red galaxies are slightly more clustered on small scales (about 100--300 kpc/h) than the null hypothesis predicts, possibly indicating that substructure within groups has some importance. At large scales (> 1 Mpc/h), the actual projected correlation functions of galaxies are biased at less than the 5% level with respect to the null hypothesis predictions. We exclude strongly the converse null hypothesis, that only the surrounding density (on any scale) matters. These results generally encourage the use of the halo model description of galaxy bias, which models the galaxy distribution as a function of host halo mass alone. We compare these results to proposed galaxy formation scenarios within the Cold Dark Matter cosmological model.

We identify galaxy groups and clusters in volume-limited samples of the SDSS redshift survey, using a redshift-space friends-of-friends algorithm. We optimize the friends-of-friends linking lengths to recover galaxy systems that occupy the same dark matter halos, using a set of mock catalogs created by populating halos of N-body simulations with galaxies. Extensive tests with these mock catalogs show that no combination of perpendicular and line-of-sight linking lengths is able to yield groups and clusters that simultaneously recover the true halo multiplicity function, projected size distribution, and velocity dispersion. We adopt a linking length combination that yields, for galaxy groups with ten or more members: a group multiplicity function that is unbiased with respect to the true halo multiplicity function; an unbiased median relation between the multiplicities of groups and their associated halos; a spurious group fraction of less than ~1%; a halo completeness of more than ~97%; the correct projected size distribution as a function of multiplicity; and a velocity dispersion distribution that is ~20% too low at all multiplicities. These results hold over a range of mock catalogs that use different input recipes of populating halos with galaxies. We apply our group-finding algorithm to the SDSS data and obtain three group and cluster catalogs for three volume-limited samples that cover 3495.1 square degrees on the sky. We correct for incompleteness caused by fiber collisions and survey edges, and obtain measurements of the group multiplicity function, with errors calculated from realistic mock catalogs. These multiplicity function measurements provide a key constraint on the relation between galaxy populations and dark matter halos.

(abridged) There are good observational reasons to believe that the progenitors of red galaxies have undergone starbursts, followed by a post-starburst phase. We investigate the environments of post-starburst galaxies by measuring \textsl{(1)} number densities in $8 h^{-1} \mathrm{Mpc}$ radius comoving spheres, \textsl{(2)} transverse distances to nearest Virgo-like galaxy clusters, and \textsl{(3)} transverse distances to nearest luminous-galaxy neighbors. We compare the post-starburst galaxies to currently star-forming galaxies identified solely by A-star excess or $\Halpha$ emission. We find that post-starburst galaxies are in the same kinds of environments as star-forming galaxies; this is our ``null hypothesis''. More importantly, we find that at each value of the A-star excess, the star-forming and post-starburst galaxies lie in very similar distributions of environment. The only deviations from our null hypothesis are barely significant: a slight deficit of post-starburst galaxies (relative to the star-forming population) in very low-density regions, a small excess inside the virial radii of clusters, and a slight excess with nearby neighbors. None of these effects is strong enough to make the post-starburst galaxies a high-density phenomenon, or to argue that the starburst events are primarily triggered by external tidal impulses (e.g., from close passages of massive galaxies). The small excess inside cluster virial radii suggests that some post-starbursts are triggered by interactions with the intracluster medium, but this represents a very small fraction of all post-starburst galaxies.

We investigate the extent to which galaxies' star-formation histories and morphologies are determined by their clustocentric distance and their nearest cluster's richness. We look at three tracers of star-formation history (i-band absolute magnitude M_i, [g-r] color, and Halpha emission line equivalent width) and two indicators of galaxy morphology (surface brightness and radial concentration) for 52,569 galaxies in the redshift range of 0.015<z<0.068. We find that our morphology indicators relate to the clustocentric distance only indirectly through their relationships with stellar population and starformation rate. Galaxies that are near the cluster center tend to be more luminous, redder and have lower Halpha EW (ie, lower star-formation rates) than those that lie near or outside the virial radius of the cluster. The detailed relationships between these galaxy properties and clustocentric distance depend on cluster richness. For richer clusters, we find that (i) the transition in color and Halpha EW from cluster center to field values is more abrupt and occurs closer to the cluster virial radius, and (ii) the color and Halpha EW distributions are overall narrower than in less rich clusters. We also find that the radial gradient seen in the luminosity distribution is strongest around the smaller clusters and decreases as cluster richness increases. We find there is a `characteristic distance' at around one virial radius (the infall region) where the change with radius of galaxy property distributions is most dramatic, but we find no evidence for infall-triggered star bursts. These results suggest that galaxies `know' the distance to, and the size of, their nearest cluster and they express this information in their star-formation histories.

We present the small-scale (0.01<r<8 h^{-1} Mpc) projected correlation function w_p(r_p) and real space correlation function xi(r) of 24520 luminous early-type galaxies from the Sloan Digital Sky Survey Luminous Red Galaxy (LRG) sample (0.16<z<0.36). ``Fiber collision'' incompleteness of the SDSS spectroscopic sample at scales smaller than 55 arcseconds prevents measurements of the correlation function for LRGs on scales smaller than ~0.3 Mpc by the usual methods. In this work, we cross-correlate the spectroscopic sample with the imaging sample, with a weighting scheme to account for the collisions, extensively tested against mock catalogs. We correct for photometric biases in the SDSS imaging of close galaxy pairs. We find that the correlation function xi(r) is surprisingly close to a r^{-2} power law over more than 4 orders of magnitude in separation r. This result is too steep at small scales to be explained in current versions of the halo model for galaxy clustering. We infer an LRG-LRG merger rate of less than or similar to 0.6 Gyr^{-1} Gpc^{-3} for this sample. This result suggests that the LRG-LRG mergers are not the main mode of mass growth for LRGs at z<0.36.

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 projected correlation function w_p(r_p) from the Sloan Digital Sky Survey for a flux-limited sample of 118,000 galaxies and for a volume limited subset of 22,000 galaxies with absolute magnitude M_r<-21. Both correlation functions show subtle but systematic departures from the best-fit power law, in particular a change in slope at r_p~1-2 Mpc/h. These departures are stronger for the volume-limited sample, which is restricted to relatively luminous galaxies. We show that the inflection point in w_p(r_p) can be naturally explained by contemporary models of galaxy clustering, according to which it marks the transition from a large scale regime dominated by galaxy pairs in separate dark matter halos to a small scale regime dominated by galaxy pairs in the same dark matter halo. For example, given the dark halo population predicted by an inflationary cold dark matter scenario, the projected correlation function of the volume-limited sample can be well reproduced by a model in which the mean number of M_r<-21 galaxies in a halo of mass M>M_1=4.74 X 10^{13}\msun/h is <N>_M=(M/M_1)^{0.89}, with 75% of the galaxies residing in less massive, single-galaxy halos. This physically motivated model has the same number of free parameters as a power law, and it fits the w_p(r_p) data better, with a \chi^2/d.o.f.=0.93 compared to 6.12 (for 10 degrees of freedom, incorporating the covariance of the correlation function errors). Departures from a power-law correlation function encode information about the relation between galaxies and dark matter halos. Higher precision measurements of these departures for multiple classes of galaxies will constrain galaxy bias and provide new tests of the theory of galaxy formation.

We investigate the origin of a nearly pointlike cluster of 5 ultrahigh energy cosmic rays at RA ~169.2deg and dec ~56.8deg, using Sloan Digital Sky Survey and other data. No particular source candidates are found near the estimated source direction, but the direction is exceptional in having a likely merging pair of galaxy clusters at 140/h Mpc, with an unusually low foreground density. Large scale shocks or another product of the merging galaxy clusters may accelerate the UHECRs, or the merging galaxy clusters may be coincidental and the UHECRs may be accelerated in a rare event of an unexceptional progenitor. Low magnetic deflections in the foreground void may explain why this is the only identified pointlike cluster of so many UHECRs.

The halo occupation distribution (HOD) describes the relation between galaxies and dark matter at the level of individual dark matter halos. The properties of galaxies residing at the centers of halos differ from those of satellite galaxies because of differences in their formation histories. Using a smoothed particle hydrodynamics (SPH) simulation and a semi-analytic (SA) galaxy formation model, we examine the separate contributions of central and satellite galaxies to the HOD, more specifically to the probability P(N|M) that a halo of virial mass M contains N galaxies of a particular class. In agreement with earlier results for dark matter subhalos, we find that the mean occupation function <N> for galaxies above a baryonic mass threshold can be approximated by a step function for central galaxies plus a power law for satellites, and that the distribution of satellite numbers is close to Poisson at fixed halo mass. For galaxy samples defined by different baryonic mass thresholds, there is a nearly linear relation between the minimum halo mass Mmin required to host a central galaxy and the mass M1 at which an average halo hosts one satellite, with M1 ~ 14 Mmin (SPH) or M1 ~ 18 Mmin (SA). The mean occupation number of young galaxies exhibits a local minimum at M ~ 10 Mmin where halos are too massive to host a young central galaxy but not massive enough to host satellites. We show that the conditional galaxy mass function at fixed halo mass cannot be described by a Schechter function because central galaxies produce a "bump" at high masses. We suggest parameterizations for the HOD and the conditional luminosity function that can be used to model observed galaxy clustering. Many of our predictions are in good agreement with recent results inferred from clustering in the Sloan Digital Sky Survey.

We use the projected correlation function w_p(r_p) of a volume-limited subsample of the Sloan Digital Sky Survey (SDSS) main galaxy redshift catalogue to measure the halo occupation distribution (HOD) of the galaxies of the sample. Simultaneously, we allow the cosmology to vary within cosmological constraints imposed by cosmic microwave background experiments in a Lambda-CDM model. We find that combining w_p(r_p) for this sample alone with the observations by WMAP, ACBAR, CBI and VSA can provide one of the most precise techniques available to measure cosmological parameters. For a minimal flat six-parameter Lambda-CDM model with an HOD with three free parameters, we find Omega_m=0.278^{+0.027}_{-0.027}, sigma_8=0.812^{+0.028}_{-0.027}, and H_0=69.8^{+2.6}_{-2.6}km s^{-1} Mpc^{-1}; these errors are significantly smaller than from CMB alone and similar to those obtained by combining CMB with the large-scale galaxy power spectrum assuming scale-independent bias. The corresponding HOD parameters describing the minimum halo mass and the normalization and cut-off of the satellite mean occupation are M_min=(3.03^{+0.36}_{-0.36})x 10^{12} h^{-1} M_sun, M_1 = (4.58^{+0.60}_{-0.60})x 10^{13} h^{-1} M_sun, and kappa=4.44^{+0.51}_{-0.69}. When more parameters are added to the HOD model, the error bars on the HOD parameters increase because of degeneracies, but the error bars on the cosmological parameters do not increase greatly. Similar modeling for other galaxy samples could reduce the statistical errors on these results, while more thorough investigations of the cosmology dependence of nonlinear halo bias and halo mass functions are needed to eliminate remaining systematic uncertainties, which may be comparable to statistical uncertainties.

We study the relationship between galaxy luminosity, color, and environment in a cosmological simulation of galaxy formation. We compare the predicted relationship with that found for SDSS galaxies and find that the model successfully predicts most of the qualitative features seen in the data, but also shows some interesting differences. Specifically, the simulation predicts that the local density around bright red galaxies is a strong increasing function of luminosity, but does not depend much on color at fixed luminosity. Moreover, we show that this is due to central galaxies in dark matter halos whose baryonic masses correlate strongly with halo mass. The simulation also predicts that the local density around blue galaxies is a strong increasing function of color, but does not depend much on luminosity at fixed color. We show that this is due to satellite galaxies in halos whose stellar ages correlate with halo mass. Finally, the simulation fails to predict the luminosity dependence of environment observed around low luminosity red galaxies. However, we show that this is most likely due to the simulation's limited resolution. A study of a higher resolution, smaller volume simulation suggests that this dependence is caused by the fact that all low luminosity red galaxies are satellites in massive halos, whereas intermediate luminosity red galaxies are a mixture of satellites in massive halos and central galaxies in less massive halos.

We analyze the halo occupation distribution (HOD), the probability for a halo of mass M to host a number of subhalos N, and two-point correlation function of galaxy-size dark matter halos using high-resolution dissipationless simulations of the concordance flat LCDM model. The halo samples include both the host halos and the subhalos, distinct gravitationally-bound halos within the virialized regions of larger host systems. We find that the first moment of the HOD, <N>(M), has a complicated shape consisting of a step, a shoulder, and a power law high-mass tail. The HOD can be described by a Poisson statistics at high halo masses but becomes sub-Poisson for <N><4. We show that the HOD can be understood as a combination of the probability for a halo of mass M to host a central galaxy and the probability to host a given number Ns of satellite galaxies. The former can be approximated by a step-like function, while the latter can be well approximated by a Poisson distribution, fully specified by its first moment <Ns>(M). We find that <Ns>~M^b with b~1 for a wide range of number densities, redshifts, and different power spectrum normalizations. This formulation provides a simple but accurate model for the halo occupation distribution found in simulations. At z=0, the two-point correlation function (CF) of galactic halos can be well fit by a power law down to ~100/h kpc with an amplitude and slope similar to those of observed galaxies. At redshifts z>~1, we find significant departures from the power-law shape of the CF at small scales. If the deviations are as strong as indicated by our results, the assumption of the single power law often used in observational analyses of high-redshift clustering is likely to bias the estimates of the correlation length and slope of the correlation 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.

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.

The halo occupation distribution (HOD) describes the bias between galaxies and dark matter by specifying (a) the probability P(N|M) that a halo of virial mass M contains N galaxies of a particular class and (b) the relative distributions of galaxies and dark matter within halos. We calculate predicted HODs for a Lambda-CDM cosmological model using an SPH hydrodynamic simulation and a semi-analytic (SA) galaxy formation model. Although the two methods predict different galaxy mass functions, their HOD predictions agree remarkably well. For mass-selected samples, the mean occupation <N(M)> exhibits a sharp cutoff at low halo masses, a slowly rising plateau for <N>~1-2, and a more steeply rising high occupancy regime. At low <N>, the mean pair and triple counts are well below Poisson expectations, with important consequences for small scale behavior of 2- and 3-point correlation functions. The HOD depends strongly on galaxy age, with high mass halos populated mainly by old galaxies and low mass halos by young galaxies. The SPH simulation supports several simplifying assumptions about HOD bias: the most massive galaxy in a halo usually lies close to the center and moves near the halo's mean velocity; satellite galaxies have the same radial profile and velocity dispersion as the dark matter; and the mean occupation at fixed halo mass is independent of the halo's larger scale environment. By applying the SPH and SA HODs to a large volume N-body simulation, we show that both methods predict slight, observable departures from a power-law galaxy correlation function. The predicted HODs are closely tied to the underlying galaxy formation physics, they offer useful guidance to theoretical models of galaxy clustering, and they will be tested empirically by ongoing analyses of galaxy redshift surveys. (Shortened)

We announce the initial release of data from the Ohio State University Bright Spiral Galaxy Survey, a BVRJHK imaging survey of a well-defined sample of 205 bright, nearby spiral galaxies. We present H-band morphological classification on the Hubble sequence for the OSU Survey sample. We compare the H-band classification to B-band classification from our own images and from standard galaxy catalogs. Our B-band classifications match well with those of the standard catalogs. On average, galaxies with optical classifications from Sa through Scd appear about one T-type earlier in the H-band than in the B-band, but with large scatter. This result does not support recent claims made in the literature that the optical and near-IR morphologies of spiral galaxies are uncorrelated. We present detailed descriptions of the H-band morphologies of our entire sample, as well as B- and H-band images for a set of 17 galaxies chosen as type examples, and BRH color-composite images of six galaxies chosen to demonstrate the range in morphological variation as a function of wavelength. Data from the survey are accessible at http://www.astronomy.ohio-state.edu/~survey/

Using cosmological N-body simulations, we investigate the influence of the matter density parameter Omega_m and the linear theory power spectrum P(k) on statistical properties of the dark matter halo population -- the mass function n(M), the halo correlation function, and halo pairwise velocity statistics. For fixed P(k), the effect of changing Omega_m is simple: the halo mass scale M_* shifts in proportion to Omega_m, pairwise velocities (at fixed M/M_*) are proportional to Omega_m^{0.6}, and halo clustering at fixed M/M_* is unchanged. If one simultaneously changes the P(k) amplitude sigma_8 to maintain a "cluster normalization" condition sigma_8 x Omega_m^{0.5} = const, then n(M) stays approximately constant near M ~ 5e14 Msun, and halo clustering and pairwise velocities are similar at fixed M. However, the shape of n(M) changes, with a decrease of Omega_m from 0.3 to 0.2 producing a ~30% drop in the number of low mass halos. One can preserve the shape of n(M) over a large dynamic range by changing the shape of P(k), but the required changes are substantial, and they significantly alter the halo clustering and halo velocities. The sensitivity of the dark halo population to cosmological model parameters has encouraging implications for efforts to constrain cosmology and galaxy bias with observed galaxy clustering, since the predicted changes in the halo population cannot easily be masked by altering the way that galaxies occupy halos. A shift in Omega_m alone would be detected by any dynamically sensitive clustering statistic; a cluster normalized change to sigma_8 and Omega_m would require a change in galaxy occupation as a function of M/M_*, which would alter galaxy clustering; and a simultaneous change to P(k) that preserves the halo mass function would change the clustering of the halos themselves.

We investigate galaxy bias in the framework of the ``Halo Occupation Distribution'' (HOD), which defines the bias of a population of galaxies by the conditional probability P(N|M) that a dark matter halo of virial mass M contains N galaxies, together with prescriptions that specify the relative spatial and velocity distributions of galaxies and dark matter within halos. By populating the halos of a cosmological N-body simulation using a variety of HOD models, we examine the sensitivity of different galaxy clustering statistics to properties of the HOD. The galaxy correlation function responds to different aspects of P(N|M) on different scales. Obtaining the observed power-law form of xi(r) requires rather specific combinations of HOD parameters, implying a strong constraint on the physics of galaxy formation; the success of numerical and semi-analytic models in reproducing this form is entirely non-trivial. Other clustering statistics such as the galaxy-mass correlation function, the bispectrum, the void probability function, the pairwise velocity dispersion, and the group multiplicity function are sensitive to different combinations of HOD parameters and thus provide complementary information about galaxy bias. We outline a strategy for determining the HOD empirically from redshift survey data. This method starts from an assumed cosmological model, but we argue that cosmological and HOD parameters will have non-degenerate effects on galaxy clustering, so that a substantially incorrect cosmological model will not reproduce the observations for any choice of HOD. Empirical determinations of the HOD as a function of galaxy type from the 2dF and SDSS redshift surveys will provide a detailed target for theories of galaxy formation, insight into the origin of galaxy properties, and sharper tests of cosmological models.

Measurements of the cosmological density parameter Omega using techniques that exploit the gravity-induced motions of galaxies constrain, in linear perturbation theory, the degenerate parameter combination beta = Omega^{0.6}/b, where the linear bias parameter b is the ratio of the fluctuation amplitudes of the galaxy and mass distributions. However, the relation between the galaxy and mass density fields depends on the complex physics of galaxy formation, and it can in general be non-linear, stochastic, and perhaps non-local. The one-parameter linear bias model is almost certainly oversimplified, which leads to the obvious question: What is the quantity beta that is actually measured by different techniques? To address this question, we estimate beta from galaxy distributions that are constructed by applying a variety of locally biased galaxy formation models to cosmological N-body simulations. We compare the values of beta estimated using three different techniques: a density-density comparison similar to the POTENT analysis, a velocity-velocity comparison similar to the VELMOD analysis, and an anisotropy analysis of the redshift-space power spectrum. In most cases, we find that beta estimated using all three methods is similar to the asymptotic value of Omega^{0.6}/b_{sigma}(R) at large R, where b_{sigma}(R) is the ratio of rms galaxy fluctuations to rms mass fluctuations on scale R. Thus, something close to the conventional interpretation of beta continues to hold even for complex bias models. Moreover, we find that beta estimates made using these three methods should, in principle, agree with each other. It is thus unlikely that non-linear or scale-dependent bias is responsible for the discrepancies that exist among current measurements of beta from different techniques.

We measure the topology (genus curve) of the galaxy distribution in a mock redshift catalog designed to resemble the upcoming Sloan Digital Sky Survey (SDSS). The catalog, drawn from a large N-body simulation of a Lambda-CDM cos- mological model, mimics the anticipated spectroscopic selection procedures of the SDSS in some detail. Sky maps, redshift slices, and 3-D contour maps of the mock survey reveal a rich and complex structure, including networks of voids and superclusters that resemble the patterns seen in the CfA redshift survey and the Las Campanas Redshift Survey (LCRS). The 3-D genus curve can be measured from the simulated catalog with superb precision; this curve has the general shape predicted for Gaussian, random phase initial conditions, but the error bars are small enough to demonstrate with high significance the subtle departures from this shape caused by non-linear gravitational evolution. These distortions have the form predicted by Matsubara's (1994) perturbative anal- ysis, but they are much smaller in amplitude. We also measure the 3-D genus curve of the radial peculiar velocity field measured by applying distance- indicator relations (with realistic errors) to the mock catalog. This genus curve is consistent with the Gaussian random phase prediction, though it is of relatively low precision because of the large smoothing length required to overcome noise in the measured velocity field. Finally, we measure the 2-D topology in redshift slices, similar to early slices from the SDSS and to slices already observed in the LCRS. The genus curves of these slices are consistent with the observed genus curves of the LCRS, providing further evidence in favor of the inflationary CDM model with Omega_M~0.4. The catalog is publicly available at http://www.astronomy.ohio-state.edu/~dhw/sdss.html.

We show that a regression of unsmoothed peculiar velocity measurements against peculiar velocities predicted from a smoothed galaxy density field leads to a biased estimate of the cosmological density parameter Omega, even when galaxies trace the underlying mass distribution and galaxy positions and velocities are known perfectly. The bias arises because the errors in the predicted velocities are correlated with the predicted velocities themselves. We investigate this bias using cosmological N-body simulations and analytic arguments. In linear perturbation theory, for cold dark matter power spectra and Gaussian or top hat smoothing filters, the bias in Omega is always positive, and its magnitude increases with increasing smoothing scale. This linear calculation reproduces the N-body results for Gaussian smoothing radii R_s > 10 Mpc/h, while non-linear effects lower the bias on smaller smoothing scales, and for R_s < 3 Mpc/h Omega is underestimated rather than overestimated. The net bias in Omega for a given smoothing filter depends on the underlying cosmological model. The effect on current estimates of Omega from velocity-velocity comparisons is probably small relative to other uncertainties, but taking full advantage of the statistical precision of future peculiar velocity data sets will require either equal smoothing of the predicted and measured velocity fields or careful accounting for the biases discussed here.

Measurements of the cosmological density parameter (Omega) using techniques that exploit the gravity-induced motions of galaxies yield, in linear perturbation theory, the degenerate parameter combination beta=Omega^{0.6}/b, where the linear bias parameter b is the ratio of the fluctuation amplitudes of the galaxy and mass distributions. However, the relation between the mass and the galaxy density fields depends on the complex physics of galaxy formation, and it can in general be non-linear, non-local, and stochastic. There is a growing consensus that the one-parameter model for bias is oversimplified. This leaves us with the following question: What is the quantity that is actually being measured by the different techniques? In order to address this question, we present estimates of beta from galaxy distributions constructed by applying a variety of biased galaxy formation models to cosmological N-body simulations. We compare the values of beta estimated using two different techniques: the anisotropy of the redshift-space power spectrum and an idealized version of the POTENT method. In most cases, we find that the bias factor b=Omega^{0.6}/beta derived from redshift-space anisotropy or POTENT is similar to the large-scale value of b_{sigma}(R) defined by the ratio of rms fluctuation amplitudes of the galaxy and mass distributions. However, non-linearity of bias and residual effects of non-linear gravitational evolution both influence beta estimates at the 10-20% level.

We examine the influence of the morphology-density(MD) relation and a wide range of simple models for biased galaxy formation on statistical measures of large scale structure. We contrast the behavior of local biasing models, in which the efficiency of galaxy formation is determined by density, geometry, or velocity dispersion of the local mass distribution, with that of non-local biasing models, in which galaxy formation is modulated coherently over scales larger than the galaxy correlation length. If morphological segregation of galaxies is governed by a local MD relation, then the correlation function of E/S0 galaxies should be steeper and stronger than that of spiral galaxies on small scales, as observed, while on large scales the correlation functions of E/S0 and spiral galaxies should have the same shape but different amplitudes. Similarly, all of our local bias models produce scale-independent amplification of the correlation function and power spectrum in the linear and mildly non-linear regimes; only a non-local biasing mechanism can alter the shape of the power spectrum on large scales. Moments of the biased galaxy distribution retain the hierarchical pattern of the mass moments, but biasing alters the values and scale-dependence of the hierarchical amplitudes S3 and S4. Pair-weighted moments of the galaxy velocity distribution are sensitive to the details of the biasing prescription. The non-linearity of the relation between galaxy density and mass density depends on the biasing prescription and the smoothing scale, and the scatter in this relation is a useful diagnostic of the physical parameters that determine the bias. Although the sensitivity of galaxy clustering statistics to the details of biasing is an obstacle to testing cosmological models, it is an asset for testing galaxy formation theories.

We measure the extinction law in a galaxy's spiral arm and interarm regions using a visual and infrared (BVRJHK) imaging study of the interacting galaxies NGC 2207 and IC 2163. This is an overlapping spiral galaxy pair in which NGC 2207 partially occults IC 2163. This geometry enables us to directly measure the extinction of light from the background galaxy as it passes through the disk of the foreground galaxy. We measure the extinction as a function of wavelength, and find that there is less extinction in the optical bands than expected from a normal Galactic extinction law. This deviation is significantly larger in the interarm region than in the spiral arm. The extinction curve in the spiral arm resembles a Milky Way $R_V=5.0$ dust model and the interarm extinction curve is flatter (``greyer'') still. We examine the effect of scattering of background galaxy light into the line of sight and find that it is negligible. We also examine the effect of an unresolved patchy dust distribution using a simple two-component dust model as well as the clumpy dust model of Witt & Gordon (1996). Both models clearly demonstrate that an unresolved patchy dust distribution can flatten the extinction curve significantly. When fit to the data, both models suggest that the observed difference between the arm and interarm extinction curves is caused by the interarm region of NGC 2207 having a higher degree of dust patchiness (density ratio between high-density and low-density phases) than the spiral arm region. We note that an unresolved patchy dust distribution will cause us to underestimate the average column depth of gas in a galaxy if based solely on the visual extinction. It is much better to use the infrared extinction for this purpose.