The Evolution of Dark Matter Halo Properties in Clusters, Filaments, Sheets and Voids
Abstract
We use a series of highresolution Nbody simulations of the concordance cosmology to investigate the redshift evolution since of the properties and alignment with the largescale structure of haloes in clusters, filaments, sheets and voids. We find that: Once a rescaling of the halo mass with , the typical mass scale collapsing at redshift , is performed, there is no further significant redshift dependence in the halo properties; The environment influences the halo shape and formation time at all investigated redshifts for haloes with masses ; and There is a significant alignment of both spin and shape of haloes with filaments and sheets. In detail, at all redshifts up to : a) Haloes with masses below tend to be more oblate when located in clusters than in the other environments; this trend is reversed at higher masses: above about , halos in clusters are typically more prolate than similar massive halos in sheets, filaments and voids. b) The haloes with in filaments spin more rapidly than similar mass haloes in clusters; haloes in voids have the lowest median spin parameters; c) Haloes with tend to be younger in voids and older in clusters; d) In sheets, halo spin vectors tend to lie preferentially within the sheet plane independent of halo mass; in filaments, instead, haloes with tend to spin parallel to the filament and higher mass haloes perpendicular to it. For halo masses , the major axis of haloes in filaments and sheets is strongly aligned with the host filament or the sheet plane, respectively. Such haloLSS alignments may be of importance in weak lensing analyses of cosmic shear. A question that is opened by our study is why, in the redshift regime that we have investigated, the mass scale for gravitational collapse, , sets roughly the threshold below which the largescale structure environment either begins to affect, or reverses, fundamental properties of dark matter haloes.
keywords:
cosmology: theory, dark matter, largescale structure of Universe – galaxies: haloes – methods: Nbody simulations1 Introduction
Numerical simulations of concordance cosmology have shown that properties of dark matter haloes do not depend only on the mass of the halo, as suggested by prior analytical work based on the excursion set theory for structure evolution (e.g. Bond et al., 1991; Lacey & Cole, 1993). Rather, they also depend on the environment in which the halo resides (e.g. Gao et al., 2005; Wechsler et al., 2006; Harker et al., 2006; Wang et al., 2007; Hahn et al., 2007, Paper I in the following). This dependence on environment is quite significant at redshift zero for low mass haloes, typically with masses . In detail, using marked statistics, Sheth & Tormen (2004) find evidence that haloes of a given mass form earlier in dense regions. High resolution simulations confirm this finding and quantify it as a function of halo mass (Gao et al., 2005; Croton et al., 2007; Harker et al., 2006; Reed et al., 2007; Maulbetsch et al., 2007). At the same time it has become clear that also other halo properties as concentration and spin correlate with local environment (AvilaReese et al., 2005; Wechsler et al., 2006; Bett et al., 2007; Macciò et al., 2006; Wetzel et al., 2007; Hahn et al., 2007). Gao & White (2007) find that haloes with, e.g., high spinparameter or formation time, tend to be more strongly clustered than younger and lowspin haloes. It is possible that this environmental dependence of halo properties has also an impact on the baryonic galaxies. Galaxy properties in the local Universe are known to vary systematically with environment (e.g. Dressler, 1980; Kauffmann et al., 2004; Blanton et al., 2005).
In this paper, we extend the investigation of the properties of dark matter haloes as a function of environment to high redshifts. In particular, we follow the definition of environment that we presented in Paper I, which associates haloes to four classes with different dynamical properties: voids, sheets, filaments and clusters. These four environments are identified on the basis of a tidal stability criterion for test particles which is inspired by the Zel’dovich approximation (Zel’dovich, 1970). We find that, at each redshift, all investigated properties of haloes show some correlation with mass and environment, and that the redshift dependence of halo properties with mass is removed when such properties are investigated as a function of the rescaled mass , where is the typical massscale collapsing at each epoch.
The paper is organised as follows. In Section 2, we briefly summarise the specifics of our Nbody simulations, the definitions of the halo properties that we study and the definitions of the four environments  clusters, sheets, filaments and voids. We present the results in Section 3 and summarise our conclusions in Section 4.
2 Numerical Simulations and Definitions
We use the three highresolution cosmological Nbody simulations described in more detail in Paper I, which were obtained with the treePM code Gadget2 (Springel, 2005). These simulations are used to follow the formation and evolution of largescale structure in a flat CDM cosmology with matter density parameter , baryonic contribution and a presentday value of the Hubble constant with with an initial power spectrum normalised to have . Each simulation consists of collisionless dark matter particles in periodic boxes of sizes , and , respectively. The corresponding particle masses are , and for the three boxes. Initial conditions were generated using the Grafic2 tool (Bertschinger, 2001). Particle positions and velocities were saved for 30 timesteps, logarithmically spaced in expansion parameter between and . The massrange of these three simulations allows us to resolve haloes with masses below up to redshifts .
The halo properties that we investigate are formation redshift, shape, and spin parameter. The formation redshift and shape parameters are defined as in Paper I; we adopt however a slightly different approach to measure the halo spin parameter than in our previous work. We summarise our definitions below.
2.1 Halo Catalogues
Haloes were identified in each snapshot using the standard friendsoffriends (FOF, Davis et al., 1985) algorithm with a linking length equal to times the mean interparticle distance. Haloes that are wellresolved in each of the three simulations are then combined into one single catalogue. Unrelaxed systems were identified and deleted from the halo catalogues. These unrelaxed systems are mainly contributed by closepair halo configurations which are spuriously linked into one single halo. To identify them, we follow Bett et al. (2007) and define the virialization parameter
(1) 
where is the total kinetic energy including the Hubble flow with respect to the centre of mass and is the total potential energy of the isolated FOF halo. The potential is computed using a tree for groups with more than particles, and via direct summation for smaller halos. The virial theorem states that the time average of vanishes for any isolated relaxed object. However, infalling material exerts a surface pressure such that (Hetznecker & Burkert, 2006). In addition, structures that are gravitationally bound have . In order to exclude accidentally linked unvirialised structures or haloes that are just about to merge, it suffices to fix bounds on . In order to directly compare with Bett et al. (2007), we make the same choice of that was adopted by those authors to set the threshold between virialized and nonvirialized structures.
Finally, we also exclude from our halo catalogues all those structures where the distance between the centre of mass and the most bound particle exceeds a fixed fraction of the largest distance between a particle in the halo and the centre of mass , i.e. .
The cleaning of the halo catalogues has a strong effect on the spin parameter distribution but only a minor influence on the other quantities that we study in this paper.
2.2 Formation Redshift
For each halo at redshift , we identify a progenitor at by identifying particles that are contained in both haloes. The main progenitor is then chosen to be the most massive halo at each redshift that contributes at least 50 per cent of its particles to the final halo. We then define the formation redshift to be the epoch at which a main progenitor which has at least half of the final mass first appears in the simulation; specifically, is found by linearly interpolating between simulation snapshots in to find the point where half of the given halo mass is accumulated.
2.3 Halo Shape
In order to determine the shape of haloes, we use the moment of inertia tensor
(2) 
where is the particle mass, is the distance of the th particle from the centre of mass of the halo and denotes the Kronecker symbol. Given the lengths of the principal axes of inertia , we then use the following definitions of sphericity and triaxiality :
(3) 
We find that a minimum of 500 particles per halo guarantees numerically reliable estimates of the shape parameters.
2.4 Halo Spin Parameter
We estimate the spin parameter (Peebles, 1969) of a halo using the simplified form (Bullock et al., 2001)
(4) 
Here all quantities with the subscript “vir” (angular momentum, mass and circular velocity) are computed within a sphere of radius around the most bound particle enclosing a mean density of , where is the critical density, and the density parameter according to the spherical collapse model. This density parameter can be approximated by (Bryan & Norman, 1998):
(5) 
with
(6) 
in a flat cosmology. Under the assumption that the halo is in dynamical equilibrium, , the spin parameter can be rewritten as
(7) 
Systematic numerical artefacts were found to be negligible for haloes consisting of at least 300 particles. The cleaning of the halo catalogues, as described in Section 2.1, has a strong influence on the distribution of . The probability for two haloes of similar mass to be erroneously linked by the halo finder grows both with the environmental density and decreasing mass of the haloes. This leads to an increasing component of unrelaxed structures of low mass in the spin distribution for which the virialisation conditions are not fulfilled. Furthermore, the angular momentum is dominated by the orbital angular momentum of the pair rather than the intrinsic spin of either one of them. Exclusion of unrelaxed objects removes the tail of these apparent highspin haloes with .
2.5 Environmental Classification
We employ the definition of environment that was described in detail in Paper I. This classification is based on a local stability criterion for the orbits of test particles in the matter distribution at a fixed epoch. A series expansion of the equation of motion for a test particle in the smoothed matter distribution gives a zero order acceleration and a first order deformation term that is determined by the tidal field tensor, i.e. the Hessian of the peculiar gravitational potential. We then classify the contraction/expansion behaviour of this first order term by the number of its positive/negative eigenvalues. In analogy with Zel’dovich theory (Zel’dovich, 1970), the application of this criterion provides four environmental categories defined by the signs of the three eigenvalues , namely:

