We consider the non-radial motions, originating in the outskirts of clusters of galaxies and we show how it may reduce the discrepancy between Cold Dark Matter (CDM) predicted X-ray temperature distribution function of clusters of galaxies and that really observed, so resolving the problem of the X-ray clusters abundance over-production predicted by the CDM model. We construct the X-ray temperature distribution function using Press-Schechter's (1974) theory and Evrard's (1990) prescriptions for the mass-temperature relation, taking also account of the non-radial motions originating from the gravitational interaction of the quadrupole moment of the protocluster. Then, we compare our X-ray temperature distribution function with the tidal field of the matter of the neighboring protostructures. We find that the model produces a reasonable cluster temperature distribution.

We calculate the density profiles of virialized halos in the case of structure evolving hierarchically from a scale-free Gaussian delta-field having a power spectrum $P(k)\propto k^n$ in a Omega=1 Universe; we suppose that the initial density contrast profile around local maxima is given by the mean peak profile introduced by Bardeen et al. (1986 hereafter BBKS). We show both that the density profiles are not power-laws but have a logarithmic slope that increases from the inner halo to its outer parts and for n \geq -1 are well approximated by Navarro et al. (1995, 1996, 1997) profile and the radius a, at which the slope alpha=-2, is a function of the mass of the halo and of the spectral index n.

In this paper, we study the role of shear fields on the evolution of density perturbations embedded in a Friedmann flat background universe, by studying the evolution of a homogeneous ellipsoid model. In this context, we show that while the effect of the shear is that of increasing the growth rate of the density contrast of a mass element, the angular momentum acquired by the ellipsoid has the right magnitude to counterbalance the shear. Finally, our result show that initial asphericities and tidal interaction induce a slowing down of the collapse after the system has broken away from the general expansion, in perfect agreement with the previrialization conjecture (Peebles & Groth 1976; Davis & Peebles 1977).

