format: text

text:
# Measuring the Geometry and Topology of Large Scale Structure using SURFGEN: Methodology and Preliminary Results

## 1 Introduction

The existence of the supercluster-void network is one of the most intriguing observational features of our universe. Redshift surveys of galaxies confirm that on very large scales the bulk of matter in the universe is concentrated in clusters and superclusters of galaxies which are separated by large almost empty regions, appropriately called voids. At moderate density thresholds $`\delta 1`$, the supercluster network percolates even though it occupies a tiny fraction of the total volume, and so has a small filling fraction.

The morphology of the supercluster-void network is quite complex and has inspired evocative descriptions such as being ‘honey-comb-like’, ‘bubble-like’, ‘a filamentary network’, ‘swiss cheese’, ‘cosmic web’ etc. Indeed the intricate weave of large scale structure is becoming increasingly evident as results from progressively larger galaxy redshift surveys (CfA, IRAS, PSCz, LCRS, 2dFGRS, SDSS) demonstrate. Quantifying the properties of the supercluster-void network is clearly one of the cardinal tasks facing cosmology today. Since the network evolved from an almost featureless Gaussian random field (as suggested by observations of the cosmic microwave background) it is important to understand how such a rich and complex cosmic tapestry – characterized by prominent non-Gaussian features (clusters, superclusters, voids) – could have arisen from small and statistically random initial conditions. Furthermore, since the network is highly evolved it is unlikely that its principle features can be described within a strictly perturbative framework which breaks down when $`\delta \text{ }\stackrel{>}{}1`$. It is therefore both interesting and revealing that approximations which successfully describe non-linear dynamics (such as the Zeldovich approximation and the Adhesion model) generically predict the formation of filaments and pancakes, which interweave and percolate to form the large scale structure of the universe (Shandarin & Zeldovich, 1989; Gurbatov, Saichev & Shandarin, 1985, 1989). Furthermore, this interweaving network of filaments, sheets and voids – predicted both by the Zeldovich approximation and the adhesion model – is readily seen in N-body simulations of gravitational clustering for a wide variety of initial power spectra (Klypin & Shandarin, 1993; Sathyaprakash, Sahni & Shandarin, 1996). It therefore appears that a filamentary distribution of large scale structure is an almost generic outcome of pressureless (dark) matter clustering from Gaussian initial conditions. (A theoretical model for understanding the morphology of large scale structure can be found in the “cosmic web” hypothesis of Bond, Kofman & Pogosyan (1996). Insightful overviews of the statistical and dynamical techniques used in studies of large scale structure can be found in: Sahni & Coles (1995); Martinez & Saar (2001); van de Weygaert (2002).)

A number of statistical measures have been advanced to quantify the pattern made by galaxies as they cluster in our universe. Prominent among them are Percolation Analysis (Zeldovich, Einasto & Shandarin, 1982; Shandarin & Zeldovich, 1983), Counts in Cells (Janes & Demarque, 1983; de Lapparent, Geller & Huchra, 1991), Minimal Spanning Trees (Barrow, Sonoda & Bhavsar, 1985), the Genus measure (Gott, Melott & Dickinson, 1986) etc. (To some extent these methods complement the traditional approach to quantify clustering using the hierarchy of correlation functions (Peebles, 1980), which become cumbersome to evaluate beyond the three point function.) Recently Mecke, Buchert and Wagner (1994) have introduced the Minkowski functionals (hereafter MFs) to cosmology. MFs contain information pertaining to geometry, connectivity (percolation) and topology (genus). In addition, the ratio’s of Minkowski functionals quantify the morphology of large scale structure. This has led to the construction of a set of measures called Shapefinders which tell us whether the distribution of matter in superclusters/voids is spherical, planar, filamentary etc. (Sahni, Sathyaprakash & Shandarin, 1998; Sathyaprakash, Sahni & Shandarin, 1998). (Applications of Shapefinders to galaxy catalogues and N-body simulations has been discussed in Basilakos, Plionis & Rowan-Robinson (2000); Kolokotronis, Basilakos & Plionis (2001). Moment-based studies of morphology predating the Shapefinders can be found in Babul & Starkman (1992); Luo & Vishniac (1995).) Thus between them, the four Minkowski functionals contain valuable information regarding both the geometrical as well as topological distribution of matter in the universe. Since Minkowski functionals are additive in nature one can glean information regarding both individual objects (galaxies, clusters, voids) as well as describe the supercluster-void network in its totality.

There have been three major attempts made to study the morphology of the LSS using MFs; these efforts differ in their approach of evaluating the MFs: (1) Firstly, Boolean grain models study the MFs of surfaces which result due to intersecting spheres decorating the input point-set (Mecke, Buchert $`\&`$ Wagner 1994). (2) Secondly, Krofton’s formulae make it possible to calculate MFs on a density field defined on a grid (Schmalzing $`\&`$ Buchert 1997). In this case the MFs are calculated by using the information of the number of cells in 1D (vertices), 2D (faces) and 3D (cuboids). (3) Finally, an alternative, resolution-dependent approach consists in employing the Koenderink invariants (Schmalzing $`\&`$ Buchert 1997, Schmalzing et al. 1999).

This paper introduces a radically new and powerful approach for computing the Minkowski functionals. This approach consists of constructing iso-density surfaces using an elaborate surface modelling scheme defined in terms of excursion sets of a density field. These surfaces are triangulated and the MFs are evaluated for the resulting closed polyhedral surface. A similar algorithm has been used in studies of two-dimensional CMB maps (Novikov, Feldman & Shandarin, 1999; Shandarin et al., 2002). However, developing a three-dimensional algorithm requires a new component that builds a surface of approximately constant (to linear order) density. It has been implemented and reported in this paper.

The MFs so constructed are tested against known results for simply connected surfaces (triaxial ellipsoids) as well as multiply connected surfaces (triaxial tori). In all cases we find that our ansatz for determining MFs on polyhedra gives results which are in excellent agreement with exact continuum values for these quantities. We have further tested the performance of SURFGEN against Gaussian random fields for which the analytic prediction for MFs is available and find that our calculations reproduce the analytic results to a remarkable precision. This gives us confidence that our ansatz is well suited for application to large N-body simulations and three-dimensional galaxy redshift surveys. Encouraged by the success on these preliminary tests, we apply our method to simulations of three cosmological models due to the Virgo consortium.

The rest of this paper is arranged as follows: In Section 2, we provide a brief overview of Minkowski functionals and present an algorithm which computes MFs for a triangulated polyhedral surface. In Section 3, we describe the SURFGEN algorithm which triangulates a general class of surfaces including iso-density surfaces describing superclusters or voids. In Section 4 we test our method against standard eikonal surfaces $``$ spheres, ellipsoids, tori $``$ and demonstrate the accuracy of our algorithm. Special boundary conditions are needed to reproduce the analytic expressions for MFs of a Gaussian random field. The related discussion along with the results for Gaussian random fields is presented at the end of the paper in an appendix. Section 5 is devoted to a detailed morphological study of cosmological N-body simulations. The models investigated are $`\mathrm{\Lambda }`$CDM, $`\tau `$CDM and SCDM at the present cosmological epoch. Our conclusions are presented in Section 6 which also discusses possible further applications of SURFGEN and the MFs.

## 2 Minkowski Functionals and Triangulated Surfaces

The large scale structure of the universe can be studied on various scales by considering the geometry and topology of excursion sets of the density contrast $`\delta (𝐱)`$ defined as $`_{TH}^+=\{𝐱|\delta (𝐱)\delta _{TH}\}`$, for overdense regions (clusters, superclusters) and $`_{TH}^{}=\{𝐱|\delta (𝐱)\delta _{TH}\}`$, for underdense regions (voids). By specifying a given density threshold one effectively defines an iso-density surface over which the Minkowski functionals should be evaluated. Depending upon whether the surface encloses overdense region of the space or underdense region, the surface refers to a cluster/supercluster or a void. Note that the term ‘cluster’ is used for any region in the excursion set ($`\delta >\delta _{TH}`$) which in general is different from the astronomical definition. The terminology reflects the cluster analysis jargon. The following four Minkowski functionals describe the morphological properties of an iso-density surface in three dimensions 111There are ($`n+1`$) Minkowski functionals in $`n`$ dimensions each defined for an ($`n1`$) dimensional hypersurface. Two dimensional MFs have been used to analyze the LCRS slices by Bharadwaj et al. (2000) and the anisotropy of the cosmic microwave background by Schmalzing & Gorski (1998) and by Novikov, Feldman & Shandarin (1999).

* Area $`S`$ of the surface,
* Volume $`V`$ enclosed by the surface,
* Integrated mean curvature $`C`$ of the surface (or integrated extrinsic curvature),

$$C=\frac{1}{2}\left(\frac{1}{R_1}+\frac{1}{R_2}\right)𝑑S,$$
(1)
where $`R_1`$ and $`R_2`$ are the principal radii of curvature at a given point on the surface.
*  Integrated intrinsic (or Gaussian) curvature $`\chi `$ of the surface – also called the Euler characteristic

$$\chi =\frac{1}{2\pi }\left(\frac{1}{R_1R_2}\right)𝑑S.$$
(2)

A related measure of topology is the genus $`G=1\chi /2`$. 222Genus measures the number of handles in excess of the number of isolated underdense regions (voids) that the surface of a cluster exhibits. Multiply connected surfaces have $`G>0`$ while $`G=0`$ for a simply connected surface such as a sphere. (The topological properties of all orientable surfaces are equivalent to those of a sphere with one or more ‘handles’. Thus a torus is homeomorphic to a sphere with one handle, while a pretzel is homeomorphic to a sphere with two handles etc.) While the genus provides information about the connectivity of a surface, the remaining three MFs are sensitive to local surface deformations and hence characterize the geometry and shape of large scale structure at varying thresholds of the density (Sahni, Sathyaprakash & Shandarin, 1998).

In nature one seldom comes across surfaces that are perfectly smooth and differentiable. As a result expressions  -  which make perfect sense for manifolds $`𝒞^n`$, $`n2`$ are woefully inadequate when it comes to determining the Minkowski functionals for real data sets which are grainy and quite often sparse. Below we describe how one can determine the Minkowski functionals (and derived quantities) for surfaces derived from real data by interpolation. The interpolation scheme which describes iso-density contours as triangulated surfaces will be presented in detail in the next Section.

We construct a polyhedral surface using an assembly of triangles in which every triangle shares its sides with each of its three neighboring triangles.

* The area of such a triangulated surface is

$$S=\underset{i=1}{\overset{N_T}{}}S_i,$$
(3)
where $`S_i`$ is the area of the $`i`$th triangle and $`N_T`$ is the total number of triangles which compose a given surface.
* The volume enclosed by this polyhedral surface is the summed contribution from $`N_T`$ tetrahedra

