Halo Occupation Distribution Parameters Research Articles

Global analysis of luminosity- and colour-dependent galaxy clustering in the Sloan Digital Sky Survey

ABSTRACT We present a halo occupation distribution (HOD) analysis of the luminosity- and colour-dependent galaxy clustering in the Sloan Digital Sky Survey. A novelty of our technique is that it uses a combination of clustering measurements in luminosity bins to perform a global likelihood analysis, simultaneously constraining the HOD parameters for a range of luminosity thresholds. We present simple, smooth fitting functions which accurately describe the resulting luminosity dependence of the best-fitting HOD parameters. To minimize systematic halo modelling effects, we use theoretical halo 2-point correlation functions directly measured and tabulated from a suite of N-body simulations spanning a large enough dynamic range in halo mass and spatial separation. Thus, our modelling correctly accounts for non-linear and scale-dependent halo bias as well as any departure of halo profiles from universality, and we additionally account for halo exclusion using the hard sphere approximation. Using colour-dependent clustering information, we constrain the satellite galaxy red fraction in a model-independent manner which does not rely on any group-finding algorithm. We find that the resulting luminosity dependence of the satellite red fraction is significantly shallower than corresponding measurements from galaxy group catalogues, and we provide a simple fitting function to describe this dependence. Our fitting functions are readily usable in generating low-redshift mock galaxy catalogues, and we discuss some potentially interesting applications as well as possible extensions of our technique.

Large-scale redshift space distortions in modified gravity theories

Measurements of redshift space distortions (RSD) provide a means to test models of gravity on large-scales. We use mock galaxy catalogues constructed from large N-body simulations of standard and modified gravity models to measure galaxy clustering in redshift space. We focus our attention on two of the most representative and popular families of modified gravity models: the Hu \& Sawicki $f(R)$ gravity and the normal branch of the DGP model. The galaxy catalogues are built using a halo occupation distribution (HOD) prescription with the HOD parameters in the modified gravity models tuned to match with the number density and the real-space clustering of {\sc boss-cmass} galaxies. We employ two approaches to model RSD: the first is based on linear perturbation theory and the second models non-linear effects on small-scales by assuming standard gravity and including biasing and RSD effects. We measure the monopole to real-space correlation function ratio, the quadrupole to monopole ratio, clustering wedges and multipoles of the correlation function and use these statistics to find the constraints on the distortion parameter, $\beta$. We find that the linear model fails to reproduce the N-body simulation results and the true value of $\beta$ on scales $s < 40\Mpch$, while the non-linear modelling of RSD recovers the value of $\beta$ on the scales of interest for all models. RSD on large scales ($s\gtrsim20$-$40\Mpch$) have been found to show significant deviations from the prediction of standard gravity in the DGP models. However, the potential to use RSD to constrain $f(R)$ models is less promising, due to the different screening mechanism in this model.

Emulating galaxy clustering and galaxy–galaxy lensing into the deeply non-linear regime: methodology, information, and forecasts

The combination of galaxy-galaxy lensing (GGL) with galaxy clustering is one of the most promising routes to determining the amplitude of matter clustering at low redshifts. We show that extending clustering+GGL analyses from the linear regime down to $\sim 0.5 \, h^{-1}$ Mpc scales increases their constraining power considerably, even after marginalizing over a flexible model of non-linear galaxy bias. Using a grid of cosmological N-body simulations, we construct a Taylor-expansion emulator that predicts the galaxy autocorrelation $\xi_{\text{gg}}(r)$ and galaxy-matter cross-correlation $\xi_{\text{gm}}(r)$ as a function of $\sigma_8$, $\Omega_m$, and halo occupation distribution (HOD) parameters, which are allowed to vary with large scale environment to represent possible effects of galaxy assembly bias. We present forecasts for a fiducial case that corresponds to BOSS LOWZ galaxy clustering and SDSS-depth weak lensing (effective source density $\sim 0.3$ arcmin$^{-2}$). Using tangential shear and projected correlation function measurements over $0.5 \leq r_p \leq 30 \, h^{-1}$ Mpc yields a 1.8% constraint on the parameter combination $\sigma_8\Omega_m^{0.58}$, a factor of two better than a constraint that excludes non-linear scales ($r_p > 2 \, h^{-1}$ Mpc, $4 \, h^{-1}$ Mpc for $\gamma_t,w_p$). Much of this improvement comes from the non-linear clustering information, which breaks degeneracies among HOD parameters that would otherwise degrade the inference of matter clustering from GGL. Increasing the effective source density to $3$ arcmin$^{-2}$ sharpens the constraint on $\sigma_8\Omega_m^{0.58}$ by a further factor of two. With robust modeling into the non-linear regime, low-redshift measurements of matter clustering at the 1-percent level with clustering+GGL alone are well within reach of current data sets such as those provided by the Dark Energy Survey.