We study the effect of the interaction between small scale clumps having velocities larger than the average escape velocity v_{\rm esc} and bound clumps (having v<v_{\rm esc}), on clusters of galaxies evolution. We find that the collapse of the bound clumps is accelerated with respect to the homogeneous case (Gunn & Gott's model - Gunn & Gott 1972). A comparison between the acceleration produced by this effect to the slowing down effect studied in Del Popolo & Gambera (1998) shows that the magnitude of the slowing down effect is larger than the acceleration produced by the effect studyed in this paper, only in the outskirts of the cluster.

In the Cold Dark Matter (hereafter CDM) scenario even isolated density peaks contain a high fraction of small scale clumps having velocities larger than the average escape velocity from the structure. These clumps populate protoclusters, especially in the peripheral regions, r > R_f (where R_f is the filtering scale). During the cluster collapse and the subsequent secondary infall, collapsing or infalling clumps (having v<v_esc) interact with the quoted unbound clumps (or high speed clumps, as we also call them) having v > v_esc. We study the interaction between these two kind of clumps by means of the impulse approximation and we find that the collapse of bound clumps is accelerated with respect to the homogeneous case (Gunn & Gott's model). The acceleration of the collapse increases with decreasing height of the peak, \nu. We finally compare the acceleration produced by this effect to the slowing down effect produced by the gravitational interaction of the quadrupole moment of the system with the tidal field of the matter of the neighboring proto-clusters studied in Del Popolo & Gambera. We find that the magnitude of the slowing down effect is larger than the acceleration produced by the effect studied in this paper, only in the outskirts of the cluster. We want to stress that the one which we study in this paper is also present in an isolated protocluster, being produced by the interaction of the collapsing clumps with the unbound substructure internal to the collapsing clumps itself while that studied in Del Popolo & Gambera (1998) is produced by substructure external to the density peak.

We use the results of a high resolution N-body simulation to investigate the role of the environment on the formation and evolution of galaxy-sized halos. Starting from a set of constrained initial conditions, we have produced a final configuration hosting a double cluster in one octant and a large void extending over two octants of the simulation box. We present results for two statistics: the relationship between 1-D velocity dispersion and mass and the probability distribution of the spin parameter $P(\lambda)$. The \svm relationship is well reproduced by the Truncated Isothermal Sphere (TIS) model introduced by Shapiro et al. (1999), although the slope is different from the original prediction. A series of \svm relationships for different values of the anisotropy parameter $\beta$, obtained using the theoretical predictions by Lokas and Mamon (2001) for NFW density profiles are found to be only marginally consistent with the data. Using some properties of the equilibrium TIS models, we construct subsamples of {\em fiducial} equilibrium TIS halos from each of the three subregions, and we study their properties. For these halos, we do find an environmental dependence of their properties, in particular of the spin parameter distribution $P(\lambda)$. We study in more detail the TIS model, and we find new relationships between the truncation radius and other structural parameters. No gravitationally bound halo is found having a radius larger than the critical value for gravithermal instability for TIS halos (\rt $\ge 34.2 r_{0}$, where $r_{0}$ is the core radius of the TIS solution). We do however find a dependence of this relationship on the environment, like for the $P(\lambda)$ statistics. These facts hint at a possible r\^{o}le of tidal fields at determining the statistical properties of halos.

We associated, in a classical formulation of "strong gravity", hadron constituents with suitable stationary, axisymmetric solutions of some new Einstein-type equations supposed to describe the strong field inside hadrons. These new equations can be obtained by the Einstein equations with cosmological term Lambda. As a consequence, Lambda and the masses M result in our theory to be scaled up, and transformed into a "hadronic constant" and into "strong masses", respectively. Due to the unusual range of Lambda and M values considered, we met a series of solutions of the Kerr-Newman-de Sitter (hereafter KNdS) type with rather interesting properties. The requirement that those solutions be stable, i.e., that their temperature (or surface gravity) be vanishingly small, implies the coincidence of at least two of their (in general, three) horizons. Imposing the stability condition of a certain horizon does yield (once chosen the values of J, q and Lambda) mass and radius of the associated black-hole (hereafter BH). In the case of ordinary Einstein equations and for stable BHs of the KNdS type, we get in particular Regge-like (hereafter RL) relations among mass M, angular momentum J, charge q and cosmological constant Lambda; which did not receive enough attention in the previous literature. Besides, we show some particular and interesting cases of these relations. Another interesting point is that, with few exceptions, all such relations (among M, J, q, Lambda) lead to solutions that can be regarded as (stable) cosmological models.

We study the joint effect of tidal torques and dynamical friction on the collapse of density peaks solving numerically the equations of motion of a shell of barionic matter falling into the central regions of a cluster of galaxies. We calculate the evolution of the expansion parameter, a(t), of the perturbation using a coefficient of dynamical friction eta_{cl} obtained from a clustered system and taking into account the gravitational interaction of the quadrupole moment of the system with the tidal field of the matter of the neighboring proto-galaxies. We show that within high-density environments, such as rich clusters of galaxies, tidal torques and dynamical friction slow down the collapse of low-nu peaks producing an observable variation of the parameter of expansion of the shell. As a consequence a bias of dynamical nature arises because high-density peaks preferentially collapse to form halos within which visible objects eventually will condense. For a standard Cold Dark Matter model this dynamical bias can account for a substantial part of the total bias required by observations on cluster scales.

We study the effect of non-radial motions on the mass function, the VDF and on the shape of clusters of galaxies using the model introduced in Del Popolo & Gambera (1998a,b; 1999). The mass function of clusters, obtained using the quoted model, is compared with the statistical data by Bahcall & Cen (1992) and Girardi et al. (1998), while the VDF is compared with the CfA data by Zabludoff et al. (1993) for local clusters. In both cases the model predictions are in good agreement with the observational data showing once more how non-radial motions can reduce many of the discrepancies between CDM model predictions and observational data. Besides we study the effect of non-radial motions on the intrinsic shape of clusters of galaxies showing that non-radial motions produce clusters less elongated with respect to CDM model.

Planets orbiting a planetesimal circumstellar disc can migrate inward from their initial positions because of dynamical friction between planets and planetesimals. The migration rate depends on the disc mass and on its time evolution. Planets that are embedded in long-lived planetesimal discs, having total mass of $10^{-4}-0.01 M_{\odot}$, can migrate inward a large distance and can survive only if the inner disc is truncated or because of tidal interaction with the star. In this case the semi-major axis, a, of the planetary orbit is less than 0.1 AU. Orbits with larger $a$ are obtained for smaller value of the disc mass or for a rapid evolution (depletion) of the disc. This model may explain several of the orbital features of the giant planets that were discovered in last years orbiting nearby stars as well as the metallicity enhancement found in several stars associated with short-period planets.

N-body codes to perform simulations of the origin and evolution of the Large Scale Structure of the Universe have improved significantly over the past decade both in terms of the resolution achieved and of reduction of the CPU time. However, state-of-the-art N-body codes hardly allow one to deal with particle numbers larger than a few 10^7, even on the largest parallel systems. In order to allow simulations with larger resolution, we have first re-considered the grouping strategy as described in Barnes (1990) (hereafter B90) and applied it with some modifications to our WDSH-PT (Work and Data SHaring - Parallel Tree) code. In the first part of this paper we will give a short description of the code adopting the Barnes and Hut algorithm \cite{barh86} (hereafter BH), and in particular of the memory and work distribution strategy applied to describe the {\it data distribution} on a CC-NUMA machine like the CRAY-T3E system. In the second part of the paper we describe the modification to the Barnes grouping strategy we have devised to improve the performance of the WDSH-PT code. We will use the property that nearby particles have similar interaction list. This idea has been checked in B90, where an interaction list is builded which applies everywhere within a cell C_{group} containing a little number of particles N_{crit}. B90 reuses this interaction list for each particle $ p \in C_{group}$ in the cell in turn. We will assume each particle p to have the same interaction list. Thus it has been possible to reduce the CPU time increasing the performances. This leads us to run simulations with a large number of particles (N ~ 10^7/10^9) in non-prohibitive times.

We study the effect of non-radial motions on the mass function, the velocity dispersion function (hereafter VDF) and on the shape of clusters of galaxies using the model introduced in Del Popolo & Gambera (1998a,b; 1999). The mass function of clusters, obtained using the quoted model, is compared with the statistical data by Bahcall & Cen (1992a,b) and Girardi et al. (1998), while the VDF is compared with the Center for Astrophysics (hereafter CfA) data by Zabludoff et al. (1993) for local clusters and those of Mazure et al. (1996) and Fadda et al. (1996). In both cases the model predictions are in good agreement with the observational data showing once more how non-radial motions can reduce many of the discrepancies between Cold Dark Matter (hereafter CDM) model predictions and observational data. Finally we study the effect of non-radial motions on the intrinsic shape of clusters of galaxies showing that non-radial motions produce clusters less elongated with respect to CDM model in agreement with de Theije et al. (1995, 1997) results.

We present and discuss a method to identify substructures in combined angular-redshift samples of galaxies within Clusters. The method relies on the use of Discrete Wavelet Transform (hereafter DWT) and has already been applied to the analysis of the Coma cluster (Gambera et al. 1997). The main new ingredient of our method with respect to previous studies lies in the fact that we make use of a 3D data set rather than a 2D. We test the method on mock cluster catalogs with spatially localized substructures and on a N-body simulation. Our main conclusion is that our method is able to identify the existing substructures provided that: a) the subclumps are detached in part or all of the phase space, b) one has a statistically significant number of redshifts, increasing as the distance decreases due to redshift distortions; c) one knows {\it a priori} the scale on which substructures are to be expected. We have found that to allow an accurate recovery we must have both a significant number of galaxies ($\approx 200$ for clusters at z$\geq 0.4$ or about 800 at z$\leq$ 0.4) and a limiting magnitude for completeness $m_B=16$. The only true limitation to our method seems to be the necessity of knowing {\it a priori} the scale on which the substructure is to be found. This is an intrinsic drawback of the method and no improvement in numerical codes based on this technique could make up for it.