$`V`$ $`=`$ $`{\displaystyle \underset{i=1}{\overset{N_T}{}}}V_i,`$
$`V_i`$ $`=`$ $`{\displaystyle \frac{1}{3}}S_i(n_j\overline{P}^j)_i.`$ (4)
Here $`V_i`$ is the volume of an individual tetrahedron whose base is a triangle on the surface. $`(n_j\overline{P}^j)`$ is the scalar product between the outward pointing normal $`\widehat{n}`$ to this triangle and the mean position vector of the three triangle vertices, for which the $`j^{th}`$ component is given by

$$\overline{𝐏}^𝐣=\frac{1}{3}(P_1^j+P_2^j+P_3^j).$$
(5)
The subscript $`i`$ in (2) refers to the $`i`$th tetrahedron, while the vectors $`P_1,P_2,P_3`$ in (5) define the location of each of three triangle vertices defining the base of the tetrahedron relative to an (arbitrarily chosen) origin (Fig. 1). The vertices are ordered anticlockwise. This ensures that the contribution to the volume from tetrahedra whose base triangles lie on opposite sides of the origin add, while the volumes of tetrahedra with base triangles lying on the same side of the origin subtract out. Thus for both possibilities we get the correct value for the enclosed volume (Fig. 1).
* The extrinsic curvature of a triangulated surface is localized in the triangle edges. As a result the integrated mean curvature $`C`$ is determined by the formula

$$C=\frac{1}{2}\underset{i,j}{}\mathrm{}_{ij}.\varphi _{ij}.ϵ,$$
(6)
where $`\mathrm{}_{ij}`$ is the edge common to adjacent triangles $`i`$ and $`j`$ and $`\varphi _{ij}`$ is the angle between the normals to these triangles $`\widehat{n}_i`$ & $`\widehat{n}_j`$.

$$\mathrm{cos}\varphi _{ij}=\widehat{n}_i.\widehat{n}_j.$$
(7)
The summation in (6) is carried out over all pairs of adjacent triangles. It should be noted that for a completely general surface, the extrinsic curvature can be positive at some (convex) points and negative at some other (concave) points on the surface. To accommodate this fact one associates a number $`ϵ=\pm 1`$ with every triangle edge in (6). $`ϵ=1`$ if the normals on adjacent triangles diverge away from the surface, indicating a locally convex surface, while $`ϵ=1`$ if the normals converge towards each other outside the surface, which is indicative of a concave surface. In the former case the centre of curvature of the surface lies within the surface-body whereas in the latter case the centre of curvature lies outside of the surface-body (Fig. 2).
* The genus of a triangulated closed polyhedral surface is given by the convenient expression

$`G`$ $`=`$ $`1{\displaystyle \frac{\chi }{2}},`$
$`\chi `$ $`=`$ $`N_TN_E+N_V,`$ (8)
where $`\chi `$ is the Euler characteristic of the triangulated surface. $`N_T,N_E,N_V`$ are, respectively, the total number of triangles, triangle-edges, and triangle-vertices defining the surface.
* As demonstrated in Sahni, Sathyaprakash & Shandarin (1998); Sathyaprakash, Sahni & Shandarin (1998) ratio’s of Minkowski functionals called “Shapefinders” provide us with an excellent measure of morphology. Therefore, in addition to determining MFs we shall also determine the Shapefinders, $`T`$ (Thickness), $`B`$ (Breadth) and $`L`$ (Length) which are defined as follows:

$$T=\frac{3V}{S},B=\frac{S}{C},L=\frac{C}{4\pi (G+1)},$$
(9)
where $`G=0`$ for simply connected surfaces and $`G>`$ 0 for multiply connected regions. The three Shapefinders (describing an object) have dimensions of length and provide us with an estimate of the ‘extension’ of the object along each of the three spatial directions: $`T`$ is the shortest and thus describes the characteristic thickness of the object, $`L`$ is the longest and thus characterizes the length of the object; $`B`$ is intermediate and can be associated with the breadth of the object. This simple interpretation obviously is relevant only for fairly simple shapes. For example, a triaxial ellipsoid has the values of $`T`$, $`B`$ and $`L`$ closely related to the lengths of the three principal axes: shortest, intermediate and the longest respectively. It is worth stressing that $`L`$ defined by eq. (9) quantifies the characteristic length between holes which should be taken into account when interpreting the results (e.g.Table LABEL:tab:top10s123 below).
* An excellent indicator of ‘shape’ is provided by the dimensionless Shapefinder statistic

$$P=\frac{BT}{B+T};F=\frac{LB}{L+B},$$
(10)
where $`P`$ and $`F`$ are measures of Planarity and Filamentarity respectively ($`P,F1`$). A sphere has $`P=F=0`$, an ideal filament has $`P=0,F=1`$ while $`P=1,F=0`$ for an ideal pancake. Other interesting shapes include ‘ribbons’ for which $`PF1`$. When combined with the genus measure, the triplet $`\{P,F,G\}`$ provides an example of shape-space which incorporates information about topology as well as morphology of superclusters and voids.333Non-geometrical shape-statistics based on mass moments etc. can give misleading results when applied to large scale structure, as demonstrated in Sathyaprakash, Sahni & Shandarin (1998).

Having presented an ansatz to calculate MFs and Shapefinders on a triangulated surface, we now discuss the surface generating algorithm (SURFGEN) which creates triangulated surfaces corresponding to iso-density contours evaluated at any desired threshold of the density field.

## 3 Surface Generating Algorithm (SURFGEN)

SURFGEN is a versatile and powerful prescription which allows us to generate and study surfaces whose physical origin can be quite varied and different. Intensity surfaces (isophotes), isotherms and iso-density surfaces constructed from 3D data provide examples of two dimensional surfaces which can be generated and studied using SURFGEN and the Minkowski functionals discussed in the previous section444SURFGEN is a modified and extended version of the Marching Cubes Algorithm (MCA) which is used in the field of medical imaging to render high resolution images of internal organs by processing the output generated by X-ray tomography (Lorenson & Cline, 1987).. The present paper is the first of a series in which SURFGEN and MFs will be used in conjunction to make a detailed morphological study of the supercluster-void network in the universe. SURFGEN will be applied to a density field $`\rho (𝐫)`$ defined on a rectangular cubic lattice.

The SURFGEN algorithm is constructed as follows:

1. A point particle distribution (from N-body simulations, galaxy catalogues) is used to reconstruct the density field on the vertices of a cubic lattice via a cloud-in-cell (CIC) approach. The lattice itself consists of a large number of closely packed elementary cubes. Any given cube is characterized by the value of the density at its eight vertices.
2. An appropriate density threshold $`\rho _{\mathrm{TH}}`$ is chosen and lattice points at this threshold are found using linear interpolation. The lattice-cells which involve such interpolated points along their sides are marked for triangulation which is to be done in the next step.
3. SURFGEN determines a closed polyhedral surface at the threshold $`\rho _{TH}`$ by triangulating elementary cubes while simultaneously moving across the lattice. Triangulated elementary cubes are then put together to make up the full two-dimensional triangulated surface.

The process of triangulation is based on the following observation.

For any arbitrary threshold of the density, all elementary cubes of the density field fall into one of the following three classes:

(i) those which have all eight vertices below the density threshold (underdense cubes – Fig. 3$``$I),

(ii) those which have all eight vertices above the density threshold (overdense cubes – Fig. 3$``$II),

(iii) cubes having both overdense and underdense vertices (surface cubes – Fig. 3$``$III).

When modelling clusters and superclusters overdense cubes will be enclosed by our surface while underdense cubes will be excluded by it. (For voids the situation will be reversed.) Contour surfaces at a prespecified density threshold will lie entirely within surface cubes. Thus the properties of surface cubes are vitally important for this surface-reconstruction exercise.

We now describe how SURFGEN constructs a surface at the desired threshold $`\rho _{TH}`$ by triangulating surface cubes. We work under the assumption that the underlying density field $`\rho (𝐫)`$ is continuous, so that in moving from an overdense site to an underdense site (along a cube edge) we invariably encounter a point at which $`\rho =\rho _{TH}`$. SURFGEN classifies a given lattice-cube in terms of the number of points (on the edges) where $`\rho =\rho _{TH}`$. The precise location of these points (found by interpolation) tells us exactly where our iso-density surface will intersect the cube. Since there are eight vertices and each vertex can either be overdense (1) or underdense (0), the number of possible configurations of the cube are $`2^8=256`$. Clearly the triangulation must be invariant upon interchanging 1’s with 0’s, and this reduces the number of independent configurations to 128. Upon invoking rotational symmetry this number further reduces to only 14. (For instance, although there are eight cubes which have a single overdense vertex, any two members of this family are related via the three dimensional rotation group. Therefore a single scheme of triangulation suffices to describe all members of this group.)

These fourteen configurations of triangulated cubes are shown in Fig. 4 (also see Lorenson & Cline (1987)).

The surface intersecting a given surface cube forms quite clearly only a portion of the full iso-density surface which we are interested in constructing. By explicitly demanding that the density field be continuous across the iso-density surface one joins triangles across the faces of neighboring cubes, thus constructing the full continuous polyhedral surface. This then is the basic prescription which may be used to construct a closed triangulated surface at any desired value of the density threshold. However, a few comments about triangulation are in order:

(i) the triangulation of the cubes having four contrasting vertices is not unique. To illustrate this let us consider Fig. 5 where the last cube belonging to the middle row of Fig. 4 is reproduced in the left panel. We note that this cube is symmetric with respect to an interchange of underdense sites with overdense sites. As a result, we could just as well triangulate it according to the prescription shown in the right panel of Fig. 5. This degeneracy in the way a given cube can be triangulated is inherent also in the last two cubes of the last row of Fig. 4. Noticeably, at all such places where we encounter one of these cubes, one scheme of triangulation has to be preferred over the other, otherwise the surface could become discontinuous. Thus to create a unique triangulated surface one must triangulate such cubes in tandem with their neighbors, ensuring that the continuity of the surface is maintained.

(ii) There are instances when two neighboring cubes have a density configuration which does not lead to the complete closure of the given surface. The resulting surface has a hole (Fig. 6a) which must be filled by constructing two additional triangles which close the surface (Fig. 6b). The probability of hole formation is finite every time the common face shared by two cubes shows two vertices of type 1 along one diagonal and of type 0 along the other diagonal.

Thus, in order to generate a closed polyhedral surface at a given value of the density threshold $`\rho =\rho _{TH}`$ SURFGEN uses the 14 independent triangulations of an elementary lattice cube in conjunction with an algorithm which ensures that the surface is continuous. We should emphasize that surfaces generated in this manner need not be simply connected (see Fig. 7). Indeed iso-density surfaces constructed at moderate density thresholds in N-body simulations as well as 3D galaxy surveys tend to show a large positive value for the genus (Eq.2). This issue will be discussed in detail in forthcoming paper(s). To enable an online calculation of the MFs, our surface modelling scheme also adheres to the following requirements:

* The vertices of all the triangles are stored in an anticlockwise order. This enforces a uniform prescription on normals which always point outward from the surface being modelled. Information regarding the directionality of normals is of great importance since it is used both for calculating the volume as well as the mean curvature.
* For online calculation of the mean curvature, information concerning a given triangle must be supplemented with information concerning the triplet of triangles which are its neighbors. This allows us to unambiguously determine the contribution to the mean curvature from a local neighborhood.
* The total number of triangles, triangle-edges and triangle-vertices of a triangulated surface are counted which enable us to determine its genus when surface construction is complete.

We might add that, although in this paper the density field is reconstructed from N-body particle positions by means of a CIC prescription, this is not a necessary prerequisite for SURFGEN which is versatile enough to be used to triangulate iso-density contours of fields reconstructed using more elaborate techniques such as Delauney tesselations (Schaap & Weygaert, 2000; van de Weygaert, 2002), Wiener reconstruction or adaptive smoothing. As it turns out, these schemes tend to maintain the structural complexity in the system without diluting the patterns in an isotropic manner. However, since SURFGEN constructs a triangulated surface by ‘marching across neighboring lattice cubes’ (Lorenson & Cline, 1987) the density field must be specified on a cubic lattice for this approach to give meaningful results. 555We should emphasize however that the ansatz developed in section 2 for determining the Minkowski functionals for a triangulated surface can easily be adapted to other triangulation schemes such as those described in O’Rourke (1998); van de Weygaert (2002).

To summarize, our surface modelling code incorporates the triangulation of the entire set of 256 possible configurations of a cube within a single scheme and makes possible online computation of the area, the volume, the integrated mean curvature and the genus of a triangulated surface. The code has been tested on a variety of standard density distributions and eikonal morphologies as well as on Gaussian random fields. In the next section, we present results based on this analysis.

## 4 Results for Standard Eikonal Surfaces

In order to test the accuracy of our ansatz we generate triangulated surfaces whose counterpart continuum surfaces have known (analytically calculable) Minkowski functionals.

### 4.1 Spherical structures

In the first exercise, we consider spherically symmetric density distributions and generate surfaces of constant density for a variety of density-thresholds starting from a chosen maximum radius down to grid-size. A convenient density distribution for this exercise is

$`\rho (i,j,k)=\{\begin{array}{cc}\frac{\rho _0}{R},\hfill & (i,j,k)(i_0,j_0,k_0),\hfill \\ \rho _0,\hfill & (i,j,k)=(i_0,j_0,k_0),\hfill \end{array}`$ (13)
where R is the distance between (i,j,k) and the centre.

$$R=\sqrt{(ii_0)^2+(jj_0)^2+(kk_0)^2}.$$
(14)
Since the density field falls off as the inverse of the distance from the centre at ($`i_0,j_0,k_0`$), any threshold $`\rho _{TH}`$ is associated with a sphere of radius $`R=\frac{\rho _0}{\rho _{TH}}`$. Thus larger spheres correspond to surfaces of lower constant density in this model. We assume $`i_0=j_0=k_0=15`$, for the centre of the sphere in our numerical calculations.

Applying the Shapefinders (9) & (10) to a sphere of radius $`R`$ one finds the simple result $`T=B=L=R`$ and $`F=P=0`$.

Figure 8 compares the area, volume, mean curvature and Shapefinders measured for a triangulated sphere against known analytical values for these quantities. This figure clearly demonstrates that exact values, and values computed using triangulation, match exceedingly well down to the very lowest scale.

We should mention that the Genus evaluated using Eq.2 is identically zero for spheres of all radii and for all possible deformations of an ellipsoid (to be discussed next). This provides an excellent independent endorsement of our methodology by demonstrating that our triangulations of a sphere and ellipsoid indeed result in closed, continuous surfaces.

### 4.2 Triaxial ellipsoid

A triaxial ellipsoid is an excellent shape with which to test a morphological statistic. This is because, depending upon the relative scales of the three axes, a triaxial ellipsoid can be oblate, prolate or spherical. We saw previously that MFs and Shapefinders give extremely accurate results for spherical surfaces. We now demonstrate that this remains true even for surfaces which are highly planar or filamentary.

The parametric form of an ellipsoid having axes $`a,`$ $`b,`$ $`c`$ and volume $`V=\frac{4\pi }{3}abc`$ is

$$𝐫=a(\mathrm{sin}\theta \mathrm{cos}\varphi )\widehat{x}+b(\mathrm{sin}\theta \mathrm{sin}\varphi )\widehat{y}+c(\mathrm{cos}\theta )\widehat{z}$$
(15)
where $`0\varphi 2\pi `$, $`0\theta \pi `$. For the purposes of our study, we systematically deform a triaxial ellipsoid and study how its shape evolves in the process. When shrunk along a single axis a triaxial ellipsoid becomes planar. Simultaneous shrinking along a second axis makes it cigar-like or filamentary.

#### 4.2.1 Oblate spheroids

We first study the accuracy of our triangulation scheme for planar configurations by considering oblate deformations of an ellipsoid. In this case two axes $`a`$ and $`b`$ are held fixed ($`ab`$), while the third axis $`c`$ is slowly shrunk leading to an increasingly planar surface.

Our results for this case are compiled in two figures. Figure 9 shows MFs as they evolve with the dimensionless variable $`c/a`$. (The exact results for the Minkowski functionals which we quote are based on the analytical expressions for MFs given in Sahni, Sathyaprakash & Shandarin (1998), we refer the reader to that paper for more details.) Figure 9 clearly demonstrates that the values of MFs obtained using triangulation match the exact values to a remarkable degree of accuracy. Indeed, for a wide range in $`c/a`$ the two distinctly different calculational algorithms give virtually indistinguishable results, thereby indicating that the triangulation ansatz is, for all practical purposes, exact!

Figure 10 shows the evolution of all the Shapefinders, i.e., $`T,B,L,P`$ and $`F`$ together with the percentage errors in the estimation of dimensional Shapefinders. We find that $`L`$ is estimated to greatest accuracy with maximum error of $`0.4\%`$, while $`T`$ and $`B`$ can be determined to an accuracy of $`0.8\%`$. We further note (right panels of Figure 10) that Planarity grows from an initially low value $`0.0`$ to a large final value $`0.4`$ as the ellipsoid becomes increasingly oblate. Filamentarity, on the other hand remains small at $`0.09`$. Both Filamentarity and Planarity are determined to great accuracy by SURFGEN.

#### 4.2.2 Prolate spheroids

Next we study prolate deformations of our oblate ellipsoid. We start with $`ba,`$ $`ca,`$ and shrink the second axis $`b`$ while keeping $`a`$ and $`c`$ fixed, so that finally $`cba`$ and our initially oblate ellipsoid becomes prolate. Our results are again summarized in two figures.

Figure 11 reveals very good agreement between measured and true values of MFs, with the former tending to be slightly smaller than the latter.

Turning to the Shapefinders we find that, with the possible exception of extremely prolate figures, the Shapefinders are remarkably well determined. Indeed, even the extremely prolate ellipsoid with axis ratio $`b/a<0.2`$ has a largest error in $`T`$ of only $`1.7\%`$ while errors in $`B`$ and $`L`$ never exceed $`1\%`$ (Fig. 12). This figure also shows the evolution of planarity ($`P`$) and filamentarity ($`F`$) as our ellipsoid becomes increasingly more prolate. In keeping with our expectation $`F`$ steadily increases from a small initial value $`0.1`$ to $`0.44`$. Planarity drops from its large initial value to a small final value $`P0.1`$ and is slightly underestimated for $`b/a\text{ }\stackrel{>}{}0.3`$. We therefore conclude that both $`P`$ and $`F`$ are determined by SURFGEN to a sufficient accuracy for prolate ellipsoidal figures.

### 4.3 The Torus and its deformations

Next we extend our analysis to manifolds which are multiply connected by considering the deformations of a torus which has an elliptical cross section and which we shall refer to as an elliptical torus. A torus is an important surface on which to test SURFGEN for two reasons: it is multiply connected and, unlike an ellipsoid, it contains regions which are convex (on its outside) as well as concave (on its inside). The elliptical torus can be described by three parameters $`a,b,c`$. In this case the elliptical toroidal tube has diameter $`2\pi b`$ and $`a`$ & $`c`$ are its radii of curvature in two mutually orthogonal directions. The elliptical torus reduces to the more familiar circular torus when $`a=c`$. We choose to work with the elliptical torus because changing the values of $`a,b,c`$ can give rise to a large variety of surfaces all of which (by definition) are multiply connected but which have very different shapes. Thus our surface modelling scheme SURFGEN can be put to a rigorous test. The parametric form for the elliptic torus is

$$𝐫=(b+c\mathrm{sin}\varphi )\mathrm{cos}\theta \widehat{x}+(b+c\mathrm{sin}\varphi )\mathrm{sin}\theta \widehat{y}+a\mathrm{cos}\varphi \widehat{z},$$
(16)
where $`a,c<b`$, $`0\varphi <2\pi `$ and $`\theta <2\pi `$. We shall compare results due to SURFGEN with the exact analytical results for four kinds of tori - a nearly spheroidal torus, a ribbon, a pancake, and a filament.

Tables LABEL:tab:mf and LABEL:tab:s123 refer to these four surfaces; relevant figures can be found in Sahni, Sathyaprakash & Shandarin (1998). Table LABEL:tab:mf shows the estimated values of the Minkowski Functionals and the percentage error in their estimation. Table LABEL:tab:s123 shows all the Shapefinders for these surfaces666It should be noted that the definition of the third Shapefinder $`L`$ is modified from that reported in Sahni, Sathyaprakash & Shandarin (1998). New definition of $`L`$ is as given in Eq.9. The newly adopted definition of $`L`$ incorporates treatment of simply connected and multiply connected surfaces in within a single scheme.. This table also shows the percentage errors in the estimation of the Shapefinders.

We should emphasize that the genus estimated by SURFGEN for all these deformations has the correct value of unity. This gives an important independent check on the self-consistency of the triangulation scheme since, the absence of even a single triangle out of several thousand, would create an artificial ‘hole’ in our surface and give the wrong value for the genus – a situation that has never occurred for any of the several dozen deformations of either the ellipsoid or the torus.

Our results clearly demonstrate that the surface modelling scheme SURFGEN provides values for the Minkowski functionals which are in excellent agreement with exact analytical formulae. We have further tested SURFGEN against a set of Gaussian random field samples. Owing to a subtle issue of the boundary conditions involved in the corresponding test, we present the corresponding results in an appendix at the end of this paper. For the time being, it suffices to make a note that our ansatz for studying iso-density contours and their morphology works exceedingly well both for simply connected as well as multiply connected surfaces and on Gaussian random fields.