The Santiago–Harvard–Edinburgh–Durham void comparison – I. SHEDding light on chameleon gravity tests

We present a systematic comparison of several existing and new void finding algorithms, focusing on their potential power to test a particular class of modified gravity models - chameleon $f(R)$ gravity. These models deviate from standard General Relativity (GR) more strongly in low-density regions and thus voids are a promising venue to test them. We use Halo Occupation Distribution (HOD) prescriptions to populate haloes with galaxies, and tune the HOD parameters such that the galaxy two-point correlation functions are the same in both f(R) and GR models. We identify both 3D voids as well as 2D underdensities in the plane-of-the-sky to find the same void abundance and void galaxy number density profiles across all models, which suggests that they do not contain much information beyond galaxy clustering. However, the underlying void dark matter density profiles are significantly different, with f(R) voids being more underdense than GR ones, which leads to f(R) voids having a larger tangential shear signal than their GR analogues. We investigate the potential of each void finder to test f(R) models with near-future lensing surveys such as EUCLID and LSST. The 2D voids have the largest power to probe f(R) gravity, with a LSST analysis of tunnel (which is a new type of 2D underdensity introduced here) lensing distinguishing at 80 and 11$\sigma$ (statistical error) f(R) models with $|f_{R0}|=10^{-5}$ and $10^{-6}$ from GR.

Galaxy–galaxy weak gravitational lensing in f(R) gravity

We present an analysis of galaxy-galaxy weak gravitational lensing (GGL) in chameleon $f(R)$ gravity - a leading candidate of non-standard gravity models. For the analysis we have created mock galaxy catalogues based on dark matter haloes from two sets of numerical simulations, using a halo occupation distribution (HOD) prescription which allows a redshift dependence of galaxy number density. To make a fairer comparison between the $f(R)$ and $\Lambda$CDM models, their HOD parameters are tuned so that the galaxy two-point correlation functions in real space (and therefore the projected two-point correlation functions) match. While the $f(R)$ model predicts an enhancement of the convergence power spectrum by up to $\sim30\%$ compared to the standard $\Lambda$CDM model with the same parameters, the maximum enhancement of GGL is only half as large and less than 5\% on separations above $\sim1$-$2h^{-1}$Mpc, because the latter is a cross correlation of shear (or matter, which is more strongly affected by modified gravity) and galaxy (which is weakly affected given the good match between galaxy auto correlations in the two models) fields. We also study the possibility of reconstructing the matter power spectrum by combination of GGL and galaxy clustering in $f(R)$ gravity. We find that the galaxy-matter cross correlation coefficient remains at unity down to $\sim2$-$3h^{-1}$Mpc at relevant redshifts even in $f(R)$ gravity, indicating joint analysis of GGL and galaxy clustering can be a powerful probe of matter density fluctuations in chameleon gravity. The scale dependence of the model differences in their predictions of GGL can potentially allow to break the degeneracy between $f(R)$ gravity and other cosmological parameters such as $\Omega_m$ and $\sigma_8$.

Angular power spectrum of galaxies in the 2MASS Redshift Survey

We present the measurement and interpretation of the angular power spectrum of nearby galaxies in the 2MASS Redshift Survey catalog with spectroscopic redshifts up to $z\approx 0.1$. We detect the angular power spectrum up to a multipole of $\ell\approx 1000$. We find that the measured power spectrum is dominated by galaxies living inside nearby galaxy clusters and groups. We use the halo occupation distribution (HOD) formalism to model the power spectrum, obtaining a fit with reasonable parameters. These HOD parameters are in agreement with the 2MASS galaxy distribution we measure toward the known nearby galaxy clusters, confirming validity of our analysis.

Approximate Bayesian computation in large-scale structure: constraining the galaxy–halo connection