In this paper we calculate the density profiles of virialized halos both in the case of structure evolving hierarchically from a scale-free Gaussian delta-field having a power spectrum P(k)=A k^n in a Omega=1 Universe and in the case of the CDM model, by using a modified version of Hoffman & Shaham's (1985) (hereafter HS) and Hoffman's (1988) model. We suppose that the initial density contrast profile around local maxima is given by the mean peak profile introduced by Bardeen et al. (1986) (hereafter BBKS), and is not just proportional to the two-point correlation function, as assumed by HS. We show that the density profiles, both for scale-free Universes and the CDM model, are not power-laws but have a logarithmic slope that increases from the inner halo to its outer parts. Both scale-free, for n >=-1, and CDM density profiles are well approximated by Navarro et al. (1995, 1996, 1997) profile. The radius a, at which the slope alpha=-2, is a function of the mass of the halo and in the scale-free models also of the spectral index n.

We study the effect of non-radial motions, originating from the gravitational interaction of the quadrupole moment of a protogalaxy with the tidal field of the matter of the neighboring protostructures, on the angular correlation function of galaxies. We calculate the angular correlation function using a Standard Cold Dark Matter (hereafter SCDM) model (Omega=1, h=0.5, n=1) and we compare it with the angular correlation function of the APM galaxy survey (Maddox et al. 1990; Maddox et al. 1996). We find that taking account of non-radial motions in the calculation of the angular correlation function gives a better agreement of the theoretical prediction of the SCDM model to the observed estimates of large-scale power in the galaxy distribution.