In the next section, we apply SURFGEN to cosmological N-body simulations and give an example of the study of the morphology in three models of structure formation: $`\mathrm{\Lambda }`$CDM, $`\tau `$CDM and SCDM.

## 5 Applications: Morphology of the Virgo Simulations

As a first application of the SURFGEN code, we report here a morphological study of cosmological simulations. We use dark matter distributions of three cosmological models simulated by the Virgo consortium. The models are (1) a flat model with $`\mathrm{\Omega }_0`$ = 0.3 ($`\mathrm{\Lambda }`$CDM), and two models with $`\mathrm{\Omega }`$ = 1, namely (2) one with the standard CDM power spectrum (SCDM) and (3) a flat model with the same power spectrum as the $`\mathrm{\Omega }_0`$ = 0.3 $`\mathrm{\Lambda }`$CDM model (also referred to as $`\tau `$CDM). The shape parameter $`\mathrm{\Gamma }`$ (=$`\mathrm{\Omega }_0`$h) for SCDM is 0.5, whereas for the other two models $`\mathrm{\Gamma }=0.21`$. Amplitude of the power spectrum in all models is set so as to reproduce the observed abundance of rich galaxy clusters at the present epoch. A detailed discussion of the cosmological parameters and simulations can be found in Jenkins et al. (1998).

The data consist of $`256^3`$ particles in a box of size 239.5 $`h^1`$Mpc. In order to carry out a detailed morphological study, we first need to smooth the data and recover the underlying density field. The issue of density field reconstruction has been well explored by several groups studying the topology of the large scale structure. Here we follow the smoothing technique used by Springel et al. (1998) which they adopted for their preliminary topological analysis of the Virgo simulations. We fit a $`128^3`$ grid on to the box. Thus, the size of each cell is 1.875 $`h^1`$Mpc. Since there are on an average 8 particles per cell, the smoothing is not too sensitive to shot-noise. The samples are fairly rich and we can smooth the fields in two steps. In the first, we apply a Cloud in Cell (CIC) smoothing technique to construct a density field on the grid. Next we smooth this field with a Gaussian kernel which offers us an extra smoothing length-scale. There are various criteria to fix this smoothing scale which depend on the relative sparseness/richness of the samples. To avoid discreteness effect (so that there are no regions in the survey where density is undefined), the simplest criteria for sparse samples has been to set $`\lambda 2.5\mathrm{}_g`$, where $`\mathrm{}_g`$ is the grid-size. There are other sophisticated criteria to deal with sparse samples which are tested and reported in literature, which need not be reiterated here on account of our samples being rich. We are left with two choices; either to smooth the fields at a scale comparable to $`r_0`$, the correlation scale, or to probe the morphology at smaller length scales with an interest to study the regularity in the behaviour of the MFs; of which the latter is more promising. The average interparticle distance in the present case is $`0.94`$ h-1Mpc, which is far outweighed by the grid-size and the correlation length. Hence, we adopt a modest scale of smoothing, $`\lambda =2`$ h-1 Mpc, which is small compared to usual standards as well as enough to provide the necessary smoothing. The Gaussian kernel for smoothing that we adopt here is 777 Since the Gaussian kernel that we use is radially symmetric, it could, in principle, diminish the true extent of anisotropy in filaments and pancakes. This effect could be minimized by considering anisotropic kernels and/or smoothing techniques based on the wavelet transform. An even more ambitious approach is to reconstruct density fields using Delaunay tesselations (van de Weygaert, 2002). Density fields reconstructed in this manner appear to preserve anisotropic features and may therefore have some advantage over conventional ‘cloud-in-cell’ techniques followed by an isotropic smoothing (Schaap & Weygaert, 2000). Thus, as far as the goal is to utilize the geometry of the patterns to discriminate between the models, such smoothing schemes should prove more powerful. We hope to return to these issues in forthcoming paper(s).

$$W(r)=\frac{1}{\pi ^{\frac{3}{2}}\lambda ^3}exp\left(\frac{r^2}{\lambda ^2}\right).$$
(17)
In this section we first study the global MFs for all the three models. Next we investigate the morphology of the percolating supercluster network. This is followed by a statistical study of the morphology of smaller structures. We summarize our results in the next section.

### 5.1 Global Minkowski Functionals

We scan the density fields at 100 values of the density threshold $`\rho _{TH}`$, all equispaced in the filling factor

$$FF_\mathrm{V}=\mathrm{\Theta }(\rho \rho _{\mathrm{TH}})d^3x,$$
(18)
where $`\mathrm{\Theta }(x)`$ is the Heaviside Theta function. $`FF_\mathrm{V}`$ measures the volume-fraction in regions which satisfy the ‘cluster’ criterion $`\rho _{\mathrm{cluster}}\rho _{\mathrm{TH}}`$ at a given density threshold $`\rho _{\mathrm{TH}}`$. In the following, we use $`FF_\mathrm{V}`$ as a parameter to label the density contours.

At each level of the density field (labelled by $`FF_\mathrm{V}`$), we construct a catalogue of clusters (overdense regions) based on a Friends of Friends (FOF) algorithm.888The clusters/superclusters discussed in the present paper are defined as connected overdense regions lying above a prescribed density threshold. Due to the large smoothing scale adopted, the overdensity in clusters ranges from $`\delta 1`$ to $`\delta 10`$, which makes them more extended ( $`\stackrel{>}{}`$ few Mpc) and less dense than the galaxy clusters in (for instance) the Abell catalogue. Next we (i) run the SURFGEN code on each of these clusters to model surfaces for each of them and (ii) determine the Minkowski Functionals (MFs) for all clusters at the given threshold (these are referred to in the literature as partial MFs). Global MFs are partial MFs summed over all clusters. Thus, at each level of the density, we first compute the partial MFs and then the global MFs. Our plots of global MFs as functions of $`FF_\mathrm{V}`$

are shown in Figures 13. This figure shows some interesting features. The three cosmological models have appreciably different morphology and hence can be distinguished from one another on the basis of morphological measures. Similarly to the percolation statistic (Dominik & Shandarin, 1992) the difference between models is brought out much better if, instead of plotting the MFs against the volume fraction (equivalently ‘volume filling factor’), we choose the ‘mass fraction’ defined as $`FF_\mathrm{M}=M_{\mathrm{total}}^1\rho (\rho _{\mathrm{TH}})dV`$. (The ‘mass fraction’ can be thought of as the ’volume or mass filling factor’ 999Since the density in the Lagrangian space is uniform there is no difference between the volume and mass fractions. in (initial) Lagrangian space, while the ‘volume fraction’ is the filling factor in (final) Eulerian space.) By employing this parameter for studying the global MFs, we essentially probe the morphology of the iso-density contours (which may refer to different thresholds of density or density contrast, but) which enclose the same fraction of the total mass. Since clustering tends to pack up the mass into progressively smaller regions of space, such a study connects to aspects of gravitational clustering in a direct way. Thus, we see that at all the values of $`FF_\mathrm{M}`$ the volume that encloses the same fraction of the total mass is least for $`\mathrm{\Lambda }`$CDM and most for $`\tau `$CDM. A further advantage of employing this parameter is that the behaviour of the rest of the global MFs is no longer restricted to follow the same pattern (like it is in case of studying global MFs with the volume filling factor as the labelling parameter). Thus, not only the peak, but also the peak-position and the shape of the MF-curves are sensitive to the models being investigated. Thus, the mass fraction parametrisation could be more useful in discriminating the models from one another as well as comparing the models with the observations. Before we commence our analysis of individual objects (connected overdense/underdense regions) we would like to draw the attention of the reader to an issue in nomenclature. The ‘clusters’ and ‘superclusters’ referred to further in this paper are elements constituting the ‘cosmic-web’ and defined on the basis of isodensity contours in N-body simulations. They should not be mistaken for the clusters and superclusters of galaxies seen in the sky.

### 5.2 Cluster abundance and percolation

In this subsection we first discuss the percolation properties of the three density fields being studied. In this context we study how the total number of clusters and the fractional volume in the largest cluster vary as we scan through a set of density levels corresponding to equispaced fractions of volume and total mass.

Figure 15 shows the volume fraction and mass fraction as a function of the density contrast. It is interesting to note that SCDM and $`\tau `$CDM show the same pattern of behaviour which differs from $`\mathrm{\Lambda }`$CDM. It is further to note that at high density thresholds ($`\delta \text{ }\stackrel{>}{}1`$) the filling factor (at a constant density threshold) is greater in $`\mathrm{\Lambda }`$CDM than in SCDM or $`\tau `$CDM, while exactly the reverse is true for underdense regions. This figure, which relates $`FF_{\mathrm{M},\mathrm{V}}`$ to $`\delta `$, serves as a reference point to all of our subsequent percolation studies. Figure 15 also shows that, relative to other models, more mass occupies less space in $`\mathrm{\Lambda }`$CDM. Thus almost $`67\%`$ of the total mass in the $`\mathrm{\Lambda }`$CDM universe resides in just $`25\%`$ of the volume. (In case of $`\tau `$CDM and SCDM $`62\%`$ mass occupies $`34\%`$ volume.)

A statistical pair which help quantify the geometry of large scale structure and its morphology are the ‘Number of Clusters Statistic’ (NCS) and the ‘Largest Cluster Statistic’ (LCS); both are shown in Figure 16. The Number of Clusters Statistic shows cluster abundance as a function of the density contrast (lower left), mass fraction (lower middle) and volume fraction (lower right). The Largest Cluster Statistic measures the fractional cluster volume occupied by the largest cluster:

$$LCS=\frac{V_{\mathrm{LC}}}{_iV_i},$$
(19)
where $`V_i`$ is the volume of the $`i`$th cluster in the sample and $`V_{\mathrm{LC}}`$ is the volume occupied by the largest cluster. Summation is over all clusters evaluated at a given threshold of the density and includes the largest cluster. In Figure 16 LCS is shown as a function of the density contrast (upper left), mass fraction (upper middle) and volume fraction (upper right).

