Galaxy interactions can have an important effect in a galaxy's evolution. Cosmological models predict a large number of small satellites around galaxies. It is important to study the effect that these small satellites can have on the host. The present work explores the effect of small N-body spherical satellites with total mass ratios in the range approx 1:1000-1:100 in inducing gas flows to the central regions of a disc galaxy with late-type morphology resembling the Milky Way. Two model galaxies are considered: barred and non-barred models; the latter one is motivated in order to isolate and understand better the effects of the satellite. Several circular and non-circular orbits are explored, considering both prograde and retrogade orientations. We show that satellites with such small mass ratios can still produce observable distortions in the gas and stellar components of the galaxy. In terms of gas flows, the prograde circular orbits are more favourable for producing gas flows, where in some cases up to $60% of the gas of the galaxy is driven to the central region. We find, hence, that small satellites can induce significant gas flows to the central regions of a disc galaxy, which is relevant in the context of fuelling active galactic nuclei.

We present the clustering properties and halo occupation distribution (HOD) modelling of very low redshift, hard X-ray-detected active galactic nuclei (AGN) using cross-correlation function measurements with Two-Micron All Sky Survey galaxies. Spanning a redshift range of $0.007 < z < 0.037$, with a median $z=0.024$, we present a precise AGN clustering study of the most local AGN in the Universe. The AGN sample is drawn from the SWIFT/BAT 70-month and INTEGRAL/IBIS eight year all-sky X-ray surveys and contains both type I and type II AGN. We find a large-scale bias for the full AGN sample of $b=1.04^{+0.10}_{-0.11}$, which corresponds to a typical host dark matter halo mass of $M_{\rm h}^{\rm typ}=12.84^{+0.22}_{-0.30}\,h^{-1} M_{\odot}$. When split into low and high X-ray luminosity and type I and type II AGN subsamples, we detect no statistically significant differences in the large-scale bias parameters. However, there are differences in the small-scale clustering which are reflected in the full HOD model results. We find that low and high X-ray luminosity AGN, as well as type I and type II AGN, occupy dark matter haloes differently, with 3.4$\sigma$ and 4.0$\sigma$ differences in their mean halo masses, respectively, when split by luminosity and type. The latter finding contradicts a simple orientation-based AGN unification model. As a by-product of our cross-correlation approach, we also present the first HOD model of 2MASS galaxies.

(Abridged version) We explore whether a scenario that combines an origin by mergers at $z\sim$1.8-1.5 with a subsequent passive evolution of the resulting S0 remnants since $z \sim$0.8-1 is compatible with observational data of S0s in the Tully-Fisher relation (TFR). We studied a set of major and minor merger experiments from the GalMer database that generate massive S0 remnants. We analysed the location of these remnants in the photometric and stellar TFRs assuming that they correspond to $z\sim0.8$ galaxies. We then estimated their evolution in these planes over the last 7 Gyr. The results were compared with data of real S0s and spirals at different redshifts. We also tested how the use of Vcirc or Vrot,max affects the results. We found that just after $\sim$1-2 Gyr of coalescence, major mergers generate S0 remnants that are outliers of the local photometric and stellar TFRs at $z\sim0.8$. After $\sim$4-7 Gyr of passive evolution in isolation, the S0 remnants move towards the local TFR, although the initial scatter among them persists. This scatter is sensitive to the indicator used for the rotation velocity: Vcirc values yield a lower scatter than when Vrot,max values are considered instead. In the planes involving Vrot,max, a clear segregation of the S0 remnants in terms of the spin-orbit coupling of the model is observed, in which the remnants of retrograde encounters overlap with local S0s hosting counter-rotating discs. The location of the S0 remnants at $z\sim 0$ agrees well with the observed distribution of local S0 galaxies in the $\sigma_0$-$M_K$, Vcirc-$\sigma_0$ and Vrot,max-$\sigma_0$ planes. Thus, massive S0 galaxies may have been formed through major mergers that occurred at high redshift and have later evolved towards the local TFR through passive evolution in relative isolation, a mechanism that would also contribute to the scatter observed in this relation.

We study some effects the inclusion of the radiation energy component in the universe, Omega_r, can have on several quantities of interest for the large-scale structure of the universe in a LCDM cosmological simulation; started at a very high redshift (z=500). In particular we compute the power spectrum density, the halo mass function, and the concentration-mass relation for haloes. We find that Omega_r has an important contribution in the long-term nonlinear evolution of structures in the universe. For instance, a lower matter density power, by approx 50%, in all scales is obtained when compared with a simulation without the radiation term. Also, haloes formed with the Omega_r taken into account are approx 20% less concentrated than when not included in the Hubble function.

