A joint analysis of data collected by the Planck and BICEP2+Keck teams has previously given $r = 0.09^{+0.06}_{-0.04}$ for BICEP2 and $r = 0.02^{+0.04}_{-0.02}$ for Keck. Analyzing BICEP2 using its published noise estimate, we had earlier (Colley & Gott 2015) found $r = 0.09 \pm 0.04$, agreeing with the final joint results for BICEP2. With the Keck data now available, we have done something the joint analysis did not: a correlation study of the BICEP2 vs. Keck B-mode maps. Knowing the correlation coefficient between the two and their amplitudes allows us to determine the noise in each map (which we check using the E-modes). We find the noise power in the BICEP2 map to be twice the original BICEP2 published estimate, explaining the anomalously high $r$ value obtained by BICEP2. We now find $r = 0.004 \pm 0.04$ for BICEP2 and $r = -0.01 \pm 0.04$ for Keck. Since $r \ge 0$ by definition, this implies a maximum likelihood value of $r = 0$, or no evidence for gravitational waves. Starobinsky Inflation ($r = 0.0036$) is not ruled out, however. Krauss & Wilzcek (2014) have already argued that "measurement of polarization of the CMB due to a long-wavelength stochastic background of gravitational waves from Inflation in the early Universe would firmly establish the quantization of gravity," and, therefore, the existence of gravitons. We argue it would also constitute a detection of gravitational Hawking radiation (explicitly from the causal horizons due to Inflation).

We derive constraints on cosmological parameters and tests of dark energy models from the combination of baryon acoustic oscillation (BAO) measurements with cosmic microwave background (CMB) and Type Ia supernova (SN) data. We take advantage of high-precision BAO measurements from galaxy clustering and the Ly-alpha forest (LyaF) in the BOSS survey of SDSS-III. BAO data alone yield a high confidence detection of dark energy, and in combination with the CMB angular acoustic scale they further imply a nearly flat universe. Combining BAO and SN data into an "inverse distance ladder" yields a 1.7% measurement of $H_0=67.3 \pm1.1$ km/s/Mpc. This measurement assumes standard pre-recombination physics but is insensitive to assumptions about dark energy or space curvature, so agreement with CMB-based estimates that assume a flat LCDM cosmology is an important corroboration of this minimal cosmological model. For open LCDM, our BAO+SN+CMB combination yields $\Omega_m=0.301 \pm 0.008$ and curvature $\Omega_k=-0.003 \pm 0.003$. When we allow more general forms of evolving dark energy, the BAO+SN+CMB parameter constraints remain consistent with flat LCDM. While the overall $\chi^2$ of model fits is satisfactory, the LyaF BAO measurements are in moderate (2-2.5 sigma) tension with model predictions. Models with early dark energy that tracks the dominant energy component at high redshifts remain consistent with our constraints. Expansion history alone yields an upper limit of 0.56 eV on the summed mass of neutrino species, improving to 0.26 eV if we include Planck CMB lensing. Standard dark energy models constrained by our data predict a level of matter clustering that is high compared to most, but not all, observational estimates. (Abridged)

We have analyzed the genus topology of the BICEP2 B-modes and find them to be Gaussian random phase as expected if they have a cosmological origin. These BICEP2 B-modes can be produced by gravity waves in the early universe, but some question has arisen as to whether these B-modes (for 50 < l < 120) may instead be produced by foreground polarized dust emission. The dust emission at 150 GHz observed by BICEP2 should be less in magnitude but have similar structure to that at 353 GHz. We have therefore calculated and mapped the B-modes in the BICEP2 region from the publicly available Q and U 353 GHz preliminary Planck polarization maps. These have a genus curve that is different from that seen in the BICEP2 observations, with features at different locations from those in the BICEP2 map. The two maps show a positive correlation coefficient of 15.2% +/- 3.9% (1-sigma). This requires the amplitude of the Planck (50 < l <120) dust modes to be low in the BICEP2 region, and the majority of the Planck 353 GHz signal in the BICEP2 region in these modes to be noise. We can explain the observed correlation coefficient of 15.2% with a BICEP2 gravity wave signal with an rms amplitude equal to 54% of the total BICEP2 rms amplitude. The gravity wave signal corresponds to a tensor-to-scalar ratio r = 0.11 +/- 0.04 (1-sigma). This is consistent with a gravity wave signal having been detected, at a 2.5-sigma level. The Planck and BICEP2 teams have recently engaged in joint analysis of their combined data|it will be interesting to see if that collaboration reaches similar conclusions.

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.