clusters (),

filaments (),

sheets (), and

voids ().
This definition relies on one free parameter, the length scale of the Gaussian filter that is used to smooth the matter distribution before obtaining the gravitational potential on a grid. As done in Paper I, we fix the smoothing scale at redshift to be , as this value provides the best agreement between the outcome of the orbitstability criterion and our aposteriori visual classification of the different environments. The smoothing length scale is related to the mass contained in the Gaussian filter at mean density via ; thus, a smoothing scale corresponds to .
A natural mass scale at any given redshift is given by the typical mass scale for collapse defined as follows. A spherical tophat perturbation collapses when its linear overdensity exceeds a value of with only a weak dependence on cosmological parameters (e.g. Eke et al., 1996). The variance of linear density fluctuations at a given mass scale is related to the linear power spectrum at redshift by
(8) 
where is the Fourier transform of a spherical tophat window function of comoving size , and is the comoving mean mass density of the universe. At a given redshift, the typical mass scale to collapse from a fluctuation is hence given by the implicit solution of
(9) 
The mass scale at is .
When computing the environmental classification at redshifts , there are two possible approaches that can be followed: i) To keep the smoothing scale () fixed to the value: the environment is thus defined over typical separations of a few Mpc in comoving space; or ii) To vary the smoothing scale. In particular, a natural choice for a redshiftdependent smoothing scale is to maintain the ratio fixed for the Gaussian filter. The respective mass scales for the high simulation snapshots investigated in this paper, i.e., and , are and . Fixing the ratio maintains the smoothing on scales of order at all redshifts.
The resulting classifications for both and at redshifts , and are shown in Figure 1 using different colours for the cluster, sheet, filament and void environments. We observe some differences between the two smoothing approaches. With a fixed smoothing length , shown in panels b) and d), the regions classified as voids, sheets and filaments remain virtually unchanged since , while the cluster environments grow substantially in size, from to higher redshifts, and extend well into the filaments at . With the constant smoothing, shown in panels c) and e), a much larger number of individual haloes change environment with redshift: at the resolution of our simulations, very few haloes are detected at in void regions, while many haloes are associated at the same redshift to relatively small cluster environments. In Tables 1 (for fixed ) and 2 (for fixed ) we show the fraction of haloes at that change their environmental class from through to by following the main progenitors of each halo with a minimum mass of in the box. For a fixed smoothing mass scale, indeed less haloes change their environmental class as the density contrasts between the environments grow through gravitational collapse.
z=0.49  
z=0  void  sheet  filament  cluster 
void  1.0  0.094  0  0 
sheet  0.28  27  1.1  0 
filament  0  5.5  54  0.29 
cluster  0  0.088  6.4  4.2 
z=1.05  
z=0  void  sheet  filament  cluster 
void  0.95  0.18  0  0 
sheet  0.56  26  1.86  0 
filament  0.028  10.3  49  0.46 
cluster  0  0.53  7.5  2.7 
z=0.49  
z=0  void  sheet  filament  cluster 
void  0.34  0.71  0.077  0 
sheet  0.050  18  9.9  0.36 
filament  0  6.2  51  2.9 
cluster  0  0.31  8.0  2.4 
z=1.05  
z=0  void  sheet  filament  cluster 
void  0  0.22  0.79  0.12 
sheet  0  5.0  19  4.1 
filament  0  4.5  46  8.8 
cluster  0  0.72  8.0  2.0 
The fraction of volume attributed to each of the four environments as a function of redshift is shown for both smoothing approaches in Figure 2. For a fixed smoothing mass , the density field asymptotically approaches Gaussianity with increasing redshift and thus the expected volume fractions for the four environments (cf. Doroshkevich, 1970). The behaviour is very different with the constant smoothing. The volume occupied by the unstable structures (sheets and voids) increases with redshift, while the volume fractions of the stable structures (clusters and filaments) decrease compared to the values at .
The redshift evolution of the median value of the smoothed overdensity field as measured at the centres of haloes in the four environments is shown in Figure 3. For a constant smoothing scale, the median overdensities grow faster than expected in linear perturbation theory as in clusters and voids. With the smoothing, however, the median overdensities are found to be essentially constant at all redshifts. In both smoothing approaches, the median overdensity in filaments is constant. The median overdensities, smoothed on scales of at redshift zero, are , , and in voids, sheets, filaments and clusters, respectively.
We note that with both smoothing approaches, and most relevantly when adopting a constant ratio for the smoothing, the haloes in the immediate surroundings of the clusters at are classified as filaments/sheets at this redshift, but they make the transition to the cluster environment by . This allows us to rigorously identify, and thus study the properties of, the haloes in these intermediatedensity “transition regions”, before their ultimate migration into the deeper potential wells of rich clusters at . We plan to investigate these haloes in future work.
3 The Redshift Evolution of Halo Properties in Different Environments
3.1 Mass Functions
The choice of smoothing scale with redshift has an impact on the analysis of the redshift evolution of the halo properties in the different environments. Starting with the halo mass functions, shown in Figure 4 for the cluster, sheet, filament and void environments at , and , there is a substantial change in their shapes when using one or the other of the smoothing approaches. Adopting a constant for the smoothing scale leads to a substantial increase in low mass haloes that are classified to be in clusters relative to the other environments. In voids, sheets and high mass filaments, the mass functions are higher when smoothing with than when adopting a constant ratio; the trend reverses for haloes in low mass filaments and clusters, for which the mass functions are instead enhanced when using the smoothing scale. The inflexion point on scales of () in the mass function of ( clusters for the smoothing reflects the increasing abundance with redshift of isolated small clusters that we mention in Section 2.5. At all redshifts of our study, the clusters have the highest mean number density of haloes, followed by filaments, sheets and voids, respectively.
3.2 Halo Formation Redshift
A closer inspection of the dependence of halo properties on the adopted smoothing scale shows that, with both approaches, there is always a threshold massscale below which the environmental influence on halo properties either switches sign or transitions from being negligible to becoming increasingly more substantial with decreasing halo mass, down to the scales which can be probed with our simulations. Figure 5 shows, for both smoothing approaches, the dependence on environment and mass (in units of ) of the median halo formation redshift . Errorbars are estimates of the error in the median, and are computed as:
(10) 
where is the number of haloes used to sample the distribution of , and and denote the 84th and 16th percentile of the distribution. These values would correspond to errors if the underlying distribution were Gaussian. The confidence region of the median formation redshift determined from all environments is represented by the shaded regions in Figure 5. As discussed in Paper I, this overall behaviour, not split by environment, is well fit by a logarithmic relation over five decades in mass at . This relation, reflecting hierarchical assembly is also present at higher redshifts. We fit a model of the form
(11) 
A leastsquares fit to all haloes extracted from our three simulation boxes provides the following parameters at the three redshifts:
We note that, in all plots and thus at all redshifts and for both smoothing approaches, there is indeed a mass scale at which the curves for the four environments meet, indicating the vanishing of significant environmental influence above this mass threshold. Specifically, below , haloes form earlier in clusters than in filaments, followed by sheets and then voids. This difference in formation time increases with decreasing mass below the threshold .
It is clear from Figure 5 that, in the case of a smoothing scale that remains constant with redshift, the threshold coincides with the massscale for gravitational collapse at , but strongly departs (and increases relative to ) at higher redshifts. Interestingly, however, in the case of the smoothing scale, the threshold mass is easily identified to lie very close to at all redshifts. This difference is simply due to the different association of haloes to the cluster and filament environments in the two smoothing approaches that is also observed in the mass functions (cf. Figure 4). Furthermore, we note that the and relations are very similar to the one at , for which we had already provided the analytic fits with environmentdependent slopes for masses in Paper I. The environmental dependence of the halo formation redshifts below the mass scale that we have found agrees with the analysis of Gao & White (2007), who find that haloes with higher formation redshifts are more strongly clustered.
Interestingly, at the highest redshift of our study (), and for the smoothing, our simulations show that haloes with masses in filaments tend to have earlier formation times than haloes of similar masses in the cluster environment, i.e., an opposite trend than the one observed at all redshifts below .
The fact that, in the case of a fixed smoothing ratio, the environmental dependence is explicitly seen to appear around , hints at a physical relevance of this scale in establishing the onset of the environmental dependence of halo properties at all redshifts. This motivates us to identify the constant ratio as the more fundamental smoothing scale in our analysis, and thus to use this smoothing scale in the remainder of our study of the redshift evolution of halo spins, shapes and alignments as a function of environment.
3.3 Halo Spin
Figure 6 shows the median spin parameter as a function of normalised mass for haloes at , and in the cluster, sheet, filament and void environments. We also plot the results integrated over all environments (shaded region in Figure 6), for comparison with previous studies. In agreement with these (e.g. Vitvitska et al., 2002), we do not observe any significant evolution of the global spin parameter with redshift. Our results also confirm a weak mass dependence of the median spin parameter as found at by Bett et al. (2007), and furthermore extend this result to significantly lower masses. Integrated over all environments, the spinmass relationship is weakly increasing up to , while the highest masses at each redshift have again a slightly smaller median spin parameter.
We find some dependence of the halo spin on environment at all redshifts of our analysis. More specifically, highmass () haloes in filaments have a higher median spin at all redshifts than comparablemass haloes in the cluster environments; in the voids, haloes with masses substantially below spin systematically more slowly (i.e. median ) than haloes of similar masses in any other environment. We have tested that this trend of void haloes having lower median spin parameters persists and actually increases when a larger scale is adopted for the smoothing, to optimize the identification of the void regions (cf. Paper I). At the lowest masses we do not observe any significant difference between haloes in clusters, filaments or sheets.
Gao & White (2007) report that, in their simulations, the most rapidly spinning haloes are more clustered than the slowest spinning haloes, which is in agreement with our earlier findings at (Paper I). Our results, after removing unrelaxed haloes as described in Section 2.1, do not support a very strong correlation between environment and spin at low masses. Still, it is slightly more likely to find the most rapidly spinning objects in environments of higher median density.
3.4 Halo Shape
In Figure 7, we show the median sphericity and triaxiality of haloes in the four environments as a function of their mass for redshifts , and . Independent of redshift and environment, haloes tend to be more spherical with decreasing mass. Over a large range of masses, haloes at are however less spherical than haloes of similar mass at ; AvilaReese et al. (2005) and Allgood et al. (2006) find a similar result in their analyses. The massshape relations, expressed as scalefree functions of , show no significant evolution with redshift up to ; the fact that the entire redshift evolution of the shapes of haloes is driven by the evolution of the mass scale for gravitational collapse, , is also supported by the independent studies quoted above.
Similarly to the case (Bett et al., 2007, Paper I), the masssphericity and the masstriaxiality relation of the global (i.e., not split for environment) sample follow a broken logarithmic relation also at high redshifts. The change in slope in these relationships occurs around . We detect however a relatively small but systematic difference in the sphericity and triaxiality of halos in low and highdensity environments, i.e., the slope in the regime weakly depends on the environment. Specifically, the median sphericity of haloes decreases systematically from the cluster environments, to the filaments, sheets and voids. At all redshifts, a weak trend is observed for haloes with masses below to be more oblate in clusters than in filaments, and more prolate in sheets than in filaments; for haloes above , there is a stronger evidence for haloes in filaments to be systematically more oblate than in clusters. These environmental differences at low masses are observed to be already in place at ; the resolution of our simulations is not adequate to properly investigate these effects at (minimum halo mass ).
3.5 Halo Alignments
3.5.1 HaloLSS alignment
Extending the analysis of Paper I to redshifts above zero, we use the directional information derived from the eigenstructure of the tidal field tensor to estimate the alignment of halo spin with the dynamical properties of the surrounding environment. Filaments and sheets have a preferred direction given by the eigenvector corresponding the single positive or negative eigenvalue. The eigenvectors indicating the direction of the filament as determined from the tidal tensor are shown in Figure 8. Given these unit eigenvectors , we compute the alignment angle . Figure 9 shows the median alignment as a function of mass at redshifts , and . At all redshifts, there is a strong tendency for sheet haloes to have a spin vector preferentially parallel to the sheet, i.e. orthogonal to the normal vector. At redshifts up to , where the errorbars of our measurements allow us to investigate trends with halo mass, this alignment increases with increasing mass. For filament haloes, there is a clear trend with halo mass: (i) haloes with masses smaller than about have spins more likely aligned with the filament in which they reside; (ii) haloes in the range to appear to be randomly aligned with respect to the largescale structure; and (iii) For , the trend appears to reverse, and more massive haloes have a weak tendency to spin orthogonally to the direction of the filament at lower redshifts^{1}^{1}1The tendency for haloes above to spin orthogonal to the host filament, shown in Figure 9 for the smoothing case, is enhanced when the smoothing is adopted. The smoothing scale not only determines the environmental split of the halo population, it also affects the scale on which the eigenvectors of the tidal field are computed. When the smoothing is performed with , the filament direction is obtained on increasingly smaller comoving scales at higher redshifts. This partially erases the stronger correlation that is observed for the most massive haloes when the smoothing is kept at constant comoving scale for all redshifts..
To further explore possible connections between the alignment of the largescale structure and the intrinsic alignment of haloes in the different environments, we search for a correlation signal between the LSS and the axis vectors of the moment of inertia ellipsoid of the haloes. In particular, we use the major axis vector to define the alignment angle , where is again the eigenvector normal to a sheet or parallel to a filament. The resulting median correlation is shown in Figure 10. We find no alignment for halo masses ; however, in both the filaments and the sheets, the halo major axis appears to be strongly aligned with the LSS for masses above about a tenth of . The strength of the alignment grows with increasing mass. This is possibly to be expected, especially for the most massive haloes, since their shape might influence the potential from which the eigenvectors are derived. Adopting a fixed smoothing scale results merely in a shift of the relations shown in Figure 10.
Results similar to ours concerning the alignments of shapes and spins with the LSS, and the transition of alignment orientation at in the filaments, are reported by AragónCalvo et al. (2007) for haloes using a definition of environment that is based on density rather than, as in our case, on the gravitational potential, as well as for haloes in the vicinity of clusters by Basilakos et al. (2006) using the moment of inertia ellipsoid of superclusters and by RagoneFigueroa & Plionis (2007) defining environment by the distance to the nearest cluster. It is clear from our present analysis that such alignments are in place at redshifts of order one, and are maintained virtually unchanged over the last eight or more billion years of evolution of structure in the universe.
3.5.2 HaloHalo alignments
We finally compute the spinspin and spinorbit correlation functions using the definitions of Porciani et al. (2002) and Bailin & Steinmetz (2005). While we show the results for the smoothing case, we stress that qualitatively the results remain unchanged when the constant smoothing is adopted.
For the spinspin correlation we have:
(12) 
where is the intrinsic angular momentum of each halo, and the average is taken over all pairs of haloes which are separated by a distance and reside in the same environment class. Similarly, the spinorbit correlation is defined as:
(13) 
where is the relative orbital angular momentum between two haloes separated by a distance .
Figure 11 shows the spinspin alignment for haloes in clusters and filaments at the three redshifts of our study; upper and lower panels show respectively the results for haloes with masses below and above . The shaded region shows the confidence area for the total sample, split by mass but not split by environment. The correlations within either of the environmental classes is never stronger than those for the total sample and all of them are consistent with no signal within . Furthermore, we find no evidence for any significant redshift evolution of these correlations.
The spinorbit correlation function is shown in Figure 12. The strong correlation that we found at in Paper I, extending out to several Mpc, is present also out to redshift with no significant changes.
4 Summary and Conclusions
We have used three Nbody simulations, tailored to cover a range of almost five decades in mass with high resolution haloes ( particles), to investigate the dependence of halo shape, spin, formation redshift and alignment as a function of mass, environment and redshift. Using the tidal stability criterion of Paper I we have classified haloes to reside in four different environments: clusters, filaments, sheets and voids. The attribution of haloes to these environments depends on one free parameter, , the length scale used to smooth the underlying mass distribution. Relating this length scale to the mass contained in the Gaussian filter, , in Paper I we optimised by visual inspection the redshift zero value of , with the typical mass scale collapsing gravitationally at redshift . At the higher redshifts that we study in this paper, we discuss two possible choices for the smoothing mass scale: i) a smoothing scale constant with redshift; and ii) a smoothing scale that varies such that remains constant with redshift. The first approach leads to the median overdensity in each environment increasing just as expected from nonlinear enhancement of density fluctuations; the second approach maintains the median density in each environment constant with redshift.
In our analysis of the redshift evolution of the halo properties we find that, when adopting a constant ratio for the smoothing, the environmental influence is roughly invariant with redshift so that the mass scale at which the environmental influence sets in is roughly given by by the mass scale . Unveiling the importance of this physical mass scale in the onset of an environmental dependence of the halo properties is a first step towards understanding the origin of the environmental role in the evolution of dark matter haloes.
Adopting the physicallymotivated smoothing at all redshifts, we have investigated the dependence of the properties of isolated dark matter haloes of masses below and above the threshold on their environment. In general, we find that all halo properties show some dependence on environment for halo masses . The strength of the correlations, however, does not change much with redshift. There is virtually no redshift evolution of the halo properties when the correlations with halo mass are expressed in terms of the normalised mass , indicating that the strongest evolution with redshift is related to the evolution of the mass scale for collapse, . In detail, our main results are summarised as follows:

