Researcher profile

Ewa L. Lokas

Ewa L. Lokas contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
45works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

45 published item(s)

preprint2022arXiv

Lopsided galactic disks in IllustrisTNG

A significant fraction of nearby late-type galaxies are lopsided. We study the asymmetry of the stellar component in a sample of well-resolved disky galaxies selected from the last snapshot of the Illustris TNG100 simulation based on their flatness and rotational support. Among 1912 disks, we identify 161 objects with significant asymmetry in terms of the m=1 Fourier mode of the stellar component within (1-2) stellar half-mass radii and describe their properties using three representative examples. The profiles of the m=1 mode typically increase with radius, and the corresponding phase is constant in the asymmetric region, signifying a global distortion. Following the evolution of the lopsided disks over time, we find that their history is rather uneventful and the occurrence of the asymmetry is fairly recent. Only about 1/3 of the lopsided disks experienced any strong interaction recently that could have led to the distortion of their shape: 24% were affected by a more massive object and 9% underwent a gas-rich merger. Still, a majority of lopsided disks show a significant increase in their recent star formation rate. The most frequent mechanism for the formation of lopsided disks thus seems to be asymmetric star formation probably related to gas accretion, although the distortions in the gas and stars are not strongly correlated. This picture is supported by the finding that the lopsided population on average contains more gas, has higher star formation rate, lower metallicity and bluer color than the remaining disks. These correlations are similar to those seen in real galaxies, even though the fraction of simulated lopsided disks (8%) is much lower than in observations (30%). The observed correlation between the presence of the asymmetry and a bar is not reproduced either. These discrepancies may be due to overquenching or insufficient resolution of IllustrisTNG simulations.

preprint2022arXiv

Multiple stellar populations in Schwarzschild modeling and the application to the Fornax dwarf

Dwarf spheroidal (dSph) galaxies are believed to be strongly dark matter dominated and thus are considered perfect objects to study dark matter distribution and test theories of structure formation. They possess resolved, multiple stellar populations that offer new possibilities for modeling. A promising tool for the dynamical modeling of these objects is the Schwarzschild orbit superposition method. In this work we extend our previous implementation of the scheme to include more than one population of stars and a more general form of the mass-to-light ratio function. We tested the improved approach on a nearly spherical, gas-free galaxy formed in the cosmological context from the Illustris simulation. We modeled the binned velocity moments for stars split into two populations by metallicity and demonstrate that in spite of larger sampling errors the increased number of constraints leads to significantly tighter confidence regions on the recovered density and velocity anisotropy profiles. We then applied the method to the Fornax dSph galaxy with stars similarly divided into two populations. In comparison with our earlier work, we find the anisotropy parameter to be slightly increasing, rather than decreasing, with radius and more strongly constrained. We are also able to infer anisotropy for each stellar population separately and find them to be significantly different.

preprint2020arXiv

The effect of warm gas on the buckling instability in galactic bars

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.

preprint2020arXiv

Tidal evolution of galaxies in the most massive cluster of IllustrisTNG-100

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.

preprint2020arXiv

Tidally induced warps of spiral galaxies in IllustrisTNG

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&#39;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.

preprint2019arXiv

Dark matter distribution and dynamics of dwarf spheroidal galaxies

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.

preprint2016arXiv

Damping of the Milky Way bar by manifold-driven spirals