Standard approaches to Bayesian parameter inference in large scale structure assume a Gaussian functional form (chi-squared form) for the likelihood. This assumption, in detail, cannot be correct. Likelihood free inferences such as Approximate Bayesian Computation (ABC) relax these restrictions and make inference possible without making any assumptions on the likelihood. Instead ABC relies on a forward generative model of the data and a metric for measuring the distance between the model and data. In this work, we demonstrate that ABC is feasible for LSS parameter inference by using it to constrain parameters of the halo occupation distribution (HOD) model for populating dark matter halos with galaxies. Using specific implementation of ABC supplemented with Population Monte Carlo importance sampling, a generative forward model using HOD, and a distance metric based on galaxy number density, two-point correlation function, and galaxy group multiplicity function, we constrain the HOD parameters of mock observation generated from selected "true" HOD parameters. The parameter constraints we obtain from ABC are consistent with the "true" HOD parameters, demonstrating that ABC can be reliably used for parameter inference in LSS. Furthermore, we compare our ABC constraints to constraints we obtain using a pseudo-likelihood function of Gaussian form with MCMC and find consistent HOD parameter constraints. Ultimately our results suggest that ABC can and should be applied in parameter inference for LSS analyses.

Analytic halo approach to the bispectrum of galaxies in redshift space

We present an analytic formula for the galaxy bispectrum in redshift space on the basis of the halo approach description with the halo occupation distribution of central galaxies and satellite galaxies. This work is an extension of a previous work on the galaxy power spectrum, which illuminated the significant contribution of satellite galaxies to the higher multipole spectrum through the non-linear redshift space distortions of their random motions. Behaviors of the multipoles of the bispectrum are compared with results of numerical simulations assuming a halo occupation distribution of the LOWZ sample of the SDSS-III BOSS survey. Also presented are analytic approximate formulas for the multipoles of the bispectrum, which is useful to understanding their characteristic properties. We demonstrate that the Fingers of God effect is quite important for the higher multipoles of the bispectrum in redshift space, depending on the halo occupation distribution parameters.

Constraining the X-ray AGN halo occupation distribution: implications foreROSITA

The X-ray emission from active galactic nucleus (AGN) is a major component of extragalactic X-ray sky. In this paper, we use the X-ray luminosity function (XLF) and halo occupation distribution (HOD) formalism to construct a halo model for the X-ray emission from AGNs. Verifying that the two inputs (XLF and HOD) are in agreement with each other, we compute the auto-correlation power spectrum in the soft X-ray band (0.5-2 keV) due to the AGNs potentially resolved by eROSITA (extended ROentgen Survey with an Imaging Telescope Array) mission and explore the redshift and mass dependence of the power spectrum. Studying the relative contribution of the Poisson and the clustering terms to the total power, we find that at multipoles $l\lesssim 1000$ (i.e. large scales), the clustering term is larger than the Poisson term. We also forecast the potential of X-ray auto-correlation power spectrum and X-ray-lensing cross-correlation power spectrum using eROSITA and eROSITA-LSST (Large Synoptic Survey Telescope) surveys, respectively, to constrain the HOD parameters and their redshift evolution. In addition, we compute the power spectrum of the AGNs lying below the flux resolution limit of eROSITA, which is essential to understand in order to extract the X-ray signal from the hot diffuse gas present in galaxies and clusters.

Galaxy–galaxy lensing in the Dark Energy Survey Science Verification data

We present galaxy-galaxy lensing results from 139 square degrees of Dark Energy Survey (DES) Science Verification (SV) data. Our lens sample consists of red galaxies, known as redMaGiC, which are specifically selected to have a low photometric redshift error and outlier rate. The lensing measurement has a total signal-to-noise of 29 over scales $0.09 < R < 15$ Mpc/$h$, including all lenses over a wide redshift range $0.2 < z < 0.8$. Dividing the lenses into three redshift bins for this constant moving number density sample, we find no evidence for evolution in the halo mass with redshift. We obtain consistent results for the lensing measurement with two independent shear pipelines, ngmix and im3shape. We perform a number of null tests on the shear and photometric redshift catalogs and quantify resulting systematic uncertainties. Covariances from jackknife subsamples of the data are validated with a suite of 50 mock surveys. The results and systematics checks in this work provide a critical input for future cosmological and galaxy evolution studies with the DES data and redMaGiC galaxy samples. We fit a Halo Occupation Distribution (HOD) model, and demonstrate that our data constrains the mean halo mass of the lens galaxies, despite strong degeneracies between individual HOD parameters.

Joint analysis of galaxy-galaxy lensing and galaxy clustering: Methodology and forecasts for Dark Energy Survey

