We study the tidal evolution of galaxies in the most massive cluster of the IllustrisTNG-100 simulation. For the purpose of this work we select 112 galaxies with the largest stellar masses at present and follow their properties in time. Using their orbital history we divide the sample into unevolved (infalling), weakly evolved (with one pericenter passage) and strongly evolved (with multiple pericenters). The samples are clearly separated by the value of the integrated tidal force from the cluster the galaxies experienced during their entire evolution and their properties depend strongly on this quantity. As a result of tidal stripping, the galaxies of the weakly evolved sample lost between 10 and 80% of their dark mass and less than 10% of stars, while those in the strongly evolved one more than 70% of dark and between 10 and 55% of stellar mass, and are significantly less or even not dark matter dominated. While 33% of the infalling galaxies do not contain any gas, this fraction increases to 67% for the weakly evolved and to 100% for the strongly evolved sample. The strongly evolved galaxies lose their gas earlier and faster (within 2-6 Gyr) but the process can take up to 4 Gyr from the first pericenter passage. These galaxies are redder and more metal rich, and at redshift z=0.5 the population of galaxies in the cluster becomes predominantly red. As a result of tidal stirring, the morphology of the galaxies evolves from oblate to prolate and their rotation is diminished thus the morphology-density relation is reproduced in the simulated cluster. The strongly evolved sample contains at least six convincing examples of tidally induced bars and six more galaxies that had their bars enhanced by their interaction with the cluster.

Warps are common features in both stellar and gaseous disks of nearby spiral galaxies with the latter usually easier to detect. Several theories have been proposed in the literature to explain their formation and prevalence, including tidal interactions with external galaxies. Observational correlations also suggest the importance of tides for warp formation. Here, we use the TNG100 run from the magnetohydrodynamical cosmological simulation suite IllustrisTNG to investigate the connection between interactions and the formation of gas warps. We find that in the sample of well-resolved gas-rich spiral galaxies ($10^{10}\lesssim\mathrm{M_{*}/M_{\odot}}\lesssim10^{11}$ at $z=0$) from the simulation TNG100-1, about $16\%$ possess the characteristic S-shaped warp. Around one third of these objects have their vertical morphology induced by interactions with other galaxies. Half of these interactions end with the perturber absorbed by the host. Warps formed in interactions are more asymmetrical than the remaining sample, however after the interaction the asymmetry decreases with time. We find that warps induced by interactions survive on average for $<1$ Gyr. The angle between the orbital angular momentum of the perturber and the angular momentum of the host's disk that most likely leads to warp formation is around 45 degrees. While our main goal is to investigate tidally induced warps, we find that during interactions in addition to tides, new gas that is accreted from infalling satellites also can contribute to warp formation.

By using $N$-body and hydro simulations, we study the formation and evolution of bars in galaxies with significant gas content focusing on the phenomenon of the buckling instability. The galaxies are initially composed of a spherical dark matter halo and only stellar, or stellar and gaseous, disks with parameters that are similar to the Milky Way and are evolved for 10 Gyr. We consider different values of the gas fraction $f =0-0.3$ and in order to isolate the effect of the gas, we kept the fraction constant during the evolution by not allowing the gas to cool and form stars. The stellar bars that form in simulations with higher gas fractions are weaker and shorter, and they do not form at all for gas fractions that are higher than 0.3. The bar with a gas fraction of 0.1 forms sooner due to initial perturbations in the gas, but despite the longer evolution, it does not become stronger than the one in the collisionless case at the end of evolution. The bars in the gas component are weaker; they reach their maximum strength around 4 Gyr and later decline to transform into spheroidal shapes. The distortion of the stellar bar during the buckling instability is weaker for higher gas fractions and weakens the bar less significantly, but it has a similar structure both in terms of radial profiles and in face-on projections. For $f=0.2,$ the first buckling lasts significantly longer and the bar does not undergo the secondary buckling event, while for $f=0.3,$ the buckling does not occur. Despite these differences, all bars develop boxy/peanut shapes in the stellar and gas component by the end of the evolution, although their thickness is smaller for higher gas fractions.

I review the current status of dynamical modelling of dwarf spheroidal galaxies focusing on estimates of their dark matter content. Starting with the simplest methods using the velocity dispersion profiles I discuss the inherent issues of mass-anisotropy degeneracy and contamination by unbound stars. I then move on to methods of increasing complexity, aiming to break the degeneracy, up to recent applications of the Schwarzschild orbit superposition method. The dynamical modelling is placed in the context of possible scenarios for the formation of dwarf spheroidal galaxies, including the tidal stirring model and mergers of dwarf galaxies. The two scenarios are illustrated with examples from simulations: a comparison between the tidal evolution of dwarfs with cuspy and cored dark matter profiles and the formation of a dwarf spheroidal with prolate rotation.

Using $N$-body simulations we study the buckling instability in a galactic bar forming in a Milky Way-like galaxy. The galaxy is initially composed of an axisymmetric, exponential stellar disk embedded in a spherical dark matter halo. The parameters of the model are chosen so that the galaxy is mildly unstable to bar formation and the evolution is followed for 10 Gyr. A strong bar forms slowly over the first few Gyr and buckles after 4.5 Gyr from the start of the simulation becoming much weaker and developing a pronounced boxy/peanut shape. We measure the properties of the bar at the time of buckling in terms of the mean acceleration, velocity and distortion in the vertical direction. The maps of these quantities in face-on projections reveal characteristic quadrupole patterns which wind up over a short time-scale. We also detect a secondary buckling event lasting much longer and occurring only in the outer part of the bar. We then study the orbital structure of the bar in periods before and after the first buckling. We find that most of the buckling orbits originate from x1 orbits supporting the bar. During buckling the ratio of the vertical to horizontal frequency of the stellar orbits decreases dramatically and after buckling the orbits obey a very tight relation between the vertical and circular frequency: $3 \nu = 4 \Omega$. We propose that buckling is initiated by the vertical resonance of the x1 orbits creating the initial distortion of the bar that later evolves as kinematic bending waves.

Strong galactic bars produced in simulations tend to undergo a period of buckling instability that weakens and thickens them and forms a boxy/peanut structure in their central parts. This theoretical prediction has been confirmed by identifying such morphologies in real galaxies. The nature and origin of this instability remains however poorly understood with some studies claiming it to be due to fire-hose instability while others relating it to vertical instability of stellar orbits supporting the bar. One of the channels for the formation of galactic bars is via the interaction of disky galaxies with perturbers of significant mass. Tidally induced bars offer a unique possibility of studying buckling instability because their formation can be controlled by changing the strength of the interaction while keeping the initial structure of the galaxy the same. We use a set of four simulations of flyby interactions where a galaxy on a prograde orbit forms a bar, which is stronger for stronger tidal forces. We study their buckling by calculating different kinematic signatures, including profiles of the mean velocity in vertical direction, as well as distortions of the bars out of the disk plane. Although our two strongest bars buckle most strongly, there is no direct relation between the ratio of vertical to horizontal velocity dispersion and the bar's susceptibility to buckling, as required by the fire-hose instability interpretation. While our weakest bar buckles, a stronger one does not, its dispersion ratio remains low and it grows to become the strongest of all at the end of evolution. Instead, we find that during buckling the resonance between the vertical and radial orbital frequencies becomes wide and therefore able to modify stellar orbits over a significant range of radii. We conclude that the vertical orbital instability is the more plausible explanation for the origin of buckling.

We present a full dynamical model of the Fornax dwarf spheroidal galaxy obtained with the spherically symmetric Schwarzschild orbit superposition method applied to the largest kinematic data set presently available. We modelled the total mass content of the dwarf with the mass-to-light ratio $\Upsilon$ varying with radius and found that Fornax is with high probability embedded in a massive and extended dark matter halo. We estimated the total mass contained within 1$\,$kpc to be $M(<1\,$kpc$)=1.25^{+0.06}_{-0.13} \times 10^8\,$M$_{\odot}$. The data are consistent with the constant mass-to-light ratio, i.e. the mass-follows-light model, only at 3$\sigma$ level, but still require a high value of $\Upsilon \approx 11.2\,$M$_{\odot}/$L$_{\odot}$. Our results are in general agreement with previous estimates of the dynamical mass profile of Fornax. As the Schwarzschild method does not require any assumptions on the orbital anisotropy of the stars, we obtained a profile of the anisotropy parameter $\beta$ as an output of our modelling. The derived anisotropy is close to zero in the centre of the galaxy and decreases with radius, but remains consistent with isotropic orbits at all radii at 1$\sigma$ confidence level.

We present a hydrodynamical model of M33 and its recent interaction with M31. This scenario was previously proposed in the literature in order to explain the distorted gaseous and stellar disks of M33, as well as the increased star formation rate in both objects around 2 Gyr ago. We used an orbit integration scheme to find which estimate of the transverse velocity of M31 more favors the interaction scenario and then tried to reproduce it in our simulations. M33 was modeled as a stellar and gaseous disk embedded in a live dark-matter halo, while M31 was approximated only with a live dark halo. In the simulations the two galaxies passed each other with a pericenter distance of 37 kpc. Tides excited a two-armed spiral structure in the M33 disk, which is found to be the predominant spiral signal in the observed galaxy and has been long known as a feature easily induced by tidal interactions. We found that the gaseous warp produced by the interaction did not resemble enough the observed one and we performed an additional simulation including the hot gas halo of M31 to show that this feature can be properly reproduced by tidal forces and ram pressure stripping acting simultaneously on the gaseous disk. In addition to the spiral arms, tidal forces produced the stellar stream similar to the observed one and triggered a star formation burst at similar radii as it is observed.

Bars in disky galaxies can be formed by interactions with other systems, including those of comparable mass. It has long been established that the effect of such interactions on galaxy morphology depends strongly on the orbital configuration, in particular the orientation of the intrinsic spin of the galactic disk with respect to its orbital angular momentum. Prograde encounters modify the morphology strongly, including the formation of tidally induced bars, while retrograde flybys should have little effect on morphology. Recent works on the subject reached conflicting conclusions, one using the impulse approximation and claiming no dependence on this angle in the properties of tidal bars. To resolve the controversy, we performed self-consistent N-body simulations of hyperbolic encounters between two identical Milky Way-like galaxies assuming different velocities and impact parameters, with one of the galaxies on a prograde and the other on a retrograde orbit. The galaxies were initially composed of an exponential stellar disk and an NFW dark halo, and were stable against bar formation in isolation for 3 Gyr. We find that strong tidally induced bars form only in galaxies on prograde orbits. For smaller impact parameters and lower relative velocities the bars are stronger and have lower pattern speeds. Stronger bars undergo extended periods of buckling instability that thicken their vertical structure. The encounters also lead to the formation of two-armed spirals with strength inversely proportional to the strength of the bars. We conclude that proper modeling of prograde and retrograde encounters cannot rely on the simplest impulse approximation.

In our previous work we confirmed the reliability of the spherically symmetric Schwarzschild orbit-superposition method to recover the mass and velocity anisotropy profiles of spherical dwarf galaxies. Here we investigate the effect of its application to intrinsically non-spherical objects. For this purpose we use a model of a dwarf spheroidal galaxy formed in a numerical simulation of a major merger of two disky dwarfs. The shape of the stellar component of the merger remnant is axisymmetric and prolate which allows us to identify and measure the bias caused by observing the spheroidal galaxy along different directions, especially the longest and shortest principal axis. The modelling is based on mock data generated from the remnant that are observationally available for dwarfs: projected positions and line-of-sight velocities of the stars. In order to obtain a reliable tool while keeping the number of parameters low we parametrize the total mass distribution as a radius-dependent mass-to-light ratio with just two free parameters we aim to constrain. Our study shows that if the total density profile is known, the true, radially increasing anisotropy profile can be well recovered for the observations along the longest axis whereas the data along the shortest axis lead to the inference an incorrect, isotropic model. On the other hand, if the density profile is derived from the method as well, the anisotropy is always underestimated but the total mass profile is well recovered for the data along the shortest axis whereas for the longest axis the mass content is overestimated.

Tens of early type galaxies have been recently reported to possess prolate rotation, i.e. significant amount of rotation around the major axis, including two cases in the Local Group. Although expected theoretically, this phenomenon is rarely observed and remains elusive. In order to explore its origin we study the population of well-resolved galaxies in the Illustris cosmological simulation. We identify 59 convincing examples of prolate rotators at the present time, more frequently among more massive galaxies, with the number varying very little with redshift. We follow their evolution back in time using the main progenitor branch galaxies of the Illustris merger trees. We find that the emergence of prolate rotation is strongly correlated with the time of the last significant merger the galaxy experienced, although other evolutionary paths leading to prolate rotation are also possible. The transition to prolate rotation most often happens around the same time as the transition to prolate shape of the stellar component. The mergers leading to prolate rotation have slightly more radial orbits, higher mass ratios, and occur at more recent times than mergers in the reference sample of twin galaxies we construct for comparison. However, they cover a wide range of initial conditions in terms of the mass ratio, merger time, radiality of the progenitor orbits, and the relative orientations of spins of the progenitors with respect to the orbital angular momenta. About half of our sample of prolate rotators were created during gas-rich mergers and the newly formed stars usually support prolate rotation.

