Scientific Purpose      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³ = 69.9 billion particles, and covering a volume of (6.592 Gpc/h)³. 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 halos, 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 light cone 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 Mpc/h scale. We are making the simulation and mock surveys publicly available. Published in ApJ (2009, vol. 701, pp. 1547-1559) [click] For the latest preprint version (arXiv:0812.1392) [click].      Authors      Juhan Kim at CAC of KIAS (email contact: kjhan@kias.re.kr)      Changbom Park at Korea Institute for Advanced Study      J.Richard Gott III at Princeton University      John Dubinski at University of Toronto Simulation Specifications cosmological model Ω(matter) Ω(baryon) Ω(Λ) spectral index H0 (Hubble constant) σ8 box Size on a side number of grids for initial conditions number of CDM particles starting redshift number of global time steps ΛCDM of WMAP 5-year 0.26 0.044 0.74 0.96 0.72 1/1.26 6592 Mpc/h 4120³ 4120³ 23 400 Outputs from the Simulation Snapshot data at z= 0, 0.1, 0.3, 0.5, 0.7 and 1.0 Mock SDSS-III surveys: eight non-overlapping all-sky (4π) surveys out to z=0.6 A list of subhalos from all-sky (4π) survey out to z=1.29 (partly overlapping with the above 8 surveys) Past lightcone slice: a wedge of Δφ=45° and 64Mpc/h width reaching the horizon Past lightcone slice-II: a wedge of Δφ=5° and Δθ=10° reaching z=23 Pictures & Movies (downloadable)      Data Archive (downloadable) Tables (downloadable) Paper & Press Release