A&A 378, 345360 (2001)
DOI: 10.1051/00046361:20011240
M. O. Hanski^{1}  G. Theureau^{2,3}  T. Ekholm^{4,1}  P. Teerikorpi^{1}
1  Tuorla Observatory, 21500 Piikkiö, Finland
2 
Observatoire de Meudon, 92195 Meudon Cedex, France
3 
Osservatorio di Capodimonte, Via Moiariello 16, 80131 Naples, Italy
4 
Observatoire de Lyon, 69561 SaintGenis Laval Cedex, France
Received 24 June 1999 / Accepted 5 September 2001
Abstract
The matter distribution around the PerseusPisces (PP) supercluster is studied by comparing
peculiar velocities given by the Kinematics of the local universe (KLUN) galaxy sample to those predicted by TolmanBondi (TB) models. To restrict the TB solutions we first solve the mass of the densest part of PP. This part is identified as a sphere at
,
Mpc having a radius of
Mpc. This sphere surrounds the main part of the PP ridge and four most prominent clusters of the region. Using viriallike mass estimators we calculate the cluster masses and obtain the upper and lower limits for the mass inside the
Mpc sphere:

.
This corresponds to a mass overdensity
,
or

.
Mass to light ratios of the clusters are M/L = 200  600
,
giving mass density ratio
 0.3, if the value of M/L is assumed to be representative elsewhere in the universe.
