Computing the Three-Point Correlation Function of Galaxies in $\mathcal{O}(N^2)$ Time

We present an algorithm that computes the multipole coefficients of the galaxy three-point correlation function (3PCF) without explicitly considering triplets of galaxies. Rather, centering on each galaxy in the survey, it expands the radially-binned density field in spherical harmonics and combines these to form the multipoles without ever requiring the relative angle between a pair about the central. This approach scales with number and number density in the same way as the two-point correlation function, allowing runtimes that are comparable, and 500 times faster than a naive triplet count. It is exact in angle and easily handles edge correction. We demonstrate the algorithm on the LasDamas SDSS-DR7 mock catalogs, computing an edge corrected 3PCF out to $90\;{\rm Mpc}/h$ in under an hour on modest computing resources. We expect this algorithm will render it possible to obtain the large-scale 3PCF for upcoming surveys such as Euclid, LSST, and DESI.

[1]  John Couch Adams,et al.  III. On the expression of the product of any two legendre’s coefficients by means of a series of Legendre’s coefficients , 1878, Proceedings of the Royal Society of London.

[2]  Phillip James Edwin Peebles,et al.  Statistical analysis of catalogs of extragalactic objects. V. Three-point correlation function for the galaxy distribution in the Zwicky catalog. , 1975 .

[3]  S. Rey,et al.  Coupling of Modes of Cosmological Mass Density Fluctuations , 1986 .

[4]  William H. Press,et al.  The Cosmological constant , 1992 .

[5]  J. Huchra,et al.  IRAS galaxies versus potent mass: Density fields, biasing and omega , 1993 .

[6]  Bias and high-order galaxy correlation functions in the APM galaxy survey , 1994, astro-ph/9407079.

[7]  Second-Order Power Spectrum and Nonlinear Evolution at High Redshift , 1993, astro-ph/9311070.

[8]  A. Szalay,et al.  A New Class of Estimators for the N-Point Correlations , 1997, astro-ph/9704241.

[9]  Donald Hamilton,et al.  The evolving universe. Selected topics on large-scale structure and on the properties of galaxies , 1998 .

[10]  J. Frieman,et al.  The Projected Three-Point Correlation Function: Theory and Observations , 1999 .

[11]  A. Lewis,et al.  Efficient computation of CMB anisotropies in closed FRW models , 1999, astro-ph/9911177.

[12]  L. Wasserman,et al.  Fast Algorithms and Efficient Statistics: N-Point Correlation Functions , 2000, astro-ph/0012333.

[13]  J. Buchler,et al.  The onset of nonlinearity in cosmology , 2001 .

[14]  S. Colombi,et al.  Large scale structure of the universe and cosmological perturbation theory , 2001, astro-ph/0112551.

[15]  Ue-Li Pen,et al.  Fast n-point correlation functions and three-point lensing application , 2003, ArXiv.

[16]  I. Szapudi.  Three-Point Statistics from a New Perspective , 2004, astro-ph/0404476.

[17]  Y. Jing,et al.  The Three-Point Correlation Function of Galaxies Determined from the Two-Degree Field Galaxy Redshift Survey , 2003, astro-ph/0311585.

[18]  Three-Point Correlation Functions of SDSS Galaxies in Redshift Space: Morphology, Color, and Luminosity Dependence , 2004, astro-ph/0403638.

[19]  I. Szapudi,et al.  Measuring the Three-Point Correlation Function of the Cosmic Microwave Background , 2005, astro-ph/0508316.

[20]  I. Szapudi,et al.  The monopole moment of the three-point correlation function of the two-degree Field Galaxy Redshift Survey , 2005, astro-ph/0505422.

[21]  A. J. Connolly,et al.  The Effect of Large-Scale Structure on the SDSS Galaxy Three-Point Correlation Function , 2006, astro-ph/0602548.

[22]  D. Eisenstein,et al.  A robust estimator of the small-scale galaxy correlation function , 2006, astro-ph/0612103.

[23]  F. Castander,et al.  Clustering of luminous red galaxies – III. Baryon acoustic peak in the three-point correlation , 2008 .

[24]  S. White,et al.  Galaxy Formation and Evolution , 2010 .

[25]  Christopher Hirata,et al.  Relative velocity of dark matter and baryonic fluids and the formation of the first structures , 2010, 1005.2416.

[26]  Uros Seljak,et al.  Primordial non-Gaussianity from the large-scale structure , 2010, 1003.5020.

[27]  R. Mehrem,et al.  The plane wave expansion, infinite integrals and identities involving spherical Bessel functions , 2009, Appl. Math. Comput..

[28]  Berkeley,et al.  Supersonic Relative Velocity Effect on the Baryonic Acoustic Oscillation Measurements , 2011, 1105.3732.

[29]  Davis,et al.  THREE-POINT CORRELATION FUNCTIONS OF SDSS GALAXIES: CONSTRAINING GALAXY-MASS BIAS , 2010, 1012.3462.

[30]  A. Connolly,et al.  THREE-POINT CORRELATION FUNCTIONS OF SDSS GALAXIES: LUMINOSITY AND COLOR DEPENDENCE IN REDSHIFT AND PROJECTED SPACE , 2010, 1007.2414.

[31]  V. Fabrikant.  Elementary exact evaluation of infinite integrals of the product of several spherical Bessel functions, power and exponential , 2012 .

[32]  G. W. Pratt,et al.  Planck 2015 results - XVII. Constraints on primordial non-Gaussianity , 2014 .

[33]  R. Nichol,et al.  Modelling the redshift-space three-point correlation function in SDSS-III , 2014, 1409.7389.

[34]  D. Eisenstein,et al.  On the signature of the baryon–dark matter relative velocity in the two- and three-point galaxy correlation functions , 2014, 1411.4052.

[35]  N. Gehrels,et al.  The Whole is Greater than the Sum of the Parts: Optimizing the Joint Science Return from LSST, Euclid and WFIRST , 2015, 1501.07897.

[36]  D. Eisenstein,et al.  Accelerating the two-point and three-point galaxy correlation functions using Fourier transforms , 2015, 1506.04746.