The square, triangle and circle in NCS denote critical values of the parameters $`\delta =\delta _{\mathrm{cluster}\mathrm{max}}`$, $`FF_\mathrm{V}=FF_{\mathrm{cluster}\mathrm{max}}`$ and $`FF_\mathrm{M}`$ at which the cluster abundance peaks in a give model. Similar etchings on the LCS curves denote values of $`\delta =\delta _{\mathrm{perc}}`$, $`FF_\mathrm{V}=FF_{\mathrm{perc}}`$ and $`FF_\mathrm{M}`$ at which the largest cluster first spans across the box in at least one direction. This is commonly referred to as the ‘percolation threshold’. It is important to note that the number of clusters in SCDM is considerably greater than the number of clusters in $`\mathrm{\Lambda }`$CDM at most density thresholds. It is also interesting that percolation takes place at higher values of the density contrast (and correspondingly lower values of the volume filling factor) in the case of $`\mathrm{\Lambda }`$CDM. ($`\delta _{\mathrm{perc}}2.3`$ for $`\mathrm{\Lambda }`$CDM, $`\delta _{\mathrm{perc}}1`$ for SCDM and $`\tau `$CDM.) Furthermore, while in the case of $`\mathrm{\Lambda }`$CDM $`\delta _{\mathrm{cluster}\mathrm{max}}\delta _{\mathrm{perc}}`$, in the case of SCDM and $`\tau `$CDM, $`\delta _{\mathrm{cluster}\mathrm{max}}\text{ }\stackrel{<}{}\delta _{\mathrm{perc}}`$. Thus in the latter two models, as the density threshold is lowered, the cluster abundance initially peaks, then, as the threshold is lowered further, neighboring clusters merge to form the percolating supercluster. In $`\mathrm{\Lambda }`$CDM on the other hand, the threshold at which the cluster abundance peaks also signals the formation of the percolating supercluster. Having said this we would like to add a word of caution: the analysis of one realization alone does not allow us to asses reliably the statistical fluctuations in the threshold of percolation. Therefore the final conclusion can be drawn only after a study of many realizations of each model.

Clearly both the LCS threshold $`\delta _{\mathrm{perc}}`$ and the NCS threshold $`\delta _{\mathrm{cluster}\mathrm{max}}`$ contain important information and most of our subsequent description of supercluster morphology will be carried out at one of these two thresholds. Table LABEL:tab:top10s123 shows the Minkowski functionals and associated Shapefinders for the 10 largest (most voluminous) superclusters compiled at the percolation threshold for the three cosmological models $`\mathrm{\Lambda }`$CDM, $`\tau `$CDM and SCDM while Table LABEL:tab:top10s123 lists values for the Planarity and Filamentarity for these superclusters along with the mass that these enclose and their genus value.

Our $`(239.5`$h-1Mpc$`)^3`$ $`\mathrm{\Lambda }`$CDM universe contains 1334 clusters & superclusters at the percolation threshold.

Of these the smallest is quasi-spherical with a radius of a few Megaparsec while the largest is extremely filamentary and percolates through the entire simulation box (see figure 17). From Figures 17 & 16 we find that the percolating $`\mathrm{\Lambda }`$CDM supercluster is a slim but massive object.101010Surface visualization is a difficult task especially if we wish to follow tunnels through superclusters to check whether our visual impression of the genus agrees with that calculated using the Euler formula (2). MATLAB has been used for surface plotting and the ‘reality’ of tunnels is ascertained by rotating surfaces and by viewing the supercluster at various angles. It contains $`4.5\%`$ of the total mass in the universe yet occupies only $`0.6\%`$ of the total volume. In SCDM the percolating supercluster contains $`4.4\%`$ of the total mass and occupies $`1.2\%`$ of the total volume, i.e., a similar mass occupies two times greater volume than in the $`\mathrm{\Lambda }`$CDM. The $`10^3`$ $`\mathrm{\Lambda }`$CDM clusters with $`\delta \delta _{\mathrm{perc}}`$ occupy $`4.4\%`$ of the total volume and contain $`33\%`$ of the total mass in the universe. Gravitational clustering thus ensures that most of the mass in the $`\mathrm{\Lambda }`$CDM universe is distributed in coherent filamentary regions which occupy very little volume but contain much of the mass. (Of the 1334 $`\mathrm{\Lambda }`$CDM clusters and superclusters identified at the percolation threshold, the ten most voluminous are extremely filamentary and contain close to $`40\%`$ of the total cluster mass (see table LABEL:tab:top10s123). The remaining $`60\%`$ of the mass is distributed amongst other 1324 objects !)

The percolating supercluster is more filamentary and less planar in $`\mathrm{\Lambda }`$CDM \[$`(P,F)_{\mathrm{\Lambda }CDM}`$=(0.13,0.81)\] than in SCDM \[$`(P,F)_{SCDM}`$=(0.14,0.78)\]. The percolating supercluster of the $`\tau `$CDM is least filamentary and most planar with $`(P,F)_{\tau CDM}`$=(0.15,0.70). It is interesting that superclusters evaluated at $`\delta _{\mathrm{perc}}`$ are more ‘porous’ (i.e. have a larger genus) in SCDM than their counterparts in $`\mathrm{\Lambda }`$CDM. For instance the percolating supercluster in $`\mathrm{\Lambda }`$CDM has only 6 ‘holes’ which tunnel all the way through it, while 20 tunnels pass through the percolating supercluster in SCDM; see Table LABEL:tab:top10s123. (This number increases as the density threshold is lowered as demonstrated in Figure 18.) Another interesting feature of gravitational clustering appears to be that although the volume occupied by clusters at $`\delta _{\mathrm{cluster}\mathrm{max}}`$ does not vary much between models, the fractional mass contained in $`\mathrm{\Lambda }`$CDM clusters is almost twice that in either $`\tau `$CDM or SCDM; see Figure 16.

### 5.3 Cluster Morphology

Cosmic density fields contain a wealth of information. As demonstrated in Fig. 16, at density thresholds significantly lower than $`\delta _{\mathrm{perc}}`$ most clusters merge to form a single percolating supercluster. On the other hand, the slightly larger NCS density contrast, $`\delta _{\mathrm{cluster}\mathrm{max}}`$, provides an excellent threshold at which to study the morphology of individual objects since it is precisely at $`\delta _{\mathrm{cluster}\mathrm{max}}`$ that the cluster abundance peaks. We study the morphology and topology of large scale structure in a two-fold manner: (i) properties of all superclusters are analysed at one of the two thresholds: $`\delta _{\mathrm{perc}},\delta _{\mathrm{cluster}\mathrm{max}}`$; (ii) the morphology of the largest supercluster is extensively probed as a function of the density contrast.

Figure 18 shows the evolution of the morphology and topology of the largest cluster as we scan through different threshold levels of the $`\mathrm{\Lambda }`$CDM density field. The density contrast, volume fraction and genus for a few of these levels have been labelled and are given in the right corner of the figure. We note that at high density thresholds the largest cluster occupies a very small fraction of the total volume and its shape is characterized by significant filamentarity ($``$0.67) and negligible planarity ($``$ 0). We note a sharp increase in the filamentarity of the largest cluster as we approach the percolation threshold (shown in the figure as a solid triangle). At smaller values of the density contrast, $`\delta \text{ }\stackrel{<}{}3`$, supercluster filamentarity rapidly drops while its planarity considerably increases. The drop in filamentarity of the supercluster is accompanied by a growth in its complexity with the result that, at moderately low thresholds ($`\delta 0.5`$), the percolating supercluster can contain several hundred tunnels. At very low thresholds ($`\delta \text{ }\stackrel{<}{}0.5`$) the percolating supercluster is an isotropic object possessing negligible values of both planarity (P) and filamentarity (F).

Since $`\delta _{\mathrm{perc}}`$ and $`\delta _{\mathrm{cluster}\mathrm{max}}`$ contain information pertaining to morphology and connectivity, we compile partial MFs for individual clusters at both these thresholds for the three cosmological models $`\mathrm{\Lambda }`$CDM, SCDM and $`\tau `$CDM. We find it convenient to work with the cumulative probability function (CPF) which we define as the normalized count of clusters with the value of a quantity Q to be greater than q at a given value ’q’. The quantity Q could be one of the MFs or one of the Shapefinders. We shall study the dependence of CPF with Q on a log-log scale.

Our results are presented in Figures 19, and 20 for clusters compiled at the NCS threshold $`\delta _{\mathrm{cluster}\mathrm{max}}`$. For large values of the Minkowski functionals ($`V,A,C`$) the CPF declines rapidly. Although the curves may look similar, some of them are statistically very different from their peers. The first four columns of Table LABEL:tab:ks show the results of the Kolmogorov-Smirnov test applied to the distribution of the MFs. (We do not carry out a similar exercise for the genus since the vast majority of clusters at the NCS threshold are simply connected.) In particular the CPF of the masses clearly distinguish the $`\mathrm{\Lambda }`$CDM model from SCDM & $`\tau `$CDM, in agreement with Fig. 19.

Figure 20 shows the Cumulative Probability Function for planarity and filamentarity in our cluster sample. The last two columns of Table LABEL:tab:ks show the results of the KS test for these statistics. For all three pairs they are among three best discriminators of the models. The value of KS statistic $`d`$ suggests that the signal comes from relatively small values of P ($`<0.1`$) and F ($`<0.25`$). The more conspicuous differences in the tail of the distribution is not statistically significant due to poor statistics. Clusters in $`\mathrm{\Lambda }`$CDM are the least anisotropic, a fact which corresponds to their relatively early formation and therefore longer evolution. From Figure 20 we also find that clusters at the NCS threshold are significantly more filamentary than they are planar. This appears to be a generic prediction of gravitational clustering as demonstrated by (1982); Klypin & Shandarin (1983); Sathyaprakash, Sahni & Shandarin (1996) and Bond, Kofman & Pogosyan (1996).

Figure 21 is a scatter plot of Shapefinders $`T,B,L`$ for clusters in $`\mathrm{\Lambda }`$CDM defined at the percolation threshold. The strong correlation between $`T`$ and $`B`$ in the left panel indicates that two (of three) dimensions defining any given cluster assume similar values and are of the same order as the correlation length. (Note that $`TB5`$h-1 Mpc for the largest superclusters in Table LABEL:tab:top10s123.) The clustering of objects near $`TB`$ ($`P0`$) in this panel suggests that our clusters/superclusters are either quasi-spherical or filamentary. The scatter plot between $`B`$ & $`L`$ in the right panel of Figure 21 breaks the degeneracy between spheres and filaments. The mass-dependence of morphology is highlighted in this panel in which larger dots denote more massive objects. This figure clearly reveals that more massive clusters/superclusters are, as a rule, also more filamentary, while smaller, less massive objects, are more nearly spherical. The concentration of points near the ‘edge’ of the scatter plot in the right panel, arises due to an abundance of low mass compact quasi-spherical objects with $`LB`$ ($`F0`$). The large number of massive superclusters with $`F>0.5`$ serves to highlight the important fact that the larger and more massive elements of the supercluster chain consist of highly elongated filaments as much as $`100`$h-1 Mpc in length, with mean diameter $`5`$h-1 Mpc (see Table LABEL:tab:top10s123). It is important to reiterate that almost 40% of the total overdense mass resides in the ten largest objects listed in table LABEL:tab:top10s123, while the remaining 60% is distributed among 1324 clusters.

We also find short filaments to be generally thinner than longer ones in agreement with predictions made by the adhesion model with regard to the formation of hierarchical filamentary structure during gravitational clustering (, 1992).