There is a strong environmental dependence of halo formation times with environment for haloes with masses . At any given mass in this mass regime, haloes in clusters tend to be older than haloes in the other environments, and haloes in voids form much more recently than in any other environment. The global halo population (with no splitting for environment) is well described by a logarithmic relation between mass and formation redshift with a roughly constant slope with redshift.

The median spin parameter of the total halo sample, not split by environment, increases weakly with mass up to around at all redshifts. There is no significant residual dependence on redshift besides the mass rescaling with . There is an additional tendency for high mass haloes in filaments to spin more rapidly than haloes of the same mass in clusters. Haloes in voids have the lowest median spin parameters.

Haloes of a given ratio have very similar median shape parameters independent of redshift in the mass and redshift ranges that we have investigated. Independent of environment, haloes are increasingly more spherical and less triaxial the lower their mass. Haloes with masses are slightly more spherical and more oblate in clusters than in filaments, and there is a hint that the situation reverses for , i.e., that haloes in filaments are more oblate than cluster halos at high masses. Lowmass haloes in voids have systematically the lowest median sphericity of similar mass halos in denser environments.

In the mass regime, haloes in sheets tend to have spin vectors in the plane of the sheets, and haloes in filaments tend to have spin vectors pointing along the filaments; above the mass scale, there is evidence that haloes in filaments reverse the previous trend and tend to have spins orthogonal to the filaments. Furthermore, haloes with masses tend to have their major axis parallel to their host sheets or filaments, with the strength of the alignment increasing with increasing mass. This may reflect the fact that, for massive haloes, the gravitational potential field is substantially influenced by their shape and thus leads to an aligned tidal field. The alignment of halo spins and major axes with the largescale structure that we have unveiled up to redshifts of order should be taken into account in studies of weak lensing maps of cosmic shear (eg. Catelan et al., 2001), especially in sheets and thus in regions surrounding voids.