In this paper we study the role of dynamical friction on the evolution of a population of large objects ($m>10^{22}$ g) at heliocentric distances $>70$ AU in the Kuiper Belt. We show that the already flat distribution of these objects must flatten further due to non-spherically symmetric distribution of matter in the Kuiper Belt. Moreover the dynamical drag, produced by dynamical friction, causes objects of masses $\geq 10^{24} g$ to lose angular momentum and to fall through more central regions in a timescale $\approx 10^9 yr$. This mechanism is able to transport inwards objects of the size of Pluto, supposing it was created beyond 50 AU, according to a Stern & Colwell's (1997b) suggestion.

In this paper we show how non-radial motions, originating from the tidal interaction of the irregular mass distribution within and around protoclusters, can solve some of the problems of the CDM model. Firstly the discrepancy between the CDM predicted two-points correlation function of clusters and the observed one. We compare the two-points correlation function, that we obtain taking account of non-radial motions, with that obtained by Sutherland & Efstathiou (1991) from the analysis of Geller & Hucra's (1988) deep redshift survey and with the data points for the APM clusters obtained by Efstathiou et al. (1992). Secondly the problem of the X-ray clusters abundance over-production predicted by the CDM model. In this case we compare the X-ray temperature distribution function, calculated using Press-Schechter theory and Evrard's (1990) prescriptions for the mass-temperature relation, taking also account of the non-radial motions, with Henry & Arnaud (1991) and Edge et al. (1990) X-ray temperature distributions for local clusters. We find that in both cases the model is in good agreement with experimental data. Finally we calculate the bias coefficient using a selection function that takes into account the effects of non-radial motions, and we show that the bias so obtained can account for a substantial part of the total bias required by observations on cluster scales.

In this paper we extend Chandrasekhar and von Neumann's analysis of the statistics of the gravitational field to systems in which particles (e.g. stars, galaxies) are not homogeneously distributed. We derive a distribution function W({F},dF/dt) giving the joint probability that a test particle is subject to a force F and an associated rate of change of F given by dF/dt. We calculate the first moment of dF/dt to study the effects of inhomogenity on dynamical friction.