We present a new study concerning the application of the Schwarzschild orbit superposition method to model spherical galaxies. The method aims to recover the mass and the orbit anisotropy parameter profiles of the objects using measurements of positions and line-of-sight velocities usually available for resolved stellar populations of dwarf galaxies in the Local Group. To test the reliability of the method, we used different sets of mock data extracted from four numerical realizations of dark matter haloes. The models shared the same density profile but differed in anisotropy profiles, covering a wide range of possibilities, from constant to increasing and decreasing with radius. The tests were done in two steps, first assuming that the mass profile of the dwarf is known and employing the method to retrieve the anisotropy only, and then varying also the mass distribution. We used two kinds of data samples: unrealistically large ones based on over 270 000 particles from the numerical realizations and small ones matching the amount of data available for the Fornax dwarf. For the large data samples we recover both the mass and the anisotropy profiles with very high accuracy. For the realistically small ones we also find a reasonably good agreement between the fitted and the input anisotropies, however the total density profiles can be significantly biased as a result of their oversensitivity to the available data. Our results therefore provide convincing evidence in favour of the applicability of the Schwarzschild method to break the mass-anisotropy degeneracy in dwarf galaxies.

Bars in galaxies may develop through a global instability or due to an interaction with another system. We study bar formation in disky dwarf galaxies orbiting a Milky Way-like galaxy. We employ $N$-body simulations to study the impact of initial orbital parameters: the size of the dwarf galaxy orbit and the inclination of its disc with respect to the orbital plane. In all cases a bar develops in the center of the dwarf during the first pericenter on its orbit around the host. Between subsequent pericenter passages the bars are stable, but at the pericenters they are usually weakened and shortened. The initial properties and details of the further evolution of the bars depend heavily on the orbital configuration. We find that for the exactly prograde orientation, the strongest bar is formed for the intermediate-size orbit. On the tighter orbit, the disc is too disturbed and stripped to form a strong bar. On the wider orbit, the tidal interaction is too weak. The dependence on the disc inclination is such that weaker bars form in more inclined discs. The bars experience either a very weak buckling or none at all. We do not observe any secular evolution, possibly because the dwarfs are perturbed at each pericenter passage. The rotation speed of the bars can be classified as slow ($R_\mathrm{CR}/l_\mathrm{bar}\sim2-3$). We attribute this to the loss of a significant fraction of the disc's rotation during the encounter with the host galaxy.

We analyze in detail the spatial distribution and kinematic properties of two different stellar populations in Andromeda II (And II) dwarf spheroidal galaxy. We obtained their detailed surface density maps, together with their radial density profiles. The two populations differ not only in age and metallicity, but also in their spatial distribution and kinematics. Old stars ($\gtrsim 11$ Gyr) follow a round distribution well fitted by truncated density profiles. These stars rotate around the projected optical major axis of the galaxy with line-of-sight velocities $v_{los}(r_h) = 16 \pm 3$ km s$^{-1}$ and a velocity gradient of $2.06 \pm 0.21$ km s$^{-1}$ arcmin$^{-1}$. Intermediate-age stars ($\lesssim 9$ Gyr) concentrate in the centre of the galaxy and form an elongated structure extending along the projected optical major axis. This structure appears to rotate with a velocity gradient of $2.24 \pm 0.22$ km s$^{-1}$ arcmin$^{-1}$, and around the optical minor axis. The centres of rotation and kinetic position angles (PA$_{\rm kin}$) of both populations differ. For intermediate-age stars we obtained PA$_{\rm kin} = 18^\circ \pm 2^\circ$ and for the old ones PA$_{\rm kin} = 63^\circ \pm 3^\circ$ in good agreement with photometric PA measured from isopleths fitted to the photometry. We conclude that the two stellar populations may not be fully relaxed and thus be the result of a merger occurred at redshift $z\sim 1.75$.

We present a novel analysis of the internal kinematics of the Fornax dwarf spheroidal galaxy. Our results are based on the largest sample of spectroscopic data for Fornax stars presently available ($> 2500$ stars), for which we have chemical and kinematic information. We introduce new software, Beacon, designed to detect chemo-kinematic patterns among stars of different stellar populations using their metallicity and velocity along the line of sight. Applying Beacon to Fornax we have detected non-negligible rotation signals around main optical axes of the galaxy, characteristic for a triaxial system partially supported by rotation. The dominant rotation pattern is relatively strong ($\sim 12$ km s$^{-1}$), but the galaxy also shows additional weaker albeit complex rotation patterns. Using the information available from the star formation history of Fornax we have also derived the average age of the different chemo-kinematic components found by Beacon, which has allowed us to obtain its kinematic history. Our results point to a possible major merger suffered by Fornax at redshift $z\sim1$, in agreement with previous works.

One of the scenarios for the formation of grand-design spiral arms in disky galaxies involves their interactions with a satellite or another galaxy. Here we consider another possibility, where the perturbation is instead due to the potential of a galaxy cluster. Using $N$-body simulations we investigate the formation and evolution of spiral arms in a Milky Way-like galaxy orbiting a Virgo-like cluster. The galaxy is placed on a few orbits of different size but similar eccentricity and its evolution is followed for 10 Gyr. The tidally induced, two-armed, approximately logarithmic spiral structure forms on each of them during the pericenter passages. The spiral arms dissipate and wind up with time, to be triggered again at the next pericenter passage. We confirm this transient and recurrent nature of the arms by analyzing the time evolution of the pitch angle and the arm strength. We find that the strongest arms are formed on the tightest orbit, however they wind up rather quickly and are disturbed by another pericenter passage. The arms on the most extended orbit, which we analyze in more detail, wind up slowly and survive for the longest time. Measurements of the pattern speed of the arms indicate that they are kinematic density waves. We attempt a comparison with observations by selecting grand-design spiral galaxies in the Virgo cluster. Among those, we find nine examples bearing no signs of recent interactions or the presence of companions. For three of them we present close structural analogues among our simulated spiral galaxies.