The joint analysis of galaxy-galaxy lensing and galaxy clustering is a promising method for inferring the growth function of large scale structure. This analysis will be carried out on data from the Dark Energy Survey (DES), with its measurements of both the distribution of galaxies and the tangential shears of background galaxies induced by these foreground lenses. We develop a practical approach to modeling the assumptions and systematic effects affecting small scale lensing, which provides halo masses, and large scale galaxy clustering. Introducing parameters that characterize the halo occupation distribution (HOD), photometric redshift uncertainties, and shear measurement errors, we study how external priors on different subsets of these parameters affect our growth constraints. Degeneracies within the HOD model, as well as between the HOD and the growth function, are identified as the dominant source of complication, with other systematic effects sub-dominant. The impact of HOD parameters and their degeneracies necessitate the detailed joint modeling of the galaxy sample that we employ. We conclude that DES data will provide powerful constraints on the evolution of structure growth in the universe, conservatively/optimistically constraining the growth function to 7.9\%/4.8\% with its first-year data that covered over 1000 square degrees, and to 3.9\%/2.3\% with its full five-year data that will survey 5000 square degrees, including both statistical and systematic uncertainties.

A pseudo-spectrum analysis of galaxy–galaxy lensing

We present the application of the pseudo-spectrum method to galaxy-galaxy lensing. We derive explicit expressions for the pseudo-spectrum analysis of the galaxy-shear cross spectrum, which is the Fourier space counterpart of the stacked galaxy-galaxy lensing profile. The pseudo-spectrum method corrects observational issues such as the survey geometry, masks of bright stars and their spikes, and inhomogeneous noise, which distort the spectrum and also mix the E-mode and the B-mode signals. Using ray-tracing simulations in N-body simulations including realistic masks, we confirm that the pseudo-spectrum method successfully recovers the input galaxy-shear cross spectrum. We also investigate the covariance of the galaxy-shear cross spectrum using the ray-tracing simulations to show that there is an excess covariance relative to the Gaussian covariance at small scales where the shot noise is dominated in the Gaussian approximation. We find that the excess of the covariance is consistent with the expectation from the halo sample variance (HSV), which originates from the matter fluctuations at scales larger than the survey area. We apply the pseudo-spectrum method to the observational data of Canada-France-Hawaii Telescope Lensing survey (CFHTLenS) shear catalogue and three different spectroscopic samples of Sloan Digital Sky Survey Luminous Red Galaxy (SDSS LRG), and Baryon Oscillation Spectroscopic Survey (BOSS) CMASS and LOWZ galaxies. The galaxy-shear cross spectra are significantly detected at the level of 7-10$\sigma$ using the analytic covariance with the HSV contribution included. We also confirm that the observed spectra are consistent with the halo model predictions with the halo occupation distribution parameters estimated from previous work. This work demonstrates the viability of galaxy-galaxy lensing analysis in the Fourier space.

Introducing decorated HODs: modelling assembly bias in the galaxy–halo connection