To further probe the morphology of clusters and superclusters we define the notion of shape-space in Figure 22. Shape-space is two dimensional with the planarity (of a cluster) plotted along the x-axis and its filamentarity along the y-axis. (One can also incorporate a third dimension showing the genus.) The first panel in Figure 22 is a scatter plot of $`P`$ and $`F`$ for clusters in $`\mathrm{\Lambda }`$CDM. The sizes of dots in the middle panel are proportional to cluster mass. We note that the most massive structures are also very filamentary. In the right panel we try to relate the shape of the structures with their topology by scaling the size of the dots with the genus value of clusters having a given morphology ($`P,F`$). As shown here, clusters which are multiply connected (larger dots are indicative of more complicated topologies) are also more filamentary. Together, the three panels show us that more massive superclusters are frequently very filamentary and often also topologically quite complex. We also see that a large number of less massive superclusters are simply connected and prolate. These structures are a few Mpc across along their two shorter dimensions and $`20`$ Mpc along the third, and therefore have appreciable filamentarity ($`F0.3`$).

## 6 Discussion and Conclusions

This paper presents a new technique for studying the geometrical and topological properties of large scale structure using the Minkowski functionals (MFs). Given a density field reconstructed by appropriately smoothing a point data set consisting (say) of a distribution of galaxies, our ansatz: (i) constructs closed polyhedral surfaces of constant density corresponding to excursion sets of a density field, using for this purpose a surface generating triangulation scheme (SURFGEN); (ii) SURFGEN evaluates the Minkowski functionals (volume, surface area, extrinsic curvature and genus) thus providing full morphological and topological information of three dimensional iso-density contour surfaces corresponding to a given data set. Evaluated in this manner, the Minkowski functionals can be used to study the properties of individual objects lying either above (clusters, superclusters) or below (voids) a given density threshold. They can also be used to study the morphological properties of the full supercluster-void network at (say) the percolation threshold. The ratio’s of Minkowski functionals (Shapefinders) are used to probe the shape of iso-density surfaces which sample the distribution of large scale structure at different thresholds of the density. (The highest density thresholds correspond to galaxies and clusters of galaxies, moderate thresholds correspond to superclusters while the lowest density thresholds characterize voids.) The performance of our ansatz has been tested against both simply and multiply connected surfaces such as the sphere, triaxial ellipsoid and triaxial torus. These three eikonal bodies can be smoothly deformed to give surfaces which are spheroidal, pancake-like (oblate) and filament-like (prolate) etc. Analytically known values of the Minkowski functionals for these surfaces allow us to test both our surface modelling scheme and our evaluation of Minkowski functionals from triangulation. Remarkably, we find that volume, area and integrated extrinsic curvature are determined to an accuracy of better than $`1\%`$ on all length-scales for both simply and multiply connected surfaces. The integrated intrinsic curvature (genus) is evaluated exactly. We also find that SURFGEN is remarkably accurate at calculating the MFs for Gaussian random fields (see the appendix). Having validated the performance of the code against eikonal surfaces and Gaussian random fields, we apply our method to cosmological simulations of large scale structure performed by the Virgo consortium. We study the geometry and topology of large scale structure in three cosmological scenarios - $`\mathrm{\Lambda }`$CDM, $`\tau `$CDM and SCDM. All three cosmologies are analysed at the present epoch (z=0) using global MFs, partial MFs and Shapefinders. (The redshift evolution of geometry and topology will be discussed in a companion paper.) Our main conclusions are summarized below:

* Using the Minkowski functionals we show that, like other diagnostics of clustering, supercluster morphology too is sensitive to the underlying cosmological parameter set characterizing our universe. Although the three cosmological models considered by us, $`\mathrm{\Lambda }`$CDM, $`\tau `$CDM and SCDM display features which are qualitatively similar, the geometrical and morphological properties of these models are sufficiently distinct to permit differentiation using MFs. We demonstrate that studying Minkowski Functionals using mass-parameterization significantly enhances their discriminatory power.
* SURFGEN employs a new contour-based method for determining topology. This method evaluates the genus for an iso-density contour surface directly, by applying the Euler formula Eq.(2) to the triangulated surface, and presents a significant improvement over earlier grid-based methods of determining topology (Melott, 1990).
* A study of percolation and cluster abundance reveals several interesting aspects of the gravitational clustering process. For all models, on decreasing the density threshold to progressively smaller values one reaches the critical percolation threshold at which the largest supercluster runs through (percolates) the simulation box. Percolation is reached at moderate values of the density contrast ranging from $`\delta _{\mathrm{perc}}2.3`$ for $`\mathrm{\Lambda }`$CDM to $`\delta _{\mathrm{perc}}1.2`$ for SCDM. The abundance of clusters reaches a maximum value at (or very near) the percolation threshold and the percolating supercluster occupies a rather small amount of space in all three cosmological models. Thus the fraction of total simulation-box volume contained in the percolating supercluster is least in $`\mathrm{\Lambda }`$CDM ($`0.6\%`$) and greatest in SCDM ($`1.2\%`$). When taken together, all overdense objects at the percolation threshold occupy $`4.4\%`$ of the total volume in the $`\mathrm{\Lambda }`$CDM model. For comparison, the volume fraction in overdense regions at the percolation threshold is $`16\%`$ in an idealised, continuous Gaussian random field (Shandarin & Zeldovich, 1989). This fraction can increase upto $`30\%`$ for Gaussian fields generated on a grid (Yess & Shandarin, 1996; Sahni, Sathyaprakash & Shandarin, 1997). The fact that clusters & supercluster occupy a very small fraction of the total volume appears to be a hallmark of the gravitational clustering process which succeeds in placing a large amount of mass ($`30\%`$ of the total, in the case of $`\mathrm{\Lambda }`$CDM) in a small region of space ($`4\%`$). The low filling fraction of the percolating supercluster in $`\mathrm{\Lambda }`$CDM ($`0.006`$) strongly suggests that this object is either planar or filamentary (Sahni, Sathyaprakash & Shandarin, 1997) and a definitive answer to this issue is provided by the Shapefinder statistic.
* Shapefinders were introduced to quantify the visual impression one has of the supercluster-void network of being a cosmic web of filaments interspersed with large voids (Sahni, Sathyaprakash & Shandarin, 1998). By applying the Shapefinder statistics via SURFGEN to realistic N-body simulations we have demonstrated that: (i) Most of the mass in the universe is contained in large superclusters which are also extremely filamentary; the vast abundance of smaller clusters and superclusters tends to be prolate or quasi-spherical. (ii) Of the three cosmological models the percolating supercluster in $`\mathrm{\Lambda }`$CDM is the most filamentary ($`F0.81`$) and the supercluster in $`\tau `$CDM the least ($`F0.7`$). (iii) The percolating supercluster in $`\mathrm{\Lambda }`$CDM is topologically a much simpler object than its counterpart in $`\tau `$CDM, with the former having only $`6`$ tunnels compared to $`19`$ in the latter. Other differences between the models are quantified in Figures 13 $``$ 22. We also show that among various morphological parameters, the Planarity and Filamentarity of clusters and superclusters is one of the most powerful statistics to discriminate between the models (see Table LABEL:tab:ks).

To summarize, this paper has demonstrated that Minkowski functionals and Shapefinders evaluated using SURFGEN provide sensitive probes of the geometry, topology and shape of large scale structure and help in distinguishing between rival cosmological models. Having established the strength and versatility of SURFGEN it would clearly be interesting to apply it to other important issues, related to those addressed in this paper, including: (i) geometry and morphology of underdense regions (voids); (ii) geometry and morphology of strongly overdense regions (clusters); (iii) redshift-space distortions of supercluster/void morphology; (iv) the time-evolution of the supercluster-void network; (v) analysis of superclusters and voids in fully three dimensional surveys such as 2dFGRS and SDSS; (vi) SURFGEN is also likely to provide useful insights in other physical and astrophysical situations in which matter is distributed anisotropically such as in the interstellar medium (Lazarian & Pogosyan (1996); Boumis, et. al. (2002)). We hope to return to some of these subjects in the near future.

## 7 Acknowledgements

We acknowledge useful discussions with Andrew Benson, Shaun Cole, Suman Datta, Carlos Frenk, Adrian Jenkins, Madan Rao, C.D.Ravikumar, Tarun Deep Saini, Yuri Shtanov, Volker Springel and David Weinberg. Our insights into programming techniques of computational geometry were honed by the excellent monograph on the subject by Joseph O’Rourke. Part of this work was done during a visit of JVS, VS and SS to the Department of Physics $`\&`$ Astronomy, Cardiff University, Cardiff under a PPARC Grant Ref. PPA/V/S/1998/00041. JVS is supported by the research fellowship of the Council of Scientific and Industrial Research (CSIR), India. VS and SS acknowledge support from the Indo-US collaborative project DST/NSF/US (NSF-RP087)/2001. The simulations studied in this paper were carried out by the Virgo Supercomputing Consortium using computers based at Computing Centre of the Max-Planck Society in Garching and at the Edinburgh Parallel Computing Centre. The data are publicly available at www.mpa-garching.mpg.de/NumCos.

## Appendix A Gaussian Random Fields and the Role of Boundary Conditions

This appendix demonstrates the great accuracy with which SURFGEN determines Minkowski functionals for Gaussian random fields (hereafter GRFs).

Before embarking on our discussion, we briefly describe the samples that we use and summarise the analytical results against which the performance of SURFGEN is to be tested.

We work with three realizations of a gaussian random field with a power law power spectrum (n = $``$1) on a $`128^3`$ grid. Each realization of the field is smoothed with a gaussian kernel of length $`\lambda =`$ 2.5 grid-units. All the fields are normalised by the standard deviation. This leaves them with zero mean and unit variance. The MFs are evaluated at a set of equispaced levels of the density field which coincide with the parameter $`\nu `$ on account of $`\sigma `$ being unity ($`\rho _{\mathrm{th}}=\nu \sigma =\nu `$). $`\nu `$ is used to label the levels and is related to the volume filling fraction through the following equation:

$$FF_V(\nu )=\frac{1}{\sqrt{2\pi }}_\nu ^{\mathrm{}}e^{\frac{t^2}{2}}𝑑t.$$
(20)
The samples are finally used for ensemble averaging.

The Minkowski Functionals (MFs) of a GRF are fully specified in terms of a length-scale $`\lambda _c`$,

$$\lambda _c=\sqrt{\frac{2\pi \xi (0)}{|\xi ^{^{\prime \prime }}(0)|}};\sigma ^2=\xi (0).$$
(21)
$`\lambda _c`$ can be analytically derived from a knowledge of the power spectrum. It can also be estimated numerically by evaluating the variance $`\xi `$(0) of the field and the variance $`\xi ^{^{\prime \prime }}`$(0) of any of its first spatial derivatives (for more details see Matsubara 2003).