We describe a new phenomenon of `bar damping&#39; 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.

preprint2016arXiv

Formation of Andromeda II via a gas-rich major merger and an interaction with M31

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.

preprint2016arXiv

On the tidal radius of satellites on prograde and retrograde orbits

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.

preprint2016arXiv

Rotating stellar populations in the Fornax dSph galaxy

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.

preprint2016arXiv

The orbital structure of a tidally induced bar

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.

preprint2016arXiv

Tidal origin of spiral arms in galaxies orbiting a cluster

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.

preprint2016arXiv

Tidally induced bars of galaxies in clusters

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.

preprint2015arXiv

Simulations of tidally induced spiral arms

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.

preprint2015arXiv

Spectral analysis of stellar orbits in a tidally induced bar

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.

preprint2015arXiv

The origin of prolate rotation in dwarf spheroidal galaxies formed by mergers of disky dwarfs

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.

preprint2015arXiv

The resonant nature of tidal stirring of disky dwarf galaxies orbiting the Milky Way

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&#39;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.

preprint2014arXiv

Andromeda II as a merger remnant

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&#39;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.

preprint2014arXiv

Dark matter in elliptical galaxies: I. Is the total mass density profile of the NFW form or even steeper?

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.

preprint2014arXiv

Dark matter in elliptical galaxies: II. Estimating the mass within the virial radius

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.

preprint2014arXiv

Estimating masses of dwarf spheroidal galaxies

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.

preprint2014arXiv

Radial orbit instability in dwarf dark matter haloes

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.

preprint2014arXiv

Tidal evolution of disky dwarf galaxies: prograde versus retrograde orbits

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&#39;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.

preprint2013arXiv

Can we measure the slopes of density profiles in dwarf spheroidal galaxies?

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.

preprint2013arXiv

Tidal Stirring of Disky Dwarfs with Shallow Dark Matter Density Profiles: Enhanced Transformation into Dwarf Spheroidals

(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.

preprint2013arXiv

Tidal tails of dwarf galaxies on different orbits around the Milky Way

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&#39;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.

preprint2012arXiv

A 2MASS All-Sky View of the Sagittarius Dwarf Galaxy: VII. Kinematics of the Main Body of the Sagittarius dSph

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&#39;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.

preprint2012arXiv

How to make an ultra-faint dwarf spheroidal galaxy: tidal stirring of disky dwarfs with shallow dark matter density profiles

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 ρ\propto r^{-α} for the dwarf DM halos, ranging from core-like (α= 0.2) to cuspy (α= 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 (ρ\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.

preprint2012arXiv

Stellar populations in tidally stirred dwarf galaxies

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.

preprint2012arXiv

The shapes of Milky Way satellites: looking for signatures of tidal stirring

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.

preprint2011arXiv

Dark matter in dwarf galaxies of the Local Group

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.

preprint2011arXiv

Evolutionary tracks of tidally stirred disky dwarf galaxies

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.

preprint2011arXiv

Formation of Dwarf Spheroidal Galaxies Via Mergers of Disky Dwarfs

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.

preprint2010arXiv

Biases in mass estimates of dSph 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.

preprint2010arXiv

Mass profiles and galaxy orbits in nearby galaxy clusters from the analysis of the projected phase-space

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.

preprint2010arXiv

On the Efficiency of the Tidal Stirring Mechanism for the Origin of Dwarf Spheroidals: Dependence on the Orbital and Structural Parameters of the Progenitor Disky Dwarfs

(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.

preprint2010arXiv

The grouping, merging and survival of subhaloes in the simulated Local Group

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.

preprint2010arXiv

The inner structure and kinematics of the Sagittarius dwarf galaxy as a product of tidal stirring

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.

preprint2010arXiv

Tidal stirring of Milky Way satellites: a simple picture with the integrated tidal force

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.

preprint2009arXiv

The mass and anisotropy profiles of galaxy clusters from the projected phase space density: testing the method on simulated data

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.

preprint2009arXiv

The orientation and kinematics of inner tidal tails around dwarf galaxies orbiting the Milky Way

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&#39;: 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.

preprint2009arXiv

The stellar structure and kinematics of dwarf spheroidal galaxies formed by tidal stirring

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.

preprint2009arXiv

Tidal evolution of disky dwarf galaxies in the Milky Way potential: the formation of dwarf spheroidals

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, σ, 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} σ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.

preprint1997arXiv

Weakly Nonlinear Density-Velocity Relation

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(Ω)^{-1} H_0^{-1} \nabla\cdot\bfv$, where $f(Ω) \simeq Ω^{0.6}$, $Ω$ 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 $Ω$ 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 $Ω$-bias degeneracy in comparisons of cosmic density and velocity fields such as IRAS-POTENT.