We show how non-radial motions, originating in the outskirts of clusters of galaxies, may reduce the discrepancy between the Cold Dark Matter (CDM) predicted X-ray temperature distribution function of clusters of galaxies and the observed one and also the discrepancy between the CDM predicted two-point correlation function of clusters of galaxies and that observed. We compare Edge et al. (1990) and Henry & Arnaud (1991) data with the distribution function of X-ray temperature, calculated using Press- Schechter's (1974 - hereafter PS) theory and Evrard's (1990) prescriptions for the mass-temperature relation and taking account of the non-radial motions originating from the gravitational interaction of the quadrupole moment of the protocluster with the tidal field of the matter of the neighboring protostructures. We find that the model produces a reasonable clusters temperature distribution. We compare the two-point cluster correlation function which takes account of the non-radial motions both with that obtained by Sutherland & Efstathiou (1991), from the analysis of Huchra's et al. (1990) deep redshift survey, and with the data points for the Automatic Plate Measuring (APM) clusters, computed by Efstathiou et al. (1992a), showing that non-radial motions reduce the discrepancy between the theoretical and the obsservational correlation function.

We study the effect of tidal torques on the collapse of density peaks through the equations of motion of a shell of barionic matter falling into the central regions of a cluster of galaxies. We calculate the time of collapse of the perturbation taking into account the gravitational interaction of the quadrupole moment of the system with the tidal field of the matter of the neighbouring proto-clusters. We show that within high-density environments, such as rich clusters of galaxies, tidal torques slow down the collapse of low \nu peaks producing an observable variation in the time of collapse of the shell and, as a consequence, a reduction in the mass bound to the collapsed perturbation. Moreover, the delay of the collapse produces a tendency for less dense regions to accrete less mass, with respect to a classical spherical model, inducing a biasing of over-dense regions toward higher mass. Finally we calculate the bias coefficient using a selection function properly defined showing that for a Standard Cold Dark Matter (SCDM) model this bias can account for a substantial part of the total bias required by observations on cluster scales.

N-body algorithms for long-range unscreened interactions like gravity belong to a class of highly irregular problems whose optimal solution is a challenging task for present-day massively parallel computers. In this paper we describe a strategy for optimal memory and work distribution which we have applied to our parallel implementation of the Barnes & Hut (1986) recursive tree scheme on a Cray T3D using the CRAFT programming environment. We have performed a series of tests to find an " optimal data distribution " in the T3D memory, and to identify a strategy for the " Dynamic Load Balance " in order to obtain good performances when running large simulations (more than 10 million particles). The results of tests show that the step duration depends on two main factors: the data locality and the T3D network contention. Increasing data locality we are able to minimize the step duration if the closest bodies (direct interaction) tend to be located in the same PE local memory (contiguous block subdivison, high granularity), whereas the tree properties have a fine grain distribution. In a very large simulation, due to network contention, an unbalanced load arises. To remedy this we have devised an automatic work redistribution mechanism which provided a good Dynamic Load Balance at the price of an insignificant overhead.

By means of numerical methods, we solve the equations of motion for the collapse of a shell of baryonic matter, made of galaxies and substructure falling into the central regions of a cluster of galaxies, taking into account the effect of the dynamical friction. The parameters on which the dynamical friction mainly depends are: the peaks' height, the number of peaks inside a protocluster multiplied by the correlation function evaluated at the origin, the filtering radius and the nucleus radius of the protocluster of galaxies. We show how the collapse time (Tau) of the shell depends on these parameters. We give a formula that links the dynamical friction coefficient (Eta) o the parameters mentioned above and an analytic relation between the collapse time and (Eta). Finally, we obtain an analytical relation between (Tau) and the mean overdensity (mean Delta) within the shell. All the analytical relations that we find are in excellent agreement with the numerical integration.