Numerous upcoming observations, such as WFIRST, BOSS, BigBOSS, LSST, Euclid, and Planck, will constrain dark energy (DE)'s equation of state with great precision. They may well find the ratio of pressure to energy density, $w$, is -1, meaning DE is equivalent to a cosmological constant. However, many time-varying DE models have also been proposed. A single parametrization to test a broad class of them and that is itself motivated by a physical picture is therefore desirable. We suggest the simplest model of DE has the same mechanism as inflation, likely a scalar field slowly rolling down its potential. If this is so, DE will have a generic equation of state and the Universe will have a generic dependence of the Hubble constant on redshift independent of the potential's starting value and shape. This equation of state and expression for the Hubble constant offer the desired model-independent but physically motivated parametrization, because they will hold for most of the standard scalar-field models of DE such as quintessence and phantom DE. Up until now two-parameter descriptions of $w$ have been available, but this work finds an additional approximation that leads to a single-parameter model. Using it, we conduct a $\chi^2$ analysis and find that experiments in the next seven years should be able to distinguish any of these time-varying DE models on the one hand from a cosmological constant on the other to 73% confidence if $w$ today differs from -1 by 3.5%. In the limit of perfectly accurate measurements of $\Omega_m$ and $H_0$, this confidence would rise to 96%. We also include discussion of the current status of DE experiment, a table compiling the techniques each will use, and tables of the precisions of the experiments for which this information was available at the time of publication.

We study the Physically Self Bound Cold Dark Matter Halo distribution which we associate with the massive galaxies within the Horizon Run 3 to estimate the accuracy in determination of the cosmological distance scale measured by the topology analysis. We apply the routine "Contour 3D" to 108 Mock Survey of $\pi$ steradians out to redshift z = 0.6, which effectively correspond to the SDSS-III BOSS survey, and compare the topology with that of a Gaussian Random Phase Field. We find that given three separate smoothing lengths $\lambda =$ 15, 21, and 34 $h^{-1}{\rm Mpc}$, the least $\chi^2$ fit genus per unit volume g yields a 1.7 % fractional uncertainty in smoothing length and angular diameter distance to $z = 0.6$. This is an improvement upon former calibrations of and presents a competitive error estimate with next BAO scale techniques. We also present three dimensional graphics of the Horizon Run 3 spherical mock survey to show a wealth of large-scale structures of the universe that are predicted in surveys like BOSS.

Large galaxy redshift surveys have long been used to constrain cosmological models and structure formation scenarios. In particular, the largest structures discovered observationally are thought to carry critical information on the amplitude of large-scale density fluctuations or homogeneity of the universe, and have often challenged the standard cosmological framework. The Sloan Great Wall (SGW) recently found in the Sloan Digital Sky Survey (SDSS) region casts doubt on the concordance cosmological model with a cosmological constant (i.e. the flat LCDM model). Here we show that the existence of the SGW is perfectly consistent with the LCDM model, a result that only our very large cosmological N-body simulation (the Horizon Run 2, HR2) could supply. In addition, we report on the discovery of a void complex in the SDSS much larger than the SGW, and show that such size of the largest void is also predicted in the LCDM paradigm. Our results demonstrate that an initially homogeneous isotropic universe with primordial Gaussian random phase density fluctuations growing in accordance with the General Relativity, can explain the richness and size of the observed large-scale structures in the SDSS. Using the HR2 simulation we predict that a future galaxy redshift survey about four times deeper or with 3 magnitude fainter limit than the SDSS should reveal a largest structure of bright galaxies about twice as big as the SGW.