We describe a new phenomenon of `bar damping' that may have played an important role in shaping the Milky Way bar and bulge as well as its spiral structure. We use a collisionless N-body simulation of a Milky Way-like galaxy initially composed of a dark matter halo and an exponential disk with Toomre parameter slightly above unity. In this configuration, dominated by the disk in the center, a bar forms relatively quickly, after 1 Gyr of evolution. This is immediately followed by the formation of two manifold-driven spiral arms and the outflow of stars that modifies the potential in the vicinity of the bar, apparently shifting the position of the L_1/L_2 Lagrange points. This modification leads to the shortening of the bar and the creation of a next generation of manifold-driven spiral arms at a smaller radius. The process repeats itself a few times over the next 0.5 Gyr resulting in further substantial weakening and shortening of the bar. The time when the damping comes to an end coincides with the first buckling episode in the bar which rebuilds the orbital structure so that no more new spiral arms are formed. The morphology of the bar and the spiral structure at this time show remarkable similarity to the present properties of the Milky Way. Later on, the bar starts to grow rather steadily again, weakened only by subsequent buckling episodes occurring at more distant parts of the disk.

Andromeda II (And II) has been known for a few decades but only recently observations have unveiled new properties of this dwarf spheroidal galaxy. The presence of two stellar populations, the bimodal star formation history (SFH) and an unusual rotation velocity of And II put strong constrains on its formation and evolution. Following Lokas et al. (2014), we propose a detailed model to explain the main properties of And II involving (1) a gas-rich major merger between two dwarf galaxies at high redshift in the field and (2) a close interaction with M31 about 5 Gyr ago. The model is based on N-body/hydrodynamical simulations including gas dynamics, star formation and feedback. One simulation is designed to reproduce the gas-rich major merger explaining the origin of stellar populations and the SFH. Other simulations are used to study the effects of tidal forces and the ram pressure stripping during the interaction between And II and M31. The model successfully reproduces the SFH of And II including the properties of stellar populations, its morphology, kinematics and the lack of gas. Further improvements to the model are possible via joint modelling of all processes and better treatment of baryonic physics.

Orbits are the key building blocks of any density distribution and their study helps us understand the kinematical structure and the evolution of galaxies. Here we investigate orbits in a tidally induced bar of a dwarf galaxy, using an $N$-body simulation of an initially disky dwarf galaxy orbiting a Milky Way-like host. After the first pericenter passage, a tidally induced bar forms in the stellar component of the dwarf. The bar evolution is different than in isolated galaxies and our analysis focuses on the period before it buckles. We study the orbits in terms of their dominant frequencies, which we calculate in a Cartesian coordinate frame rotating with the bar. Apart from the well-known x$_1$ orbits we find many other types, mostly with boxy shapes of various degree of elongation. Some of them are also near-periodic, admitting frequency ratios of 4/3, 3/2 and 5/3. The box orbits have various degrees of vertical thickness but only a relatively small fraction of those have banana (i.e. smile/frown) or infinity-symbol shapes in the edge-on view. In the very center we also find orbits known from the potential of triaxial ellipsoids. The elongation of the orbits grows with distance from the center of the bar in agreement with the variation of the shape of the density distribution. Our classification of orbits leads to the conclusion that more than $80 \%$ of them have boxy shapes, while only $8 \%$ have shapes of classical x$_1$ orbits.

We present new radial velocity (RV) measurements of old (horizontal branch) and intermediate-age (red clump) stellar tracers in the Carina dwarf spheroidal. They are based on more than 2,200 low-resolution spectra collected with VIMOS at VLT. The targets are faint (20<V<21.5 mag), but the accuracy at the faintest limit is <9 kms-1. These data were complemented with RV measurements either based on spectra collected with FORS2 and FLAMES/GIRAFFE at VLT or available in the literature. We ended up with a sample of 2748 stars and among them 1389 are candidate Carina stars. We found that the intermediate-age stellar component shows a well defined rotational pattern around the minor axis. The western and the eastern side of the galaxy differ by +5 and -4 km s-1 when compared with the main RV peak. The old stellar component is characterized by a larger RV dispersion and does not show evidence of RV pattern. We compared the observed RV distribution with N-body simulations for a former disky dwarf galaxy orbiting a giant MilkyWay-like galaxy (Lokas et al. 2015). We rotated the simulated galaxy by 60 degrees with respect to the major axis, we kept the observer on orbital plane of the dwarf and extracted a sample of stars similar to the observed one. Observed and predicted Vrot/{\sigma} ratios across the central regions are in remarkable agreement. This evidence indicates that Carina was a disky dwarf galaxy that experienced several strong tidal interactions with the Milky Way. Owing to these interactions, Carina transformed from a disky to a prolate spheroid and the rotational velocity transformed into random motions.

Using N-body simulations we study the formation and evolution of tidally induced bars in disky galaxies in clusters. Our progenitor is a massive, late-type galaxy similar to the Milky Way, composed of an exponential disk and an NFW dark matter halo. We place the galaxy on four different orbits in a Virgo-like cluster and evolve it for 10 Gyr. As a reference case we also evolve the same model in isolation. Tidally induced bars form on all orbits soon after the first pericenter passage and survive until the end of the evolution. They appear earlier, are stronger, longer and have lower pattern speeds for tighter orbits. Only for the tightest orbit the properties of the bar are controlled by the orientation of the tidal torque from the cluster at pericenters. The mechanism behind the formation of the bars is the angular momentum transfer from the galaxy stellar component to its halo. All bars undergo extended periods of buckling instability that occur earlier and lead to more pronounced boxy/peanut shapes when the tidal forces are stronger. Using all simulation outputs of galaxies at different evolutionary stages we construct a toy model of the galaxy population in the cluster and measure the average bar strength and bar fraction as a function of clustercentric radius. Both are found to be mildly decreasing functions of radius. We conclude that tidal forces can trigger bar formation in cluster cores, but not in the outskirts, and thus cause larger concentrations of barred galaxies towards cluster center.

A tidal radius is a distance from a satellite orbiting in a host potential beyond which its material is stripped by the tidal force. We derive a revised expression for the tidal radius of a rotating satellite which properly takes into account the possibility of prograde and retrograde orbits of stars. Besides the eccentricity of the satellite orbit, the tidal radius depends also on the ratio of the satellite internal angular velocity to the orbital angular velocity. We compare our formula to the results of two $N$-body simulations of dwarf galaxies orbiting a Milky Way-like host on a prograde and retrograde orbit. The tidal radius for the retrograde case is larger than for the prograde. We introduce a kinematic radius separating stars still orbiting the dwarf galaxy from those already stripped and following the potential of the host galaxy. We find that the tidal radius matches very well the kinematic radius. Our results provide a connection between the formalism of the tidal radius derivation and the theory of resonant stripping.

We apply the Schwarzschild orbit superposition method to mock data in order to investigate the accuracy of recovering the profile of the orbit anisotropy. The mock data come from four numerical realizations of dark matter haloes with well defined anisotropy profiles. We show that when assuming a correct mass distribution we are able to determine the anisotropy with high precision and clearly distinguish between the models.

The origin of grand design spiral structure in galaxies is still under debate but one of promising scenarios involves tidal interactions. We use N-body simulations to study the evolution of a Milky Way-size galaxy in a Virgo-like cluster. The galaxy is placed on a typical eccentric orbit and evolved for 10 Gyr. We find that grand design spiral arms are triggered by pericenter passages and later on they wind up and dissipate. The arms formed in the simulations are approximately logarithmic, but are also dynamic, transient and recurrent.

Using numerical analysis of fundamental frequencies we study the orbital structure of a tidally induced bar formed in a simulated dwarf galaxy orbiting a Milky Way-like host. We find that only about 10% of stars have frequencies compatible with x1 orbits, the classical periodic orbits in a barred potential. The rest of the stars follows box orbits parallel to the bar, with varying degree of elongation.

Motivated by the discovery of prolate rotation of stars in Andromeda II, a dwarf spheroidal companion of M31, we study its origin via mergers of disky dwarf galaxies. We simulate merger events between two identical dwarfs changing the initial inclination of their disks with respect to the orbit and the amount of orbital angular momentum. On radial orbits the amount of prolate rotation in the merger remnants correlates strongly with the inclination of the disks and is well understood as due to the conservation of the angular momentum component of the disks along the merger axis. For non-radial orbits prolate rotation may still be produced if the orbital angular momentum is initially not much larger than the intrinsic angular momentum of the disks. The orbital structure of the remnants with significant rotation is dominated by box orbits in the center and long-axis tubes in the outer parts. The frequency analysis of stellar orbits in the plane perpendicular to the major axis reveals the presence of two families roughly corresponding to inner and outer long-axis tubes. The fraction of inner tubes is largest in the remnant forming from disks oriented most vertically initially and is responsible for the boxy shape of the galaxy. We conclude that prolate rotation results from mergers with a variety of initial conditions and no fine tuning is necessary to reproduce this feature. We compare the properties of our merger remnants to those of dwarfs resulting from the tidal stirring scenario and the data for Andromeda II.

Using N-body simulations we study the tidal evolution of initially disky dwarf galaxies orbiting a Milky Way-like host, a process known to lead to the formation of dwarf spheroidal galaxies. We focus on the effect of the orientation of the dwarf galaxy disk's angular momentum with respect to the orbital one and find very strong dependence of the evolution on this parameter. We consider four different orientations: the exactly prograde, the exactly retrograde and two intermediate ones. Tidal evolution is strongest for the exactly prograde and weakest for the exactly retrograde orbit. In the prograde case the stellar component forms a strong bar and remains prolate until the end of the simulation, while its rotation is very quickly replaced by random motions of the stars. In the retrograde case the dwarf remains oblate, does not form a bar and loses rotation very slowly. In the two cases of intermediate orientation of the disk, the evolution is between the two extremes, suggesting a monotonic dependence on the inclination. We interpret the results in terms of the resonance between the angular velocity of the stars in the dwarf and its orbital motion by comparing the measurements from simulations to semi-analytic predictions. We conclude that resonant effects are the most important mechanism underlying the tidal evolution of disky dwarf galaxies.

Elliptical galaxies are modelled as Sersic luminosity distributions with density profiles (DPs) for the total mass adopted from the DPs of haloes within dissipationless LambdaCDM N-body simulations. Ellipticals turn out to be inconsistent with cuspy low-concentration NFW models representing the total mass, nor are they consistent with a steeper -1.5 inner slope, nor with the shallower models proposed by Navarro et al. 04, nor with NFW models 10 times more concentrated than predicted, as deduced from several X-ray observations: the mass models, extrapolated inwards, lead to local mass-to-light ratios that are smaller than the stellar value inside an effective radius (R_e), and to central aperture velocity dispersions that are much smaller than observed. This conclusion remains true as long as there is no sharp steepening (slope < -2) of the dark matter (DM) DPs just inside 0.01 virial radii. The too low total mass and velocity dispersion produced within R_e by an NFW-like total mass profile suggests that the stellar component should dominate the DM one out to at least R_e. It should then be difficult to kinematically constrain the inner slope of the dark matter DP of ellipticals. The high concentration parameters deduced from X-ray observations appear to be a consequence of fitting an NFW model to the total mass DP made up of a stellar component that dominates inside and a DM component that dominates outwards. An appendix gives the virial mass dependence of the concentration parameter, central density, and total mass of the Navarro et al. model. In a 2nd appendix are given single integral expressions for the velocity dispersions averaged along the line-of-sight, in circular apertures and in thin slits, for general luminosity density and mass distributions, with isotropic orbits.

Elliptical galaxies are modelled with a a 4-component model: Sersic stars, LCDM dark matter (DM), hot gas and central black hole. DM is negligible in the inner regions, which are dominated by stars and the central black hole. This prevents any kinematical estimate (using a Jeans analysis) of the inner slope of the DM density profile. The gas fraction rises, but the baryon fraction decreases with radius, at least out to 10 effective radii (R_e). Even with line-of-sight velocity dispersion (VD) measurements at 4 to 6 R_e with 20 km/s accuracy and perfectly known velocity anisotropy, the total mass within the virial radius (r_v) is uncertain by a factor over 3. The DM distributions found in LCDM simulations are consistent with the stellar VD profiles, but appear inconsistent with the low VDs measured by Romanowsky et al. (2003) of planetary nebulae between 2 and 5 R_e, which imply such low M/Ls that the baryon fraction within r_v must be greater than the universal value. Replacing the NFW DM model by the new model of Navarro et al. (2004) decreases slightly the VD at a given radius. So, given the observed VD measured at 5 R_e, the inferred M/L within r_v is 40% larger than predicted with the NFW model. Folding in the slight (strong) radial anisotropy found in LCDM (merger) simulations, which is well modelled (much better than with the Osipkov-Merritt formula) with beta(r) = 1/2 r/(r+a), the inferred M/L within r_v is another 1.6 (2.4) times higher than for the isotropic NFW model. Thus, the DM model and radial anisotropy can partly explain the low PN VDs, but not in full. In an appendix, single integral expressions are derived for the VDs in terms of the tracer density and total mass profiles, for 3 anisotropic models: radial, Osipkov-Merritt, and the model above, for general radial profiles of luminosity density and mass.

Using N-body simulations we study the phenomenon of radial orbit instability occurring in dark matter haloes of the size of a dwarf galaxy. We carried out simulations of seven spherical models, with the same standard NFW density profile but different anisotropy profiles of particle orbits. Four of them underwent instability: two with a constant positive anisotropy, one with an anisotropic core and an isotropic envelope and one with a very small isotropic core and an anisotropic envelope. Haloes affected by the instability become approximately axisymmetric and prolate, with the profile of the shortest-to-longest axis ratio increasing with radius. The lower limit for the central value of this axis ratio is 0.3 for an NFW halo. The density profiles of the haloes did not change significantly, whereas the velocity distributions became axisymmetric. The total angular momentum rose due to large-amplitude oscillations of its components perpendicular to the symmetry axis of the halo. We also studied orbits of individual particles assigning them to classical orbit families in triaxial potentials. We find that the membership of a given particle in a family depends on its initial total angular momentum and its components along the principal axes of matter distribution.

Using N-body simulations, we study the properties of a bar induced in a dwarf galaxy as a result of tidal interaction with the Milky Way. The dwarf is initially composed of a disk embedded in a dark matter halo and we follow its evolution on a typical orbit for 10 Gyr. It undergoes an evolution typical of tidally stirred dwarf galaxies: it loses mass, the stellar component transforms from a disk to a spheroid and the rotation of the stars is partially replaced by random motions. A tidally induced bar forms at the first pericentre passage and survives until the end of the evolution. Fourier decomposition of the distribution of stars reveals that only even modes are significant and preserve a hierarchy so that the bar mode is always the strongest. They show a characteristic profile with a maximum, similar to simulated bars forming in isolated galaxies and observed bars in real galaxies. We adopt the maximum of the bar mode as a measure of the bar strength and we estimate the bar length by comparing the density profiles along the bar and perpendicular to it. Both the bar strength and the bar length decrease with time, mainly at pericentres, as a result of tidal torques acting at those times and not to secular evolution. The pattern speed of the bar varies significantly on a time scale of 1 Gyr and seems to be fully controlled by the orientation of the tidal torque from the Milky Way. The bar is never tidally locked but we discover a hint of a 5/2 orbital resonance between the third and fourth pericentre passage. The speed of the bar expressed in terms of the ratio of the bar length to the corotation radius decreases in the long run so that the bar changes from initially rather fast to slow in the later stages. The boxy/peanut shape is clearly visible in the edge-on view between the second and third pericentre and its occurrence is preceded by a short period of buckling instability.

Precise measurements of mass in dark matter dominated dwarf spheroidal galaxies are of great importance for testing the theories of structure formation. We use $N$-body simulations of the tidal evolution of a dwarf galaxy orbiting the Milky Way to generate mock kinematical data sets and use them to test the reliability of a simple mass estimator proposed by Wolf et al. The evolution of the initially disky dwarf galaxy embedded in a dark matter halo was traced for 10 Gyr on a rather tight orbit. After about half of the time a dwarf spheroidal galaxy is formed that retains some remnant rotation and a non-spherical shape. Observing the triaxial galaxy along each of its principal axes we measure its half-light radius and the line-of-sight velocity dispersion and use them to estimate the mass. We find that the mass is significantly overestimated when the dwarf is seen along the longest axis of the stellar component and underestimated when observed along the shortest axis. We provide a formula that quantifies the systematic error in the estimated mass with respect to the true one as a function of the galaxy shape and line of sight.

The formation of dwarf spheroidal galaxies in the Local Group from disky progenitors via tidal interaction with a bigger host is one of the most promising scenarios of their origin. Using N-body simulations we study the process by following the evolution of a disky dwarf orbiting a Milky Way-like host. We focus on the effect of the orientation of the dwarf galaxy disk's angular momentum with respect to the orbital one. We find a strong dependence of the efficiency of the transformation from a disk to a spheroid on the disk orientation. The effect is strongest for the exactly prograde and weakest for the exactly retrograde orbit. In the prograde case the stellar component forms a strong bar and remains prolate until the end of the evolution, while its rotation is very quickly replaced by random motions of the stars. In the retrograde case the dwarf remains oblate, does not form a bar and loses rotation very slowly. Our results suggest that resonant effects are the most important mechanism underlying the evolution while tidal shocking plays only a minor role.

Using N-body simulations we study the origin of prolate rotation recently detected in the kinematic data for And II, a dSph satellite of M31. We propose an evolutionary model for the origin of And II involving a merger between two disky dwarf galaxies whose structural parameters differ only in their disk scale lengths. The dwarfs are placed on a radial orbit towards each other with their angular momenta inclined by 45 deg to the orbital plane and by 90 deg with respect to each other. After 5 Gyr of evolution the merger remnant forms a stable triaxial galaxy with rotation only around the longest axis. The origin of this rotation is naturally explained as due to the symmetry of the initial configuration which leads to the conservation of angular momentum components along the direction of the merger. The stars originating from the two dwarfs show significantly different surface density profiles while having very similar kinematics in agreement with the properties of separate stellar populations in And II. We also study an alternative scenario for the formation of And II, via tidal stirring of a disky dwarf galaxy. While intrinsic rotation occurs naturally in this model as a remnant of the initial rotation of the disk, it is mostly around the shortest axis of the stellar component. The rotation around the longest axis is induced only occasionally and remains much smaller that the system's velocity dispersion. We conclude that although tidal origin of the velocity distribution in And II cannot be excluded, it is much more naturally explained within the scenario involving a past merger event. Thus, in principle, the presence of prolate rotation in dSph galaxies of the Local Group and beyond may be used as an indicator of major mergers in their history or even as a way to distinguish between the two scenarios of their formation.

We present a phenomenological description of the properties of tidal tails forming around dwarf galaxies orbiting the Milky Way. For this purpose we use collisionless N-body simulations of dwarfs initially composed of a disk embedded in an NFW dark matter halo. The dwarfs are placed on seven orbits around the Milky Way-like host, differing in size and eccentricity, and their evolution is followed for 10 Gyr. In addition to the well-studied morphological and dynamical transformation of the dwarf's main body, the tidal stripping causes them to lose a substantial fraction of mass both in dark matter and stars which form pronounced tidal tails. We focus on the properties of the stellar component of the tidal tails thus formed. We first discuss the break radii in the stellar density profile defining the transition to tidal tails as the radii where the profile becomes shallower and relate them to the classically defined tidal radii. We then calculate the relative density and velocity of the tails at a few break radii as a function of the orbital phase. Next, we measure the orientation of the tails with respect to an observer placed at the centre of the Milky Way. The tails are perpendicular to this line of sight only for a short period of time near the pericentre. For most of the time the angles between the tails and this line of sight are low, with orbit-averaged medians below 42 degrees for all, even the almost circular orbit. The median angle is typically lower while the maximum relative density higher for more eccentric orbits. The combined effects of relative density and orientation of the tails suggest that they should be easiest to detect for dwarf galaxies soon after their pericentre passage.

Using collisionless N-body simulations of dwarf galaxies orbiting the Milky Way (MW) we construct realistic models of dwarf spheroidal (dSph) galaxies of the Local Group. The dwarfs are initially composed of stellar disks embedded in dark matter haloes with different inner density slopes and are placed on an eccentric orbit typical for MW subhaloes. After a few Gyr of evolution the stellar component is triaxial as a result of bar instability induced by tidal forces. Observing the simulated dwarfs along three principal axes of the stellar component we create mock data sets and determine the their half-light radii and line-of-sight velocity dispersions. Using the estimator proposed by Wolf et al. we calculate masses within half-light radii. The masses obtained this way are over(under)estimated by up to a factor of two when the line of sight is along the longest (shortest) axis of the stellar component. We then divide the initial stellar distribution into an inner and outer population and trace their evolution in time. The two populations, although affected by tidal forces, retain different density profiles even after a few Gyr. We measure the half-light radii and velocity dispersions of the stars in the two populations along different lines of sight and use them to estimate the slope of the mass distribution in the dwarfs following the method proposed by Walker & Penarrubia. The inferred slopes are systematically over- or underestimated, depending on the line of sight. In particular, when the dwarf is seen along the longest axis of the stellar component, a significantly shallower density profile is inferred than the real one measured from the simulations. Since most dSphs are non-spherical and their orientation with respect to our line of sight is unknown, the method can be reliably applied only to a large sample of dwarfs when these systematic errors are expected to be diminished.

(Abridged) The origin of dSphs in the Local Group (LG) remains an enigma. The tidal stirring model posits that late-type, rotationally-supported dwarfs resembling present-day dwarf irregular (dIrr) galaxies can transform into dSphs via interactions with Milky Way-sized hosts. Using collisionless N-body simulations, we investigate for the first time how tidal stirring depends on the dark matter (DM) density distribution in the central stellar region of the progenitor disky dwarf. Specifically, we explore various asymptotic inner slopes gamma of the dwarf DM density profiles (rho \propto r^{-gamma} as r -> 0). For a given orbit inside the primary, rotationally-supported dwarfs embedded in DM halos with core-like density distributions (gamma = 0.2) and mild density cusps (gamma = 0.6) demonstrate a substantially enhanced likelihood and efficiency of transformation into dSphs compared to their counterparts with steeper DM density profiles (gamma = 1). Such shallow DM distributions are akin to those of observed dIrrs, highlighting tidal stirring as a plausible model for the LG morphology-density relation. When gamma <1, a single pericentric passage can induce dSph formation and disky dwarfs on low-eccentricity or large-pericenter orbits are able to transform into dSphs; these new results allow the tidal stirring model to explain the existence of virtually all known dSphs across a wide range of distances from their hosts. A subset of rotationally-supported dwarfs with gamma <1 are eventually disrupted by the primary; those that survive as dSphs are generally on orbits that are biased towards lower eccentricities and/or larger pericenters relative to those of typical CDM satellites. The latter could explain the rather peculiar orbits of several classic LG dSphs such as Fornax, Leo I, Tucana, and Cetus.

Using N-body simulations we study the evolution of separate stellar populations in dwarf galaxies in the context of the tidal stirring scenario for the formation of dwarf spheroidal (dSph) galaxies in the Local Group. The dwarf galaxies, initially composed of a stellar disk and a dark matter halo, are placed on seven different orbits around the Milky Way. The stars are divided into two populations, within and outside the half-light radius, and their positions are followed for 10 Gyr. We find that the populations retain different density distributions even over such long timescales. Some of the stars of the outer population migrate to the central part of the dwarf forming an extended core while the stars of the inner population develop a tail in the outer parts. In addition, the outer population is more heavily stripped by tidal forces from the Milky Way and may become subdominant at all radii on tight enough orbits. We conclude that the tidal stirring model is fully compatible with the presence of multiple stellar populations in dSph galaxies.

We have assembled a large-area spectroscopic survey of giant stars in the Sagittarius (Sgr) dwarf galaxy core. Using medium resolution (R ~15,000), multifiber spectroscopy we have measured velocities of these stars, which extend up to 12 degrees from the galaxy's center (3.7 core radii or 0.4 times the King limiting radius). From these high quality spectra we identify 1310 Sgr members out of 2296 stars surveyed distributed across 24 different fields across the Sgr core. Additional slit spectra were obtained of stars bridging from the Sgr core to its trailing tail. Our systematic, large area sample shows no evidence for significant rotation, a result at odds with the ~20 km/s rotation required as an explanation for the bifurcation seen in the Sgr tidal stream; the observed small (<= 4 km/s) velocity trend along primarily the major axis is consistent with models of the projected motion of an extended body on the sky with no need for intrinsic rotation. The Sgr core is found to have a flat velocity dispersion (except for a kinematically colder center point) across its surveyed extent and into its tidal tails, a property that matches the velocity dispersion profiles measured for other Milky Way dwarf spheroidal (dSph) galaxies. We comment on the possible significance of this observed kinematical similarity for the dynamical state of the other classical Milky Way dSphs in light of the fact that Sgr is clearly a strongly tidally disrupted system.

In recent years the Sloan Digital Sky Survey has unraveled a new population of ultra-faint dwarf galaxies (UFDs) in the vicinity of the Milky Way (MW) whose origin remains a puzzle. Using a suite of collisionless N-body simulations, we investigate the formation of UFDs in the context of the tidal stirring model for the formation of dwarf spheroidal galaxies in the Local Group (LG). Our simulations are designed to reproduce the tidal interactions between MW-sized host galaxies and rotationally supported dwarfs embedded in 10^9 M_sun dark matter (DM) halos. We explore a variety of inner density slopes \rho \propto r^{-\alpha} for the dwarf DM halos, ranging from core-like (\alpha = 0.2) to cuspy (\alpha = 1), and different dwarf orbital configurations. Our experiments demonstrate that UFDs can be produced via tidal stirring of disky dwarfs on relatively tight orbits, consistent with a redshift of accretion by the host galaxy of z \sim 1, and with intermediate values for the halo inner density slopes (\rho \propto r^{-0.6}). The inferred slopes are in excellent agreement with those resulting from both the modeling of the rotation curves of dwarf galaxies and recent cosmological simulations of dwarf galaxy formation. Comparing the properties of observed UFDs with those of their simulated counterparts, we find remarkable similarities in terms of basic observational parameters. We conclude that tidal stirring of rotationally supported dwarfs represents a viable mechanism for the formation of UFDs in the LG environment.

We study the shapes of Milky Way satellites in the context of the tidal stirring scenario for the formation of dwarf spheroidal galaxies. The standard procedures used to measure shapes involve smoothing and binning of data and thus may not be sufficient to detect structural properties like bars, which are usually subtle in low surface brightness systems. Taking advantage of the fact that in nearby dwarfs photometry of individual stars is available we introduce discrete measures of shape based on the two-dimensional inertia tensor and the Fourier bar mode. We apply these measures of shape first to a variety of simulated dwarf galaxies formed via tidal stirring of disks embedded in dark matter halos and orbiting the Milky Way. In addition to strong mass loss and randomization of stellar orbits, the disks undergo morphological transformation that typically involves the formation of a triaxial bar after the first pericenter passage. These tidally induced bars persist for a few Gyr before being shortened towards a more spherical shape if the tidal force is strong enough. We test this prediction by measuring in a similar way the shape of nearby dwarf galaxies, satellites of the Milky Way. We detect inner bars in Ursa Minor, Sagittarius, LMC and possibly Carina. In addition, six out of eleven studied dwarfs show elongated stellar distributions in the outer parts that may signify transition to tidal tails. We thus find the shapes of Milky Way satellites to be consistent with the predictions of the tidal stirring model.

We review basic properties of the population of dwarf galaxies in the Local Group focusing on dwarf spheroidal galaxies found in the immediate vicinity of the Milky Way. The evidence for dark matter in these objects is critically assessed. We describe the methods of dynamical modelling of such objects, using a few examples of the best-studied dwarfs and discuss the sources of uncertainties in mass estimates. We conclude with perspectives for dwarf galaxies as targets for dark matter detection experiments.

We perform collisionless N-body simulations to investigate whether binary mergers between rotationally-supported dwarfs can lead to the formation of dwarf spheroidal galaxies (dSphs). Our simulation campaign is based on a hybrid approach combining cosmological simulations and controlled numerical experiments. We select merger events from a Constrained Local UniversE (CLUES) simulation of the Local Group (LG) and record the properties of the interacting dwarf-sized halos. This information is subsequently used to seed controlled experiments of binary encounters between dwarf galaxies consisting of exponential stellar disks embedded in cosmologically-motivated dark matter halos. These simulations are designed to reproduce eight cosmological merger events, with initial masses of the interacting systems in the range ~ (5-60) x 10^7 Mo, occurring quite early in the history of the LG, more than 10 Gyr ago. We compute the properties of the merger remnants as a distant observer would and demonstrate that at least three of the simulated encounters produce systems with kinematic and structural properties akin to those of the classic dSphs in the LG. Tracing the history of the remnants in the cosmological simulation to z=0, we find that two dSph-like objects remain isolated at distances larger than 800 kpc from either the Milky Way or M31. These systems constitute plausible counterparts of the remote dSphs Cetus and Tucana which reside in the LG outskirts, far from the tidal influence of the primary galaxies. We conclude that merging of rotationally-supported dwarfs represents a viable mechanism for the formation of dSphs in the LG and similar environments.

Using collisionless N-body simulations, we investigate the tidal evolution of late-type, rotationally supported dwarfs inside Milky Way-sized host galaxies. Our study focuses on a wide variety of dwarf orbital configurations and initial structures. During the evolution, the disky dwarfs undergo strong mass loss, the stellar disks are transformed into spheroids, and rotation is replaced by random motions of the stars. Thus, the late-type progenitors are transformed into early-type dwarfs as envisioned by the tidal stirring model for the formation of dwarf spheroidal (dSph) galaxies in the Local Group. We determine the photometric properties of the dwarfs, including the total visual magnitude, the half-light radius and the central surface brightness as they would be measured by an observer near the galactic center. Special emphasis is also placed on studying their kinematics and shapes. We demonstrate that the measured values are biased by a number of observational effects including the increasing angle of the observation cone near the orbital pericenter, the fact that away from the pericenter the tidal tails are typically oriented along the line of sight, and the fact that for most of the evolution the stellar components of the dwarfs are triaxial ellipsoids whose major axis tumbles with respect to the line of sight. Finally, we compare the measured properties of the simulated dwarfs to those of dwarf galaxies in the Local Group. The evolutionary tracks of the dwarfs in different parameter planes and the correlations between their different properties, especially the total magnitude and the surface brightness, strongly suggest that present-day dSph galaxies may have indeed formed from late-type progenitors as proposed by the tidal stirring scenario.

The tidal stirring model envisions the formation of dwarf spheroidal (dSph) galaxies in the Local Group via the tidal interaction of disky dwarf systems with a larger host galaxy like the Milky Way. These progenitor disks are embedded in extended dark halos and during the evolution both components suffer strong mass loss. In addition, the disks undergo the morphological transformation into spheroids and the transition from ordered to random motion of their stars. Using collisionless N-body simulations we construct a model for the nearby and highly elongated Sagittarius (Sgr) dSph galaxy within the framework of the tidal stirring scenario. Constrained by the present known orbit of the dwarf, the model suggests that in order to produce the majority of tidal debris observed as the Sgr stream, but not yet transform the core of the dwarf into a spherical shape, Sgr must have just passed the second pericenter of its current orbit around the Milky Way. In the model, the stellar component of Sgr is still very elongated after the second pericenter and morphologically intermediate between the strong bar formed at the first pericenter and the almost spherical shape existing after the third pericenter. This is thus the first model of the evolution of the Sgr dwarf that accounts for its observed very elliptical shape. At the present time there is very little intrinsic rotation left and the velocity gradient detected along the major axis is almost entirely of tidal origin. We model the recently measured velocity dispersion profile for Sgr assuming that mass traces light and estimate its current total mass within 5 kpc to be 5.2 x 10^8 M_sun. To have this mass at present, the model requires that the initial virial mass of Sgr must have been as high as 1.6 x 10^10 M_sun, comparable to that of the Large Magellanic Cloud, which may serve as a suitable analog for the pre-interaction, Sgr progenitor.

Most of dwarf spheroidal galaxies in the Local Group were probably formed via environmental processes like the tidal interaction with the Milky Way. We study this process via N-body simulations of dwarf galaxies evolving on seven different orbits around the Galaxy. The dwarf galaxy is initially composed of a rotating stellar disk and a dark matter halo. Due to the action of tidal forces it loses mass and the disk gradually transforms into a spheroid while stellar motions become increasingly random. We measure the characteristic scale-length of the dwarf, its maximum circular velocity, mass, shape and kinematics as a function of the integrated tidal force along the orbit. The final properties of the evolved dwarfs are remarkably similar if the total tidal force they experienced was the same, independently of the actual size and eccentricity of the orbit.

(Abridged) The tidal stirring model posits the formation of dSph galaxies via the tidal interactions between rotationally-supported dwarfs and MW-sized host galaxies. Using a set of collisionless N-body simulations, we investigate the efficiency of the tidal stirring mechanism. We explore a wide variety of dwarf orbital configurations and initial structures and demonstrate that in most cases the disky dwarfs experience significant mass loss and their stellar components undergo a dramatic morphological and dynamical transformation: from disks to bars and finally to pressure-supported spheroidal systems with kinematic and structural properties akin to those of the classic dSphs in the Local Group (LG). Our results suggest that such tidal transformations should be common occurrences within the currently favored cosmological paradigm and highlight the key factor responsible for an effective metamorphosis to be the strength of the tidal shocks at the pericenters of the orbit. We demonstrate that the combination of short orbital times and small pericenters, characteristic of dwarfs being accreted at high redshift, induces the strongest transformations. Our models also indicate that the transformation efficiency is affected significantly by the structure of the progenitor disky dwarfs. Lastly, we find that the dwarf remnants satisfy the relation Vmax = \sqrt{3} * sigma, where sigma is the 1D, central stellar velocity dispersion and Vmax is the maximum halo circular velocity, with intriguing implications for the missing satellites problem. Overall, we conclude that the action of tidal forces from the hosts constitutes a crucial evolutionary mechanism for shaping the nature of dwarf galaxies in environments such as that of the LG. Environmental processes of this type should thus be included as ingredients in models of dwarf galaxy formation and evolution.

We analyze kinematic data of 41 nearby (z<0.1) relaxed galaxy clusters in terms of the projected phase-space density using a phenomenological, fully anisotropic model of the distribution function. We apply the Markov Chain Monte Carlo approach to place constraints on total mass distribution approximated by the universal NFW profile and the profile of the anisotropy of galaxy orbits. We find the normalization of the mean mass-concentration relation is c=6.9_{-0.7}^{+0.6} at the virial mass M_v=5x10^{14}M_sun. Assuming a one-to-one correspondence between sigma_8 and the normalization of the mass-concentration relation in the framework of the concordance model we estimate the normalization of the linear power spectrum to be sigma_8=0.91_{-0.08}^{+0.07}. Our constraints on the parameters of the mass profile are compared with estimates from other methods. We show that galaxy orbits are isotropic at the cluster centres (with the mean ratio of the radial-to-tangential velocity dispersions sigma_r/sigma_theta=0.97+/-0.04) and radially anisotropic at the virial sphere (with the mean ratio sigma_r/sigma_theta=1.75^{+0.23}_{-0.19}). Although the value of the central anisotropy appears to be universal, the anisotropy at the virial radius differs between clusters within the range 1<(sigma_r/sigma_theta)<2. Utilizing the Bautz-Morgan morphological classification and information on the prominence of a cool core we select two subsamples of galaxy clusters corresponding to less and more advanced evolutionary states. It is demonstrated that less evolved clusters have shallower mass profiles and their galaxy orbits are more radially biased at the virial sphere. This property is consistent with the expected evolution of the mass profiles as well as with the observed orbital segregation of late and early type galaxies.

Using a high resolution N-body simulation of a two-component dwarf galaxy orbiting in the potential of the Milky Way, we study two effects that lead to significant biases in mass estimates of dwarf spheroidal galaxies. Both are due to the strong tidal interaction of initially disky dwarfs with their host. The tidal stripping of dwarf stars leads to the formation of strong tidal tails that are typically aligned with the line of sight of an observer positioned close to the host center. The stars from the tails contaminate the kinematic samples leading to a velocity dispersion profile increasing with the projected radius and resulting in an overestimate of mass. The tidal stirring of the dwarf leads to the morphological transformation of the initial stellar disk into a bar and then a spheroid. The distribution of stars in the dwarf remains non-spherical for a long time leading to an overestimate of its mass if it is observed along the major axis and an underestimate if it seen in the perpendicular direction.

We use a simulation performed within the Constrained Local UniversE Simulation (CLUES) project to study a realistic Local Group-like object. We employ this group as a numerical laboratory for studying the evolution of the population of its subhaloes from the point of view of the effects it may have on the origin of different types of dwarf galaxies. We focus on the processes of tidal stripping of the satellites, their interaction, merging and grouping before infall. The tidal stripping manifests itself in the transition between the phase of mass accretion and mass loss seen in most subhaloes, which occurs at the moment of infall on to the host halo, and the change of the shape of their mass function with redshift. Although the satellites often form groups, they are loosely bound within them and do not interact with each other. The infall of a large group could however explain the observed peculiar distribution of the Local Group satellites, but only if it occurred recently. Mergers between prospective subhaloes are significant only during an early stage of evolution, i.e. more than 7 Gyr ago, when they are still outside the host haloes. Such events could thus contribute to the formation of more distant early type Milky Way companions. Once the subhaloes enter the host halo the mergers become very rare.

Using high-resolution N-body simulations we study the stellar properties of dwarf spheroidal galaxies resulting from the tidally induced morphological transformation of disky dwarfs on a cosmologically motivated eccentric orbit around the Milky Way. Dwarf galaxy models initially consist of an exponential stellar disk embedded in an extended spherical dark matter halo. Depending on the initial orientation of the disk with respect to the orbital plane, different final configurations are obtained. The least evolved dwarf is triaxial and retains a significant amount of rotation. The more evolved dwarfs are prolate spheroids with little rotation. We show that the final density distribution of stars can be approximated by a simple modification of the Plummer law. The kinematics of the dwarfs is significantly different depending on the line of sight which has important implications for mapping the observed stellar velocity dispersions of dwarfs to subhalo circular velocities. When the dwarfs are observed along the long axis, the measured velocity dispersion is higher and decreases faster with radius. In the case where rotation is significant, when viewed perpendicular to the long axis, the effect of minor axis rotation is detected, as expected for triaxial systems. We model the velocity dispersion profiles and rotation curves of the dwarfs by solving the Jeans equations for spherical and axisymmetric systems and adjusting different sets of free parameters. We find that the mass is typically overestimated when the dwarf is seen along the long axis and underestimated when the observation is along the short or intermediate axis. The effect of non-sphericity cannot however bias the inferred mass by more than 60 percent in either direction, even for the most strongly stripped dwarf which is close to disruption.

Using high-resolution collisionless N-body simulations we study the properties of tidal tails formed in the immediate vicinity of a two-component dwarf galaxy evolving in a static potential of the Milky Way (MW). The stellar component of the dwarf is initially in the form of a disk and the galaxy is placed on an eccentric orbit motivated by CDM-based cosmological simulations. We measure the orientation, density and velocity distribution of the stars in the tails. Due to the geometry of the orbit, in the vicinity of the dwarf, where the tails are densest and therefore most likely to be detectable, they are typically oriented towards the MW and not along the orbit. We report on an interesting phenomenon of `tidal tail flipping': on the way from the pericentre to the apocentre the old tails following the orbit are dissolved and new ones pointing towards the MW are formed over a short timescale. We also find a tight linear relation between the velocity of stars in the tidal tails and their distance from the dwarf. Using mock data sets we demonstrate that if dwarf spheroidal (dSph) galaxies in the vicinity of the MW are tidally affected their kinematic samples are very likely contaminated by tidally stripped stars which tend to artificially inflate the measured velocity dispersion. The effect is stronger for dwarfs on their way from the peri- to the apocentre due to the formation of new tidal tails after pericentre. Realistic mass estimates of dSph galaxies thus require removal of these stars from kinematic samples.