Evidence for clustering within the Coma cluster is found by means of a multiscale analysis of the combined angular-redshift distribution. We have compiled a catalogue of 798 galaxy redshifts from published surveys from the region of the Coma cluster. We examine the presence of substructure and of voids at different scales ranging from $\sim 1 to \sim 16 h^{-1}$ Mpc, using subsamples of the catalogue, ranging from $cz=3000$ km/s to $cz=28000$ km/s. Our substructure detection method is based on the wavelet transform and on the segmentation analysis. The wavelet transform allows us to find out structures at different scales and the segmentation method allows us a quantitative statistical and morphological analysis of the sample. From the whole catalogue we select a subset of 320 galaxies, with redshifts between cz=5858 km/s and cz=8168 km/s that we identify as belonging to the central region of Coma and on which we have performed a deeper analysis, on scales ranging from $180 h^{-1}$ kpc to $1.44 h^{-1}$ Mpc. Our results are expressed in terms of the number of structures or voids and their sphericity for different values of the threshold detection and at all the scales investigated. According to our analysis, there is strong evidence for multiple hierarchical substructure, on scales ranging from a few hundreds of kpc to about $4 h^{-1}$ Mpc. The morphology of these substructures is rather spherical. On the scale of $720 h^{-1}$ kpc we find two main subclusters which where also found before, but our wavelet analysis shows even more substructures, whose redshift position is approximatively marked by these bright galaxies: NGC 4934 & 4840, 4889, 4898 & 4864, 4874 & 4839, 4927, 4875.

We solve numerically the equations of motion for the collapse of a shell of baryonic matter, made of galaxies and substructure of $ 10^{6} M_{\odot} \div 10^{9} M_{\odot}$, taking into account the dynamical friction and the parameters on which it depends: the peaks' height $\nu_{c}$, the number of peaks inside a protocluster $N_{tot}$ multiplied by the correlation function $\xi (r)$ evaluated at r=0, the filtering radius $R_{f}$ and the core radius of the cluster of galaxies, r_0. We show how the collapse time of a shell of baryonic matter depends on those parameters.

We solve numerically the equations of motion for the collapse of a shell of baryonic matter falling into the central regions of a cluster of galaxies, taking into account of the presence of the substructure inducing dynamical friction. The evolution of the expansion parameter a(t) of the perturbation is calculated in spherical systems. The effect of dynamical friction is to reduce the binding radius and the total mass accreted by the central regions. Using a peak density profile given by Bardeen et al. (1986) we show how the binding radius of the perturbation is modified by dinamical friction. We show how dynamical friction modifies the collapse parameter of the perturbation slowing down the collapse.

We solve numerically the equations of motion for the collapse of a shell of baryonic matter, made of galaxies and substructure of $ 10^{6} M_{\odot}- 10^{9} M_{\odot}$, falling into the central regions of a cluster of galaxies taking account of dynamical friction. We calculate the evolution of the expansion parameter, a(t), of the perturbation using a coefficient of dynamical friction, $ \eta_{0}$, calculated for a perturbation in which clustering is absent and a coefficient $ \eta_{cl}$ obtained from a clustered one. The effect of the dynamical friction is to slow down the collapse (V. Antonuccio-Delogu \& S. Colafrancesco 1994, hereafter AC) producing an observable variation of the parameter of expansion of the shell. The effect increases with increasing $ \eta $ and with the increasing of clustering. Finally, we show how the collapse time depends on $ \eta_{0}$ and $ \eta_{cl}$. keywords: {cosmology: theory- galaxies: formation}

We show how the results given by several authors relatively to the mass of a density peak are changed when small scale substructure induced by dynamical friction are taken into account. The peak mass obtained is compared to the result of Peacock \& Heavens (1990) and to the peak mass when dynamical friction is absent to show how these effects conspire to reduce the mass accreted by the peak.