The connection between galaxies and dark matter halos is often inferred from data using probabilistic models, such as the Halo Occupation Distribution (HOD). Conventional HOD formulations assume that only halo mass governs the galaxy-halo connection. Violations of this assumption, known as galaxy assembly bias, threaten the HOD program. We introduce decorated HODs, a new, flexible class of models designed to account for assembly bias. Decorated HODs minimally expand the parameter space and maximize the independence between traditional and novel HOD parameters. We use decorated HODs to quantify the influence of assembly bias on clustering and lensing statistics. For SDSS-like samples, the impact of assembly bias on galaxy clustering can be as large as a factor of two on r ~ 200 kpc scales and ~15% in the linear regime. Assembly bias can either enhance or diminish clustering on large scales, but generally increases clustering on scales r <~ 1 Mpc. We performed our calculations with Halotools, an open-source, community-driven python package for studying the galaxy-halo connection (http://halotools.readthedocs.org). We conclude by describing the use of decorated HODs to treat assembly bias in otherwise conventional likelihood analyses.

COSMIC EMULATION: FAST PREDICTIONS FOR THE GALAXY POWER SPECTRUM

The halo occupation distribution (HOD) approach has proven to be an effective method for modeling galaxy clustering and bias. In this approach, galaxies of a given type are probabilistically assigned to individual halos in N-body simulations. In this paper, we present a fast emulator for predicting the fully nonlinear galaxy power spectrum over a range of freely specifiable HOD modeling parameters. The emulator is constructed using results from 100 HOD models run on a large LCDM N-body simulation, with Gaussian Process interpolation applied to a PCA-based representation of the galaxy power spectrum. The total error is currently ~3% (~2% in the simulation and ~1% in the emulation process) from z=1 to z=0, over the considered parameter range. We use the emulator to investigate parametric dependencies in the HOD model, as well as the behavior of galaxy bias as a function of HOD parameters. The emulator is publicly available at http://www.hep.anl.gov/cosmology/CosmicEmu/emu.html.

THE WEAK LENSING SIGNAL AND THE CLUSTERING OF BOSS GALAXIES. II. ASTROPHYSICAL AND COSMOLOGICAL CONSTRAINTS

We perform a joint analysis of the abundance, the clustering and the galaxy-galaxy lensing signal of galaxies measured from Data Release 11 of the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey (SDSS III-BOSS) in our companion paper, Miyatake et al. (2014). The lensing signal was obtained by using the shape catalog of background galaxies from the Canada France Hawaii Telescope Legacy Survey, which was made publicly available by the CFHTLenS collaboration, with an area overlap of about 105 deg$^2$. We analyse the data in the framework of the halo model in order to fit halo occupation parameters and cosmological parameters ($\Omega_{\rm m}$ and $\sigma_8$) to these observables simultaneously, and thus break the degeneracy between galaxy bias and cosmology. Adopting a flat $\Lambda$CDM cosmology with priors on $\Omega_b h^2$, $n_{\rm s}$ and $h$ from the analysis of WMAP 9-year data, we obtain constraints on the stellar mass-halo mass relation of galaxies in our sample. Marginalizing over the halo occupation distribution parameters and a number of other nuisance parameters in our model, we obtain $\Omega_{\rm m}=0.310^{+0.019}_{-0.020}$ and $\sigma_8=0.785^{+0.044}_{-0.044}$ (68% confidence). We demonstrate the robustness of our results with respect to sample selection and a variety of systematics such as the halo off-centering effect and possible incompleteness in our sample. Our constraints are consistent, complementary and competitive with those obtained using other independent probes of these cosmological parameters. The cosmological analysis is the first of its kind to be performed at a redshift as high as $0.53$.

Constraining halo occupation distribution and cosmic growth rate using multipole power spectrum

Abstract We propose a new method of measuring halo occupation distribution (HOD) together with cosmic growth rate using multipole components of galaxy power spectrum Pl(k). The non-linear redshift-space distortion due to the random motion of satellite galaxies, i.e. Fingers-of-God, generates high-l multipole anisotropy in galaxy clustering, such as the hexadecapole (l = 4) and tetra-hexadecapole (l = 6), which are sensitive to the fraction and the velocity dispersion of satellite galaxies. Using simulated samples following the HOD of luminous red galaxies, we find that the input HOD parameters are successfully reproduced from Pl(k), and that high-l multipole information help to break the degeneracy among HOD parameters. We also show that the measurements of the cosmic growth rate as well as the satellite fraction and velocity dispersions are significantly improved by adding the small-scale information of high-l multipoles.

Cosmology with massive neutrinos I: towards a realistic modeling of the relation between matter, haloes and galaxies

By using a suite of large box-size N-body simulations that incorporate massive neutrinos as an extra set of particles, with total masses of 0.15, 0.30, and 0.60 eV, we investigate the impact of neutrino masses on the spatial distribution of dark matter haloes and on the distribution of galaxies within the haloes. We compute the bias between the spatial distribution of dark matter haloes and the overall matter and cold dark matter distributions using statistical tools such as the power spectrum and the two-point correlation function. Overall we find a scale-dependent bias on large scales for the cosmologies with massive neutrinos. In particular, we find that the bias decreases with the scale, being this effect more important for higher neutrino masses and at high redshift. However, our results indicate that the scale-dependence in the bias is reduced if the latter is computed with respect to the cold dark matter distribution only. We find that the value of the bias on large scales is reasonably well reproduced by the Tinker fitting formula once the linear cold dark matter power spectrum is used, instead of the total matter power spectrum. We also investigate whether scale-dependent bias really comes from purely neutrino's effect or from nonlinear gravitational collapse of haloes. For this purpose, we address the Ων-σ8 degeneracy and find that such degeneracy is not perfect, implying that neutrinos imprint a slight scale dependence on the large-scale bias. Finally, by using a simple halo occupation distribution (HOD) model, we investigate the impact of massive neutrinos on the distribution of galaxies within dark matter haloes. We use the main galaxy sample in the Sloan Digital Sky Survey (SDSS) II Data Release 7 to investigate if the small-scale galaxy clustering alone can be used to discriminate among different cosmological models with different neutrino masses. Our results suggest that different choices of the HOD parameters can reproduce the observational measurements relatively well, and we quantify the difference between the values of the HOD parameters between massless and massive neutrino cosmologies.

Non-Gaussianity of the cosmic infrared background anisotropies – I. Diagrammatic formalism and application to the angular bispectrum

We present the first halo model based description of the Cosmic Infrared Background (CIB) non-Gaussianity (NG) that is fully parametric. To this end, we introduce, for the first time, a diagrammatic method to compute high order polyspectra of the 3D galaxy density field. It allows an easy derivation and visualisation of the different terms of the polyspectrum. We apply this framework to the power spectrum and bispectrum, and we show how to project them on the celestial sphere in the purpose of the application to the CIB angular anisotropies. Furthermore, we show how to take into account the particular case of the shot noise terms in that framework. Eventually, we compute the CIB angular bispectrum at 857 GHz and study its scale and configuration dependencies, as well as its variations with the halo occupation distribution parameters. Compared to a previously proposed empirical prescription, such physically motivated model is required to describe fully the CIB anisotropies bispectrum. Finally, we compare the CIB bispectrum with the bispectra of other signals potentially present at microwave frequencies, which hints that detection of CIB NG should be possible above 220 GHz.

THE HALO OCCUPATION DISTRIBUTION OF X-RAY-BRIGHT ACTIVE GALACTIC NUCLEI: A COMPARISON WITH LUMINOUS QUASARS

We perform halo occupation distribution (HOD) modeling of the projected two-point correlation function (2PCF) of high-redshift (z~1.2) X-ray-bright active galactic nuclei (AGN) in the XMM-COSMOS field measured by Allevato et al. The HOD parameterization is based on low-luminosity AGN in cosmological simulations. At the median redshift of z~1.2, we derive a median mass of (1.02+0.21/-0.23)x10^{13} Msun/h for halos hosting central AGN and an upper limit of ~10% on the AGN satellite fraction. Our modeling results indicate (at the 2.5-sigma level) that X-ray AGN reside in more massive halos compared to more bolometrically luminous, optically-selected quasars at similar redshift. The modeling also yields constraints on the duty cycle of the X-ray AGN, and we find that at z~1.2 the average duration of the X-ray AGN phase is two orders of magnitude longer than that of the quasar phase. Our inferred mean occupation function of X-ray AGN is similar to recent empirical measurements with a group catalog and suggests that AGN halo occupancy increases with increasing halo mass. We project the XMM-COSMOS 2PCF measurements to forecast the required survey parameters needed in future AGN clustering studies to enable higher precision HOD constraints and determinations of key physical parameters like the satellite fraction and duty cycle. We find that N^{2}/A~5x10^{6} deg^{-2} (with N the number of AGN in a survey area of A deg^{2}) is sufficient to constrain the HOD parameters at the 10% level, which is easily achievable by upcoming and proposed X-ray surveys.

THE LARGE-SCALE THREE-POINT CORRELATION FUNCTION OF SLOAN DIGITAL SKY SURVEY LUMINOUS RED GALAXIES

We present new measurements of the redshift-space three-point correlation function (3PCF) of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS). Using the largest dataset to date, the Data Release 7 (DR7) LRGs, and an improved binning scheme compared to previous measurements, we measure the LRG 3PCF on large scales up to ~90 Mpc/h, from the mildly non-linear to quasi-linear regimes. Comparing the LRG correlations to the dark matter two- and three-point correlation functions, obtained from N-body simulations we infer linear and non-linear bias parameters. As expected, LRGs are highly biased tracers of large scale structure, with a linear bias b1 ~ 2; the LRGs also have a large positive non-linear bias parameter, in agreement with predictions of galaxy population models. The use of the 3PCF to estimate biasing helps to also make estimates of the cosmological parameter {\sigma_8}, as well as to infer best-fit parameters of the Halo Occupation Distribution parameters for LRGs. We also use a large suite of public mock catalogs to characterize the error covariance matrix for the 3PCF and compare the variance among simulation results with jackknife error estimates.