We present a new method of constraining the mass and velocity anisotropy profiles of galaxy clusters from kinematic data. The method is based on a model of the phase space density which allows the anisotropy to vary with radius between two asymptotic values. The characteristic scale of transition between these asymptotes is fixed and tuned to a typical anisotropy profile resulting from cosmological simulations. The model is parametrized by two values of anisotropy, at the centre of the cluster and at infinity, and two parameters of the NFW density profile, the scale radius and the scale mass. In order to test the performance of the method in reconstructing the true cluster parameters we analyze mock kinematic data for 20 relaxed galaxy clusters generated from a cosmological simulation of the standard LCDM model. We use Bayesian methods of inference and the analysis is carried out following the Markov Chain Monte Carlo approach. The parameters of the mass profile are reproduced quite well, but we note that the mass is typically underestimated by 15 percent, probably due to the presence of small velocity substructures. The constraints on the anisotropy profile for a single cluster are in general barely conclusive. Although the central asymptotic value is determined accurately, the outer one is subject to significant systematic errors caused by substructures at large clustercentric distance. The anisotropy profile is much better constrained if one performs joint analysis of at least a few clusters. In this case it is possible to reproduce the radial variation of the anisotropy over two decades in radius inside the virial sphere.

We conduct high-resolution collisionless N-body simulations to investigate the tidal evolution of dwarf galaxies on an eccentric orbit in the Milky Way (MW) potential. The dwarfs originally consist of a low surface brightness stellar disk embedded in a cosmologically motivated dark matter halo. During 10 Gyr of dynamical evolution and after 5 pericentre passages the dwarfs suffer substantial mass loss and their stellar component undergoes a major morphological transformation from a disk to a bar and finally to a spheroid. The bar is preserved for most of the time as the angular momentum is transferred outside the galaxy. A dwarf spheroidal (dSph) galaxy is formed via gradual shortening of the bar. This work thus provides a comprehensive quantitative explanation of a potentially crucial morphological transformation mechanism for dwarf galaxies that operates in groups as well as in clusters. We compare three cases with different initial inclinations of the disk and find that the evolution is fastest when the disk is coplanar with the orbit. Despite the strong tidal perturbations and mass loss the dwarfs remain dark matter dominated. For most of the time the 1D stellar velocity dispersion, \sigma, follows the maximum circular velocity, V_{\rm max}, and they are both good tracers of the bound mass. Specifically, we find that M_{\rm bound} \propto V_{\rm max}^{3.5} and V_{\rm max} \sim \sqrt{3} \sigma in agreement with earlier studies based on pure dark matter simulations. The latter relation is based on directly measuring the stellar kinematics of the simulated dwarf and may thus be reliably used to map the observed stellar velocity dispersions of dSphs to halo circular velocities when addressing the missing satellites problem.