For a GRF in three dimensions the four MFs (per unit volume) are (Tomita, H., 1990; Matsubara, 2003):

$$V(\nu )=\frac{1}{2}\frac{1}{2}\mathrm{\Phi }(\frac{\nu }{\sqrt{2}});\mathrm{\Phi }(x)=\frac{2}{\sqrt{\pi }}_0^x^{}\mathrm{exp}\left(\frac{x^2}{2}\right)𝑑x^{},$$
(22)
where $`\mathrm{\Phi }`$(x) is the error function.

$$S(\nu )=\frac{2}{\lambda _c}\sqrt{\frac{2}{\pi }}\mathrm{exp}\left(\frac{\nu ^2}{2}\right),$$
(23)
$$C(\nu )=\frac{\sqrt{2\pi }}{\lambda _c^2}\nu \mathrm{exp}\left(\frac{\nu ^2}{2}\right),$$
(24)
$$G(\nu )=\frac{1}{\lambda _c^3\sqrt{2\pi }}(1\nu ^2)\mathrm{exp}\left(\frac{\nu ^2}{2}\right).$$
(25)
Let us now describe the boundary conditions which we adopt and which incorporate the periodic nature of the GRFs generated using a FFT routine. GRFs defined on a grid are, by construction, periodic in nature. As a result, two overdense excursion sets visibly separated from each other inside the box and touching two opposite sides of the box would constitute a single cluster. To make this more clear, let us turn to Fig. 23 which shows four clusters (labelled 1, 2, 3 and 4 respectively) in the left panel. The boundaries of these clusters $``$ as drawn by an analog of SURFGEN in 2-D $``$ are shown in bold, and we note that the sides of the box (bold dashed lines) are incorporated when we define the clusters. However, due to the periodicity in the density field, these four clusters are in fact one single cluster (right panel). We note that the portion of the contours due to the boundaries of the box occurs in the interior of the actual cluster (shown with a bold dashed line).

It is easy to extrapolate this situation to 3$``$D. Note that the interpolation method employed by SURFGEN will correctly estimate the global volume of the excursion set even if we worked with the field as it is. In 2$``$D this is same as estimating the area enclosed by the full cluster without having to piece the four “clusters” (labelled by 1, 2, 3 and 4) together. Our estimation of global genus will also turn out to be exact, and as it can be shown, this will hold true even if either of these four constituent clusters exhibits a nontrivial topology. On the other hand, if we work with the same set of contours for estimating global area and mean curvature, we will wrongly be adding an excess contribution due to the sides of the box. (In 2$``$D, this translates to adding the contribution of the dashed line segments to the perimeter of the contour in the right panel of Figure 23). The actual area and mean curvature contribution comes from the interface between overdense and underdense regions (such as the 2$``$D contours drawn in solid bold line in both the panels) as opposed to the artificial boundaries of the clusters (shown dashed).

We note that incorporating periodic boundary conditions in order to construct a total contour as shown in the right panel of Fig. 23 is nontrivial from the numerical point of view. When the excursion set touches all the boundaries (as, for example, at a very low thresold of density, below the percolation threshold), the corresponding algorithm could further run into an infinite loop, which makes it impractical for implementation. In the light of this, it would be more manegeable if we restricted ourselves to using contours which are closed at the boundaries of the box (like the ones in the left panel of Fig. 23), but devise a way to exclude the excess contribution to surface area and mean curvature from the boundaries of the box. A uniform prescription of this type, when applied to all the appropriate clusters at all the thresholds of density will lead to correct estimation of all the global MFs. This then, is the choice of boundary conditions incorporated in our calculation of a GRF. A futher cautionary note concerns the familiar W-shape of the genus-curve. Recall that the genus-curve peaks at $`\nu =0`$ signifying a sponge-like topology of the medium at the mean density threshold. Note further that G($`\nu =\pm \sqrt{3}`$)= $``$2G(0)$`\mathrm{exp}\left(\frac{3}{2}\right)`$ are the two negative minima of the genus curve. In fact for $`|\nu |>`$ 1, the genus is always negative, and approaches zero from below when $`|\nu |1`$. For $`\nu <1`$, this negative genus refers to isolated underdense regions (bubbles) in the overdense excursion set. Similarly, the negative genus for $`\nu >1`$ should be interpreted as being caused by isolated overdense regions (meatballs) in the underdense excursion set. To summarize, in order to determine the genus for a GRF we need to estimate the genus of the overdense excursion set for $`\nu `$0 and the genus of the underdense excursion set for $`\nu >`$ 0, thereby making use of the symmetry property of the GRFs. (Note however that in our earlier determination of the genus curve for N-body simulations in this paper, we have dealt with the topology of overdense regions only, because of which the genus curve is asymmetric with respect to $`\nu =0`$ and is always positive when $`\nu >1`$.)

Having elaborated on the nature of boundary conditions used in our evaluation of MF’s for a GRF, we now present our results.

Figure 24 shows the global MFs averaged over 3 realizations of the density field (solid lines with 1$`\sigma `$ error bars) along with the exact analytical results (shown dotted). There are in all 40 levels for every realization. Figure 24 clearly shows the remarkable agreement between exact theoretical results (22) - (25) and numerical estimates obtained using SURFGEN. We therefore conclude that SURFGEN determines the MF’s of a GRF to great precision.

The above discussion demonstrated the important role played by boundary conditions in reproducing the analytical predictions for GRFs. Having tested the performance of SURFGEN against GRFs, we are faced with the following choice of boundary conditions when dealing with N$``$body simulations and mock/real galaxy catalogues: (i) we could either correct for the boundaries when dealing with the overdense regions which encounter the faces of the survey-volume, or, (ii) we could avoid doing this and treat all measurements with a tacit understanding that area and mean curvature contain an excess contribution which arises because of boundary effects. Of the two options, the former is computationally more demanding. In addition, its effect is appreciable only significantly after the system has percolated, i.e., when there are a large number of clusters touching the boundaries of the box. The corresponding value of the filling factor (FFV $`>`$ FFperc) is useful only when one wishes to compare two samples through the trends in their global MFs. Since the structural elements of the cosmic web (superclusters, voids) are usually identified near the percolation threshold, and are therefore not very sensitive to boundary effects, we have chosen option (ii) in our analysis of N-body simulations in this paper. Option (ii) is also more suited for dealing with real galaxy catalogues which do not satisfy periodic boundary conditions. We should also point out that since we compare two samples under the same conditions $``$ the same volume and dimensions of the box, same resolution, etc.$``$ the excess contribution to the area and mean curvature caused by box boundaries enters in identical fashion for both the samples. We therefore prefer option (ii) to (i) and all our simulations in Section 5 are analysed without correcting for contributions from the boundaries. (It may also be noted that the contribution from the boundaries becomes less important as we deal with larger survey volumes.)


text:
# Deformations of plane algebraic curves and integrable systems of hydrodynamic type

## 1 Deformation of plane algebraic curves

Let $`𝒞`$ be a plane algebraic curve on $`^2=\{(k,p)\}`$:

$$𝒞:p^N=\underset{j=1}{\overset{N}{}}u_j(k)p^{Nj},u(k)[k].$$
(1)
Then consider a deformation of the curve $`𝒞`$ as the following equation of conservation law;

$$\frac{p}{t}=\frac{}{x}Q,\mathrm{with}Q\frac{[k,p]}{𝒞}=\underset{k=1}{\overset{N}{}}[k]p^{Nk},$$
(2)
where $`x,t`$ are the deformation parameters. The equation (2) describes a (quasi-linear) system of equations of hydrodynamic type for $`u_j=u_j(x,t;k)`$ in the curve $`𝒞`$. Then from a direct computation we find:

###### Lemma 1.1.

Let $`Q`$ be given by

$$Q=\underset{j=1}{\overset{N}{}}\alpha _j(k)p^{Nj},\alpha _j(k)[k].$$

Then the functions $`u_j`$ in the curve $`𝒞`$ satisfy

$$\frac{u_l}{t}=\underset{j=1}{\overset{N}{}}N\gamma _l^{(j)}\frac{\alpha _j}{x}+\underset{i=1}{\overset{N}{}}\underset{j=1}{\overset{N}{}}(Nj)\gamma _l^{(i+j)}\left(\alpha _j\frac{u_i}{x}u_j\frac{\alpha _i}{x}\right),$$

where $`\gamma _l^{(j)}`$ are defined by

$$p^{2Nj1}=\underset{l=1}{\overset{N}{}}\gamma _l^{(j)}p^{Nk},1j2N1.$$

The last equation leads to the recursion relation for $`\gamma _l^{(i)}`$ for $`l=1,\mathrm{},N`$ with $`\gamma _{N+1}^{(i)}=0`$,