This is the fourth paper in a series that reports on our investigation of the clustering properties of active galactic nuclei (AGN) identified in the ROSAT All-Sky Survey (RASS) and Sloan Digital Sky Survey (SDSS). In this paper we investigate the cause of the X-ray luminosity dependence of the clustering of broad-line, luminous AGN at 0.16<z<0.36. We fit the H-alpha line profile in the SDSS spectra for all X-ray and optically-selected broad-line AGN, determine the mass of the super-massive black hole (SMBH), M_BH, and infer the accretion rate relative to Eddington (L/L_EDD). Since M_BH and L/L_EDD are correlated, we create AGN subsamples in one parameter while maintaining the same distribution in the other parameter. In both the X-ray and optically-selected AGN samples we detect a weak clustering dependence with M_BH and no statistically significant dependence on L/L_EDD. We find a difference of up to 2.7sigma when comparing the objects that belong to the 30% least and 30% most massive M_BH subsamples, in that luminous broad-line AGN with more massive black holes reside in more massive parent dark matter halos at these redshifts. These results provide evidence that higher accretion rates in AGN do not necessarily require dense galaxy environments in which more galaxy mergers and interactions are expected to channel large amounts of gas onto the SMBH. We also present semi-analytic models which predict a positive M_DMH dependence on M_BH, which is most prominent at M_BH ~ 10^{8-9} M_SUN.

Motivated by recent inferred form of the halo occupation distribution (HOD) of X-ray selected AGNs, in the COSMOS field by Allevato et al. (2012), we investigate the HOD properties of moderate X-ray luminosity Active Galactic Nuclei (mXAGNs) using a simple model based on merging activity between dark matter halos (DMHs) in a $\Lambda$-CDM cosmology. The HODs and number densities of the simulated mXAGNs at $z=0.5$, under the above scenarios to compare with Allevato et al. (2012) results. We find that the simulated HODs of major and minor mergers, and the observed for mXAGNs are consistent among them. Our main result is that minor mergers, contrary to what one might expect, can play an important role in activity mAGNs.

To date, 17 circumbinary planets have been discovered. In this paper, we focus our attention on the stability of the Kepler circumbinary planetary systems with only one planet, i.e. Kepler-16, Kepler-34, Kepler-35, Kepler-38, Kepler-64 and Kepler-413. In addition to their intrinsic interest, the study of such systems is an opportunity to test our understanding of planetary system formation and evolution around binaries. The investigation is done by means of numerical simulations. We perform numerical integrations of the full equations of motion of each system with the aim of checking the stability of the planetary orbit. The investigation of the stability of the above systems consists of three numerical experiments. In the first one we perform a long term (1Gyr) numerical integration of the nominal solution of the six Kepler systems under investigation. In the second experiment, we look for the critical semimajor axis of the six planetary orbits, and finally, in the third experiment, we construct two dimensional stability maps on the eccentricity-pericentre distance plane. Additionally, using numerical integrations of the nominal solutions we checked if this solutions were close to the exact resonance.

We study the distribution of intragroup dark matter in small groups of dark matter halos of galaxy-like size in a LCDM cosmology. These groups are identified using a physical criterion and may be an appropriate representation of small galaxy groups. We quantify the amount of intragroup dark matter and characterize its distribution. We find that compact associations of halos, as well as those intermediate and loose groups, have rather flat intragroup dark matter profiles with logarithmic slopes of gamma ~0 and ~-0.2$, respectively. Hence, the intragroup dark matter of these halo systems does not follow the same cuspy tendency that halos of galaxies have. In intermediate and loose galaxy-size halo associations the intragroup matter tends to be <50% that of the total mass of the group, and in compact associations is <20% within their group radius.

We study the effect of a massive planetesimal disk on the dynamical stability of the outer planets assuming, as has been suggested recently, that these were initially locked in a compact and multiresonant configuration as a result of gas-driven migration in a protoplanetary disk. The gravitational interaction among all bodies in our simulations is included self-consistently using the Mercury6.5 code. Several initial multiresonant configurations and planetesimal disk models are considered. Under such conditions a strong dynamical instability, manifested as a rapid giant planet migration and planetesimal disk dispersal, develops on a timescale of less than 40 Myr in most cases. Dynamical disk heating due to the gravitational interactions among planetesimals leads to more frequent interactions between the planetesimals and the ice giants Uranus and Neptune, in comparison to models in which planetesitmal-planetesimal interactions are neglected. On account of the rapid evolution of the multiresonant configurations obtained with fully self-consistent simulations, our results are inconsistent with the dynamical instability origin of the Late Heavy Bombardment as currently considered by the Nice model for the Solar System.