We model the large kinematic data sets for the four Milky Way dwarf spheroidal (dSph) satellites: Carina, Fornax, Sculptor and Sextans, recently published by Walker et al. The member stars are selected using a reliable dynamical interloper removal scheme tested on simulated data. Our member selection is more restrictive than the one based on metallicity indicators as it removes not only contamination due to Milky Way stars but also the unbound stars from the tidal tails. We model the cleaned data sets by adjusting the solutions of the Jeans equations to the profiles of the projected velocity dispersion and kurtosis. The data are well reproduced by models where mass follows light and the best-fitting stellar orbits are isotropic to weakly tangential, as expected from the tidal stirring scenario. The Fornax dwarf, with more than 2400 member stars, is a dSph galaxy with the most accurately determined mass to date: its 1 sigma error following from the sampling errors of the velocity moments is below 5 percent. With mass-to-light ratio of 97 solar units, Sextans seems to be the most dark matter dominated of the four dSph galaxies.

We model the recently published kinematic data set for Leo I dSph galaxy by fitting the solutions of the Jeans equations to the velocity dispersion and kurtosis profiles measured from the data. We demonstrate that when the sample is cleaned of interlopers the data are consistent with the assumption that mass follows light and isotropic stellar orbits with no need for an extended dark matter halo. Our interloper removal scheme does not clean the data of contamination completely, as demonstrated by the rotation curve of Leo I. When moving away from the centre of the dwarf, the rotation appears to be reversed. We interpret this behaviour using the results of an N-body simulation of a dwarf galaxy possessing some intrinsic rotation, orbiting in the Milky Way potential and show that it can be reproduced if the galaxy is viewed almost along the tidal tails so that the leading (background) tail contaminates the western part of Leo I while the trailing (foreground) tail the eastern one. We show that this configuration leads to a symmetric and Gaussian distribution of line-of-sight velocities. The simulation is also applied to test our modelling method on mock data sets. We demonstrate that when the data are cleaned of interlopers and the fourth velocity moment is used the true parameters of the dwarf are typically within 1 \sigma errors of the best-fitting parameters. Restricting the fitting to the inner part of Leo I our best estimate for the anisotropy is \beta = -0.2^{+0.3}_{-0.4} and the total mass M = (4.5 +/- 0.7) x 10^7 M_sun. The mass-to-light ratio including the errors in mass, brightness and distance is M/L_V = 8.2 +/- 4.5 solar units.