We obtain constraints on cosmological parameters from the spherically averaged redshift-space correlation function of the CMASS Data Release 9 (DR9) sample of the Baryonic Oscillation Spectroscopic Survey (BOSS). We combine this information with additional data from recent CMB, SN and BAO measurements. Our results show no significant evidence of deviations from the standard flat-Lambda CDM model, whose basic parameters can be specified by Omega_m = 0.285 +- 0.009, 100 Omega_b = 4.59 +- 0.09, n_s = 0.96 +- 0.009, H_0 = 69.4 +- 0.8 km/s/Mpc and sigma_8 = 0.80 +- 0.02. The CMB+CMASS combination sets tight constraints on the curvature of the Universe, with Omega_k = -0.0043 +- 0.0049, and the tensor-to-scalar amplitude ratio, for which we find r < 0.16 at the 95 per cent confidence level (CL). These data show a clear signature of a deviation from scale-invariance also in the presence of tensor modes, with n_s <1 at the 99.7 per cent CL. We derive constraints on the fraction of massive neutrinos of f_nu < 0.049 (95 per cent CL), implying a limit of sum m_nu < 0.51 eV. We find no signature of a deviation from a cosmological constant from the combination of all datasets, with a constraint of w_DE = -1.033 +- 0.073 when this parameter is assumed time-independent, and no evidence of a departure from this value when it is allowed to evolve as w_DE(a) = w_0 + w_a (1 - a). The achieved accuracy on our cosmological constraints is a clear demonstration of the constraining power of current cosmological observations.

We present two large cosmological N-body simulations, called Horizon Run 2 (HR2) and Horizon Run 3 (HR3), made using 6000^3 = 216 billions and 7210^3 = 374 billion particles, spanning a volume of (7.200 Gpc/h)^3 and (10.815 Gpc/h)^3, respectively. These simulations improve on our previous Horizon Run 1 (HR1) up to a factor of 4.4 in volume, and range from 2600 to over 8800 times the volume of the Millennium Run. In addition, they achieve a considerably finer mass resolution, down to 1.25x10^11 M_sun/h, allowing to resolve galaxy-size halos with mean particle separations of 1.2 Mpc/h and 1.5 Mpc/h, respectively. We have measured the power spectrum, correlation function, mass function and basic halo properties with percent level accuracy, and verified that they correctly reproduce the LCDM theoretical expectations, in excellent agreement with linear perturbation theory. Our unprecedentedly large-volume N-body simulations can be used for a variety of studies in cosmology and astrophysics, ranging from large-scale structure topology, baryon acoustic oscillations, dark energy and the characterization of the expansion history of the Universe, till galaxy formation science - in connection with the new SDSS-III. To this end, we made a total of 35 all-sky mock surveys along the past light cone out to z=0.7 (8 from the HR2 and 27 from the HR3), to simulate the BOSS geometry. The simulations and mock surveys are already publicly available at http://astro.kias.re.kr/Horizon-Run23/.

We measure the acoustic scale from the angular power spectra of the Sloan Digital Sky Survey III (SDSS-III) Data Release 8 imaging catalog that includes 872,921 galaxies over ~ 10,000 deg^2 between 0.45<z<0.65. The extensive spectroscopic training set of the Baryon Oscillation Spectroscopic Survey (BOSS) luminous galaxies allows precise estimates of the true redshift distributions of galaxies in our imaging catalog. Utilizing the redshift distribution information, we build templates and fit to the power spectra of the data, which are measured in our companion paper, Ho et al. 2011, to derive the location of Baryon acoustic oscillations (BAO) while marginalizing over many free parameters to exclude nearly all of the non-BAO signal. We derive the ratio of the angular diameter distance to the sound horizon scale D_A/r_s= 9.212 + 0.416 -0.404 at z=0.54, and therefore, D_A= 1411+- 65 Mpc at z=0.54; the result is fairly independent of assumptions on the underlying cosmology. Our measurement of angular diameter distance D_A is 1.4 \sigma higher than what is expected for the concordance LCDM (Komatsu et al. 2011), in accordance to the trend of other spectroscopic BAO measurements for z >~ 0.35. We report constraints on cosmological parameters from our measurement in combination with the WMAP7 data and the previous spectroscopic BAO measurements of SDSS (Percival et al. 2010) and WiggleZ (Blake et al. 2011). We refer to our companion papers (Ho et al. 2011; de Putter et al. 2011) for investigations on information of the full power spectrum.