We estimate a radial density distribution around the PP core using two toy models and a smoothed density distribution observed for IRAS galaxies. The cosmological density parameters
and the PP mass are free parameters in the TB calculations. The KLUN velocities, obtained by TullyFisher relation and the normalized distance method, are adjusted by the Local Group (LG) infall velocity towards PP. Comparison of the TB velocities to KLUN data points indicates that the infall velocity
km s^{1}. Allowing
to vary within the limits given above we get constraints for the value of ;
0.4 are prefered to the more extreme values,
or 1. A choice of either
or 0 do not cause any significant changes in the results.
The validity of the TB model in complex environments is studied with an Nbody simulation. There we see that the radially averaged velocity fields around simulated clusters are compatible with the corresponding TB velocities. This confirms the applicability of the TB model around large galaxy concentrations, providing that smoothed density and radially averaged velocity fields are used.
Key words: galaxies: clusters: general  galaxies: clusters: PersPisc supercluster  cosmology: theory  largescale structure of the universe
TB results are valid with appropriate smoothing; the observed velocity field must be averaged over concentric shells around a galaxy concentration, when compared to the TB results. Thus the TB model does not give the results in as much details as the various more commonly used gravitational instability  reconstruction methods. On the other hand, it does not require any assumptions for the breakdown of the theory in the nonlinear region. In order to test whether the TB model may be applied to an apparently still less favourable situation than Virgo, we study the PerseusPisces supercluster, where a TBlike deviation in the Hubble flow has already been observed (Theureau et al. 1998a). Our criterion for the success of the TB method is whether the mass inferred using TB is comparable to virial and luminosity M/L estimate for the region, and whether the peculiar velocity field is adequately described by the model.
The local d<70 h^{1} Mpc^{} universe is dominated by two large mass concentrations, the Great Attractor and the PerseusPisces (PP) complexes. According to Sigad et al. (1998) and Dekel et al. (1999) both regions have overdensities 1.5, where , is the local, and the average density. Both studies used a Gaussian smoothing of 12 h^{1} Mpc and assumed . Thus the velocity field surrounding PP should reveal a strong infall towards the supercluster's center of mass. For a galaxy at around half way from PP to us the infall velocity should be 500 km s^{1}. However, the infall signature is difficult to reveal because of large uncertainties in redshift independent distance measurements, e.g. those provided by TullyFisher (TF) or Fundamental Plane (FP) methods. Uncertainties are caused by statistical biases difficult to correct for in a heterogeneous environment. Additional inaccuracies are raised by the filamentary shape of the structure, consisting of an alignment of several subclusters, which complicates considerably any attempt to model the velocity field.
In the past ten years PP has been investigated by numerous authors using various distance criteria and statistical methods. The region is now well covered by great amount of data. From a large reshift sample of 3311 galaxies in the region bounded by 22 RA , +20 Dec. , Wegner et al. (1993) described the PerseusPisces complex as follows: the supercluster consists of a continuous arrangement of high density clusters and groups, among which the Pisces, A262, and Perseus (A426) clusters are the most prominent features. Its main ridge extends at least 50 h^{1} Mpc from Pegasus eastward to the Perseus cluster, with a width in the plane of the sky of about 510 h^{1} Mpc, and a typical depth in redshift of 250500 km s^{1}. All the objects in the ridge are lying roughly at the same distance from us, at 5000 km s^{1}.
First kinematical studies (Willick 1990, 1991; Courteau et al. 1993) based on an rband TF sample of 355 galaxies, suggested that the whole structure takes part on a large scale flow defined by the Great Attractor (see e.g. LyndenBell et al. 1988; Dekel 1994), with a mean velocity of 350 km s^{1} towards the Local Group (LG). Han & Mould (1992) confirmed this result independently using a whole sky TF sample of 21 clusters. Hudson et al. (1997) combined their inverse FP relation measurements of seven clusters in the PP region with data of nine other clusters from the literature and concluded that the whole region has a mean CMB frame bulk motion of km s^{1} towards galactic coordinates . POTENT reconstruction method of the density and velocity field for Mark III galaxies (Dekel et al. 1999) gave different results. There it was found that the center of PP is roughly at rest in the LG frame and the regions between the LG and PP are infalling towards PP.
Most studies of the velocity field of galaxies suggest that the kinematics of the local universe can be described by a model in which galaxies are infalling to two mass concentrations, one in PP and one in HydraCentaurus (e.g. Han & Mould 1992; Dekel et al. 1999). However, Freudling et al. (1995) and da Costa et al. (1996) used a method of analysis, based on MonteCarlo simulations, in order to control biases and gain a better sampling of the sky. Their studies revealed a clear infall motion towards PP, while the GA was found less dominant than originally believed.
The differences in the results stated above may be caused by uncertainties in the reconstruction method, in the corrections for the selection effects, and/or in the distance indicator calibration. For example, an error in the TF zero point would cause a systematic error that resembles a Hubble flow in the peculiar velocity field. Dekel et al. (1999) also note that the results in the PP region are more sensitive to noise than the other regions.
In a previous study using the KLUN galaxies, a full infall pattern, both in the front and the backside of the PP supercluster was shown, for the first time. A new Malmquist bias correction method for the TF distance data was used (Theureau et al. 1998a). In the present paper we investigate the dynamics of galaxies in the PP surroundings in more detail. We study the spatial distribution and peculiar velocities of the galaxies using the Malmquist corrected TF distances. The observed velocity field is averaged over spherical shells around the PP center and compared to the results given by TB calculations.
The outline of the paper is as follows: in Sect. 2 we briefly present the data and summarize the normalized distance method (NDM) that has been used for deriving unbiased average distances and peculiar velocities. In Sect. 3 we calculate the virial masses of the most prominent clusters in the PP region and estimate the mass of the whole central region of PP from these masses. Section 4 includes a description of the TolmanBondi model and its application to PP. The model is also tested on an Nbody simulation of galaxy clusters (Sect. 4.2). In Sect. 5 we compare the TB model to the infall pattern observed with the KLUN data. Finally, the main results are summarized in Sect. 6.
The main difficulties in velocity field studies are caused by the large uncertainties in distance determination and by the resulting biases. Attempts to overcome these problems require a few things. Firstly, one must know thoroughly the sample, how it was collected, and to which extent it can be assumed complete. Secondly, systematic effects in the distance criterion have to be modeled, so that one can derive for a given application an unbiased set of distances. In particular, biases will be different if one estimates peculiar velocities at a given derived distance, at a given true distance, or at a given observed radial velocity. Thirdly, the sample has to be deep enough and should cover a sufficiently large solid angle in the sky for the average relations, needed for the bias correction, to be free from local perturbations in the density or the velocity field. This last requirement is similar as needed for a satisfactory determination of the Hubble constant.
Figure 1: The vs. diagram of KLUN galaxies. The line, indicating the average , is (approximately) horizontal in the unbiased region at small . Malmquist bias makes the curve to rise at larger .  
Open with DEXTER 
We use the KLUN (Kinematics of the local universe) sample consisting of 6600 spiral galaxies having measurements of isophotal diameter D_{25}, HI line width, radial velocity, and Bmagnitude. This sample has been fully described in previous papers of this series, e.g. by Theureau et al. (1997), so we only summarize here the main characteristics.
KLUN is selected according to apparent diameter, and is complete down to (Paturel et al. 1994), covering the type range SaSdm (T=18). The data are extracted from LEDA (LyonMeudon Extragalactic Database, Paturel et al. 1997) and complemented with optical redshifts measured with ESO and OHP telescopes, and HI spectra obtained with Nançay and Parkes radio telescopes (Bottinelli et al. 1992, 1993; di Nella et al. 1996; Theureau et al. 1998b). According to LEDA precepts, all astrophysical parameters have been homogenized and reduced to a standard and common system (Paturel et al. 1997). Isophotal D_{25} diameters and apparent Bmagnitudes are corrected for galactic extinction according to Fouqué & Paturel (1985), and for inclination, i.e. the opacity effect, as explained in Bottinelli et al. (1995). HI line widths, reduced to the standard levels of 20% and 50%, are corrected for internal velocity dispersion according to Tully & Fouqué (1985). Heliocentric radial velocities are changed to the rest frame of the LG centroid following Yahil et al. (1977) and corrected for a Virgocentric component, assuming an infall of the LG towards the Virgo cluster v_{0} = 200 km s^{1} (Theureau et al. 1997; Jerjen & Tammann 1993) and a radial velocity of Virgo km s^{1} (Mould et al. 1980), within a linear infall model (Peebles 1976).
In the TF relation studies we need to make some restrictions to the data. Galaxies close to the Galactic plane ( ) are excluded due to the large uncertainties in the Galactic extinction correction. Faceon galaxies, , are excluded to prevent large errors on calculation of the rotational velocity parameter , where is the maximal rotational velocity of a galaxy, measured from the HI line width.
Absolute TF distances are obtained using the calibration given by Theureau (1998). TF scatter was there reduced by 30% by taking into account the dependence of the TF zero point on mean surface brightnesses of the galaxies. The NDM first order Malmquist corrections are calculated on the basis of the whole sky sample, according to Theureau et al. (1998a). The sample with thus corrected TF distances, restricted to the strictly diameter complete part, contains finally 2800 spiral galaxies.
The concept of normalized distances was introduced in a theoretical discussion by Teerikorpi (1984). The method was first applied by Bottinelli et al. (1986) and improved in a series of papers, most recently by Theureau et al. (1997). A theoretical description in terms of probability densities was given by Theureau et al. (1998a), where the authors showed that the method can be extended to make Malmquist bias corrections to galaxy distances even beyond the socalled unbiased plateau. The unbiased plateau has been profitably used for determination of the unbiased value of H_{0} and has the advantage that one does not need modeldependent bias corrections.
In order to understand how the NDM works, let us define classes of galaxies as follows: objects belong to a given class if they have identical values of log rotational velocity p, Galactic extinction correction ,
inclination R_{25}, morphological type T and mean surface brightness .
The main idea of the NDM is that for a given magnitude or diameter limited sample, and for a given TF scatter the behaviour of the Malmquist bias^{} as a function of true distance is similar for all galaxies belonging to the same class. The class is assumed to have a Gaussian distribution of absolute magnitudes. Then, for each class, there is a unique curve of evolution of the bias, which is conveniently expressed as the curve of observed average value of the Hubble parameter
vs. distance (d), where
is the distance given by TF relation,
(2) 
Figure 2: LEDA (o) and KLUN galaxies (+, * and ) in the PP region. Large circle is the 15 h^{1} Mpc sphere at containing the main concentration. Small circles are the four densest clusters, Perseus, A262, 0122+3305 and Pisces, from up left to down right.  
Open with DEXTER 
Essentially, the normalized distance technique means superimposing all the bias curves
of the different classes of objects by using an appropriate scaling on the distance:
(3) 
The average bias correction is valid only if the following assumptions are fulfilled:
As for assumption H4, a list of recent measurements of largescale bulk motions is given by Strauss (2000). The list is somewhat divergent; roughly half of the surveys find no bulk flow at scales 6000 km s^{1}, the rest claim bulk flows of 300600 km s^{1}. We will return to the question of largescale flow within the KLUN sample in the future, for now we assume that assumption H4 is true for our data.
Finally, the corrected distance is given by
(4) 
We emphasize that this correction is statistical in nature. When applied to an individual galaxy, it gives an accurate average correction only if the velocity field model (providing and d) is correct. When the velocity field is distorted, say, in the direction of a supercluster, it still gives a first order correction, but a complete correction would require an iterative procedure (see Sect. 5.3). Even the first order correction essentially improves the derived velocity field, as was shown by Theureau et al. (1998a) for the PP supercluster, but it is notable that the peculiar velocities will be somewhat suppressed if only a first order correction is done.
Figure 3: The Hubble diagram of KLUN galaxies in the PP region. is the NDM corrected TF distance.  
Open with DEXTER 
Figure 2 shows a map of the PP region. Circles stand for galaxies from LEDA, having radial velocities between 3000 and 7000 km s^{1}. This limitation allows us to look at the main part of the PP concentration, lying at km s^{1}, practically free from any contamination from fore/background objects. LEDA is regularly updated for all published data and thus contains all the recent PPsurveys.
Symbols other than circles are KLUN galaxies. +:s are objects within the velocity interval stated above, *:s have V>7000 km s^{1} and :s V<3000 km s^{1}. Only the KLUN galaxies satisfying the completeness limit, , and other criteria demanded by our TFmethod (see Sect. 2.2), are plotted. In total, Fig. 2 contains 2016 LEDA galaxies and 266 objects from KLUN.
l  b  N  V  M  M/L  
1  150.5  13.5  152  5144  2.0  1.8  812  310490 
2  136.7  25.0  104  4878  2.0  1.7  46  210370 
3  130.5  29.0  67  5085  1.5  1.3  46  260510 
4  127.1  30.0  63  5096  1.5  1.3  35  270550 
We seek a rough estimate for the amount of matter in the PP region. First we determine the masses of the densest clusters in the region. We define the densest consentrations as the ones having the projected galaxy overdensity . Here , the two dimensional sky densities are calculated in 1 degree^{2} regions, and is the mean value of the whole region shown in Fig. 2. The limiting factor 10 has no special physical meaning. It was chosen to give a reasonable number of dense galaxy concentrations. There are four clusters satisfying this criterium. The characteristics of these clusters are listed in Table 1. Masses are estimated using the virial theorem and the alternative methods presented by Heisler et al. (1985), see Table 2 below. The effects of possible outsiders and interlopers are tested by using various radii and redshift intervals in defining cluster members. Table 2 shows some of the results. For the mass evaluations in Table 1 we adopt Vinterval of 30007000 km s^{1}, and use virial radii of for Perseus and A262, and for 0122+3305 and Pisces. These radii are drawn in Fig. 2.
N  M/L  
Perseus  
115  5107  5.5  6.6  4.9  5.5  244362  
131  5131  6.4  7.9  6.1  6.9  275393  
152  5144  8.1  11.8  8.1  10.1  309493  
172  5119  9.4  14.3  10.0  12.5  314525  
183  5148  10.3  16.1  11.2  14.2  320546  
A262  
62  4965  2.8  3.4  2.3  2.7  195342  
93  4915  3.6  4.1  2.9  3.5  168284  
104  4878  4.3  5.6  3.8  4.6  206369  
116  4892  5.1  7.0  4.7  5.6  220405  
132  4924  5.9  8.1  5.6  6.6  233410  
0122+3305  
49  5036  2.7  2.9  2.5  2.5  236333  
67  5085  4.5  5.9  3.6  4.7  260511  
88  5063  5.6  7.1  4.5  5.8  261502  
98  5037  6.2  7.9  5.1  6.5  266499  
125  5052  7.1  9.1  5.5  7.5  220440  
Pisces  
47  5073  1.9  2.9  1.7  2.1  203409  
63  5096  2.9  4.6  2.7  3.3  270550  
77  5105  4.2  7.1  3.7  5.1  308713  
95  5107  5.2  8.3  4.7  6.1  329702  
122  5187  9.0  16.2  7.1  11.6  3851074 
Figure 4: a) shows the observed absolute magnitude histograms for the Perseus galaxies (gray bins) and the three other clusters combined (hatched bins). Two Schechter function curves are drawn on both of these histograms. Lower lines are fitted to the histograms at , upper lines are calculated by assuming completeness of for Perseus and for the other clusters, following the latitude dependent completeness estimate by Marinoni et al. (1999). b) shows the corresponding total luminosities in 0.5 magnitude bins. The lines are the values predicted by the upper Schechter curves of a) By estimating the ratio of the gray/hatched areas to total areas under the Schechter curves in b) we can calculate the amount of luminosity of nondetected galaxies. This is for Perseus and of the total luminosity for the others.  
Open with DEXTER 
For the M/L ratio derivations in Tables 1 and 2 we calculate the total luminosities of the clusters from LEDA Bband brightnesses and estimate selection effects
using the Schechter function
Marinoni et al. (1999) observed for galaxies in large groups and . The difference between our value for Perseus and the one given by Marinoni et al. is probably due to environmental differences between Perseus and other clusters and/or uncertainties in the brightness measurements caused by the Galactic absorption. Comparison of the total observed luminosities of the clusters to the values obtained by integrating Eq. (5) indicates that the contribution of nondetected galaxies is for Perseus and for the others, with bootstrap errors. The limits for the M/L ratios are then calculated with the upper and lower mass values from Table 1 and by estimating the true total luminosities to be , i.e. 1.641.80 (Perseus) and 1.601.95 (others) the observed values. Resulting M/L ratios are all in the range of 210550 , in good agreement with other recent estimations of M/L ratios of rich clusters (e.g. Bahcall et al. 1995; Carlberg et al. 1997).
Adding up the masses of the four clusters in Table 1, we get 23 . The center of mass of PP, using these clusters, is at km s^{1}, . Distances of the clusters from this point are 11.5, 3.9, 9.9 and 12.5 h^{1} Mpc. A sphere at , with a radius of 15 h^{1} Mpc surrounds all these main concentrations. In addition to the four clusters, this sphere contains 513 other galaxies. If we assume that the M/L ratio does not grow in scales larger than the rich clusters, as suggested by Bahcall et al. (1995, 2000), we can then estimate 47 , where is the mass inside the 15 h^{1} Mpc sphere. The selection effects are assumed to be the same for the galaxies in the four clusters and for the other galaxies in the 15 h^{1} Mpc sphere.
However, there is some evidence that M/L increases on scales larger than clusters. For example Small et al. (1998) found the mass to light ratio to be 564 for the Corona Borealis supercluster and 726 for a background supercluster. Small et al. used the virial and projected mass estimators for the whole bound supercluster. The TolmanBondi calculations (next section) show that the PP is also bound, so we may apply the virial estimators on the PP also. We get 2.6, 4.2, 2.8, and 3.3 for the 836 galaxies in the central 15 h^{1} Mpc sphere in PP, using the virial, projected, medium, and average mass estimators, respectively. The corresponding sum of the galaxy luminosities gives 4.4 . Taking from above the more conservative limits for the number of nondetected galaxies, this number is then increased by a factor 1.601.95, giving the total luminosity  . Finally, the mass to light ratio, using these values, is 330600 . However, since this mass for the whole PP core is about equal to the sum of the four clusters derived above, we doubt the accuracy of the virial and other estimators in this case, and adopt the values 47 for .
These upper and lower limits for give a mean density 24 12 , where is the critical density of the universe. In terms of mass overdensity this is 4, assuming 0.20.4. The lower value of here corresponds naturally to the lower limit of , and high to the higher . The result is in agreement with the value ( ) by Sigad et al. (1998) and Dekel et al. (1999).
TolmanBondi (TB) equations determine the kinematics of a spherically symmetric mass concentration in an expanding universe (see e.g. Olson & Silk 1979; Ekholm 1996 for main features of TB, and Teerikorpi et al. 1992; Ekholm & Teerikorpi 1994; Ekholm et al. 1999, 2000 for some recent applications). In practice we study the behaviour of a spherical layer of matter at distance r from the center of symmetry. Each of these shells have two dynamical constants. In Newtonian terms we can think of these constants as the mass inside the radius of the shell,
(6) 
(7) 
We attempt to solve the present day velocities of the mass shells, v_{r}(r), for the spherically symmetric mass distribution, with the initial condition that at some point in the past these mass shells have all been very close together. That is, we assume the big bang to have happened at t=0. This moment is acquired from the underlying Friedmann model, by calculating the age of the Universe
for the set of cosmological parameters H_{0},
,
and
,
where
,
is the mass density, and
Then we find the value of the energy term E(r) for each present day radius r_{0}, by demanding
r(T_{0},E)=r_{0} and r(0,E)=0, and integrating backwards in time. The integration is performed numerically for both r and v_{r} with
For , where the velocities grow considerably, the numerical methods confront some difficulties. Even with adapted stepsize algorithms we can not avoid all the problems. Using a fourth order RungeKutta method, we obtain inconsistent results, , when . Fortunately our main concern is to calculate the dynamics at larger t, where the method proved to be robust. A comparison to the parameterized TB solution showed that the problems at small t did not have any effect on our resulting values of E(r). With the E(r) thus obtained, the present day velocities v_{r}(r_{0})_{E=E(r0)} can be calculated.
Figure 5: The (141 h^{1} Mpc)^{3} cube of GIF Nbody simulation galaxies. The simulation started with 256^{3} dark halos, producing the 15 000 galaxies seen in this figure. The four circles indicate the clusters for which we studied the applicability of the TB calculations.  
Open with DEXTER 
Figure 6: The four GIF galaxy clusters extracted for further study. Figures show XYplanes from the previous figure, centered to the cluster centers. The width in Z direction of each plot is 30 Mpc. Circles with radii are drawn. It is notable that these complexes are not spherically symmetric, resembling the situation with PP.  
Open with DEXTER 
The TB model is a rough approximation of the dynamics of galaxy clustering. The reader may suspect that the assumption of spherical matter distribution around a galaxy cluster is perilous. In reality galaxies seem to accumulate into filamentary and walllike structures. However, an averaged galaxy density field may be well described with a spherical model; the IRAS map, using a 12 h^{1} Mpc smoothing, shows a roughly spherical density field around PP (Sect. 4.3). And the distribution of dark matter may be even more spherical.
How about the peculiar velocities? The manybody gravitational interplay determining the dynamics of galaxy clustering is a complicated issue, can it really be described by smooth radial TB velocities? That may seem doubtful, but, again, we rely on the power of averaging. The TB results are compared only to the average radial velocities around the center of mass.
The correspondence between the TB solutions and the averaged observational quantities is here studied with a cosmological Nbody simulation. In a simulated galaxy cluster we have the advantage of knowing the true distribution of matter and the true velocities, not distorted by selection effects or observational uncertainties. And even better, we can study the evolution of the cluster and compare the matter and velocity distribution of the galaxies to those predicted by the TB model at different epochs.
In what follows we use the Nbody simulation data provided by the GIF group, a GermanIsraeli collaboration working with the Virgo consortium project (http://starwww.dur.ac.uk/~frazerp/virgo/). The GIF team investigates the clustering of galaxies in a hierarchical universe (Kauffmann et al. 1999a, 1999b; Diaferio et al. 1999). They simulate structure formation in a periodical (141 h^{1} Mpc)^{3} cube using a cold dark matter, flat curvature, nonzero cosmological constant model with , , and h=0.7. The simulation starts with 256^{3} dark matter halos with masses . The galaxies are formed at concentrations of halos, and the effects of gas cooling, star formation, and supernova feedback are taken into account. The data files at http://www.mpagarching.mpg.de/Virgo/data_download.html have the positions, velocities, and luminosities of galaxies and distribution and masses of halos at six redshifts, z=0.00, 0.20, 0.42, 0.62, 0.82, and 1.05, corresponding to epochs t/T_{0}= 1.00, 0.97, 0.79, 0.67, 0.57, and 0.49. There are two sets of files, one for the galaxies (luminous matter) and one for the halos (luminous + dark matter).
Figure 5 shows the galaxies in the whole (141 h^{1} Mpc)^{3} cube at z=0.
We study the evolution of four large clusters in the simulation. These clusters are circled in Fig. 5 and shown individually in Fig. 6. For each of these clusters we approximate the central radius
by eye, and calculate the masses and the radially averaged ,
see Table 3. These values are similar to what we
observe in the PP supercluster.
M  
1  10  4.5  8.0 
2  14  7.2  4.3 
3  8  3.0  10.3 
4  10  4.4  7.9 
Figure 7: The density fields around the four simulated clusters at the six epochs between the present and z=1.05. Points show the averaged of the simulation, calculated using a 3 h^{1} Mpc smoothing. Error bars are the standard deviations. Gray line is the observed density field at z=0 evolved backwards in time via the TB model. Thick black line is the TB evolved "step function'' toy model. are the adopted cluster radii and are the average overdensities within these radii at z=0.  
Open with DEXTER 
Figure 8: The evolution of the velocities v_{r}(r) around the four clusters. Crosses and error bars are the averaged v_{r}(r) and their standard deviations observed in the simulation. Gray lines correspond to the observed density distribution and black lines to the toy model density distributions, evolved within the TB models.  
Open with DEXTER 
The overdensities are calculated from the halo number distribution, with a Mpc smoothing. At this smoothing length the matter distribution in the simulation is unbiased according to Kauffmann et al. (1999a). Then we determine the center of mass and calculate the radially averaged around each of the four clusters. For the TB calculations we use either this observed (at z=0) or a "step function'' toy model where , and . Figure 7 shows the evolution of in the simulation and by means of the TB model. The toy model is very instructive in showing how the density enhancements grow in amplitude and converge in comoving radius with time.
As expected, in the simulation the structures form at z>1 and there is not much evolution in from z=1.05 to z=0. Contrarily, in the TB evolved clusters the overdensities are low at z=1.05 but grow quite steadily to the present day values. At larger r the given by TB model follow the corresponding values of the simulation quite well at each epoch.
Figure 8 shows the comparison of the observed averaged peculiar velocities in the simulated clusters to the TB results. A similar time dependent factor as in the evolution of the density field is seen here. At earlier epochs the peculiar velocities given by the TB model are somewhat smaller than the v_{r}(r) in the simulation. At z=0 the TB values are higher, although the differences are not at a significant level.
Our conclusion is that concerning the evolution of density and velocity fields around clusters there is a marginal difference between the simulations and the TB calculations. The overdensities are quite well approximated at large r at each epoch, and the present day peculiar velocities v_{r}(r) are sufficiently well derived, that is, with < deviation, using the observed present day .
Now, how to use TB method for PerseusPisces where the distribution of luminuous matter apparently is not spherically symmetric? Certain restrictions and assumptions must be applied. Firstly, we isolate a sphere containing the main part of the mass, as was done in Sect. 3. Then we assume that any substructures or dynamical features inside the sphere do not affect a galaxy outside it; the sphere can effectively be replaced by a point mass located at the center. Outside the sphere, we assume the density to decrease to the average density of the universe following spherical symmetry:
The assumptions may seem doubtful. Is it not oversimplified to approximate the filamentary structure seen in the PP region with a spherical model? The filaments must produce observable nonradial motions. The main justification for our method lies in the power of averaging. Small scale nonradial peculiar velocities are smoothed away in a similar manner as a smoothed density map blurs away small scale substructure. Compare the smoothed IRAS map by Sigad et al. (1998) to the galaxy map in Fig. 2. The IRAS map shows a spherical density distribution centered in the PP, quite unlike the first glance on the individual galaxy distribution suggests. Nbody simulations indicate (Sect. 4.2) that the same thing happens for the velocity fields, justifying the usage of the TB models. In any case, over 90 percent of the matter is dark and may be more symmetrically distributed. The final criterion is how well the model describes the infall pattern and predicts the total mass.
For the density distribution outside the sphere we consider three models: two toy models, a step function
Figure 10 shows schematically the system under study. A galaxy is located at distance d from us. The lines of sight towards PP and the galaxy form an angle .
The distance from PP to the galaxy is then
Figure 9: Assumed spherically symmetric density distributions , and , and the corresponding masses inside radius r from the PP center. Overdensity and mass M_{1}(r) are plotted with continuous, and M_{2}(r) with dashed, and and with dasheddotted lines. See text for details.  
Open with DEXTER 
Figure 10: Galaxy distances from the PP and LG, and observed and PP rest frame velocities are plotted in this schematic picture. Distances are normalized to the PP distance; , .  
Open with DEXTER 
Figure 11: LEDA galaxies within angular distance of from the PP center. LG is at the apex of the cone, the gray sphere indicates the central PP region. The figure does not give a truthful representation of the galaxy distribution: distances based on redshift velocities are distorted, selection effects are strong in larger distances and projection effects in the regions close to the cone axis. The figure is a projection perpendicular to the Galactic plane, which is marked with a dashed line.  
Open with DEXTER 
Peculiar velocities v_{r} of the KLUN galaxies are obtained in the following manner:
To avoid any difficulties in combining possibly incompatible individual distance measurements we do not use ordinary weighted averages in the bin distance determinations. The distances are deduced with more sophisticated Bayesian methods (Press 1997): probability distributions of "good'' and "bad'' measurements are calculated as
(17) 
The probability that the binned set of data (
)
have a distance
is given by
Actual bin distance is adopted from the peak of the vs. curve, and 95% confidence levels are deduced from calculating the area below the curve (Fig. 12). Figure 13 shows the individual and binned values for all the data points in the KLUNPP sample. For comparison the binned distances are calculated also using the simple weighted mean, see Fig. 12. There are no significant differences to the values obtained with the bayesian method, the two methods deviate <1 for each .
Figure 12: An example of using the vs. curve in the bin distance determination. Individual galaxies are plotted on the axis. Vertical lines show the adopted and the 95% confidence levels. For comparison, the simple weighted mean bin distance is also drawn with error bars. KLUNPP galaxies with 2500 km s^{1}<V< 3000 km s^{1} are shown in this example.  
Open with DEXTER 
Figure 13: KLUNPP galaxies, bin distances and 95% error bars in a V vs. d diagram.  
Open with DEXTER 
The peculiar velocity field around PP is determined similarly. Firstly each galaxy is assigned its bin distance ()
and 95% limit errors (denoted by
). PP distances ()
are calculated by Eq. (12), and errors by
Peculiar velocities are obtained with Eq. (14) and errors by
Galaxies are then binned into PP distance intervals of . For each bin we calculate PP distances and peculiar velocities with 95% confidence levels using the Bayesian methods stated above. Notice that for the peculiar velocity calculations we need to assume a value for our infall velocity towards PP. That is, we have in Eq. (14), as a free parameter.
Figure 14: Binned KLUNPP data points and TB v_{r}(r) curves. KLUN points are fixed with . The error bars are confidence intervals. TB curves are calculated for the three density models, , , and . For each density, there is a set of curves, corresponding to the seven different cosmological models, (0.1,0.9), (0.1,0.0), (0.2,0.8), (0.2,0.0), (0.4,0.6), (0.4,0.0), and (1.0,0.0).  
Open with DEXTER 
We study how well our TB models produce the velocity field of the KLUN galaxies when the models are varied by
To put this into numbers, we describe the goodness of fit by
0.1  0.9  2.7_{ 0.6}^{+ 1.4}  0.47_{ 0.14}^{+ 0.08} 
0.1  0.0  3.1_{ 1.0}^{+ 1.4}  0.41_{ 0.11}^{+ 0.14} 
0.2  0.8  3.9_{ 1.4}^{+ 1.8}  0.49_{ 0.15}^{+ 0.24} 
0.2  0.0  3.9_{ 1.4}^{+ 1.8}  0.49_{ 0.15}^{+ 0.24} 
0.4  0.6  5.5_{ 1.4}^{+ 2.2}  0.53_{ 0.16}^{+ 0.20} 
0.4  0.0  5.5_{ 1.4}^{+ 2.2}  0.53_{ 0.16}^{+ 0.20} 
1.0  0.0  9.4_{ 1.8}^{+ 2.9}  0.55_{ 0.17}^{+ 0.20} 
Figure 15 shows the formal confidence regions for the values of and obtained from fitting the TB curves to the KLUN data points. The confidence regions are defined in the following manner: the value of is calculated for each pair of parameters ( , ). Then a point in the (v, M) space belongs to the confidence region if , where is the minimum value of in the whole region of v=0600 km s^{1} and M=2 , and is the average of the errors of the binned KLUN velocities. This definition thus gives all the points in ( , ) giving within observational errors.
The confidence regions in Fig. 15 suggest that the LG infall towards PP is less than about 100 km s^{1}. In each case the best fit, plotted in the figure as a black triangle for and a white circle for , is obtained with . The best fit values of are coupled to the value of ; the correspondence with the range  , which was obtained in Sect. 3, is better for 0.4 than for or 1.
As mentioned above, the mass is varied by adjusting a bias parameter for the values given by Sigad et al. (1998). Let be the overdensity at different distances around the PP center, obtained using the contour plot in Fig. 6 of Sigad et al. (1998). Since Sigad et al. use in their figure, must be multiplied by for a general case of . Then we assume a bias of the form , where gives the density which we use in the TB calculations. Increasing would mean decreasing the value of , for example when , requires and means . The best fit values for and , and their formal confidence regions, derived from the regions in Fig. 15, are listed in Table 4.
Our conclusion from comparing the KLUNPP peculiar velocity data to the TBmodeled velocity field is that we get a good correspondence between the two when
Figure 15: The formal confidence regions (see text) for and . The best fit values are marked with triangles for and circles for . The dotted lines show the limits obtained in Sect. 3.  
Open with DEXTER 
The first results can be used to correct the observed velocities in our NDM distance determination. With the best fit TB model shown in Fig. 14 we calculate peculiar velocities at different distances r, which are used to correct the input kinematical distance scale. The whole process of binning and solving for PP distances and peculiar velocities is then repeated with these new kinematical distances as a reference scale.
A big disadvantage in this iterative procedure is that at each iteration we lose about half of our data points due to the constraints of Eq. (20) and . Figure 16 shows how we have only 6 binned data points containing 90 individual galaxies after the first iteration. Our conclusion is that the iterative process does not in this case improve the results.
Ekholm & Teerikorpi (1994) used in their study of the peculiar velocities in the Great Attractor region a socalled pclass analysis. There the possibility of the not being normally distributed is taken into account by taking only the galaxies lying in the unbiased plateau of the diagram. Applying the same procedure here did not give disagreeing results with the velocity patterns shown e.g. in Fig. 14, but, again, the number of points given by the pclass analysis was small.
Figure 16: Binned data points in v_{r} vs. r diagram after the iterative process for correcting distance determinations. Each time the PP distances are calculated, approximately half of the data are lost due to restrictions (Eq. (20)) and . Only six binned data points are now present, rather than the 13 shown in Fig. 14. In this figure we have adjusted the KLUN points with , and the curves show the best fitted TB solutions using . The best fit values are 2.9, 3.7, 5.5, and 9.4 , for 0.1, 0.2, 0.4, and 1.0, respectively, corresponding to bias values 0.45, 0.52, 0.53, and 0.55.  
Open with DEXTER 
We have studied the velocity field around the PerseusPisces supercluster using the KLUN galaxy sample and TolmanBondi modeling. In the TB calculations we used a smoothed density field of IRAS galaxies, as given by Sigad et al. (1998). There is a good correspondence with the peculiar velocities given by the KLUN data and the TB models, if the infall velocity of the LG towards the PP center is assumed to be less than 100 km s^{1} and if we add a bias in the Sigad et al. density field, . The resulting small infall velocity of the LG towards the PP region, km s^{1} disagrees with the km s^{1} of the earlier kinematical studies (see Sect. 1), but is consistent with the recent POTENT analysis (Dekel 1999).
The goodness of fit between the KLUN and TB velocities is not affected by varying the cosmological parameters and within reasonable limits. However, requiring the mass of the central region of PP to be in the limits obtained with virial calculations,  , we get constraints for the value of . Then the intermediate values 0.4 are preferred to the more extreme or . The value of , assumed to be either or 0, does not have a significant effect on the results. In fact, the infall velocities around growing structures being little influenced by the presence of is not so unexpected (Peebles 1984; Lahav et al. 1991).
Calculation of the virial masses of the most prominent clusters in the PP region resulted in cluster mass to luminosity ratios M/L = 200  . Combining this with a luminosity density of the universe W Hz^{1} Mpc^{3} (Lilly et al. 1996), and assuming that the M/L of the PP region is valid in the rest of the universe, we obtain 0.3. Allowing M/L to grow at larger scales would mean that this is a lower limit. The result agrees with recent estimates of . In comparison, in the study of the Hubble flow around the Virgo cluster, Teerikorpi et al. (1992) found that the constraint M/L= constant is inconsistent with or . A smaller was possible if M/L in Virgo is larger than in the field, contrary to the usual bias scenario. The same is valid here.
One of our aims was to study the applicability of the TB method in a complex environment like the PP supercluster. TB solutions are for spherically symmetric mass distributions only, and looking at the galaxy distribution in Fig. 2, PP hardly can be taken as such. However, a smoothed density map, like the IRAS map by Sigad et al. (1998), has a spherical structure. Similarly, the small scale dynamics are smoothed away, when we use an averaged velocity map, where the radial parts of the peculiar velocities are averaged over a shell around the main galaxy concentration. This is shown not only in the compatibility of the TB results and the KLUN data, but also in a study of an Nbody simulation. There it was confirmed that the radially averaged velocity field can be calculated by the TB method, using the observed smoothed density field. There were marginal differences between the TB solutions and the simulation when we looked at the peculiar velocities around large clusters at different epochs. However, at the present epoch these deviations are negligible. We argue that the TolmanBondi method can be used for a reliable deduction of a radially averaged velocity field around large galaxy concentrations.
Acknowledgements
We have made use of the LyonMeudon Extragalactic Database (LEDA) supplied by the LEDA team at the CRALObservatoire de Lyon (France). M.H.'s work has been supported by Finnish Cultural Foundation. We acknowledge the support by Academy of Finland projects "Cosmology in the local galaxy universe'' and "Galaxy streams and structures in the local universe''. The work was initiated while M.H. was visiting Osservatorio di Capodimonte in Naples. The staff in the observatory was very kind and warmhearted and deserves to be mentioned here. We thank Dr. Yu. Baryshev for his advice concerning the TB model, and the referee for useful suggestions concerning various parts of the article.