We study the distribution function (DF) of dark matter particles in haloes of mass range 10^{14}--10^{15}\Msun. In the numerical part of this work we measure the DF for a sample of relaxed haloes formed in the simulation of a standard \LambdaCDM model. The DF is expressed as a function of energy E and the absolute value of the angular momentum L, a form suitable for comparison with theoretical models. By proper scaling we obtain the results that do not depend on the virial mass of the haloes. We demonstrate that the DF can be separated into energy and angular momentum components and propose a phenomenological model of the DF in the form f_{E}(E)[1+L^{2}/(2L_{0}^{2})]^{-\beta_{\infty}+\beta_{0}}L^{-2\beta_{0}}. This formulation involves three parameters describing the anisotropy profile in terms of its asymptotic values (\beta_{0} and \beta_{\infty}) and the scale of transition between them (L_{0}). The energy part f_{E}(E) is obtained via inversion of the integral for spatial density. We provide a straightforward numerical scheme for this procedure as well as a simple analytical approximation for a typical halo formed in the simulation. The DF model is extensively compared with the simulations: using the model parameters obtained from fitting the anisotropy profile, we recover the DF from the simulation as well as the profiles of the dispersion and kurtosis of radial and tangential velocities. Finally, we show that our DF model reproduces the power-law behaviour of phase space density Q=\rho(r)/\sigma^{3}(r).

NGC 5128 is a well-studied elliptical galaxy with an excellent kinematic data set for planetary nebulae which has been modelled up till now only by solving the Jeans equation for spherical systems. As a first approximation beyond spherical symmetry we model the galaxy as an axisymmetric system flattened by rotation with isotropic velocity distribution. We propose a new version of such an isotropic rotator having cuspy density profile with NFW-like behavior. The solutions of the Jeans equations for a single component of such form do not reproduce the data well: the rotation curve rises too slowly with radius and the velocity dispersion profile drops too fast. The data are well fitted however by a system built with two components: a more compact, flattened, less massive, fast-rotating and cold `bulge' and an extended, almost spherical, more massive, slow-rotating and hot `halo'. This picture agrees well with the results of recent N-body simulations of galaxy mergers which tend to produce oblate rotating spheroids. The total mass of the system is estimated to be (9.1 +/- 3.5) x 10^11 solar masses and the mass-to-light ratio is only 26 +/- 17 solar units.

Using an example of the Abell 2634 galaxy cluster we discuss the effect of contamination of kinematic data by interlopers and its impact on mass modelling. The cluster data reveal rich substructure along the line of sight. We demonstrate that it is necessary to apply a few independent methods of interloper removal in order to obtain a reliable sample of cluster members. We present results of three such methods which are commonly used in the literature and have been recently extensively tested on simulated data. Only two of them lead to consistent and reliable samples of cluster galaxies. For both of them we provide parameters of the best-fitting NFW density profile by fitting an isotropic solution of the Jeans equation to the velocity dispersion profiles.

We discuss how different approaches to selecting member stars in kinematic samples of dwarf spheroidal galaxies affect the estimates of their mass and anisotropy of stellar orbits. We demonstrate that the selection of members is an additional source of error compared to the usual uncertainties due to the sampling of velocity moments. As an example we use the kinematic data set for 202 stars in the Fornax dwarf galaxy for which we model the velocity dispersion profile and estimate the mass-to-light ratio and anisotropy assuming that mass follows light. We also show that stronger constraints on these parameters can be obtained if kurtosis of the velocity distribution is included in the analysis. Using the Besancon model of the Milky Way we demonstrate that the majority of contamination in Fornax probably comes from the Milky Way stars.

We summarize the method of mass modelling of galaxy clusters based on reproducing the dispersion and kurtosis of the projected velocity distribution of galaxies. The models are parametrized within the framework of the NFW density profile, characterized by the virial mass and concentration, together with the constant anisotropy of galaxy orbits. The use of velocity dispersion alone does not allow to constrain all the three parameters from kinematic data due to the mass-anisotropy degeneracy. The degeneracy is broken by introducing the fourth velocity moment, the kurtosis. We tested the method based on fitting both moments on mock data sets drawn from simulated dark matter haloes and showed it to reproduce reliably the properties of the haloes. The method has been applied to estimate the mass, concentration and anisotropy of more than 20 clusters which allowed us to confirm, for the first time using kinematic data, the mass-concentration relation found in N-body simulations.

A new revision of the gamma flux that we expect to detect in Imaging Atmospheric Cherenkov Telescopes (IACTs) from neutralino annihilation in the Draco dSph is presented in the context of the minimal supersymmetric standard models (MSSM) compatible with the present phenomenological and cosmological constraints, and using the dark matter (DM) density profiles compatible with the latest observations. This revision takes also into account the important effect of the Point Spread Function (PSF) of the telescope, and is valid not only for Draco but also for any other DM target. We show that this effect is crucial in the way we will observe and interpret a possible signal detection. Finally, we discuss the prospects to detect a possible gamma signal from Draco for current or planned gamma-ray experiments, i.e. MAGIC, GLAST and GAW. Even with the large astrophysical and particle physics uncertainties we find that the chances to detect a neutralino annihilation signal in Draco seem to be very scarce for current experiments. However, the prospects for future IACTs with upgraded performances (especially lower threshold energies and higher sensitivities) such as those offered by the CTA project, might be substantially better.

We study the origin and properties of the population of unbound stars in the kinematic samples of dwarf spheroidal galaxies. For this purpose we have run a high resolution N-body simulation of a two-component dwarf galaxy orbiting in a Milky Way potential. The dwarf is placed on a highly eccentric orbit, its initial stellar component is in the form of an exponential disk and it has a NFW-like dark matter halo. After 10 Gyrs of evolution the dwarf produces a spheroidal stellar component and is strongly tidally stripped so that mass follows light and the stars are on almost isotropic orbits. From this final state, we create mock kinematic data sets for 200 stars by observing the dwarf in different directions. We find that when the dwarf is observed along the tidal tails the kinematic samples are strongly contaminated by unbound stars from the tails. We also study another source of possible contamination by adding stars from the Milky Way. We demonstrate that most of the unbound stars can be removed by the method of interloper rejection proposed by den Hartog & Katgert. We model the cleaned-up kinematic samples using solutions of the Jeans equation with constant mass-to-light ratio and velocity anisotropy parameter. We show that the Jeans analysis, when applied to cleaned samples, allows us to reproduce the mass and mass-to-light ratio of the dwarf with accuracy typically better than 25 percent. The analysis was applied to the new data for the Fornax dSph galaxy. We show that after careful removal of interlopers the velocity dispersion profile of Fornax can be reproduced by a model in which mass traces light with a mass-to-light ratio of 11 solar units and isotropic orbits. We demonstrate that most of the contamination in the kinematic sample of Fornax probably originates from the Milky Way. (Abridged)

We study the effect of contamination by interlopers in kinematic samples of galaxy clusters. We demonstrate that without the proper removal of interlopers the inferred parameters of the mass distribution in the cluster are strongly biased towards higher mass and lower concentration. The interlopers are removed using two procedures previously shown to work most efficiently on simulated data. One is based on using the virial mass estimator and calculating the maximum velocity available to cluster members and the other relies on the ratio of the virial and projected mass estimators. We illustrate the performance of the methods in detail using the example of A576, a cluster with a strong uniform background contamination, and compare the case of A576 to 15 other clusters with different degree of contamination. We model the velocity dispersion and kurtosis profiles obtained for the cleaned data samples of these clusters solving the Jeans equations to estimate the mass, concentration and anisotropy parameter. We present the mass-concentration relation for the total sample of 22 clusters.

The aim of this paper is to study the efficiency of different approaches to interloper treatment in dynamical modelling of galaxy clusters. Using cosmological N-body simulation of standard LCDM model, we select 10 massive dark matter haloes and use their particles to emulate mock kinematic data in terms of projected galaxy positions and velocities as they would be measured by a distant observer. Taking advantage of the full 3D information available from the simulation, we select samples of interlopers defined with different criteria. The interlopers thus selected provide means to assess the efficiency of different interloper removal schemes found in the literature. We study direct methods of interloper removal based on dynamical or statistical restrictions imposed on ranges of positions and velocities available to cluster members. In determining these ranges, we use either the velocity dispersion criterion or a maximum velocity profile. We also generalize the common approaches taking into account both the position and velocity information. Another criterion is based on the dependence of the commonly used virial mass and projected mass estimators on the presence of interlopers. We find that the direct methods exclude on average 60-70 percent of unbound particles producing a sample with contamination as low as 2-4 percent. Next, we consider indirect methods of interloper treatment which are applied to the data stacked from many objects. In these approaches, interlopers are treated in a statistical way as a uniform background which modifies the distribution of cluster members. Using a Bayesian approach, we reproduce the properties of composite clusters and estimate the probability of finding an interloper as a function of distance from the object centre.

A new revision of the gamma flux that we expect to detect in Imaging Atmospheric Cherenkov Telescopes (IACTs) from SUSY dark matter annihilation in the Draco dSph is presented using the dark matter density profiles compatible with the latest observations. This revision takes also into account the important effect of the Point Spread Function (PSF) of the Cherenkov telescope. We show that this effect is crucial in the way we will observe and interpret a possible signal profile in the telescope. Given these new considerations, some light can be shed on the recent signal excess reported by the CACTUS experiment.

Elements of kinematical and dynamical modeling of elliptical galaxies (Es) are presented. In projection, NFW models resemble Sersic models, but with a very narrow range of shapes (m=3+/-1). The total density profile of Es cannot be NFW-like because the predicted local M/L and aperture velocity dispersion within an effective radius (Re) are much lower than observed. Stars must then dominate Es out to a few Re. Fitting an NFW model to the total density profile of Sersic+NFW (stars+dark matter [DM]) Es results in very high concentration parameters, as found by X-ray observers. Kinematical modeling of Es assuming an isotropic NFW DM model underestimates M/L at the virial radius by a factor of 1.6 to 2.4, because dissipationless LCDM halos have slightly different density profiles and slightly radial velocity anisotropy. In N-body+gas simulations of Es as merger remnants of spirals embedded in DM halos, the slope of the DM density profile is steeper when the initial spiral galaxies are gas-rich. The Hansen & Moore (2006) relation between anisotropy and the slope of the density profile breaks down for gas and DM, but the stars follow an analogous relation with slightly less radial anisotropies for a given density slope. Using kurtosis (h_4) to infer anisotropy in Es is dangerous, as h_4 is also sensitive to small levels of rotation. The stationary Jeans equation provides accurate masses out to 8 Re. The discrepancy between the modeling of Romanowsky et al. (2003), indicating a dearth of DM in Es, and the simulations analyzed by Dekel et al. (2005), which match the spectroscopic observations of Es, is partly due to radial anisotropy and to observing oblate Es face-on. However, one of the 15 solutions to the orbit modeling of Romanowsky et al. is found to have an amount and concentration of DM consistent with LCDM predictions.

We study the mass distribution in six nearby (z<0.06) relaxed Abell clusters of galaxies A0262, A0496, A1060, A2199, A3158 and A3558. Given the dominance of dark matter in galaxy clusters we approximate their total density distribution by the NFW formula characterized by virial mass and concentration. We also assume that the anisotropy of galactic orbits is reasonably well described by a constant and that galaxy distribution traces that of the total density. Using the velocity and position data for 120-420 galaxies per cluster we calculate, after removal of interlopers, the profiles of the lowest-order even velocity moments, dispersion and kurtosis. We then reproduce the velocity moments by jointly fitting the moments to the solutions of the Jeans equations. Including the kurtosis in the analysis allows us to break the degeneracy between the mass distribution and anisotropy and constrain the anisotropy as well as the virial mass and concentration. The method is tested in detail on mock data extracted from N-body simulations of dark matter haloes. We find that the best-fitting galactic orbits are remarkably close to isotropic in most clusters. Using the fitted pairs of mass and concentration parameters for the six clusters we conclude that the trend of decreasing concentration for higher masses found in cosmological N-body simulations is consistent with the data. By scaling the individual cluster data by mass we combine them to create a composite cluster with 1465 galaxies and perform a similar analysis on such sample. The estimated concentration parameter then lies in the range 1.5 < c < 14 and the anisotropy parameter in the range -1.1 < \beta < 0.5 at the 95 percent confidence level.