The Sloan Digital Sky Survey (SDSS) surveyed 14,555 square degrees, and delivered over a trillion pixels of imaging data. We present a study of galaxy clustering using 900,000 luminous galaxies with photometric redshifts, spanning between $z=0.45$ and $z=0.65$, constructed from the SDSS using methods described in Ross et al. (2011). This data-set spans 11,000 square degrees and probes a volume of $3h^{-3} \rm{Gpc}^3$, making it the largest volume ever used for galaxy clustering measurements. We present a novel treatment of the observational systematics and its applications to the clustering signals from the data set. In this paper, we measure the angular clustering using an optimal quadratic estimator at 4 redshift slices with an accuracy of ~15% with bin size of delta_l = 10 on scales of the Baryon Acoustic Oscillations (BAO) (at l~40-400). We derive cosmological constraints using the full-shape of the power-spectra. For a flat Lambda CDM model, when combined with Cosmic Microwave Background Wilkinson Microwave Anisotropy Probe 7 (WMAP7) and H_0 constraints from 600 Cepheids observed by HST, we find \Omega_\Lambda = 0.73 +/- 0.019 and H_0 to be 70.5 +/- 1.6 km/s/Mpc. For an open Lambda CDM model, when combined with WMAP7 + HST, we find $\Omega_K = 0.0035 +/- 0.0054, improved over WMAP7+HST alone by 40%. For a wCDM model, when combined with WMAP7+HST+SN, we find w = -1.071 +/- 0.078, and H_0 to be 71.3 +/- 1.7 km/s/Mpc, which is competitive with the latest large scale structure constraints from large spectroscopic surveys such as SDSS Data Release 7 (DR7) (Reid et al. 2010, Percival et al. 2010, Montesano et al. 2011) and WiggleZ (Blake et al. 2011). The SDSS-III Data Release 8 (SDSS-III DR8) Angular Clustering Data allows a wide range of investigations into the cosmological model, cosmic expansion (via BAO), Gaussianity of initial conditions and neutrino masses. (abridged)

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)

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

We present the second report of our systematic search for strongly lensed quasars from the data of the Sloan Digital Sky Survey (SDSS). From extensive follow-up observations of 136 candidate objects, we find 36 lenses in the full sample of 77,429 spectroscopically confirmed quasars in the SDSS Data Release 5. We then define a complete sample of 19 lenses, including 11 from our previous search in the SDSS Data Release 3, from the sample of 36,287 quasars with i<19.1 in the redshift range 0.6<z<2.2, where we require the lenses to have image separations of 1"<\theta<20" and i-band magnitude differences between the two images smaller than 1.25 mag. Among the 19 lensed quasars, 3 have quadruple-image configurations, while the remaining 16 show double images. This lens sample constrains the cosmological constant to be \Omega_\Lambda=0.84^{+0.06}_{-0.08}(stat.)^{+0.09}_{-0.07}(syst.) assuming a flat universe, which is in good agreement with other cosmological observations. We also report the discoveries of 7 binary quasars with separations ranging from 1.1" to 16.6", which are identified in the course of our lens survey. This study concludes the construction of our statistical lens sample in the full SDSS-I data set.