There is no evidence for a significant spinspin correlation between neighbouring haloes. There is in contrast a substantial halo spinorbit alignment, whose strength appears to be independent of mass, environment and redshift up to : haloes in close pairs tend to spin preferentially parallel to the orbital angular momentum of the pair.
An important conclusion that we draw from our study is that the environmental influence on halo properties shows an intriguing dependence on the halo mass, and appears to be essentially modulated by the typical mass scale of gravitational collapse at each redshift. Our data suggests that the mass scale might indeed play the role of a bifurcation point below which many of the median properties of dark matter haloes either begin to feel the influence of their largescale environment, or show an opposite response to their largescale environment relative to the more massive haloes. The existence of such a thresholding mass scale in the environmenthalo relationship is yet to be understood.
Acknowledgements
OH acknowledges support from the Swiss National Science Foundation. All simulations were performed on the Gonzales cluster at ETH Zurich, Switzerland.
References
 Allgood et al. (2006) Allgood B., Flores R. A., Primack J. R., Kravtsov A. V., Wechsler R. H., Faltenbacher A., Bullock J. S., 2006, MNRAS, 367, 1781
 AragónCalvo et al. (2007) AragónCalvo M. A., van de Weygaert R., Jones B. J. T., van der Hulst J. M., 2007, ApJ, 655, L5
 AvilaReese et al. (2005) AvilaReese V., Colín P., Gottlöber S., Firmani C., Maulbetsch C., 2005, ApJ, 634, 51
 Bailin & Steinmetz (2005) Bailin J., Steinmetz M., 2005, ApJ, 627, 647
 Basilakos et al. (2006) Basilakos S., Plionis M., Yepes G., Gottlöber S., Turchaninov V., 2006, MNRAS, 365, 539
 Bertschinger (2001) Bertschinger E., 2001, ApJS, 137, 1
 Bett et al. (2007) Bett P., Eke V., Frenk C. S., Jenkins A., Helly J., Navarro J., 2007, MNRAS, 376, 215
 Blanton et al. (2005) Blanton M. R., Eisenstein D., Hogg D. W., Schlegel D. J., Brinkmann J., 2005, ApJ, 629, 143
 Bond et al. (1991) Bond J. R., Cole S., Efstathiou G., Kaiser N., 1991, ApJ, 379, 440
 Bryan & Norman (1998) Bryan G. L., Norman M. L., 1998, ApJ, 495, 80
 Bullock et al. (2001) Bullock J. S., Dekel A., Kolatt T. S., Kravtsov A. V., Klypin A. A., Porciani C., Primack J. R., 2001, ApJ, 555, 240
 Catelan et al. (2001) Catelan P., Kamionkowski M., Blandford R. D., 2001, MNRAS, 320, L7
 Croton et al. (2007) Croton D. J., Gao L., White S. D. M., 2007, MNRAS, 374, 1303
 Davis et al. (1985) Davis M., Efstathiou G., Frenk C. S., White S. D. M., 1985, ApJ, 292, 371
 Doroshkevich (1970) Doroshkevich A. G., 1970, Astrophysics, 6, 320
 Dressler (1980) Dressler A., 1980, ApJ, 236, 351
 Eke et al. (1996) Eke V. R., Cole S., Frenk C. S., 1996, MNRAS, 282, 263
 Gao et al. (2005) Gao L., Springel V., White S. D. M., 2005, MNRAS, 363, L66
 Gao & White (2007) Gao L., White S. D. M., 2007, MNRAS, 377, L5
 Hahn et al. (2007) Hahn O., Porciani C., Carollo C. M., Dekel A., 2007, MNRAS, 375, 489 (Paper I)
 Harker et al. (2006) Harker G., Cole S., Helly J., Frenk C., Jenkins A., 2006, MNRAS, 367, 1039
 Hetznecker & Burkert (2006) Hetznecker H., Burkert A., 2006, MNRAS, 370, 1905
 Kauffmann et al. (2004) Kauffmann G., White S. D. M., Heckman T. M., Ménard B., Brinchmann J., Charlot S., Tremonti C., Brinkmann J., 2004, MNRAS, 353, 713
 Lacey & Cole (1993) Lacey C., Cole S., 1993, MNRAS, 262, 627
 Macciò et al. (2006) Macciò A. V., Dutton A. A., van den Bosch F. C., Moore B., Potter D., Stadel J., 2007, MNRAS, 378, 55
 Maulbetsch et al. (2007) Maulbetsch C., AvilaReese V., Colín P., Gottlöber S., Khalatyan A., Steinmetz M., 2007, ApJ, 654, 53
 Peebles (1969) Peebles P. J. E., 1969, ApJ, 155, 393
 Porciani et al. (2002) Porciani C., Dekel A., Hoffman Y., 2002, MNRAS, 332, 325
 RagoneFigueroa & Plionis (2007) RagoneFigueroa C., Plionis M., 2007, MNRAS, 377, 1785
 Reed et al. (2007) Reed D. S., Governato F., Quinn T., Stadel J., Lake G., 2006, MNRAS, 378, 777
 Sheth & Tormen (2004) Sheth R. K., Tormen G., 2004, MNRAS, 350, 1385
 Springel (2005) Springel V., 2005, MNRAS, 364, 1105
 Vitvitska et al. (2002) Vitvitska M., Klypin A. A., Kravtsov A. V., Wechsler R. H., Primack J. R., Bullock J. S., 2002, ApJ, 581, 799
 Wang et al. (2007) Wang H. Y., Mo H. J., Jing Y. P., 2007, MNRAS, 375, 633
 Wechsler et al. (2006) Wechsler R. H., Zentner A. R., Bullock J. S., Kravtsov A. V., Allgood B., 2006, ApJ, 652, 71
 Wetzel et al. (2007) Wetzel A. R., Cohn J. D., White M., Holz D. E., Warren M. S., 2007, ApJ, 656, 139
 Zel’dovich (1970) Zel’dovich Y. B., 1970, A&A, 5, 84