This is the third paper in a series that reports on our investigation of the clustering properties of AGNs identified in the ROSAT All-Sky Survey (RASS) and Sloan Digital Sky Survey (SDSS). In this paper, we extend the redshift range to 0.07<z<0.50 and measure the clustering amplitudes of both X-ray and optically-selected SDSS broad-line AGNs with and without radio detections as well as for X-ray selected narrow-line RASS/SDSS AGNs. We measure the clustering amplitude through cross-correlation functions (CCFs) with SDSS galaxies and derive the bias by applying a halo occupation distribution (HOD) model directly to the CCFs. We find no statistically convincing difference in the clustering of X-ray and optically-selected broad-line AGNs, as well as with samples in which radio-detected AGNs are excluded. This is in contrast to low redshift optically-selected narrow-line AGNs, where radio-loud AGNs are found in more massive halos than optical AGNs without a radio-detection. The typical dark matter halo masses of our broad-line AGNs are log M_DMH/[h^(-1) M_SUN] ~ 12.4-13.4, consistent with the halo mass range of typical non-AGN galaxies at low redshifts. We find no significant difference between the clustering of X-ray selected narrow-line AGNs and broad-line AGNs. We confirm the weak dependence of the clustering strength on AGN X-ray luminosity at a ~2 sigma level. Finally, we summarize the current picture of AGN clustering to z~1.5 based on three dimensional clustering measurements.

It has been suggested that the ejection to interplanetary space of terrestrial crustal material, accelerated in a large impact, may result in the interchange of biological material between Earth and other Solar System bodies. In this paper, we analyze the fate of debris ejected from Earth by means of numerical simulations of the dynamics of a large collection of test particles. This allows us to determine the probability and conditions for the collision of ejecta with other planets of the Solar System. We also estimate the amount of particles falling-back to Earth as a function of time after being ejected. We find that, in general, the collision rates of Earth ejecta with Venus and the Moon, as well as the fall-back rates, are consistent with results reported in the literature. By considering a larger number of particles than in all previous calculations we have also determined directly the collision probability with Mars and, for the first time, computed collision probabilities with Jupiter. We find that the collision probability with Mars is greater than values determined from collision cross section estimations previously reported.

We study some of the effects that magnetic turbulent fluctuations have on the dynamics of pickup O+ ions in the magnetic polar regions of the Mars ionosheath. In particular we study their effect on the bulk velocity profiles of ions as a function of altitude over the magnetic poles, in order to compare them with recent Mars Express data; that indicate that their average velocity is very low and essentially in the anti-sunward direction. We find that, while magnetic field fluctuations do give rise to deviations from simple ExB-drift gyromotion, even fluctuation amplitudes much greater than those of in situ measurements are {\it not} able to reproduce the vertical velocity profile of O+ ions. We conclude that other physical mechanisms, different from a pure charged particle dynamics, are acting on pickup ions at the Martian terminator. A possibility being a viscous-like interaction between the Solar Wind and the Martian ionosphere at low altitudes.

We study the dynamics of the upsilon-Andromedae planetary system proposed by Curiel et al. (2011). We focus on the study of the 3:1 Mean Motion Resonance between upsilon-Andromedae-d and the recently discovered upsilon-Andromedae-e (hereafter upsilon-And-d and upsilon-And-e). Numerical simulations of the dynamics of the four planet system are conducted. The previously reported apsidal resonance between upsilon-And-c and upsilon-And-d is confirmed. In addition, we find that upsilon-And-d and upsilon-And-e are also in an apsidal resonance condition. Our results further indicate that the upsilon-Andromedae planetary system configuration is in the middle of a stability island in the semimajor axis-eccentricity domain. Additionally, we performed numerical integrations of the planetary configuration over 500 Myr and found it to be stable. We conclude that, within the uncertainties in the value of the orbital parameters, it is likely that upsilon-Andromedae planetary system will remain stable for a long timescale.