The use of standard rulers, such as the scale of the Baryonic Acoustic oscillations (BAO), has become one of the more powerful techniques employed in cosmology to probe the entity driving the accelerating expansion of the Universe. In this paper, the topology of large scale structure (LSS) is used as one such standard ruler to study this mysterious `dark energy'. By following the redshift evolution of the clustering of luminous red galaxies (LRGs) as measured by their 3D topology (counting structures in the cosmic web), we can chart the expansion rate and extract information about the equation of state of dark energy. Using the technique first introduced in (Park & Kim, 2009), we evaluate the constraints that can be achieved using 3D topology measurements from next-generation LSS surveys such as the Baryonic Oscillation Spectroscopic Survey (BOSS). In conjunction with the information that will be available from the Planck satellite, we find a single topology measurement on 3 different scales is capable of constraining a single dark energy parameter to within 5% and 10% when dynamics are permitted. This offers an alternative use of the data available from redshift surveys and serves as a cross-check for BAO studies.

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

In support of the new Sloan III survey, which will measure the baryon oscillation scale using the luminous red galaxies (LRGs), we have run the largest N-body simulation to date using $4120^3 = 69.9$ billion particles, and covering a volume of $(6.592 h^{-1} {\rm Gpc})^3$. This is over 2000 times the volume of the Millennium Run, and corner-to-corner stretches all the way to the horizon of the visible universe. LRG galaxies are selected by finding the most massive gravitationally bound, cold dark matter subhalos, not subject to tidal disruption, a technique that correctly reproduces the 3D topology of the LRG galaxies in the Sloan Survey.We have measured the covariance function, power spectrum, and the 3D topology of the LRG galaxy distribution in our simulation and made 32 mock surveys along the past lightcone to simulate the Sloan III survey. Our large N-body simulation is used to accurately measure the non-linear systematic effects such as gravitational evolution, redshift space distortion, past light cone space gradient, and galaxy biasing, and to calibrate the baryon oscillation scale and the genus topology. For example, we predict from our mock surveys that the baryon acoustic oscillation peak scale can be measured with the cosmic variance-dominated uncertainty of about 5% when the SDSS-III sample is divided into three equal volume shells, or about 2.6% when a thicker shell with $0.4<z<0.6$ is used. We find that one needs to correct the scale for the systematic effects amounting up to 5.2% to use it to constrain the linear theories. And the uncertainty in the amplitude of the genus curve is expected to be about 1% at 15 $h^{-1}$Mpc scale. We are making the simulation and mock surveys publicly available.

We measure the 3D genus topology of large scale structure using Luminous Red Galaxies (LRGs) in the Sloan Digital Sky Survey and find it consistent with the Gaussian random phase initial conditions expected from the simplest scenarios of inflation. This studies 3D topology on the largest scales ever obtained. The topology is sponge-like. We measure topology in two volume-limited samples: a dense shallow sample studied with smoothing length of 21h^{-1}Mpc, and a sparse deep sample studied with a smoothing length of 34h^{-1}Mpc. The amplitude of the genus curve is measured with 4% uncertainty. Small distortions in the genus curve expected from non-linear biasing and gravitational effects are well explained (to about 1-\sigma accuracy) by N-body simulations using a subhalo-finding technique to locate LRGs. This suggests the formation of LRGs is a clean problem that can be modeled well without any free fitting parameters. This bodes well for using LRGs to measure the characteristic scales such as the baryon oscillation scale in future deep redshift surveys.

Using the largest cosmological reionization simulation to-date (~24 billion particles), we use the genus curve to quantify the topology of neutral hydrogen distribution on scales > 1 Mpc as it evolves during cosmological reionization. We find that the reionization process proceeds primarily in an inside-out fashion, where higher density regions become ionized earlier than lower density regions. There are four distinct topological phases: (1) Pre-reionization at z ~ 15, when the genus curve is consistent with a Gaussian density distribution. (2) Pre-overlap at 10 < z < 15, during which the number of HII bubbles increases gradually with time, until percolation of HII bubbles starts to take effect, characterized by a very flat genus curve at high volume fractions. (3) Overlap at 8 < z < 10, when large HII bubbles rapidly merge, manifested by a precipitous drop in the amplitude of the genus curve. (4) Post-overlap at 6 < z < 8, when HII bubbles have mostly overlapped and the genus curve is consistent with a diminishing number of isolated neutral islands. After the end of reionization (z < 6), the genus of neutral hydrogen is consistent with Gaussian random phase, in agreement with observations.

We present a unified picture on the evolution of galaxy luminosity and morphology. Galaxy morphology is found to depend critically on the local environment set up by the nearest neighbor galaxy in addition to luminosity and the large scale density. When a galaxy is located farther than the virial radius from its closest neighbor, the probability for the galaxy to have an early morphological type is an increasing function only of luminosity and the local density due to the nearest neighbor ($\rho_n$). The tide produced by the nearest neighbor is thought to be responsible for the morphology transformation toward the early type at these separations. When the separation is less than the virial radius, i.e. when $\rho_n > \rho_{\rm virial}$, its morphology depends also on the neighbor's morphology and the large-scale background density over a few Mpc scales ($\rho_{20}$) in addition to luminosity and $\rho_n$. The early type probability keeps increasing as $\rho_n$ increases if its neighbor is an early type. But the probability decreases as $\rho_n$ increases when the neighbor is a late type. The cold gas streaming from the late type neighbor can be the reason for the morphology transformation toward late type. The overall early-type fraction increases as $\rho_{20}$ increases when $\rho_n > \rho_{\rm virial}$. This can be attributed to the hot halo gas of the neighbor which is confined by the pressure of the ambient medium held by the background mass. We have also found that galaxy luminosity depends on $\rho_n$, and that the isolated bright galaxies are more likely to be recent merger products. We propose a scenario that a series of morphology and luminosity transformation occur through distant interactions and mergers, which results in the morphology--luminosity--local density relation.

Tissot indicatrices have provided visual measures of local area and isotropy distortions. Here we show how large scale distortions of flexion (bending) and skewness (lopsidedness) can be measured. Area and isotropy distortions depend on the map projection metric, flexion and skewness, which manifest themselves on continental scales, depend on the first derivatives of the metric. We introduce new indicatrices that show not only area and isotropy distortions but flexion and skewness as well. We present a table showing error measures for area, isotropy, flexion, skewness, distances, and boundary cuts allowing us to compare different world map projections. We find that the Winkel-Tripel projection (already adopted for world maps by the National Geographic), has low distortion on most measures and excellent quality overall.

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

We have independently measured the genus topology of the temperature fluctuations in the cosmic microwave background seen in the Wilkinson Microwave Anisotropy Probe (WMAP) 3-year data. A genus analysis of the WMAP data indicates consistency with Gaussian random-phase initial conditions, as predicted by standard inflation. We set 95% confidence limits on non-linearities of -101 < f_{nl} < 107. We also find that the observed low l (l <= 8) modes show a slight anti-correlation with the Galactic foreground, but not exceeding 95% confidence, and that the topology defined by these modes is consistent with that of a Gaussian random-phase distribution (within 95% confidence).

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

It is useful to have mathematical criteria for evaluating errors in map projections. The Chebyshev criterion for minimizing rms (root mean square) local scale factor errors for conformal maps has been useful in developing conformal map projections of continents. Any local error criterion will be minimized ultimately by map projections with multiple interruptions, on which some pairs of points that are close on the globe are far apart on the map. Since it is as bad to have two points on the map at two times their proper separation as to have them at half their proper separation, it is the rms logarithmic distance, s, between random points in the mapped region that we will minimize. The best previously known projection of the entire sphere for distances is the Lambert equal-area azimuthal with an rms logarithmic distance error of s=0.343. For comparison, the Mercator has s=0.444, and the Mollweide has s=0.390. We present new projections: the "Gott equal-area elliptical" with perfect shapes on the central meridian, the "Gott-Mugnolo equal-area elliptical" and the "Gott-Mugnolo azimuthal" with rms logarithmic distance errors of s=0.365, s=0.348, and s=0.341 respectively, which improve on previous projections of their type. The "Gott-Mugnolo azimuthal" has the lowest distance errors of any map and is produced by a new technique using "forces" between pairs of points on a map which make them move so as to minimize s. The "Gott equal-area elliptical" projection produces a particularly attractive map of Mars, and the "Gott-Mugnolo azimuthal" projection produces an interesting map of the moon.

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] We have produced a cleaned map of the Wilkinson Microwave Anisotropy Probe (WMAP) 3-year data using an improved foreground subtraction technique. We perform an internal linear combination (ILC) to subtract the Galactic foreground emission from the temperature fluctuations observed by the WMAP. We divide the whole sky into hundreds of pixel groups with similar foreground spectral indices over a range of WMAP frequencies, apply the ILC for each group, and obtain a CMB map with foreground emission effectively reduced. With the resulting foreground-reduced ILC map, we have investigated the known anomalies in CMB maps at large scales. The quadrupole and the octopole powers measured from our ILC map are \delta T_2^2 = 276_{-126}^{+94} uK^2 and \delta T_3^2 = 952_{-83}^{+64} uK^2, respectively. The measured quadrupole power is lower than the value expected in the concordance LCDM model (1250 uK^2), in which the probability of finding a quadrupole power lower than the measured value is 5.7%. We have confirmed that the quadrupole and the octopole are strongly aligned with angle \theta_{23} = 11.8_{-8.0}^{+6.4} degree, and are planar with high planarity parameters t=0.98_{-0.02}^{+0.02} for l=2 and t=0.91_{-0.03}^{+0.02} for l=3. The observed angular separation \theta_{23} is marginally statistically significant because the probability of finding the angular separation as low as the observed value is 4.3%. However, the observed planarity is not statistically significant. The probability of observing such a planarity as high as the measured t values is over 18%. The ILC simulations show that the residual foreground emission in the ILC map does not affect the estimated values significantly.

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

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

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

We have studied the dependence of topology of large scale structure on tracer, gravitational evolution, redshift space distortion, and cosmology. A series of large N-body simulations of the $\Lambda$CDM and SCDM models that have evolved 1.1 or 8.6 billion particles, are used in the study. Evolution of the genus statistic, used as a topology measure, from redshift 8 to 0 is accurately calculated over a wide range of smoothing scales using the simulations. The tracers of large scale structure considered are the CDM matter, biased peaks in the initial density field, dark halos, and `galaxies' populating the dark halos in accordance with a Halo Occupation Distribution (HOD) model. We have found that the effects of biasing, gravitational evolution, and initial conditions on topology of large scale structure are all comparable. The redshift space distortion effects are relatively small down to about 5 $h^{-1}$Mpc for all tracers except for the high threshold part of the genus curve. The gravitational effects are found to be well-modeled by the analytic perturbation theory when the CDM matter distribution is considered. But the direction of gravitational evolution of topology can be even reversed for different tracers. For example, the shift parameter of the genus curve evolves in opposite directions for matter and HOD `galaxies' at large scales. At small scales there are interesting deviations of the genus curve of dark halos and `galaxies' from that of matter in our initially Gaussian simulations. The deviations should be understood as due to combined effects of gravitational evolution and biasing.

The current standard theory of the origin of the Moon is that the Earth was hit by a giant impactor the size of Mars causing ejection of iron poor impactor mantle debris that coalesced to form the Moon. But where did this Mars-sized impactor come from? Isotopic evidence suggests that it came from 1AU radius in the solar nebula and computer simulations are consistent with it approaching Earth on a zero-energy parabolic trajectory. But how could such a large object form in the disk of planetesimals at 1AU without colliding with the Earth early-on before having a chance to grow large or before its or the Earth's iron core had formed? We propose that the giant impactor could have formed in a stable orbit among debris at the Earth's Lagrange point $L_4$ (or $L_5$). We show such a configuration is stable, even for a Mars-sized impactor. It could grow gradually by accretion at $L_4$ (or $L_5$), but eventually gravitational interactions with other growing planetesimals could kick it out into a chaotic creeping orbit which we show would likely cause it to hit the Earth on a zero-energy parabolic trajectory. This paper argues that this scenario is possible and should be further studied.

Since no one lives forever, all a life preserver can really do is prolong life for longer than would have otherwise been the case. With this rather limited definition in mind we explore in this paper whether in principle you can take a life preserver with you to protect you (for a while at least) against the tidal forces encountered on a trip inside a black hole.

We construct the error distribution of Hubble constant ($H_0$) measurements from Huchra's compilation of 461 measurements of $H_0$ and the WMAP experiment central value $H_0$ = 71 km s$^{-1}$ Mpc$^{-1}$. This error distribution is non-Gaussian, with significantly larger probability in the tails of the distribution than predicted by a Gaussian distribution. The 95.4 % confidence limits are 7.0 $\sigma$ in terms of the quoted errors. It is remarkably well described by either a widened $n = 2$ Student's $t$ distribution or a widened double exponential distribution. These conclusions are unchanged if we use instead the central value $H_0$ = 67 km s$^{-1}$ Mpc$^{-1}$ found from a median statistics analysis of a major subset of $H_0$ measurements used here.

We have independently measured the genus topology of the temperature fluctuations in the cosmic microwave background seen by the Wilkinson Microwave Anisotropy Probe (WMAP). A genus analysis of the WMAP data indicates consistency with Gaussian random-phase initial conditions, as predicted by standard inflation.

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

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

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

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

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

We have studied topology of the distribution of the high redshift galaxies identified in the Hubble Deep Field (HDF) North and South. The two-dimensional genus is measured from the projected distributions of the HDF galaxies at angular scales from $3.8''$ to $ 6.1''$. We have also divided the samples into three redshift slices with roughly equal number of galaxies using photometric redshifts to see possible evolutionary effects on the topology. The genus curve of the HDF North clearly indicates clustering of galaxies over the Poisson distribution while the clustering is somewhat weaker in the HDF South. This clustering is mainly due to the nearer galaxies in the samples. We have also found that the genus curve of galaxies in the HDF is consistent with the Gaussian random phase distribution with no significant redshift dependence.

We have simulated cosmic microwave background (CMB) anisotropy maps for several COBE-DMR-normalized cold dark matter (CDM) cosmogonies, to make predictions for the upcoming MAP experiment. We have studied the sensitivity of the simulated MAP data to cosmology, sky coverage, and instrumental noise. With accurate knowledge of instrumental noise, MAP data will discriminate among the cosmogonies considered, and determine the topology of the initial fluctuations. A correlation function analysis of the simulated MAP data results in a very accurate measurement of the acoustic Hubble radius at decoupling. A low-density open CDM model with Omega_0=0.4 can be distinguished from the Omega_0=1 fiducial CDM model or a Lambda CDM model with > 99% confidence from the location of the acoustic "valley" in the correlation function. A genus analysis of the simulated MAP data indicates that in cosmogonies with Gaussian random-phase initial conditions, a shift of the zero-crossing point of the genus curve near the mean temperature threshold level nu=0 should not exceed |Delta nu| = 0.01 (0.04) when the total effective FWHM smoothing is 0.3 deg (1.0 deg). The asymmetry of the genus curve at the positive and negative threshold levels should not exceed |Delta g/g(nu=1,-1)| = 0.8% (4%) at 0.3 deg (1.0 deg) FWHM smoothing. Deviations of the observed MAP data in excess of these small values will be evidence for non-Gaussian behavior. The amplitude of the genus curve is a measure of the shape of the power spectrum at the smoothing scale. Even with the expected amount of instrumental noise and partial sky coverage (due to the Galaxy), the MAP data should allow discrimination amongst the cosmogonies considered at more than 99% confidence solely from a genus amplitude analysis.

In a flat, k=0 cosmology with galaxies that approximate singular isothermal spheres, gravitational lens image separations should be uncorrelated with source redshift. But in an open k=-1 cosmology such gravitational lens image separations become smaller with increasing source redshift. The observed separations do become smaller with increasing source redshift but the effect is even stronger than that expected in an Omega=0 cosmology. The observations are thus not compatible with the "standard" gravitational lensing statistics model in a flat universe. We try various open and flat cosmologies, galaxy mass profiles, galaxy merging and evolution models, and lensing aided by clusters to explain the correlation. We find the data is not compatible with any of these possibilities within the 95% confidence limit, leaving us with a puzzle. If we regard the observed result as a statistical fluke, it is worth noting that we are about twice as likely to observe it in an open universe (with 0<Omega<0.4) as we are to observe it in a flat one. Finally, the existence of an observed multiple image lens system with a source at z=4.5 places a lower limit on the deceleration parameter: q_0 > -2.0.