The Abell 1689 galaxy cluster has recently become a subject of intensive study. Thanks to its intermediate redshift (z=0.183) its mass distribution can be reconstructed using numerous methods including gravitational lensing, galaxy kinematics and X-ray imaging. The methods used to yield conflicting mass estimates in the past and recently the cluster mass distribution has been claimed to be in conflict with standard CDM scenarios due to rather large concentration and steep mass profile obtained from detailed studies of Broadhurst et al. using lensing. By studying in detail the kinematics of about 200 galaxies with measured redshifts in the vicinity of the cluster we show that the cluster is probably surrounded by a few structures, quite distant from each other, but aligned along the line of sight. We support our arguments by referring to cosmological N-body simulations and showing explicitly that distant, non-interacting haloes can produce entangled multi-peak line-of-sight velocity distributions similar to that in A1689. We conclude that it is difficult to estimate the cluster mass reliably from galaxy kinematics, but the value we obtain after applying a simple cut-off in velocity agrees roughly with the mass estimated from lensing. The complicated mass distribution around the cluster may however increase the uncertainty in the determination of the density profile shape obtained with weak lensing.

We present interpretations of the line-of-sight velocity distribution of stars in the Draco dwarf spheroidal galaxy in terms of CDM and MOND assuming constant mass-to-light ratio and anisotropy. We estimate the two parameters by fitting the line-of-sight velocity dispersion and kurtosis profiles for stellar samples differing by a number of stars rejected as interlopers. The results of the fitting procedure for CDM, high mass-to-light ratio (131-141 solar units in V-band) and weakly tangential orbits, are similar for different samples, but the quality of the fit is much worse when fewer interlopers are removed. For MOND, the derived mass-to-light ratio (21 solar units) is too large to be explained by the stellar content of the galaxy.

Using cosmological N-body simulations we study the line-of-sight velocity distribution of dark matter haloes focusing on the lowest-order even moments, dispersion and kurtosis, and their application to estimate the mass profiles of cosmological structures. For each of the ten massive haloes selected from the simulation box we determine the virial mass, concentration and the anisotropy parameter. In order to emulate observations from each halo we choose randomly 300 particles and project their velocities and positions along the line of sight and on the surface of the sky, respectively. After removing interlopers we calculate the profiles of the line-of-sight velocity moments and fit them with the solutions of the Jeans equations. The estimates of virial mass, concentration parameter and velocity anisotropy obtained in this way are in good agreement with the values found from the full 3D analysis.

We study the distribution of dark matter in the Draco dwarf spheroidal galaxy by modelling the moments of the line-of-sight velocity distribution of stars obtained from new velocity data of Wilkinson et al. The luminosity distribution is approximated by a Sersic profile fitted to the data by Odenkirchen et al. We assume that the dark matter density profile is given by a formula with an inner cusp and an outer exponential cut-off, as recently proposed by Kazantzidis et al. as a result of simulations of tidal stripping of dwarfs by the potential of the Milky Way. The dark matter distribution is characterized by the total dark mass and the cut-off radius. The models have arbitrary velocity anisotropy parameter assumed to be constant with radius. We estimate the three parameters by fitting both the line-of-sight velocity dispersion and kurtosis profiles, which allows us to break the degeneracy between the mass distribution and velocity anisotropy. The results of the fitting procedure turn out to be very different depending on the stellar sample considered, that is on our choice of stars with discrepant velocities to be discarded as interlopers. For our most reliable sample, the model parameters remain weakly constrained, but the robust result is the preference for weakly tangential stellar orbits and high mass-to-light ratios. The best-fitting total mass is then 7 10^7 M_sun, much lower than recent estimates, while the mass-to-light ratio is M/L_V = 300 and almost constant with radius. If the binary fraction in the stellar population of Draco turns out to be significant, the kurtosis of the global velocity distribution will be smaller and the orbits inferred will be more tangential, while the resulting mass estimate lower.

Using cosmological N-body simulations we study the radial velocity distribution in dark matter haloes focusing on the lowest-order even moments, dispersion and kurtosis. We determine the properties of ten massive haloes in the simulation box approximating their density distribution by the NFW formula characterized by the virial mass and concentration. We also calculate the velocity anisotropy parameter of the haloes and find it mildly radial and increasing with distance from the halo centre. The radial velocity dispersion of the haloes shows a characteristic profile with a maximum, while the radial kurtosis profile decreases with distance starting from a value close to Gaussian near the centre. We therefore confirm that dark matter haloes possess intrinsically non-Gaussian, flat-topped velocity distributions. We find that the radial velocity moments of the simulated haloes are quite well reproduced by the solutions of the Jeans equations obtained for the halo parameters with the anisotropy measured in the simulations. We also study the radial velocity moments for a composite cluster made of ten haloes out to ten virial radii. In this region the velocity dispersion decreases systematically to reach the value of the background, while kurtosis increases from below to above the Gaussian value of 3 signifying a transition from a flat-topped to a strongly peaked velocity distribution with respect to the Gaussian, which can be interpreted as the dominance of ordered flow with a small dispersion. We illustrate the transition by showing explicitly the velocity distribution of the composite cluster in a few radial bins.

We provide the first direct lifting of the mass/anisotropy degeneracy for a cluster of galaxies, by jointly fitting the line of sight velocity dispersion and kurtosis profiles of the Coma cluster, assuming an NFW tracer density profile, a generalized-NFW dark matter profile and a constant anisotropy profile. We find that the orbits in Coma must be quasi-isotropic, and find a mass consistent with previous analyses, but a concentration parameter 50% higher than expected in cosmological N-body simulations. We then test the accuracy of our method on realistic non-spherical systems with substructure and streaming motions, by applying it to the ten most massive structures in a cosmological N-body simulation. We find that our method yields fairly accurate results on average (within 20%), although with a wide variation (factor 1.7 at 1 sigma) for the concentration parameter, with decreased accuracy and efficiency when the projected mean velocity is not constant with radius.

We use $N$-body simulations to measure mass functions in flat cosmological models with quintessence characterized by constant $w$ with $w=-1$, -2/3 and -1/2. The results are compared to the predictions of the formula proposed by Jenkins et al. at different redshifts, in terms of FOF masses as well as Abell masses appropriate for direct comparison to observations. The formula reproduces quite well the mass functions of simulated haloes in models with quintessence. We use the cluster mass function data at a number of redshifts from Carlberg et al. to constrain $\Omega_0$, $\sigma_8$ and $w$. The best fit is obtained in the limit $w \to 0$, but none of the values of $w$ in the considered range $-1 \le w < 0$ can actually be excluded. However, the adopted value of $w$ affects significantly the constraints in the $\Omega_0-\sigma_8$ plane. Taking into account the dependence on $w$ we find $\Omega_0=0.32 \pm 0.15$ and $\sigma_8=0.85_{-0.12}^{+0.38}$ (68% c.l.). Since less negative $w$ push the confidence regions toward higher $\Omega_0$ and lower $\sigma_8$ we conclude that relaxing the assumption of $w=-1$ typically made in such comparisons may resolve the discrepancy between recent cluster mass function results (yielding rather low $\Omega_0$ and high $\sigma_8$) and most other estimates. The fact that high $w$ values are preferred may however also point towards some unknown systematics in the data or the model with constant $w$ being inadequate.

Using cosmological $N$-body simulations of dark matter haloes, we study the effects of non-sphericity, substructure and streaming motions in reproducing the structure and internal kinematics of clusters of galaxies from kinematical measurements. Fitting an NFW model to the 3D density profile, we determine the virial mass, concentration parameter and velocity anisotropy of the haloes, and then calculate the profiles of projected velocity moments, as they would be measured by a distant observer. Using these mock data, we apply a Jeans analysis for spherical objects to reproduce the line-of-sight velocity dispersion and kurtosis and fit the three parameters. We find that the line-of-sight velocity dispersion and kurtosis profiles of a given halo can vary considerably with the angle of view of the observer. We show that the virial mass, concentration parameter and velocity anisotropy of the haloes can be reproduced satisfactorily independently of the halo shape, although the virial mass tends to be underestimated, the concentration parameter overestimated, while the recovered anisotropy is typically more tangential than the true one. The mass, concentration and velocity anisotropy of haloes are recovered with better precision when their mean velocity profiles are near zero.

Using high resolution N-body simulations we address the problem of emptiness of giant 20 Mpc/h diameter voids found in the distribution of bright galaxies. Are the voids filled by dwarf galaxies? Do cosmological models predict too many small dark matter haloes inside the voids? Can the problems of cosmological models on small scales be addressed by studying the abundance of dwarf galaxies inside voids? We find that voids in the distribution of 10^12 Msun/h haloes (expected galactic magnitudes ~ M_*) are almost the same as the voids in 10^11 Msun/h haloes. Yet, much smaller haloes with masses 10^9 Msun/h and circular velocities v_circ about 20 km/s readily fill the voids: there should be almost 1000 of these haloes in a 20 Mpc/h void. A typical void of diameter 20 Mpc/h contains about 50 haloes with v_circ > 50 km/s. The haloes are arranged in a pattern, which looks like a miniature Universe: it has the same structural elements as the large-scale structure of the galactic distribution of the Universe. There are filaments and voids; larger haloes are at the intersections of filaments. The only difference is that all masses are four orders of magnitude smaller. There is severe (anti)bias in the distribution of haloes, which depends on halo mass and on the distance from the centre of the void. Large haloes are more antibiased and have a tendency to form close to void boundaries. The mass function of haloes in voids is different from the ``normal'' mass function. It is much steeper for high masses resulting in very few M33-type galaxies (v_circ about 100 km/s). We present an analytical approximation for the mass function of haloes in voids.

We study velocity moments of elliptical galaxies in the Coma cluster using Jeans equations. The dark matter distribution in the cluster is modelled by a generalised formula based upon the results of cosmological N-body simulations. Its inner slope (cuspy or flat), concentration, and mass within the virial radius are kept as free parameters, as well as the velocity anisotropy, assumed independent of position. We show that the study of line-of-sight velocity dispersion alone does not allow to constrain the parameters. By a joint analysis of the observed profiles of velocity dispersion and kurtosis we are able to break the degeneracy between the mass distribution and velocity anisotropy. We determine the dark matter distribution at radial distances larger than 3% of the virial radius and we find that the galaxy orbits are close to isotropic. Due to limited resolution, different inner slopes are found to be consistent with the data and we observe a strong degeneracy between the inner slope $\alpha$ and concentration $c$: the best-fitting profiles have the two parameters related with $c=19 - 9.6 \alpha$. Our best-fitting NFW profile has concentration $c=9$, which is 50% higher than standard values found in cosmological simulations for objects of similar mass. The total mass within the virial radius of $2.9 h_{70}^{-1}$ Mpc is $1.4 \times 10^{15} h_{70}^{-1} M_{\sun}$ (with 30% accuracy), 85% of which is dark. At this distance from the cluster centre, the mass-to-light ratio in the blue band is $351 h_{70}$ solar units. The total mass within the virial radius leads to estimates of the density parameter of the Universe, assuming that clusters trace the mass-to-light ratio and baryonic fraction of the Universe, with $\Omega_0=0.29 \pm 0.1$.

I review the main characteristics of structure formation in the quintessential Universe. Assuming equation of state w=p/\varrho=$const I provide a brief description of the background cosmology and discuss the linear growth of density perturbations, the strongly nonlinear evolution, the power spectra and rms fluctuations as well as mass functions focusing on the three values w=-1, -2/3 and -1/3. Finally I describe the presently available and future constraints on w.

We generalize the spherical collapse model for the formation of bound objects to apply in a Universe with arbitrary positive cosmological constant. We calculate the critical condition for collapse of an overdense region and give exact values of the characteristic densities and redshifts of its evolution.

Using numerical simulations and the Sheth-Tormen approximation we study the mass function of dark matter halos in voids. We find that the void mass function is significantly lower and its shape is different than that of the field halos. We predict that in the standard LCDM model a void with radius 10 h^-1 Mpc should have 50 halos with circular velocity v_c > 50 km/s and 600 halos with v_c > 20 km/s.

We study the distribution of dark matter in dwarf spheroidal galaxies by modelling the moments of their line-of-sight velocity distributions. We discuss different dark matter density profiles, both cuspy and possessing flat density cores. The predictions are made in the framework of standard dynamical theory of two-component (stars and dark matter) spherical systems with different velocity distributions. We compare the predicted velocity dispersion profiles to observations in the case of Fornax and Draco dwarfs. For isotropic models the dark haloes with cores are found to fit the data better than those with cusps. Anisotropic models are studied by fitting two parameters, dark mass and velocity anisotropy, to the data. In this case all profiles yield good fits but the steeper the cusp of the profile, the more tangential is the velocity distribution required to fit the data. To resolve this well-known degeneracy of density profile versus velocity anisotropy we obtain predictions for the kurtosis of the line-of-sight velocity distribution for models found to provide best fits to the velocity dispersion profiles. It turns out that profiles with cores typically yield higher values of kurtosis which decrease more steeply with distance than the cuspy profiles, which will allow to discriminate between the profiles once the kurtosis measurements become available. We also show that with present quality of the data the alternative explanation of velocity dispersions in terms of Modified Newtonian Dynamics cannot yet be ruled out.