$$\gamma _l^{(i1)}=\{\begin{array}{cccc}u_l\gamma _1^{(i)}+\gamma _{l+1}^{(i)},\hfill & \mathrm{for}\hfill & 1iN1\hfill & \\ \delta _{l,iN+1},\hfill & \mathrm{for}\hfill & Ni2N1\hfill & \end{array}$$

which can be used to find the explicit form of $`\gamma _l^{(i)}`$’s in terms of $`u_j`$’s.

A hierarchy associated with the system (2) may be obtained by choosing proper $`u_j(k),\alpha _j(k)[k]`$ for $`j=1,\mathrm{},N`$, so that the systems for $`n=1,2,\mathrm{}`$,

$$\frac{p}{t_n}=\frac{}{x}Q^{(n)},\mathrm{with}Q^{(n)}=\underset{j=1}{\overset{N}{}}\alpha _j^{(n)}(k)p^{Nj},$$

are compatible, i.e. $`^2p/t_nt_m=^2p/t_mt_n`$. (In this paper we choose $`Q^{(0)}=p`$, that is, $`p/t_0=Q^{(0)}/x`$, and identify $`t_0=x`$.) A general scheme to construct such $`Q^{(n)}`$ has not been found. It is however quite interesting to classify the curves which allow integrable deformations of this type. In this report, we give several examples of the hierarchy associated with the deformations of the curve $`𝒞`$ in (1).

###### Remark 1.2.

The curve (1) contains several known examples of the hierarchy, which are obtained by the reductions of the dispersionless KP (and mKP) hierarchy (see also Remark 3.2). The examples contain a model of the Landau-Ginzburg theory with rational potentials . This is simply obtained by the following choices of polynomials for $`u_i`$’s in (1): Fix an index of $`u_i`$’s, say $`n`$, and set

$$u_n=k^nv_n,\mathrm{and}u_i=v_i,\mathrm{for}in,$$
(3)
where $`v_i`$’s do not depend on $`k`$. Then we have

$$k^n=p^n+v_1p^{n1}+\mathrm{}+v_n+\frac{v_{n+1}}{p}+\mathrm{}+\frac{v_N}{p^{Nn}}.$$
(4)
which gives a rational Landau-Ginzburg potential. An interesting question is how one can extend the ring $`[k,p]/𝒞`$ to a rational ring. Also one should include a log-term to define a complete set of primary fields.

There is another interesting example, which is obtained by the dipersionless limit of the vector nonlinear Schrödinger equation (Zhakharov reduction),

$$k=p+\underset{j=1}{\overset{N1}{}}\frac{\rho _j}{pv_j}.$$

The hierarchy associated to this curve provides a model of the topological field theory with $`N1`$ punctures.

###### Remark 1.3.

The system (2) has an infinite number of conserved densities which are explicitely obtained by expressing $`p`$ in the Laurent series of $`k`$.

## 2 Deformations of hyperelliptic curves $`(N=2)`$

The curve $`𝒞`$ with $`N=2`$ corresponds to hyperelliptic curve of genus determined by the degrees of polynomials $`u_j(k),j=1,2`$;

$$p^2=u_1(k)p+u_2(k),u_1,u_2[k].$$

The deformation equation (2) then gives

$$\{\begin{array}{cccc}& \frac{u_1}{t}=\frac{}{x}\left(u_1\alpha _1+2\alpha _2\right)\hfill & & \\ & & & \\ & \frac{u_2}{x}=2u_2\frac{\alpha _1}{x}+\alpha _1\frac{u_2}{t}u_1\frac{\alpha _2}{x}.\hfill & & \end{array}$$
(5)
As a particular case with $`u_1=kv_1,u_2=v_2`$ (see Remark 1.2), we have

$$k=p+v_0+\frac{v_1}{p},$$

which corresponds to the curve of the dispersioness Toda equation discussed in . This curve is also related to the dispersionless AKNS system for the variable $`p^{}=p+v_0`$.

We also note that the variable $`u_1(k)`$ can be eliminated by the (gauge) transformation, $`pp+u_1/2`$, and then from (5) one can take $`Q`$ in the form,

$$Q=\alpha _1(k)p,(\mathrm{i}.\mathrm{e}.\alpha _2=0).$$

Thus we have the following two examples for the case $`N=2`$, where the polynomial $`u_2(k)`$ is either odd or even.

### 2.1 The $`\mathrm{BH}_m`$ hierarchy

The $`\mathrm{BH}_m`$ hierarchy is defined on the singular sectors of the Burgers-Hopf (BH) hiearchy, where the solution of the hiearachy has a shock type singularity. Here the number $`m`$ is the genus of the curve $`𝒞`$ given by

$$y^2=u_2(k)=k^{2m+1}+\underset{j=0}{\overset{2m}{}}v_jk^{2mj}.$$
(6)
In particular, the case with $`m=0`$ gives the BH hierarchy (or the dispersionless KdV hierarchy), i.e. with $`u_2=k+v_0`$,

$$p^2=k+v,\mathrm{and}Q^{(n)}=\alpha _1^{(n)}(k)p=[k^{n+\frac{1}{2}}]_{+p},$$

where $`[k^a]_{+p}`$ represents the polynomial part of $`k^a`$ in $`p`$. The BH hierarchy in terms of the variable $`v`$ is given by

$$\frac{v}{t_n}=c_nv^n\frac{v}{x},\mathrm{with}c_n=(1)^n\frac{(2n+1)!!}{2^nn!},n=0,1,2,\mathrm{}.$$

In the case with $`m0`$, we have:

###### Proposition 2.1.

The following system associated with the curve (6) forms an integrable hierachy of a deformation of the curve;

$$\frac{p}{t_n}=\frac{}{x}Q^{(n)},\mathrm{with}Q^{(n)}=\alpha _1^{(n)}p=\left[\frac{k^{m+n+\frac{1}{2}}}{p}\right]_+p,$$

where $`[]_+`$ is the polynomial part of $`k`$.

The form $`Q^{(n)}`$ can be obtained as follows: First assume $`\alpha _1^{(n)}(k)`$ be a monic polynomial of degree $`n`$ in $`k`$. Then from (2) we have

$$\alpha _1^{(n)}(k)p=k^{n+m+\frac{1}{2}}+O\left(k^{m\frac{1}{2}}\right),$$

from which we have

$$\alpha _1^{(n)}=\left[\frac{k^{n+m+\frac{1}{2}}}{p}\right]_+.$$

The compatibility between the flows in the hierarchy can be shown directly. The explicit form of the $`\mathrm{BH}_m`$ hierarchy is given by the following form with the polynomial $`\alpha _1^{(n)}(k)`$ for $`k`$ being replaced by a matrix $`K`$,

$$\frac{U}{t_n}=\alpha _1^{(n)}(K)\frac{U}{x},\mathrm{with}U=(u_0,u_1,\mathrm{},u_{2m})^T,$$

where $`K`$ is the (companion) matrix given by

$$K=\left(\begin{array}{cccccc}u_0\hfill & 1\hfill & \mathrm{}\hfill & \mathrm{}\hfill & 0\hfill & \\ u_1\hfill & 0\hfill & \mathrm{}\hfill & \mathrm{}\hfill & 0\hfill & \\ \mathrm{}\hfill & \mathrm{}\hfill & \mathrm{}\hfill & \mathrm{}\hfill & \mathrm{}\hfill & \\ u_{2m1}\hfill & 0\hfill & \mathrm{}\hfill & \mathrm{}\hfill & 1\hfill & \\ u_{2m}\hfill & 0\hfill & \mathrm{}\hfill & \mathrm{}\hfill & 0\hfill & \end{array}\right).$$

For example, $`\alpha _1^{(1)}(K)=K(1/2)u_0I`$ with $`I=(2m+1)\times (2m+1)`$ identity matrix.

It is also interesting to note that the systems in the hierarchy can be written in the Riemann invariant form,

$$\frac{\kappa _i}{t_n}=\alpha _1^{(n)}\frac{\kappa _i}{x},\mathrm{for}i=0,1,2,\mathrm{},2m,$$

where the Riemann invariants $`\kappa _i`$’s are the roots of the polynomial associated with the curve $`𝒞`$, i.e.

$$p^2=\underset{j=1}{\overset{N}{}}(k\kappa _j).$$

The $`u_j`$’s are then given by the elementary symmetric polynomials of $`\kappa _j`$’s, i.e.

$$u_j=(1)^{j+1}\underset{1i_1<\mathrm{}<i_jN}{}\kappa _{i_1}\mathrm{}\kappa _{i_j},\mathrm{for}j=0,\mathrm{},2m.$$

The $`\kappa _j`$ are also the eigenvalues of the matrix $`K`$. Then the solution can be obatined by the generalized hodograph form,

$$\mathrm{\Omega }_j(v_1,\mathrm{},v_{2m+1},x,t_1,\mathrm{})=\underset{n=0}{\overset{\mathrm{}}{}}\alpha _1^{(n)}(\kappa _j)t_n=0,j=1,2,\mathrm{},2m+1.$$

The regularity of the solution is obtained by the smoothness of the curve, which corresponds to the distinction of the roots $`\kappa _j`$. The singular structure of the solution can be described by their intersections of the functions $`\mathrm{\Omega }_j`$ in the hodograph solution.

### 2.2 The dispersionless $`\mathrm{JM}_m`$ hierarchy.

The dispersionless $`\mathrm{JM}_m`$ ($`\mathrm{dJM}_m`$) hierarchy is given by a dispersionless (classical) limit of the hidden integrable hierachy of the Jaulent-Miodic equation. The $`\mathrm{dJM}_m`$ hierarchy is given by the similar form as the $`\mathrm{BH}_m`$ hierarchy:

###### Proposition 2.2.

For the hyperelliptic curve of genus $`m`$ given by

$$p^2=k^{2m+2}+\underset{j=1}{\overset{2m+2}{}}v_jk^{2m+2j},$$

we have an integrable hierarchy of deformation of the curve,

$$\frac{p}{t_n}=\frac{}{x}Q^{(n)},\mathrm{with}Q^{(n)}=\alpha _1^{(n)}p=\left[\frac{k^{m+n+1}}{p}\right]_+p.$$

For example in the case with $`m=0,n=1`$, we have

$$\{\begin{array}{cccc}& \frac{v_1}{t_2}=\frac{}{x}\left(v_2\frac{3}{4}v_1^2\right)\hfill & & \\ & & & \\ & \frac{v_2}{t_2}=\frac{1}{2}\frac{}{x}\left(v_1v_2\right)\hfill & & \end{array}$$

which describes classical shallow water waves. Then the $`\mathrm{dJM}_m`$ hierarchy can be considered as an integrable system defined on the singular sector of the codimension $`m`$ in the solution of the $`\mathrm{dJM}_0`$ hierarchy. As in the case of the BH hierarchy, the singular sectors are given by the intersection structure of the hodograph solution $`\mathrm{\Omega }_j,j=1,2`$ on $`^{\mathrm{}}=\{(x,t_1,\mathrm{},)\}`$.

## 3 Examples of non-hyperelliptic case

A natural extension of the hyperelliptic case ($`N=2`$) may be given by the case with the curve,

$$p^N=u_N(k)=k^M+\underset{j=1}{\overset{M}{}}v_jk^{Mj},$$
(7)
where $`N,M`$ are positive integers with $`N3`$. The genus of the (irreducible projective) curve $`𝒞^2`$ for $`N>M`$ is given by Max Noether’s genus formula

$$g=\frac{1}{2}((N1)(N2)(N1)(NM1)\mathrm{g}.\mathrm{c}.\mathrm{d}.(N,NM)+1),$$

where g.c.d.($`A,B`$) implies the greatest common divisor of $`A`$ and $`B`$. If $`M>N`$, the genus is given by the same formula with the exchange $`NM`$. Also if $`M=N`$, the curve is smooth, and the genus is

$$g=\frac{(N1)(N2)}{2}.$$

For the curve (7), we have a similar result as in the case of hyperelliptic curve:

###### Proposition 3.1.

The following integrable hierarchy of hydrodynamic type provides a deformation of the curve (7);

$$\frac{p}{t_n}=\frac{}{x}Q^{(n)},\mathrm{with}Q^{(n)}=\left[k^{\frac{M}{N}+n}\right]_+p.$$

The compatibility can be shown in the same way as the hyperelliptic case and it can be found in our paper.

###### Remark 3.2.

One should note that the hierarchy in Proposition 3.1 is not the same as a dispersionless Lax reduction, i.e.

$$k^N=[k^N]_{+p}=\mathrm{polynomial}\mathrm{of}\mathrm{degree}N\mathrm{in}p.$$

However the Lax reduction is a special case of the curve (1) where all of the $`u_j(k)`$’s do not depend on $`k`$ except $`u_N(k)=k^N`$.

## Acknowledgments

YK would like to thank the organizers for a financial support for this meeting.