This is the second paper of a series that reports on our investigation of the clustering properties of AGNs in the ROSAT All-Sky Survey (RASS) through cross-correlation functions (CCFs) with Sloan Digital Sky Survey (SDSS) galaxies. In this paper, we apply the Halo Occupation Distribution (HOD) model to the CCFs between the RASS Broad-line AGNs with SDSS Luminous Red Galaxies (LRGs) in the redshift range 0.16<z<0.36 that was calculated in paper I. In our HOD modeling approach, we use the known HOD of LRGs and constrain the HOD of the AGNs by a model fit to the CCF. For the first time, we are able to go beyond quoting merely a `typical' AGN host halo mass, M_h, and model the full distribution function of AGN host dark matter halos. In addition, we are able to determine the large-scale bias and the mean M_h more accurately. We explore the behavior of three simple HOD models. Our first model (Model A) is a truncated power-law HOD model in which all AGNs are satellites. With this model, we find an upper limit to the slope (\alpha) of the AGN HOD that is far below unity. The other two models have a central component, which has a step function form, where the HOD is constant above a minimum mass, without (Model B) or with (Model C) an upper mass cutoff, in addition to the truncated power-law satellite component, similar to the HOD that is found for galaxies. In these two models we find the upper limits of \alpha < 0.95 and \alpha < 0.84 for Model B and C respectively. Our analysis suggests that the satellite AGN occupation increases slower than, or may even decrease with, M_h, in contrast to the satellite's HODs of luminosity-threshold samples of galaxies, which, in contrast, grow approximately as \propto M_h^\alpha with \alpha\approx 1. These results are consistent with observations that the AGN fraction in groups and clusters decreases with richness.

We model the interaction of the solar wind with the plasma tail of a comet by means of numerical simulations, taking into account the effects of viscous-like forces.A 2D hydrodynamical, two species, finite difference code has been developed for the solution of the time dependent continuity, momentum and energy conservation equations, as applied to the problem at hand. Velocity, density and temperature profiles across the tail are obtained. Several cases with different flow parameters are considered in order to study the relative importance of viscous-like effects and the coupling between species on the flow dynamics. Assuming a Mach number equal to 2 for the incident solar wind as it flows past the comet nucleus, the flow exhibits three transitions with location and properties depending on the Reynolds number for each species and on the ratio of the timescale for inter-species coupling to the crossing time of the free flowing solar wind. By comparing our results with the measurements taken in situ by the Giotto spacecraft during its flyby of comet Halley we constrain the flow parameters for both plasmas. In the context of our approximations, we find that our model is qualitatively consistent with the in situ measurements as long as the Reynolds number of the solar wind protons and of cometary H2O+ ions is low, less than 100, suggesting that viscous-like momentum transport processes may play an important role in the interaction of the solar wind and the plasma environment of comets.

We determine both representations of the Fundamental Plane [Re Sig^a Ie^b and Re (Sig^2/Ie)^lambda] and the luminosity-effective phase space density [L fe^(-gamma)] scaling relation for N-body remnants of binary mergers of spiral-like galaxies. The main set of merger simulations involves a mass-ratio of the progenitors in the range of about 1:1 to 1:5, harboring or not a bulge-like component, and are constructed using a cosmological motivated model. Equal-mass mergers are also considered. Remnants lead to average values for the scaling indices of <a>~1.6, <b>~0.6, <lambda>~0.7, and <gamma>~0.65. These values are consistent with those of K--band observations (Mobasher et al. 1999) of ellipticals: <a>~1.5, <b>~0.8, <lambda>~0.7, and <gamma>~0.60. The b index is, however, not well reproduced. This study does not allow us to establish a conclusive preference for models with or without a bulge as progenitors. Our results indicate that the L--fe and FP scalings might be determined to a large extent by dissipationless processes, a result that appears to be in contradiction to other dissipationless results.

Sersic parameters characterising the density profiles of remnants formed in collision-less disc galaxy mergers are obtained; no bulge is included in our simulations. For the luminous component we find that the Sersic index is n ~ (1.5,5.3) with <n> ~ 3 +/- 1 and an effective radius of R_e ~ (1.6,12.9) kpc with <R_e> ~ 5 +/- 3 kpc. A strong correlation of n with the central projected density I_0 is found [n ~ I_0^(-0.14)] which is consistent with observations. No positive linear correlation between the size (R_e) and structure (n) of our remnants is found; we do not advocate the existence of this. The photometric plane (PHP) of the luminous component [n ~ R_e^(0.05) I_0^(0.15)] agrees well, within the uncertainties and the assumption of a constant mass-to-light ratio, with those observationally determined for ellipticals. We found that the surface defined by Sersic parameters (n, R_e, mu_0) in log-space is not a true plane, but a pseudo-plane with a small curvature at low values of n owed to intrinsic properties of the Sersic model. The dark haloes of the remnants have a 3-dimensional Sersic index of <n> ~ 4 +/- 0.5 that are smaller than the ones obtained for dark haloes in LCDM cosmologies <n> ~ 6 +/- 1. A tight dark Sersic ``plane'' (DSP) is also defined by the parameters of the remnants haloes with n \~ r_e^(0.07) rho_0^(0.10). We conclude that collision-less merger remnants of pure disc galaxies have Sersic properties and correlations consistent with those of observed in early-type galaxies and local remnants. It seems that a ``primordial'' bulge in spirals is not a necessary condition to form bona fide ellipticals on grounds of the Sersic properties of remnants.

We carried out three self-consistent $N$-body simulations of galaxy mergers, with a cuspy dark matter profile, in order to study if the inner dark density slope is preserved in the remnants. In these simulations the progenitors include both a stellar disk and an intrinsic angular momentum for the halos, unlike previous similar studies. The mass-ratios of progenitor galaxies are about 1:1, 1:3 and 1:10. We find that the initial cuspy density profile of the dark halos is preserved in the remnants for the cases considered here.

We present scaling relations between the virial velocity (V) and the one-dimensional central velocity dispersion (Sig0); the gravitational radius (Rv) and the effective radius (Re); and the total mass (M) and the luminous mass (ML) found in N-body simulations of binary mergers of spiral galaxies. These scalings are of the form V^2 ~ Sig0^alpha, Rv ~ Re^beta and M ~ ML^gamma. The particlar values obtained for (alpha,beta,gamma) depend on the method of fitting used [ordinary least-squares (ols) or orthogonal distance regression (odr)], the assumed profile [de Vaucouleurs (deV) or Sersic (S)], and the size of the radial interval where the fit is done. The alpha and gamma indexes turn out more sensitive to the fitting procedure, obtaining for the ols a mean alpha_ols=1.51 and gamma_ols=0.69, while for the odr alpha_odr=2.35 and gamma_odr=0.76. The beta index depends more on the adopted type of profile, with beta_deV=0.13 and beta_S=0.27. We conclude that dissipationless formed remnants of mergers have a strong breaking of structural and kinematical homology.

We set up a series of self-consistent N-body simulations to investigate the fundamental plane of merger remnants of spiral galaxies. These last ones are obtained from a theoretical Tully-Fisher relation at z=1, assuming a constant mass-to-light ratio within the LambdaCDM cosmogony. Using a Sersic growth curve and an orthogonal fitting method, we found that the fundamental plane of our merger remnants is described by the relation Re ~ sigma^{1.48} Ie^{-0.75} which is in good agreement with that reported from the Sloan Digital Sky Survey Re ~ sigma^{1.49} Ie^{-0.75}. However, the R^{1/4}-profile leads to a fundamental plane given by Re ~ sigma^{1.79} Ie^{-0.60}. In general, the correlation found in our merger remnants arises from homology breaking (V^2 ~ sigma^nu, Rg ~ Re^eta) in combination with a mass scaling relation between the total and luminous mass, $M ~ ML^gamma. Considering an orthogonal fitting method, it is found that 1.74<nu<1.79, 0.21<eta<0.52 and 0.80<gamma<0.90 depending on the adopted profile (Sersic or R^{1/4}).

A series of N-body simulations aimed to study the dynamics of small groups of galaxies are presented. The properties of small galaxy groups are very well reproduced, and those of Hickson's compact groups are well reproduced by the most advanced stage of collapsing groups. We find no overmerging problem in our simulations. An important fraction of groups (~40%) initially in virial equilibrium can last for ~10 Gyr without complete merging. These results provide an alternative solution to the overmerging expected in Hicksons's compact groups.

By means of N-body simulations we study the global dynamics of triplets of galaxies, considering initial conditions starting from `maximum expansion' and in virial equilibrium. Unlike previous studies we treat galaxies self-consistently. We do not consider the influence of a primordial common halo of dark matter. Our results indicate that a low number of triple mergers are expected at the present epoch (~10%) for collapsing triplets. Initially virialised conditions yield ~5% of triple mergers in ~10 Gyr of evolution; hence, the 3-galaxy problem has stable states. No overmerging problem for these small groups of galaxies is found. Their geometrical properties, as reflected by the Agekyan-Anosova map, do not show an excess of extreme hierarchical triplets. The median velocity dispersion of observed compact triplets (~100 km/s) is not well reproduced in our models at the present epoch. However, about 10% of simulated triplets starting from maximum expansion have dynamical properties very similar to them. Our median values agree, however, very well with new data on triples. We find that the median mass estimators performs somewhat better than the virial mass estimate. The VME underestimates mass by \~35%. We conclude that Karachentsev's compact triplets probably represent the most advanced stage of gravitational clustering of initially diffuse triplets.