We present predictions for the line-of-sight velocity dispersion profiles of dwarf spheroidal galaxies and compare them to observations in the case of the Fornax dwarf. The predictions are made in the framework of standard dynamical theory of spherical systems with different velocity distributions. The stars are assumed to be distributed according to Sersic laws with parameters fitted to observations. We compare predictions obtained assuming the presence of dark matter haloes (with density profiles adopted from N-body simulations) versus those resulting from Modified Newtonian Dynamics (MOND). If the anisotropy of velocity distribution is treated as a free parameter observational data for Fornax are reproduced equally well by models with dark matter and with MOND. If stellar mass-to-light ratio of $1 M_{\sun}/L_{\sun}$ is assumed, the required mass of the dark halo is $1.5 \times 10^9 M_{\sun}$, two orders of magnitude bigger than the mass in stars. The derived MOND acceleration scale is $a_0 = 2.1 \times 10^{-8}$ cm/s$^2$. In both cases a certain amount of tangential anisotropy in the velocity distribution is needed to reproduce the shape of the velocity dispersion profile in Fornax.

We generalize the spherical collapse model for the formation of dark matter halos to apply in a universe with arbitrary positive cosmological constant. We calculate the critical condition for collapse of an overdense region and give exact values of the characteristic densities and redshifts of its evolution.

Using the standard dynamical theory of spherical systems, we calculate the properties of spherical galaxies and clusters whose density profiles obey the universal form first obtained in high resolution cosmological N-body simulations by Navarro, Frenk & White. We adopt three models for the internal kinematics: isotropic velocities, constant anisotropy and increasingly radial Osipkov-Merritt anisotropy. Analytical solutions are found for the radial dependence of the mass, gravitational potential, velocity dispersion, energy and virial ratio and we test their variability with the concentration parameter describing the density profile and amount of velocity anisotropy. We also compute structural parameters, such as half-mass radius, effective radius and various measures of concentration. Finally, we derive projected quantities, the surface mass density and line-of-sight as well as aperture velocity dispersion, all of which can be directly applied in observational tests of current scenarios of structure formation. On the mass scales of galaxies, if constant mass-to-light is assumed, the NFW surface density profile is found to fit well Hubble-Reynolds laws. It is also well fitted by Sersic R^(1/m) laws, for m ~ 3, but in a much narrower range of m and with much larger effective radii than are observed. Assuming in turn reasonable values of the effective radius, the mass density profiles imply a mass-to-light ratio that increases outwards at all radii.

Using the formalism of the spherical infall model the structure of collapsed and virialized dark halos is calculated for a variety of scale-free initial conditions. In spite of the scale-free cosmological nature of the problem, the collapse of individual objects is not self-similar. Unlike most of previous calculations the dynamics used here relies only on adiabatic invariants and not on self-similarity. The paper focuses on the structure of the innermost part of the collapsed halos and addresses the problem of central density cusps. The slopes of density profiles at 1% of virial radius are calculated for a variety of cosmological models and are found to vary with the mass of the halos and power spectrum of the initial conditions. The inner slopes range between r^-2.3 and r^-2 with the limiting case of r^-2 reached for the largest masses. The steep cusps found here correspond to the limiting case where all particles move on radial orbits. The introduction of angular momentum will make the density profile shallower. We expect this to resolve the discrepancy found between the calculated profiles and the ones found in high resolution N-body simulations, where the exponent ranges from -0.5 to -1.5. The robust prediction here is that collisionless gravitational collapse in an expanding universe is expected to form density cups and not halos with a core structure.

I apply the modified spherical infall model to test the dependence of the shape of density profiles of dark matter halos on redshift. The model predicts that the density profile of a halo of a given mass steepens with time and its concentration parameter changes like c~1/(1+z_0)^2.

I propose a modification of the spherical infall model for the evolution of density fluctuations with initially Gaussian probability distribution and scale-free power spectra. I introduce a generalized form of the initial density distribution around an overdense region and cut it off at half the inter-peak separation accounting in this way for the presence of neighbouring fluctuations. Contrary to the original predictions of the model, resulting density profiles within virial radii no longer have power-law dependence on the distance, but are well fitted by the universal formula of changing slope obtained as a result of N-body simulations. The profiles of halos of a given mass are in general flatter than the corresponding ones from the simulations, but the trend of steeper profiles for smaller masses is reproduced. The profiles of galaxy size objects are in general in better agreement with N-body simulations than those of larger ones.

I propose a modification of the spherical infall model for the evolution of density fluctuations with initially Gaussian probability distribution and scale-free power spectra in Einstein-de Sitter universe as developed by Hoffman & Shaham. I introduce a generalized form of the initial density distribution around an overdense region and cut it off at half the inter-peak separation accounting in this way for the presence of the neighbouring fluctuations. Contrary to the original predictions of Hoffman & Shaham the resulting density profiles within virial radii no longer have power-law shape but their steepness increases with distance. The profiles of halos of galactic mass are well fitted by the universal profile formula of changing slope obtained as a result of N-body simulations by Navarro, Frenk & White. The trend of steeper profiles for smaller masses and higher spectral indices is also reproduced. The agreement between the model and simulations is better for smaller masses and lower spectral indices which suggests that galaxies form mainly by accretion while formation of clusters involves merging.

We present results of the investigations of the statistical properties of a joint density and velocity divergence probability distribution function (PDF) in the mildly non-linear regime. For that purpose we use both perturbation theory results, extended here for a top-hat filter, and numerical simulations. In particular we derive the quantitative (complete as possible up to third order terms) and qualitative predictions for constrained averages and constrained dispersions -- which describe the nonlinearities and the stochasticity properties beyond the linear regime -- and compare them against numerical simulations. We find overall a good agreement for constrained averages; however, the agreement for constrained dispersions is only qualitative. Scaling relations for the Omega-dependence of these quantities are satisfactory reproduced. Guided by our analytical and numerical results, we finally construct a robust phenomenological description of the joint PDF in a closed analytic form. The good agreement of our formula with results of N-body simulations for a number of cosmological parameters provides a sound validation of the presented approach. Our results provide a basis for a potentially powerful tool with which it is possible to analyze galaxy survey data in order to test the gravitational instability paradigm beyond the linear regime and put useful constraints on cosmological parameters. In particular we show how the nonlinearity in the density--velocity relation can be used to break the so-called Omega-bias degeneracy in cosmic density-velocity comparisons.

Using third-order perturbation theory, we derive a relation between the mean divergence of the peculiar velocity given density and the density itself. Our calculations assume Gaussian initial conditions and are valid for Gaussian filtering of the evolved density and velocity fields. The mean velocity divergence turns out to be a third-order polynomial in the density contrast. We test the power spectrum dependence of the coefficients of the polynomial for scale-free and standard CDM spectra and find it rather weak. Over scales larger than about 5 megaparsecs, the scatter in the relation is small compared to that introduced by random errors in the observed density and velocity fields. The relation can be useful for recovering the peculiar velocity from the associated density field, and also for non-linear analyses of the anisotropies of structure in redshift surveys.

Using weakly nonlinear conditional PDF for the density field around an overdense region we find that the expected density contrast around a peak is always smaller while its rms fluctuation larger than in the linear case. We apply these results to the spherical model of collapse as developed by Hoffman & Shaham (1985) and find that in the case of open Universe the effect of weakly nonlinear interactions is to decrease the scale from which a peak gathers mass and the mass itself as well as to steepen the final profile of the virialized protoobject.

Using the two-point Edgeworth series up to second order we construct the weakly nonlinear conditional probability distribution function for the density field around an overdense region. This requires calculating the two-point analogues of the skewness parameter $S_{3}$. We test the dependence of the two-point skewness on distance from the peak for scale-free power spectra and Gaussian smoothing. The statistical features of such conditional distribution are given as the values obtained within linear theory corrected by the terms that arise due to weakly nonlinear evolution. The expected density around the peak is found to be always below the linear prediction while its rms fluctuation is always larger than in the linear case. We apply these results to the spherical model of collapse as developed by Hoffman & Shaham (1985) and find that in general the effect of weakly nonlinear interactions is to decrease the scale from which a peak gathers mass and therefore also the mass itself. In the case of open universe this results in steepening of the final profile of the virialized protoobject.

We rigorously derive weakly nonlinear relation between cosmic density and velocity fields up to third order in perturbation theory. The density field is described by the mass density contrast, $\de$. The velocity field is described by the variable $\te$ proportional to the velocity divergence, $\te = - f(\Omega)^{-1} H_0^{-1} \nabla\cdot\bfv$, where $f(\Omega) \simeq \Omega^{0.6}$, $\Omega$ is the cosmological density parameter and $H_0$ is the Hubble constant. Our calculations show that mean $\de$ given $\te$ is a third order polynomial in $\te$, $\lan \de \ran|_{\te} = a_1 \te + a_2 (\te^2 - \s_\te^2) + a_3 \te^3$. This result constitutes an extension of the formula $\lan \de \ran|_{\te} = \te + a_2 (\te^2 - \s_\te^2)$, found by Bernardeau~(1992) which involved second order perturbative solutions. Third order perturbative corrections introduce the cubic term. They also, however, cause the coefficient $a_1$ to depart from unity, in contrast with the linear theory prediction. We compute the values of the coefficients $a_p$ for scale-free power spectra, as well as for standard CDM, for Gaussian smoothing. The coefficients obey a hierarchy $a_3 \ll a_2 \ll a_1$, meaning that the perturbative series converges very fast. Their dependence on $\Omega$ is expected to be very weak. The values of the coefficients for CDM spectrum are in qualitative agreement with the results of N-body simulations by Ganon et al. (1996). The results provide a method for breaking the $\Omega$-bias degeneracy in comparisons of cosmic density and velocity fields such as IRAS-POTENT.

We present a series of N-body experiments which confirm the reality of the previrialization effect. We also use weakly nonlinear perturbative approach to study the phenomenon. These two approaches agree when the rms density contrast, $\sigma$, is small; more surprisingly, they remain in agreement when $\sigma \approx 1$. When the slope of the initial power spectrum is $n>-1$, nonlinear tidal interactions slow down the growth of density fluctuations and the magnitude of the suppression increases when $n$ (i.e. the relative amount of small scale power) is increased. For $n<-1$ we see an opposite effect: the fluctuations grow more rapidly than in linear theory. The transition occurs at $n=-1$ when the weakly nonlinear correction to $\sigma$ is close to zero and the growth rate is close to linear. Our results resolve recent controversy between two N-body studies of previrialization. Peebles (1990) assumed $n=0$ and found strong evidence in support of previrialization, while Evrard \& Crone (1992), who assumed $n=-1$, reached opposite conclusions. As we show here, the initial conditions with $n=-1$ are rather special because the nonlinear effects nearly cancel out for that particular spectrum. In addition to our calculations for scale-free initial spectra, we show results for a more realistic spectrum of Peacock \& Dodds (1994). Its slope near the scale usually adopted for normalization is close to $-1$, so $\sigma$ is close to linear. Our results retroactively justify linear normalization at 8$h^{-1}$ Mpc, while also demonstrating the danger and limitations of this practice.

An attractive and simple hypothesis for the formation of large-scale structure is that it developed by gravitational instability from primordial fluctuations with an initially Gaussian probability distribution. Non-linear gravitational evolution drives the distribution away from the Gaussian form, generating measurable skewness and kurtosis even when the variance of the fluctuations is much smaller than unity. We use perturbation theory to compute the kurtosis of the mass density field and the velocity divergence field that arises during the weakly non-linear evolution of initially Gaussian fluctuations. We adopt an Einstein--de~Sitter universe for the perturbative calculations, and we discuss the generalization to a universe of arbitrary $\Omega$. We obtain semi-analytic results for the case of scale-free, power-law spectra of the initial fluctuations and final smoothing of cosmic fields with a Gaussian filter. We also give an exact analytical formula for the dependence of the skewness of these fields on the power spectrum index. We show that the kurtosis decreases with the power spectrum index, and we compare our more accurate results for the kurtosis to previous estimates from Monte Carlo integrations. We also compare our results to values obtained from cosmological N-body simulations with power-law initial spectra. Measurements of the skewness and kurtosis parameters can be used to test the hypothesis that structure in the universe formed by gravitational instability from Gaussian initial conditions.