Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the analysis of cosmic large - scale structure in the distribution of galaxies contributed significantly to establish the current cosmological concordance model . with the increasing degree of perfection of the standard two - point correlation analysis , the related measurements reached unprecedented precision    . however , the available and upcoming galaxy redshift survey data @xcite can not merely be characterized by its two - point correlation properties . already the filamentary structure that we see in those surveys indicates that the higher - order correlations play an important role . there have been different approaches to use this information for improving our knowledge about the constituents and evolution of the universe . the most direct approach of course is to simply measure the higher - order correlations directly . this has been done for both simulated and observed data starting for early catalogs with @xcite and @xcite . especially measurements of the three - point function have been standard for recent surveys @xcite and provide useful complimentary information @xcite . for higher orders , however , there are computational and conceptual challenges . on the computational side , increasing the order of the correlation function measured unleashes the curse of combinatorics . measuring all pairs in current surveys is feasible , for all triangles there are methods to do it @xcite , but beyond that it becomes increasingly hard . and then , from the conceptual side , even if measuring the full @xmath7-point function was possible , it would require modeling and covariances for a large number of data bins . already for a modest number of bins    for each parameter of the @xmath7-point function this is a challenge . therefore , it is worthwhile to consider other approaches that contain information on certain aspects of the higher - order clustering without referring to the full correlation functions .    since the beginning of large - scale structure analysis , a large number of such methods have been used , including moments of counts in cells , void probability @xcite , structure functions on minimal spanning trees , wavelet methods , the genus , etc . recently also approaches that attempt to bring the higher - order information back into the two - point correlation function through non - linear transforms were used like    a logarithmic transformation @xcite . the method ( re-)considered in this paper is using minkowski functionals    . they also also fall into the category of methods that provide condensed access to higher - order information . this information is encoded in the morphology of the density field . minkowski functionals access it by quantifying extended regions of the field through their geometrical properties like volume and surface ( see sec .  [ sub : minkowski - functionals ] ) . for the construction of these extended regions , two different methods have been widely used in the past :    in the first application of minkowski functionals for analyzing large - scale structure in the galaxy distribution in ( see @xcite for a short motivation and @xcite for a brief overview ) , the regions considered were defined by the union of balls around every galaxy location . we shall use the same method , referred to as germ - grain model , here and explain the details in sec . [ sub : germ  grain ] . after a couple of applications in the analysis of galaxy and cluster catalogs , the germ - grain model has been less used in recent years . the second and most popular form of the application of minkowski functionals works with extended regions that are defined by iso - contours of the field under investigation . it has been used for the higher - order contributions to the dark matter overdensity field @xcite , for the weak lensing shear field @xcite , for local morphology of cosmic structure and for the morphology of reionization bubbles @xcite . another very active field of recent application is the study of iso - temperature contour maps of the cosmic microwave background ( cmb ) , where minkowski functionals helped to establish stringent bounds on the non - gaussianity of the cmb . so the iso - contour minkowski functionals are a well established tool in various domains of cosmology . however , as we shall discuss in sec . [ sub : germ  grain ] , the germ - grain model accesses different aspects of the higher - order distribution and can be better suited for the analysis of galaxy redshift surveys . we have demonstrated some of these advantages already in the analysis of the sloan digital sky survey ( sdss ) data release seven ( dr7 ) luminous red galaxies ( lrgs ) in @xcite .    in the present paper , we concentrate more explicitly on the higher - order aspects of the minkowski functionals . to this end , we exploit the analytically known connection to the integrals over higher - order correlation functions . using direct measurements of the two - point correlation function , we can isolate the contribution of higher - order correlation functions to the minkowski functionals . this allows us to quantify the evidence for non - gaussianity of the galaxy distribution ( see @xcite and @xcite for related recent studies of non - gaussianity with isodensity contour functionals ) . the analysis is demonstrated with the final data from the baryon oscillation spectroscopic survey ( boss ) @xcite that is now publicly available within the dr12 of sdss - iii . being the largest spectroscopic galaxy redshift data set yet observed , it provides an ideal testbed for the analysis and reveals the strong presence of higher - order correlations already at redshifts of @xmath8 . we proceed as follows : sec . [ sec : minkowski - functionals - in - the ] recalls the properties of minkowski functionals in general and the germ - grain model in particular . we write out the connection of the functionals to higher - order correlations explicitly in [ sub : densdef ] and lay out in [ sub : density - dependence ] , how knowing the concrete form allows us to isolate the part of the functionals that is only dependent on higher - orders . we explain the subtleties of a gaussian reference model in the discrete case in [ sub : gauss - poisson ] and recall the basic properties of our minkowski code in [ sub : samplingcode ] .    in sec . [ sec : data ] we turn to the description of the sdss dr12 data used . [ sub : cmass ] summarizes the data and mocks we are using . in [ sub : ceometry ] , we address issues related to the imperfections in the data , constructing a custom angular mask that defines a high quality region . in [ sub : radial ] we turn to the radial distribution of the data and choose the analysis volume and the radial grid . in [ sub : weights ] we finally discuss the effect of the regions that can not be accounted for by our custom mask .    sec . [ sec : higherorders ] then turns to the core of the analysis . in [ sub : mfdensities ] we present the minkowski functionals of the mocks and the data and explain how we transform the data to be able to isolate the higher - order only part . this part is then examined in [ sub : highsubtract ] . in [ sub : seriescoefficients ] we demonstrate that the contributions to this higher - order part range at least up to the @xmath9-point function . to wrap up we discuss in [ sub : voidfraction ] an interesting aspects about the galaxy distribution using the base minkowski functional to study the volume in extreme voids in the observed field . we conclude in sec . [ sec : conclusion ] . in this section , we shall provide the basic notions necessary for understanding the quantities we used in our analysis . we summarize some key points already covered in the dr7 analysis paper @xcite , but with an emphasis on the crucial points that are the source for common misconceptions . minkowski functionals ( often called mfs later on ) quantify the morphology of extended bodies . they are defined as scalar functionals that map the shape of a body to a set of real numbers . a quite surprising result of integral geometry @xcite is that any such scalar functional , defined for bodies in euclidean @xmath10-dimensional space , that fulfills some basic requirements ( motion invariance , additivity and conditional continuity ) can be expressed as a linear combination of @xmath11 base functionals . there are different choices for the normalization of the base functionals . we will use one where the four mfs that exist for @xmath12 are denoted @xmath13 and related to geometrical properties of the body as follows : @xmath14 here , @xmath15 is the volume of the body , @xmath16 is its surface area , @xmath17 the integral mean curvature of the surface , and @xmath18 the euler characteristic ( the integral gaussian curvature of the surface ) .    instead of using the functionals @xmath19 directly , we will mainly work with the mf densities @xmath20 : @xmath21 here @xmath22 is the volume within the survey boundary . in the practical calculation of the @xmath20 discussed in sec . [ sub : samplingcode ] , @xmath22 is the volume of the region that is more than two ball radii away from our custom angular mask of fig . [ fig : custommask ] . as mentioned , mfs are only defined for extended bodies . for the analysis of large - scale structure , however , the galaxies are usually treated as points . for points , all quantities in eq . ( [ eq : minkrelgeo ] ) are , a priori , trivial . for a meaningful analysis of large - scale structure , one therefore has to smooth the point distribution in some manner to define extended regions . the two most popular choices for defining extended regions are :    1 . isodensity contours 2 .   the germ - grain model    the first approach attempts to reconstruct the underlying density field from the observed galaxies and to define the extended bodies as the interior of the surfaces of constant density . the value of this density threshold can then be varied and the functionals are studied as a function of this parameter .    in the second approach , the extended bodies are constructed by surrounding every galaxy ( the germ ) by an extended body ( the grain ) , which in our case is a ball of a certain radius , and calculating the functionals for the union of all these balls . the diagnostic parameter in this case is the common radius of the balls . [ fig : boolean ] shows how a structured body emerges with increasing radius . are surrounded by balls of a common radius @xmath23 . then the mfs of the body formed by the union of all the balls are studied as a function of @xmath23 . the parts of the body outside of the survey region are not taken into account . with increasing @xmath23 , more and more balls intersect and a complex structure develops . through the single and multiple intersections of neighboring balls , it encodes information about the @xmath7-point correlation hierarchy . plot from @xcite.,scaledwidth=45.0% ]    only those parts of the structure are taken into account that lie within the survey boundaries , in fig . [ fig : boolean ] illustrated as a cube . even though the basic analysis tool is the same in both cases , the properties of the two approaches are quite different . the isodensity contours , at least in theory , are more directly linked the underlying physical density field , whereas the germ - grain radii are purely diagnostic parameters that merely specify the scale at which the correlations are probed . regarding the temporal evolution , the mfs for isodensity contours turn out to be relatively independent of the ( linear ) growth of the density field , if the value of the threshold density is adjusted accordingly . on the germ - grain side , there is no corresponding relabeling and , therefore , as we also shall see below , they depend significantly on the amplitude of the power spectrum . besides these conceptual differences , there are some technical advantages to the calculation of germ - grain mfs :    1 . the evaluation directly uses the point distribution and does not use a smoothing procedure that is an art by itself . the use of balls does not introduce additional ( smoothing ) parameters . the construction of smooth surfaces can be noisy if the density is low . the rougher smoothing procedure in the germ - grain model on the other hand always constructs some well defined surface .    for both smoothing methods there have been detailed analysis how the values of the measured functionals are related to the parameters of the underlying density distribution . for the germ - grain case we recall this relation in the following section . due to the simplicity of the geometrical construction as a union of intersecting balls , it is relatively straightforward to relate the germ - grain mfs to correlation properties of the analyzed point distribution . for the case of a poisson distribution this has been achieved in @xcite , with the only relevant parameter being the density of the random sample . in the presence of structure , the dependence on higher order correlations has first been demonstrated in . in the following , we shall use the form derived in @xcite and detailed in the appendix of @xcite . in these works it has been shown that the ensemble average of the germ - grain mf densities defined in eq . ( [ eq : minkdens ] ) , for the case of balls as grains , can be expressed in an analytic manner as a function of the quantities @xmath24 and @xmath25 . @xmath24 is the ensemble average density of the point distribution studied , and the @xmath25 are the unnormalized modified mfs from the dr7 paper @xcite . the analytic relation reads    @xmath26    and provides the implicit definition of the @xmath25 . in the rest of the paper we will drop the explicit ensemble averaging brackets and use @xmath27 as a shorthand , as we never refer to the mf density of one individual realization . we discuss this aspect in more detail in sec . [ sub : samplingcode ] .    for a poisson distribution , the @xmath25 are simply the mfs @xmath28 of the balls of common radius @xmath23 that we use to obtain extended bodies : @xmath29 in the poisson case eq . ( [ eq : minkdensdef ] ) reduces to the result of @xcite . for a point distribution with structure , the @xmath25 are modified by weighted integrals of the dimensionless higher - order connected correlation functions @xmath30 . this modification is what we are interested in in this paper . it is convenient to define it as a multiplicative dimensionless correction to the poisson value , i.e. to study the quantities @xmath31 we will refer to @xmath32 as the modified mfs . they can be easily obtained from the quantities we measure , i.e. the average densities @xmath20 , by inverting eqs . ( [ eq : minkdensdef ] ) : @xmath33 in the analysis below , we will mainly use these modified mfs instead of the densities @xmath20 , as they are more directly related to the correlation properties of the point distribution via ( [ eq : minkcorrcon ] ) . from @xcite the relation of @xmath32 to the correlations of order @xmath34 , @xmath35 , is of the form @xmath36 with @xmath37 being simply the mfs of the _ intersection _ of @xmath34 balls . @xmath38 denotes a ball around the origin and @xmath39 a ball around the point @xmath40 . in our model all intersecting balls have the same radius @xmath23 . for general germ - grain models , this is not a necessary requirement . one could imagine to increase the weight of some points by assigning them larger radii than other points . using elliptic or even more irregular grains would also be possible . however , this would complicate the interpretation as we expect that it would lead to a more complicated relation of the functionals to the correlation hierarchy than the one of eq . ( [ eq : minkcorrcon ] ) . the advantage of the @xmath32 of being more directly related to the correlation properties is counterbalanced by an increase in the error . in @xcite we found weaker constraints when using the @xmath25 as compared to the @xmath20 . the main reason for this is the additional error that is introduced by removing the exponential suppression @xmath41 from @xmath42 in ( [ eq : minkdensdef ] ) . for @xmath43 , for example , @xmath44 so the relative error for the fraction of the volume not covered by any ball enters also for the higher functionals . this spoils to some extent the good precision that can be obtained for @xmath42 , because @xmath45 has to be measured in the full volume and @xmath42 only on the surface . especially for small @xmath46 it becomes costly to determine @xmath45 well enough . we will therefore restrict our analysis of the @xmath32 to parameter regions where @xmath46 is larger than @xmath47 . the modified mfs @xmath32 of eq . ( [ eq : minkcorrcon ] ) not only depend on our diagnostic parameter , i.e. the radius @xmath23 of the balls that we are using to smooth the galaxy distribution , but also on the sample density . the rather complicated @xmath23 dependence is entirely encoded in the weight functions @xmath48 in the integrals over correlation functions . the dependence on the mean density of the galaxy sample @xmath24 , however , is of an easier form . it is a simple power series in @xmath24 @xmath49 with coefficients @xmath50 and @xmath51 that contain all of the radius and correlation function dependence . @xmath52 is the volume of a ball @xmath53 of radius @xmath23 . we choose to define the second index of @xmath54 such that it directly corresponds to the order of the integrated correlation function .    the formal split into @xmath24 dependence and correlation function dependence in eq . ( [ eq : powerseriesdecomp ] ) allows us two things : first , we can remove the poisson part and the part that depends on the two - point correlation function , if we measure the sample density and the two - point correlations independently . this leaves us with a quantity @xmath55 that only includes higher - order contributions . this will prove the significance of their contribution to the mfs in sec . [ sub : highsubtract ] . second , we can extract the integrals over correlation functions by fitting a polynomial to the density evolution of @xmath56 . this allows to define a fitting formula for the form of the mfs in the sdss by giving the values and covariances of the coefficients @xmath57 . we shall use this in sec . [ sub : seriescoefficients ] . we close with the remark that , in general , a finite number of coefficients is not sufficient to describe the full shape of a power series of the sort ( [ eq : powerseriesdecomp ] ) . for cooked up point distributions with small @xmath30 , but a single large @xmath58 at some higher @xmath59 , giving only the lower coefficients will not determine @xmath32 well . for the assumed structure of the universe as the product of standard structure formation , however , we find that the coefficients @xmath60 are well behaved and the series @xmath32 converges reasonably well ( see sec .  [ sub : seriescoefficients ] ) .      in order to illustrate the relations ( [ eq : minkcorrcon ] ) , let us consider the simplest non - trivial case . in analogy with the continuous case , one can define a gaussian point distribution by the absence of all higher - order connected correlation functions . in this case the hierarchy of correlations in ( [ eq : minkcorrcon ] ) is truncated after the linear term and we have    @xmath61    the weight functions @xmath62 are non - zero only for the range of @xmath63-values for which the balls intersect . for a coordinate system centered around the origin , i.e. the center of @xmath53 , the @xmath40-values therefore lie in a sphere of radius @xmath64 . a configuration of balls of radius @xmath23 therefore probes the correlation function up to a scale @xmath64 . we shall therefore plot the radial dependence as a function of the diameter @xmath65 instead of the radius . for isotropic correlations ( [ eq : minkgaussexp ] ) becomes @xmath66 the shape of the weights @xmath67 is shown in fig . [ fig : integrationwindow ] . and their fourier transforms can be found in @xcite ( in the fourier form of @xmath68 in eq . ( 16 ) there is a typo : the prefactor of the integral should be @xmath69 instead of @xmath70 ) . ]     of eq . ( [ eq : minkgaussexp-1 ] ) , for a ball radius of @xmath71 . all functions have been normalized to their maximal value . the functionals of higher index @xmath72 probe the correlation function at larger distances . ( same plot already in @xcite ) [ fig : integrationwindow],scaledwidth=45.0% ]    let us just recall that the weights are such that @xmath73 and @xmath74 have a particularly interesting form . @xmath73 directly contains the amplitude of fluctuations on the scale of the balls @xmath75 , @xmath76 and @xmath74 is a simple integral over the two - point function @xmath77 which could therefore be recovered from the derivative of @xmath74 . let us pause here and elaborate a bit on the gauss - poisson process as such ( see @xcite for a more complete discussion of what follows ) . it is tempting to think of the truncated quantities ( [ eq : minkgaussexp-1 ] ) as being the mfs of a point distribution that represents the gaussian limit of an actual point distribution that has higher - order correlations . however , in general there will not be a gauss - poisson point process that actually realizes these mfs , as in order for a gauss - poisson process to exist , @xmath78 and @xmath24 have to fulfill certain conditions .    already in the continuum case , the clustering amplitude is limited by the fact that the matter density can only be positive . the density distribution , therefore , is in any case only approximated by a gaussian . this approximation of course is , excellent for early times and large scales , but when matter fluctuations @xmath75 grow too large it breaks down . this is also the point where higher - order correlations become important .    in the same manner , the clustering amplitude of a gauss - poisson point process is limited . in this case , however , also the mean density of the point sample @xmath24 plays a role . for low enough density a thinning of the point distribution goes to an approximated poisson distribution . this effect can compensate larger fluctuations and it can be shown @xcite that the corresponding constraint is @xmath79 where @xmath80 are the fluctuations in the number density of cells @xmath81 stemming from structure ( @xmath82 are the full fluctuations including the poisson term ) , @xmath83 is the volume of one of these cells , and @xmath84 the mean number of particles in such a cell . this constraint has two interesting implications : first , we can already tell from the amplitude of the two - point function if the field will have a significant amount of higher order clustering . when we evaluate ( [ eq : fluctconstraint ] ) for the structure that we measure in the boss galaxies , we find that the constraint actually is violated . second , to study the gaussian limit on large scales it is not sufficient to simply increase the domain size ( the ball size in our case ) . in some cases , where the fluctuations on a region decay faster than its volume the constraint indicates that the structure for larger cells becomes more gaussian . in general , however , one will have to reduce the number density . in the continuous case , by smoothing over larger and larger regions , the amplitude of the fluctuations between these regions decreases and the gaussian approximation becomes better and better . in the discrete case , the analogous coarse graining procedure involves a subsampling of the density field . in the same way that the number of regions is reduced in the continuous case when averaging the density over larger cells , individual points of the downsampled point distribution represent larger regions . we will demonstrate this effect in sec . [ sub : mfdensities ] where we show the mfs for a couple of densities . the lower this density becomes , the better the gaussian model fits . in addition to the requirement ( [ eq : fluctconstraint ] ) , @xcite also showed that the correlation function of a gauss - poisson point process has to satisfy @xmath85 and can not be negative    @xmath86    any violation of ( [ eq : fluctconstraint ] )  ( [ eq : xiconst ] ) already implies the existence of higher - order correlations @xmath30 with @xmath87 . we note finally that the condition ( [ eq : fluctconstraint ] ) prevents ( [ eq : v0-sigma ] ) and therefore the volume fraction @xmath45 in ( [ eq : minkdensdef ] ) from becoming negative . so we could have derived a similar constraint on a gauss - poisson process from the requirement that @xmath88 . the procedure we use to measure the functionals is the same as in the dr7 analysis @xcite . as the quantities that are related to the hierarchy of correlation functions via ( [ eq : minkdensdef ] ) are , strictly speaking , the ensemble averages of the mf densities @xmath20 , we approximate these averages by creating randomly subsampled realizations of the point distribution and averaging over their mf densities . mathematically speaking there is no guarantee that these averages should converge to those of an ensemble of independent realizations ( see also @xcite ) , but practically the procedure is working quite well . for studies on small scales one could in principle use subfields of the full survey as independent realizations and perform an average that is then more likely to converge ( even though it also requires the assumption of ergodicity and is only valid approximately as neighboring regions are not completely independent ) . for the scales that we are interested in , however , there is only one realization of the survey and therefore we use the subsampling . for details on how we generate the samples see the end of sec . [ sub : radial ] . the calculation of the @xmath20 for a given downsampled realization is performed by chipmink ( code for high  speed investigation of partial minkowski functionals ) , introduced in the dr7 paper @xcite . in a first step , it generates neighbor lists for every galaxy in the sample . galaxies are considered as neighbors if they are closer than twice the maximal radius of the balls . to determine @xmath89-@xmath90 , the code considers galaxies that are more than @xmath64 away from the boundary . using the list of neighbors the code calculates the uncovered surface , integrated mean curvature and euler characteristic of the ball put around the galaxy . these partial functionals are then summed to give the global @xmath20 . the volume used for obtaining the densities from the summed @xmath91 is the effective part that is more than @xmath64 away from the boundary .    for the volume functional @xmath45 another set of neighbor lists is generated that contains all galaxies within @xmath92 around points from a random sample filling the same volume as the data , but also eventual holes in the data that are not treated as internal boundaries . around each primary random point we generate @xmath4 random points in spheres of increasing radius @xmath93 . the volume fraction covered by the balls is then estimated by counting the fraction of the secondary randoms that have a distance of less than @xmath23 to any of the galaxies in the neighbor list of the primary random point . as compared to the dr7 code , we made a couple of enhancements . besides improved memory management , the code can now handle boundaries of arbitrary shape in ra - dec . it can use data with non - constant density , uses a larger fraction of points for radii smaller than @xmath94 , and samples the @xmath95-plane more efficiently than the earlier version . we plan to make also the current version publicly available . having completed the discussion of the methods , we now turn to a description of the data . the galaxy sample we are analyzing was observed by the sloan digital sky survey ( sdss ) @xcite using the 2.5-meter sloan telescope @xcite . photometric data used for target selection was obtained using a drift - scanning mosaic ccd camera @xcite and information on the camera , photometry and photometric pipeline can be found in @xcite . all the photometry was re - processed and released in the eighth data release @xcite . based on the so obtained target catalog , the spectroscopic redshifts are taken by the spectrographs described in @xcite . the final step of spectroscopic data reduction and redshift determination is presented in @xcite .    the data analyzed here is contained in the twelfth data release ( dr12 ) @xcite of sdss data and was obtained during the third stage of sdss @xcite within the baryon oscillation spectroscopic survey ( boss ) @xcite . it is described in detail in @xcite . the full data set consists of the spectroscopic redshifts of 836,347 galaxies . the area probed is divided in a northern and a southern part and two different target redshift ranges . for our analysis we use the largest contiguous part of the survey , which is the northern cmass sample . given the technical progress in the boundary treatment of our code , mentioned in sec . [ sub : samplingcode ] , it would now also be possible to use it for the combined sample analyzed in @xcite . the mock catalogs we are using for estimating errors and covariances were constructed using two different methods . one is based on the quick particle mesh technique ( qpm ) described in @xcite and the other one is a hybrid method @xcite which employs what the authors call `` augmented '' second - order lagrangian perturbation theory . it was used for generating the multidark(md)-patchy mocks @xcite . we shall work with both mock catalogs in our analysis , but concentrate on the md - patchy mocks . , @xmath96 and @xmath97 away from the nearest boundary.[fig : custommask],scaledwidth=45.0% ]    for taking the data of the northern cmass sample , @xmath98 objects were targeted over an area of @xmath99 . out of these , @xmath100 turned out to be galaxies and obtained redshifts from good boss spectra @xcite . they are distributed over @xmath101 . their redshifts range approximately from @xmath102 to @xmath103 . the central region of redshift @xmath104-@xmath105 that is used in most of the cmass analysis ( see e.g. @xcite ) contains @xmath106 galaxies . the md - patchy and qpm mock catalogs are constructed such that they trace the number density of the data . they are available for the central range of @xmath104-@xmath105 . their angular distribution fills the full survey geometry . for matching the angular distribution of the data as well , regions of bad data quality or missing observations can be masked out . the qpm mocks contain on average @xmath107 galaxies . after application of the full veto mask the average number is @xmath108 . for the md - patchy mocks , the unmasked number of galaxies is @xmath109 and the masked one @xmath110 . our custom radial selection of a redshift range of @xmath111 , explained in section [ sub : radial ] , reduces the number of galaxies that we use to @xmath0 .      in order to compare the data to the galaxy samples created in the mock catalogs , it would be sufficient to calculate the mfs for the masked mocks and the data . however , ideally we would like to derive quantities that are independent of effects of the survey boundary . this is not necessary for the direct comparison , but desirable for later comparison of the data with other simulations without having to apply the survey mask to them . when estimating the correlation function , the effect of boundaries ( inner boundaries such as holes as well as the outer boundary ) can be taken into account by using random samples with the exact same boundaries as the data . the excess of data over random pairs in a given bin is then solely due to structure .    in the case of the mfs , unfortunately there is no simple analogue to this approach . we tried to correct the mfs of the masked mocks with the mfs of a random distribution and the same mask . a comparison with the unmasked mocks showed that this correction is not effective in the mf case . in order to minimize the influence of boundaries , we therefore choose to restrict our analysis to regions with good data quality . this is achieved by using only galaxies that are more than @xmath64 away from the closest boundary , where @xmath23 is the common radius of the balls around every galaxy ( see sec .  [ sub : germ  grain ] ) . of course , especially for large @xmath23 , this reduces the volume probed and so the number of galaxies that can be used . we therefore have to find a compromise as to which boundaries we include and which we can ignore without affecting the result too much .    as discussed in detail in @xcite , there are at least seven different effects that lead to bad or missing data . for each effect there is a veto mask that can be used in the correlation function analysis . the effects are :    * the central bolt obscures the center of the plate * ly@xmath112 targets are given priority in fiber assignment * no observations around bright stars * no observations around bright objects * non - photometric conditions * spectra with bad seeing * spectra degraded by extinction    the regions obscured by the first three effects are quite small individually , but relatively evenly distributed over the sky . including them , and especially the @xmath64 region around them , would reduce the usable area by a large amount , so we can not correct for their effect . the fourth effect is also present all over the sky , however , the obscured regions have a wide range of sizes . therefore we can exclude the largest ones .    . the points indicate the minimal value of the density of the 1000 qpm mocks in the respective bin . the reference density with respect to which we define our downsampling is @xmath113 . this combination of reference density and redshift range maximizes the number of downsampled galaxies in the boxed blue shaded region . we also use the points above the reference density by downsampling always from the full density.[fig : radialselection],scaledwidth=45.0% ]    the combination of the lower four masks is shown as the black regions in fig . [ fig : custommask ] . of those regions we select those with an area on the sky of more than about @xmath114 . we densely sample these regions and their boundaries with points that are shown in blue . these are used to determine the distance of a galaxy from the nearest boundary . the projected holes created by ball diameters @xmath64 of @xmath115 , @xmath96 and @xmath116 are shown as the regions not filled by orange , blue and yellow color . at these diameters the effective volume left is @xmath117 , @xmath118 and @xmath119 , respectively . for the largest scales probed of @xmath120 , only @xmath121 remains . as the mask we constructed does not take into account all imperfections , we have to test the extent to which the result still depends on the survey boundaries . in order to study this , we removed the galaxies from the mocks that fall into one of the seven masks described above . then we determined the distance of the remaining galaxies from our boundaries that are shown in fig . [ fig : custommask ] . we calculated their germ - grain mfs and those of the galaxy in the corresponding cubic mocks periodically wrapped to remove any boundary . the resulting comparison is shown in fig . [ fig : maskcorrection ] . it also contains the case where we use the full veto mask to remove the unobserved galaxies , but only take the outer boundaries into account in the calculation of the mfs . we find that our custom mask greatly improves the match with the boundaryless case as compared to not taking inner boundaries into account . after the correction the periodic mfs are now within the error bars of the sky mfs . the sky values of @xmath73 in fig . [ fig : maskcorrection ] have been corrected for a systematic @xmath122 offset from the periodic case . we shall apply the same correction throughout the rest of the paper . as the offset has been derived from a comparison of the sky mocks to the periodic mocks , one could think of using the deviations of fig . [ fig : maskcorrection ] to apply density and scale dependent corrections to both the sky mocks and the data , to empirically remove the effect of the boundary . we decided not to apply this correction as it would make the results for the mfs of the data completely dependent on this particular set of mocks .    from fig . [ fig : maskcorrection ] , we also see that our correction is largest for the fourth functional @xmath74 and least important for @xmath73 . this reflects the property of the integration kernel of fig . [ fig : integrationwindow ] . with increasing index @xmath72 , the functionals probe larger and larger scales and are hence more sensitive to the large scale boundary effects we removed . in addition to the angular selection , we also apply some radial cuts . as described in sec .  [ sub : density - dependence ] we want to study the mfs of the galaxy distribution as a function of ball radius , but also sample density . a controlled variation of sample density requires a roughly constant density throughout the sample . otherwise the summation of partial mfs ( see sec .  [ sub : samplingcode ] ) from different regions of the survey would combine @xmath20 values effectively measured at different densities , removing the clear distinction between the @xmath24 and the @xmath30 dependence in eq . ( [ eq : minkdensdef ] ) .    to generate a homogeneous sample we downsample the density to a common value for all radial shells . the concept is shown in fig . [ fig : radialselection ] .    -@xmath24-space . for the smaller diameters there are @xmath123 density points , for the largest diameter used in the analysis @xmath124 , there are @xmath125 . the lines of constant @xmath45 are calculated for a poisson point process and correspond to \{0.99 , 0.95 , 0.75 , 0.55 , 0.35 , 0.15 , 0.05 , 0.01 } from top to bottom.[fig : densitygrid],scaledwidth=45.0% ]    there , we plot the binned density distribution of the qpm mocks . for each bin we show the minimum of the 1000 qpm mocks . using this density distribution , we choose a reference density such that the number of points in the constant density region of the considered redshift range is maximal . this density is @xmath126 and corresponds to @xmath127 galaxies . the total number of galaxies in this redshift range is @xmath0 . as we used the minimum of all mocks , we can be sure that we can downsample from this reference density and obtain the same constant density for every single mock realization . the sampling points in radius and density that we used are shown in fig . [ fig : densitygrid ] . for ball diameters up to @xmath128 , there are @xmath123 density points ranging from @xmath129 to @xmath130 of the reference density . for diameters above @xmath128 we decrease the maximal density and therefore the number of density points . the reason for this is two - fold : first , as discussed in sec . [ sub : densdef ] , for very small fractions of uncovered volume ( i.e. large radius and/or high density ) , the precision of the determination of @xmath45 is no longer high enough to give reasonable error bars for the derived quantities @xmath32 . second , the runtime is much worse if the galaxies have a large numbers of neighbors . the density cutoff in fig . [ fig : densitygrid ] is therefore roughly chosen along the lines of constant @xmath45 . [ fig : densitygrid ] includes these lines for the poisson case , which we used for a first estimate of reasonable density ranges .    from the range of diameters shown in fig . [ fig : densitygrid ] , we choose @xmath131 to @xmath132 for our analysis . for smaller scales than @xmath131 we expect to pick up too much effects from the ignored holes in the survey and for larger scales than @xmath132 , the volume far enough from the boundary becomes too small .    for our independent measurement of the correlation functions of the data and the mocks , we chose to use those regions that are more than @xmath133 away from the closest boundary that we treat . the structure that we measured for these regions is shown in fig . [ fig : meancorrfun ] . the concrete monte carlo procedure we use is then to randomly select galaxies from the full number of galaxies in our redshift range ( and _ not _ a pre - downsampled reference density as in the dr7 paper ) for each @xmath134-point of the grid in fig .  [ fig : densitygrid ] , to obtain a sample with the constant target density in each radial bin . we then calculate the mfs for this configuration . we repeat the random selection and average the resulting mfs over these different realizations of the subsampled point distribution . the number of realizations is largest for the smallest densities and ranges from @xmath135 for a fraction of @xmath129 to @xmath136 for a fraction of @xmath130 . scatter of the mocks.[fig : meancorrfun],scaledwidth=45.0% ]      in addition to the regions that could not be observed or had bad data quality and that therefore had to be masked out , there are also imperfections in the rest of the observed survey volume . even though the completeness of the survey is high , see @xcite , there are still galaxies that have not obtained robust redshifts ( no - z ) . in addition due to the physical size of the fibers of the spectrograph , galaxies separated by less than @xmath137 can not be observed simultaneously ( close pairs ) . in the correlation function analysis , these holes are accounted for by upweighting nearby galaxies . in our case this can not be done as easily . in principle one could use imaging data for the angular position and the redshift of nearby galaxies to mimic the procedure in the correlation function case and fill in the holes . however , as we are already ignoring a large number of holes and looking at relatively large scales , we decided not to take these holes into account and to mask them in the mocks consistently . however , we include another class of weights correcting for spurious density fluctuations due to stellar density and seeing . for the correlation function analysis , these have been combined into a systematic weight ( see @xcite ) . we use this systematic weight to adjust the probability of selecting a galaxy when we downsample to a given target density . as the modifications are below @xmath138 and not factors of @xmath139 to @xmath125 as in the case of the no - z and close pairs weights , the fluctuations can be accounted for relatively easily by this procedure . in this section we first present the measured values for the observables @xmath20 introduced in sec . [ sub : densdef ] . we study their radial and density dependence based on the md - patchy mocks . then we turn to the analysis of the higher - order contributions . for this analysis , we need the @xmath32 that are obtained from the measured @xmath20 through eqs . ( [ eq : inverse - mink ] ) . from the @xmath32 we can then simply subtract the two - point part of the expansion ( [ eq : minkcorrcon ] ) to extract the @xmath140 part ( [ eq : defmodminkhigher ] ) that only depends on higher orders . we will see that it contributes significantly to the full @xmath32 . finally we analyze the coefficients @xmath60 of ( [ eq : minkcorrcon ] ) order by order and compare them to correlation function integrals .      as mentioned in sec . [ sub : samplingcode ] , the quantities calculated by chipmink are the the mf densities @xmath20 . their dependence on the ball diameter is shown in fig .  [ fig : minfundens ] and based on the mean of @xmath141 md - patchy realizations . we only used @xmath141 mocks in this case because fig . [ fig : minfundens ] is mainly for illustration and the large range of diameters is not optimal for the evaluation of the @xmath32 as discussed at the end of sec . [ sub : densdef ] .        in order to illustrate also the density dependence , we plot @xmath142 different densities ranging from @xmath129 to @xmath130 of our reference density of @xmath126 ( the same density points as in fig . [ fig : densitygrid ] , however with an increased maximal radius ) . for the most extreme cases of @xmath129 and @xmath130 , we also plot the mfs of a poisson point process as dashed lines . for the high density case , the functionals for the mocks are clearly non - poissonian , whereas in the low density case the downsampling erased already a lot of the structure and the difference to a poisson sample is less important . let us qualitatively describe the difference for the euler characteristic @xmath90 ( @xmath143 ) . the faster drop of the functional with structure indicates that the individual balls aggregate earlier , which reduces the number of independent components @xmath144 faster than for the poisson case . then the holes @xmath17 in the union of multiple balls prevail up to larger radii and in the end there are less cavities @xmath145 but at larger radii .    for our quantitative analysis in sec . [ sec : higherorders ] we will mainly use the modified mfs @xmath32 , eq . ( [ eq : minkcorrcon ] ) . therefore , we use eqs . ( [ eq : inverse - mink ] ) to transform from the initial observables @xmath146 to @xmath147 . fig . [ fig : modminkfun ] shows the dependence of the @xmath32 on the ball diameter for the same densities as in fig . [ fig : minfundens ] , but now for 900 md - patchy mocks .        in these plots the convergence of the functionals to those of the poisson point process with decreasing sample density is more visible than in fig . [ fig : minfundens ] . however , given the small error bars , even for the lowest density of @xmath148 the structure is well detected . the differing number of radial points for different densities illustrates our grid of points sampling the @xmath134 plane of fig . [ fig : densitygrid ] . for the highest density , @xmath125 radial bins of @xmath149 are possible until the survey volume is significantly filled , for the lowest density there are @xmath150 bins ( of which only the first @xmath151 are shown here ) . ccccccccc + & & & & & +   + & @xmath152 & @xmath153 & & @xmath152 & @xmath154 & & @xmath152 & @xmath154 + @xmath155 & 1810 & 41.1 & & 171 & 10.1 & & 29.8 & 1.31 + @xmath156 & 1120 & 31.7 & & 88.1 & 5.93 & & 32.1 & 1.54 + @xmath157 & 454 & 19.1 & & 46.9 & 2.93 & & 29.6 & 1.28 + @xmath158 & 182 & 10.6 & & 29.7 & 1.30 & & 23.5 & 0.694 +      we will now turn to the comparison of the density dependence of the modified mfs @xmath32 in the mocks and in the sdss data . instead of probing the grid of fig . [ fig : densitygrid ] for a fixed density as in fig . [ fig : modminkfun ] , we show in fig .  [ fig : modminkdensdep ] the @xmath32 as a function of the sampling density for a ball diameter of @xmath133 . the form of @xmath32 in this direction of parameter space is easier to describe than in the @xmath159 direction , due to the power series form of the @xmath24-dependence , see eq . ( [ eq : powerseriesdecomp ] ) . in the plots of fig . [ fig : modminkdensdep ] , the poisson - case corresponds to the horizontal line at @xmath160 , thanks to the normalization of the @xmath32 in eq . ( [ eq : minkcorrcon ] ) . the data are clearly different from @xmath160 , showing that there is structure in the galaxy distribution . the straight lines that depart from the poisson case at zero density show the @xmath32 truncated after the term linear in @xmath24 , see eq . ( [ eq : minkgaussexp ] ) . they have been calculated by measuring the two - point correlation functions for each of the @xmath161 qpm and @xmath162 md - patchy mocks and integrating them with the corresponding weight functions @xmath163 of eq . ( [ eq : minkgaussexp-1 ] ) . the plot shows the average and the @xmath164 fluctuations between the different mock realizations . the average of the correlation functions that formed the basis of these lines is shown in fig . [ fig : meancorrfun ] . the clear deviation of the full @xmath32 of the data and the mocks from these lines shows that the integrals over higher - order correlation functions in the series expansion of @xmath32 actually do play a significant role for the structure in the sdss galaxies . this becomes even more apparent in the plots of fig . [ fig : modminksubtract ] . cccccc + & & & +   + & @xmath152 & @xmath153 & & @xmath152 & @xmath154 + @xmath155 & 10300 & 97.3 & & 290 & 5.50 + @xmath156 & 11600 & 104 & & 282 & 5.22 + @xmath157 & 12300 & 107 & & 221 & 2.76 + @xmath158 & 12500 & 108 & & 252 & 4.04 +    for these we use the integrated two - point correlations and subtract them off the corresponding @xmath32 for each mock individually . from these subtracted quantities @xmath140 ( eq .  ( [ eq : defmodminkhigher ] ) ) we calculate the mean of the mocks and the covariances of the density points . the plots of the residuals show , that the higher - order contributions to the modified mfs are quite well described by the structure in the md - patchy mocks . the data happens to lie close to the @xmath164 line of the fluctuations between the @xmath162 mocks for all four mfs . it is quite encouraging that the mocks that were constructed to match the two- and three - point correlations only also capture what is going on in the higher orders quite well .    as the data points in both fig . [ fig : modminkdensdep ] and fig . [ fig : modminksubtract ] are strongly correlated , we need to take into account the full covariance matrix in order to estimate the significance of the remaining difference . the results of this calculation are shown in tab . [ tab : chisqmockdat ] . for the curves of fig . [ fig : modminkdensdep ] and fig . [ fig : modminksubtract ] with @xmath123 density points , we first evaluate the @xmath3 values and then convert their @xmath165-values into those of a two sided deviation for a gaussian , which we express as multiples of one standard deviation . we inflate the error bars by multiplying the inverse covariance matrix by a factor of @xmath166 , in order to account for the finite number of mocks we use for estimating the covariance matrix @xcite . we find that the md - patchy mocks are a much better match of the data than the qpm mocks which are about @xmath167 away from the data in this case . for the higher components only , i.e. the central column of tab . [ tab : chisqmockdat ] , the agreement is even better and does not exceed @xmath168 for any of the functionals . so also quantitatively with covariances included the impression from the residual plots of fig . [ fig : modminksubtract ] that the typical deviation is around @xmath169 is correct . when we combine all the information of the relevant @xmath170 points in our sampling grid of fig . [ fig : densitygrid ] , i.e. ball diameters from @xmath131 to @xmath132 and their complete density dependence , we find in tab . [ tab : chisqradcomb ] that the deviation of the higher - order md - patchy modified mfs is at most @xmath171 . for some diameters as in tab . [ tab : chisqmockdat ] the agreement is nearly perfect . this is very encouraging , as the mocks have only been generated to match the two- and three - pt . given the combined values of tab . [ tab : chisqradcomb ] , it is clear that there is still some room for improvement in the construction of realistic mocks to match the higher orders even better .    regarding the deviation of the @xmath140 from the case of vanishing higher - order correlations ( the zero line in fig .  [ fig : modminksubtract ] ) , the conclusion is clear . there is strong evidence for higher - order correlations in the mocks and the data as seen by the mfs . even though this is not surprising , as at least the three - point function of the sdss data is well known ( see e.g. @xcite and @xcite ) , it is still interesting to see the statistical power that we have at hand for testing models with different higher - order correlations .    for all of the quantities of tabs . [ tab : chisqmockdat ] and [ tab : chisqradcomb ] , one should keep in mind that the different mfs are not independent . so in order to derive one single value that would quantify the mismatch between mocks and data , one would also here have to take the covariances into account . for a single radius this is possible . for the full data set , we are limited by the number of mocks that we analyzed as an estimate of the @xmath172 covariance matrix would be noisy .    finally , let us illustrate the scale dependence of the higher order contributions . [ fig : modminksubtractrad ] shows the deviation of the full @xmath32 ( [ eq : minkcorrcon ] ) from the two - point part ( [ eq : minkgaussexp ] ) , i.e. @xmath140 of ( [ eq : defmodminkhigher ] ) . so it is the same quantities as in fig . [ fig : modminksubtract ] , but now as a function of ball diameter for four fixed densities @xmath173 . the higher the density , the larger is the deviation from the two - point part , because the higher - order contributions play a more important role in the @xmath24-power series . in the radial direction , the functionals probe the higher - order correlations at different scales . for the two - point part we have seen in fig . [ fig : integrationwindow ] that the integration kernel picks up correlations at larger scales for increasing functional index @xmath72 . this is most probably also the case for the higher - order correlations ( we explicitly checked it for the three - point integrals ) and so @xmath174 probes the largest scales that we accessed in this analysis . the decreasing slope of the radial dependence reflects a similar behavior in the three point integral of @xmath175 . for even larger scales , the latter integral becomes negative and the difference to the two - point case decreases again . having found a significant contribution of higher - order correlations to the mfs , we now want to study which orders actually contribute . using the relation in eq . ( [ eq : minkcorrcon ] ) we can decompose the functionals in different contributions from higher orders by fitting a polynomial of order @xmath7 to the density dependence of the @xmath32 as described in sec . [ sub : density - dependence ] . this decomposition is not necessary for using the higher - order mf information in the comparison of the model to the data , as we have shown in the previous section , but it illustrates the convergence properties of the @xmath32 series .    in fig . [ fig : higherfit ] we evaluate the convergence of the fit for increasing polynomial order @xmath7 .     of eq . ( [ eq : minkintegr ] ) . to the @xmath24-dependence of @xmath43 at a diameter of @xmath176 we fit polynomials in @xmath24 of increasing order taking the full correlations of the points into account and extract the @xmath177 from the best fit coefficients . the linear to third order polynomials in the large panel give a terrible fit already visually . the fourth and fifth order in the `` residual '' part of the plot indicate convergence at 5.-order . the corresponding @xmath3 values can be found in tab . [ tab : higherfitchi ] and confirm this fact.[fig : higherfit],scaledwidth=45.0% ]    the deviations of the fitted model from the data are quantified in tab . [ tab : higherfitchi ] .    from both presentations it is clear that the linear term that only includes the two - point coefficient @xmath178 is a very bad fit to the density evolution . the table indicates that only after the inclusion of the fifth term , which corresponds to integrals over the six - point function , the polynomial fit captures the shape of the modified mfs well . however , this only means that we have some sensitivity of the mfs of the mocks to contributions from up to the six - point function . it does not mean that these contributions are also significantly different from the ones in the data . in fact , as we have found in sec .  [ sub : highsubtract ] the higher - order contributions of the mocks and the data actually agree relatively well .    in order to evaluate how strongly the coefficients @xmath60 of eq . ( [ eq : powerseriesdecomp ] ) are biased by truncating the power series at a finite @xmath7 , we plot in fig . [ fig : evwithfitorder ] the dependence of the mean of the best fit values of @xmath179 , @xmath180 and @xmath181 of the @xmath162 md - patchy mocks on the polynomial order @xmath7 for a fit at a diameter of @xmath133 . as high values of @xmath7 also increase the chance of picking up spurious fluctuations , we decided to put a prior on the higher order coefficients . to this end we fit a three - parameter pad model to the density evolution of the mean mfs extracted from @xmath6 periodic md - patchy cubes . taylor expanding this model provides some model for the higher - order coefficients . we choose a gaussian prior around the so obtained values with a standard deviation corresponding to the value of the mean @xmath182 . we only apply this prior to coefficients beyond the four - point function coefficient @xmath183 . the first two coefficients @xmath178 and @xmath184 are already well measured without any prior , and for @xmath183 it is sufficient to regularize the higher orders .        fig . [ fig : evwithfitorder ] shows that the value of the first coefficient @xmath178 is independent of @xmath7 as soon as @xmath7 is high enough that we have a good fit , i.e. @xmath185 from fig . [ fig : higherfit ] . the second coefficient @xmath184 is also largely independent of fit order , but needs a higher fit order @xmath7 than @xmath178 to be unbiased . for @xmath183 there is still some remaining uncertainty from the fit order , which is , however , still much smaller than the error from the scatter of the mocks . [ fig : errorellipse ] shows evolution of the mean value together with the @xmath164 contours of the parameters @xmath180 and @xmath181 for different @xmath7 starting at @xmath186 , to include also the covariances between the parameters .     and @xmath181 of the fit to @xmath187 of fig . [ fig : higherfit ] for a diameter of @xmath176 . the leftmost point is for a third order polynomial and the order is increasing from left to right with a maximum order of @xmath188 . the error bars indicate @xmath164 intervals . starting from the ninth order , error bars include all higher order points.[fig : errorellipse],scaledwidth=45.0% ]    below the order @xmath189 that provides the first good fit in tab . [ tab : higherfitchi ] , the best fit values are strongly dependent on fit order . for orders starting at @xmath190 the shift becomes much smaller than the scatter of the mocks . for the following plots , we will use polynomials of order @xmath191 for the @xmath43 fit at @xmath192 . to determine the fit order for the radii with a smaller number of density points , i.e. beyond @xmath128 , we constructed and analyzed test data . these were created from a pad model fit to the cut sky mocks and provide a model with known coefficients . then we applied the fitting procedure to these test data and found that for high enough @xmath7 , the coefficients we put in were recovered within the errors estimated by the fluctuations from the mocks . from the results of these test data fits , we determined for every diameter the minimum @xmath7 needed . for the larger diameters this value is lower due to the decreasing contribution of higher orders . the reason for this decrease is the steep scale dependence for the higher - order coefficients , an example of which we will see in fig . [ fig : threeptcoeffs ] for the three - point term ( see also fig .  14 of @xcite for the scale dependence of the coefficients of @xmath73 for a log - normal model ) . the decomposition into @xmath7-point components also allows us to test to which extent the deviation of the data from the mocks in fig . [ fig : modminkdensdep ] could be simply the result of a misnormalization of the power spectrum . the green line in fig . [ fig : modminkdensdep ] therefore shows the mean of the md - patchy mocks for a reduction in the power of @xmath193 . it has been obtained by taking the coefficients of the @xmath9-th order fit of fig . [ fig : higherfit ] and multiplying the integrated correlation functions with the power of the reduction factor that corresponds to the number of power spectra involved , i.e. @xmath194 for @xmath35 . of course this is only a very simplified model , given that the fit is not unbiased , especially for higher orders , but it should still provide a useful approximation . from the result in fig . [ fig : higherfit ] we see that at least approximately it removes the offset of the data from the mocks and makes them agree even better . this seems to be consistent with the analysis of the power spectrum in @xcite . our redshift range roughly corresponds to their second redshift bin , for which their fig . 1 indicates that the data also is slightly lower than the mean of the md - patchy mocks . cccc + fit order & coeff . & @xmath152 & p - value + @xmath160 & @xmath179 & 7859 & @xmath195 + @xmath139 & @xmath180 & 1509 & @xmath196 + @xmath197 & @xmath181 & 306 & @xmath198 + @xmath199 & @xmath200 & 47 & @xmath201 + @xmath189 & @xmath202 & 11.6 & @xmath203 +    for the two cases @xmath178 and @xmath184 where the convergence in fig . [ fig : evwithfitorder ] is convincing , we compare in figs . [ fig : twoptcoeffs ] and [ fig : threeptcoeffs ] the extracted coefficients to reference models . [ fig : twoptcoeffs ] shows the two - point coefficients @xmath178 . for those , @xmath204 coincides with @xmath205 ( see eq . ( [ eq : v0-sigma ] ) ) and @xmath206 corresponds to the volume averaged correlation function ( see eq . ( [ eq : v3-xi ] ) ) . in addition to the average of the coefficients extracted from the md - patchy mocks , we show the coefficients of the data and the expected coefficients derived from integrals over the average two - point correlation function that we measured independently . the overall agreement of the two independent determinations is reassuring . it provides a confirmation that the fit to @xmath56 can extract individual terms from the power series ( [ eq : minkcorrcon ] ) . there is , however , for some of the functionals a relatively large offset . it is not clear if this offset is due to residual bias in the fitting method , or if our method to estimate the correlation function of a mock does not include the boundaries in the same way than the mfs . [ fig : threeptcoeffs ] then shows the @xmath184 . the connection to other known quantities is not as apparent as in the two - point case and we did not obtain an independent estimate of the three - point correlation integrals . what is shown in fig . [ fig : threeptcoeffs ] instead is the three - point function model of @xcite . it is based on leading - order perturbation theory and includes redshift space distortions . for the plots in fig . [ fig : threeptcoeffs ] we fixed the three bias parameters to @xmath207 , @xmath208 and @xmath209 . the @xmath210-value is the one we used in fig . [ fig : twoptcoeffs ] and the @xmath211-value is obtained from eq . ( 23 ) in @xcite and based on local lagrangian biasing @xcite . we did not attempt to fit for @xmath212 and so the lines in fig . [ fig : threeptcoeffs ] do not contain any free parameter . given that , the match of the model with the extracted coefficients is encouraging . in addition to being related to the hierarchy of @xmath7-point correlation functions , the mfs also contain information about the cosmic web . for example , the partial functionals probe the local morphology around the galaxy sites , and for an isodensity contour smoothing `` shapefinders '' have been defined to quantify the amount of filamentary vs. planarity of the local density field @xcite . for the global germ - grain mf densities that we study here , there are different quantities that can be extracted . the euler characteristic @xmath213 , for example , indicates for large balls how many extreme voids ( i.e. regions that are further away from any galaxy than @xmath92 ) there are in the galaxy distribution . the area functional @xmath89 gives an idea of their area and @xmath45 their volume . more generally , the comparison of the void fraction @xmath46 to the one of a poisson distribution of the same density indicates how inhomogeneously these void regions are distributed .    in fig . [ fig : freevolumefraction ] , we therefore show @xmath46 as a function of the ball diameter and compare it to the corresponding poisson case .     for a ball of diameter @xmath159 , where the first values correspond to the lower ( blue ) points and the values restarting at @xmath197 to the upper ( red ) points . the densities are @xmath214 and @xmath215 . especially for the higher density the solid poisson lines lie way below the values for the data , so the structure leads to a significantly higher void fraction than we would expect in an unclustered field.[fig : freevolumefraction],scaledwidth=45.0% ]    especially for the higher density of @xmath216 the difference is quite remarkable . for a diameter of @xmath128 , the volume is a factor of @xmath217 larger than in the poisson case . for the most extreme case of @xmath132 , it is nearly @xmath218 . in this latter case the expected number of galaxies in a ball of diameter @xmath132 would be @xmath219 . with the standard poisson fluctuation one would expect a fraction of @xmath220 to be uncovered . including the observed structure however , there is @xmath221 of the @xmath222 in our analyzed region that remain uncovered . of course one has to make sure that this conclusion is not significantly affected by the holes in the survey . we tested this using the periodic cubic md - patchy mocks and plot the resulting void fractions as purple and green points in fig . [ fig : freevolumefraction ] . for the region that we evaluated , the cubic mocks for the difference is small , below @xmath223 for @xmath224 and below @xmath225 for @xmath226 . however , for the more extreme diameters the deviations will grow , but we do not expect that this affects the general picture . in this paper , we perform a new analysis that combines minkowski functionals with the standard two - point correlation function to measure the amount of non - gaussianity of the density field . we find in sec . [ sub : highsubtract ] that the significance of the higher - order contribution to the mfs is of the order of @xmath227 ( see tab . [ tab : chisqradcomb ] and fig . [ fig : modminksubtract ] ) . this demonstrates that this part of the mfs can be well enough measured to give meaningful constraints on the higher - order properties of the galaxy distribution . given this constraining power it is reassuring that we find from fig . [ fig : modminksubtract ] and tab . [ tab : chisqmockdat ] that the md - patchy mocks , that were constructed to match the two- and three - point properties of the measured galaxy distribution , match fairly well to the higher - order mfs . this indicates that as far as the higher - orders accessible by the mfs are concerned , the mock galaxy catalogs are quite good . there is still some room for improvement as the combined deviation is as large as @xmath228 for some functionals ( see tab .  [ tab : chisqradcomb ] ) , but this might be addressed by a @xmath193 shift in the normalization of the mock power spectrum ( see fig .  [ fig : modminkdensdep ] ) .    on the technical side , our careful boundary treatment minimizes the effects of the survey geometry on our measurement of the modified mfs ( see fig .  [ fig : maskcorrection ] ) . the corrected values lie within the error bars of a boundaryless periodic box . for the qpm mocks , we find a much larger deviation ( beyond @xmath229 ) from the data than for md - patchy ( see tab .  [ tab : chisqmockdat ] ) . as in this case already the two - point correlations are overnormalized in the range of scales that we are interested in ( see fig . [ fig : meancorrfun ] ) , a reduction of the amplitude of the initial gaussian field could also help in this case . we show that in the the higher - order part of the mfs calculated for the sdss dr12 cmass galaxies contains at least information up to the integrated @xmath9-point function ( see tab . [ tab : higherfitchi ] ) . this underlines the power of the mfs to bring condensed higher - order information into a simple functional form . aw thanks antonio cuesta for providing support with the qpm correlation functions . we thank graziano rossi and thomas buchert for valuable comments on the manuscript . the work of aw was supported by the german research organization dfg , grant no . wi 4501/1 - 1 . dje is supported by u.s . department of energy grant de - sc0013718 and as a simons foundation investigator . funding for sdss - iii has been provided by the alfred p. sloan foundation , the participating institutions , the national science foundation , and the u.s . department of energy office of science . the sdss - iii web site is http://www.sdss3.org/.    sdss - iii is managed by the astrophysical research consortium for the participating institutions of the sdss - iii collaboration including the university of arizona , the brazilian participation group , brookhaven national laboratory , carnegie mellon university , university of florida , the french participation group , the german participation group , harvard university , the instituto de astrofisica de canarias , the michigan state / notre dame / jina participation group , johns hopkins university , lawrence berkeley national laboratory , max planck institute for astrophysics , max planck institute for extraterrestrial physics , new mexico state university , new york university , ohio state university , pennsylvania state university , university of portsmouth , princeton university , the spanish participation group , university of tokyo , university of utah , vanderbilt university , university of virginia , university of washington , and yale university . t. , 1995 , in mcket j.  p. , gottlber s. , mller v. , eds , large scale structure in the universe robust morphological measures for large scale structure . p.  156 , http://arxiv.org/abs/astro-ph/9412061[astro-ph/9412061 ] , http://adsabs.harvard.edu/abs/1995lssu.conf..156b[ads ]                                    hadwiger h. , 1957 , vorlesungen ber inhalt , oberflche und isoperimetrie . die grundlehren der mathematischen wissenschaften in einzeldarstellungen mit besonderer bercksichtigung der anwendungsgebiete ; 93 , springer , http://cosmologist.info/isbn/3-540-02151-5 , [ isbn : 3 - 540 - 02151 - 5 , ]            m. , schmalzing j. , buchert t. , wagner h. 1996 , in weiss a. , raffelt g. , hillebrandt w. , von feilitzsch f. , buchert t. , eds , astro - particle physics the significance of the fluctuations in the iras 1.2 jy galaxy catalogue .. pp 8398 , http://adsabs.harvard.edu/abs/1996app..conf...83k[ads ]                    r. et  al . , 2001 , in harnden jr . f.  r. , primini f.  a. , payne h.  e. , eds , astronomical data analysis software and systems x vol . 238 of astronomical society of the pacific conference series , the sdss imaging pipelines . p.  269 , http://arxiv.org/abs/astro-ph/0101420[astro-ph/0101420 ] , http://adsabs.harvard.edu/abs/2001aspc..238..269l[ads ]                    m.  c. et  al . , 2011 , in american astronomical society meeting abstracts # 218 the log - density as a better cosmological density variable . p. # 233.04 , http://adsabs.harvard.edu/abs/2011aas...21823304n[ads ]                    m. , buchert t. 1996 , in weiss a. , raffelt g. , hillebrandt w. , von feilitzsch f. , buchert t. , eds , astro - particle physics applications of minkowski - functionals to the statistical analysis of dark matter models . , http://arxiv.org/abs/astro-ph/9509014[astro-ph/9509014 ] , http://adsabs.harvard.edu/abs/1996app..conf..251p[ads ]                      j. , kerscher m. , buchert t. 1996 , in bonometto s. , primack j.  r. , provenzale a. , eds , dark matter in the universe minkowski functionals in cosmology . p.  281 , http://arxiv.org/abs/astro-ph/9508154[astro-ph/9508154 ] , http://adsabs.harvard.edu/abs/1996dmu..conf..281s[ads ]
Abstract:

Answer: we probe the higher - order clustering of the galaxies in the final data release ( dr12 ) of the sloan digital sky survey baryon oscillation spectroscopic survey ( boss ) using the method of germ - grain minkowski functionals ( mfs ) . our sample consists of @xmath0 boss galaxies from the northern galactic cap in the redshift range @xmath1@xmath2 . we show the mfs to be sensitive to contributions up to the six - point correlation function for this data set . we ensure with a custom angular mask that the results are more independent of boundary effects than in previous analyses of this type . we extract the higher - order part of the mfs and quantify the difference to the case without higher - order correlations . the resulting @xmath3 value of over @xmath4 for a modest number of degrees of freedom , @xmath5 , indicates a @xmath6-sigma deviation and demonstrates that we have a highly significant signal of the non - gaussian contributions to the galaxy distribution . this statistical power can be useful in testing models with differing higher - order correlations . comparing the galaxy data to the qpm and multidark - patchy mocks , we find that the latter better describes the observed structure . from an order - by - order decomposition we expect that , for example , already a reduction of the amplitude of the md - patchy mock power spectrum by 5% would remove the remaining tension . [ firstpage ]    methods : data analysis  methods : statistical  cosmology : large - scale structure of universe  cosmology : observations