We present a study of the dynamics of triplets of galaxies. In particular we focus on properties that are directly comparable with observations: e.g., velocity dispersion, mean harmonic radius, and virial mass. We find, for example, that the median velocity dispersion of our simulations is smaller (<50 km/s) than the observed value (~100 km/s). A possible solution to this discrepancy is to consider triplets as non-isolated systems; e.g. subject to a tidal perturbation from a poor cluster in a Hubble time. We suggest that the dynamics of triplets is intimately linked to the large-scale structure, and evolution, of the universe.

Optical red spectra of a set of 18 bright barred spiral galaxies are presented. The study is aimed at determining the local kinematics, and the physical conditions of ionized gas in the compact nucleus (inside a diameter of 5") and in the circumnuclear regions (inside a diameter of 20"). Only 8 galaxies showed bright emission from their east and west side of the nucleus. The spectrum of each region was analized separately. In other 10 galaxies the line emission was so weak that we were only able to obtain an average spectrum of the central emission. No emission was detected in the remaining 8 galaxies. An estimate of the dynamical mass is presented based on the observed velocities in the circumnuclear regions. In NGC 4314 and NGC 6951, that show H_alpha emission distributed in circumnuclear ring structures, we determine the [NII]/H_alpha and [SII]/H_alpha ratios for the eastern and western regions of the rings. The velocity difference for the two sides is used to derive the rotation velocity of the gas around the compact nucleus. The ratio [NII]6583/H_alpha is a factor of 2 larger in the compact nucleus of NGC 6951 than in its western side. The electron gas densities have been estimated from the [SII] lines ratio.

We present a preliminary report on the dynamical evolution of triplets of spiral galaxies that lead to merger remnants that resemble an elliptical. We find that large projected stellar tidal tails ~200 Kpc can develop during the binary formation, and up to ~350 Kpc in the third galaxy when this merges with the binary. When gas is treated as test-particles, gaseous tidal tails may reach up to ~ 1 Mpc in radius in a Hubble time; this radius is similar to the size of the common dark halo formed during the interaction. It is suggested that gaseous tails may provide a mean to determine to some extent the size of the common dark halo formed during the merging of spirals. The stellar remnant shows an elliptical-like surface density profile, but its phase-space retains more information that can be used to determine whether it was assembled from spirals or not.

The two most common mass estimators that use velocity dispersions and positions of objects in a stellar system, the virial mass estimator (VME) and the projected mass estimator (PME), are revisited and tested using N-body experiments. We consider here only spherical, isolated, and isotropic velocity dispersion systems. We have found that the PME can overestimate masses by ~ 20%, for realistic cluster mass profiles, if applied only to regions around the total system's "effective" radius. The VME can yield a correct mass at different radii provided that the total potential energy term is correctly taken into account and the system is completely sampled; otherwise, it may lead to similar errors as the PME. A surface pressure (3PV) term recently alluded to be usually neglected when using the VME and therefore required as a correction term is here found not necessarily required, although it can be used to yield a reasonable correction term. The preferred method here, however, is the virial theorem due to its simplicity and better agreement with N-body experiments. The possible reasons for the mass discrepancies found when using the PME and the VME in some N-body simulations are also briefly discussed.

The projected (PME) and virial mass estimator (VME) are revisited and tested using an N-body equilibrium system. It is found that the PME can overestimate the mass by approx 40% if a cluster is sampled only about its effective radius, Re. The exact value of this error depends on the mass profile of the cluster. The VME can yield errors of approx 20% if we are also able to sample about one Re. Both estimators yield the correct total mass for an equilibrium system. Theoretically, if we know the total extent of the system, the VME provides an accurate mass at different radii provided its gravitational potential term is correctly taken into consideration. The VME can provide acceptable masses for equilibrium systems with anisotropic velocity dispersions.