Showing 1–50 of 112 results for author: Łokas, E L

A Sagittarius-like simulated dwarf spheroidal galaxy from TNG50

Authors: Ewa L. Lokas

Abstract: The Sagittarius dwarf spheroidal (Sgr dSph) galaxy provides one of the most convincing examples of tidal interaction between satellite galaxies and the Milky Way (MW). The main body of the dwarf was recently demonstrated to have an elongated, prolate, bar-like shape and to possess some internal rotation. Whether these features are temporary results of the strong tidal interaction at the recent per… ▽ More The Sagittarius dwarf spheroidal (Sgr dSph) galaxy provides one of the most convincing examples of tidal interaction between satellite galaxies and the Milky Way (MW). The main body of the dwarf was recently demonstrated to have an elongated, prolate, bar-like shape and to possess some internal rotation. Whether these features are temporary results of the strong tidal interaction at the recent pericenter passage or are due to a disky progenitor is a matter of debate. I present an analog of Sgr selected among bar-like galaxies from the TNG50 simulation of the IllustrisTNG project. The simulated dwarf is initially a disky galaxy with mass exceeding $10^{11}$ M$_\odot$ and evolves around a MW-like host on a tight orbit with seven pericenter passages and a period of about 1 Gyr. At the second pericenter passage, the disk transforms into a bar and the bar-like shape of the stellar component is preserved until the end of the evolution. The morphological transformation is accompanied by strong mass loss, leaving a dwarf with a final mass of below $10^{9}$ M$_\odot$. The gas is lost completely and the star formation ceases at the third pericenter passage. At the last pericenters, the dwarf possesses a bar-like shape, a little remnant rotation, and the metallicity gradient, which are consistent with observations. The more concentrated metal-rich stellar population rotates faster and has a lower velocity dispersion than the more extended metal-poor one. The metallicity distribution evolves so that the most metal-poor stars are stripped first, which explains the metallicity gradient detected in the Sgr stream. This study demonstrates that a dSph galaxy with properties akin to the Sgr dwarf can form from a disky progenitor with a mass of above $10^{11}$ M$_\odot$ by tidal evolution around the MW in the cosmological context. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Comments: 10 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 687, A82 (2024)

A new tidal scenario for double bar formation

Authors: Marcin Semczuk, Ewa L. Łokas, Adriana de Lorenzo-Cáceres, E. Athanassoula

Abstract: Double bars make up a significant fraction of barred galaxies. We propose a new formation scenario for double bars that involves tidal interactions. We demonstrate the viability of this scenario using two examples of simulated galaxies from run TNG50-1 of the IllustrisTNG project. In the proposed scenario the inner bar forms first, either in isolation, via instabilities, or through previous tides.… ▽ More Double bars make up a significant fraction of barred galaxies. We propose a new formation scenario for double bars that involves tidal interactions. We demonstrate the viability of this scenario using two examples of simulated galaxies from run TNG50-1 of the IllustrisTNG project. In the proposed scenario the inner bar forms first, either in isolation, via instabilities, or through previous tides. The outer bar forms later from the material that is tidally distorted by a strong interaction. The inner and outer bars formed this way rotate with different pattern speeds and can be mistaken for a single bar when their phases align. The double-barred structure is stable and can last for at least 3 Gyr. The inner bars of the tidally induced double bars can also have big sizes, which can possibly explain the origin of sizable inner bars recently found in some galaxies. △ Less

Submitted 29 September, 2023; originally announced September 2023.

Comments: 5 pages, 6 figures, submitted to MNRAS Letters

Preprocessing in small groups: Three simulated galaxies interacting prior to cluster infall

Authors: Ewa L. Lokas

Abstract: The formation of galaxy clusters is a complicated process that probably involves the accretion of galaxies in groups, as observed in nearby clusters, such as Virgo and Fornax. The members of the groups undergo "preprocessing" prior to cluster infall, which affects their stellar populations and morphology. In this paper I present an extreme example of such an accretion event selected from the Illus… ▽ More The formation of galaxy clusters is a complicated process that probably involves the accretion of galaxies in groups, as observed in nearby clusters, such as Virgo and Fornax. The members of the groups undergo "preprocessing" prior to cluster infall, which affects their stellar populations and morphology. In this paper I present an extreme example of such an accretion event selected from the IllustrisTNG100 simulation. The group, composed of three full-sized disky galaxies and a number of smaller satellites, is accreted early, with the first pericenter around the cluster at redshift z=1.3. Before the infall, the three galaxies interact strongly in pairs within the group, which produces tidally induced bars in the two more massive ones. The interactions also lead to mass exchange and trigger some star formation activity resulting in temporary rejuvenation of their stellar populations. After infall, they all undergo seven pericenter passages around the cluster, experiencing strong mass loss in the dark matter and gas components, as well as reddening of the stellar populations. Their tidally induced bars are, however, preserved and even enhanced probably due to the loss of gas via ram-pressure stripping in the intracluster medium. The study demonstrates that group accretion can happen very early in cluster formation and proposes another scenario for the formation of tidally induced bars. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: 5 pages, 3 figures, accepted for publication in A&A

Journal ref: A&A 678, A147 (2023)

Internal kinematics of dwarf satellites of MW/M31-like galaxies in TNG50

Authors: Alberto Manuel Martínez-García, Andrés del Pino, Ewa L. Łokas, Roeland P. van der Marel, Antonio Aparicio

Abstract: We present a kinematic study of a thousand of dwarf satellites of MW/M31-like hosts from the IllustrisTNG50 simulation. Internal kinematics were derived for all the snapshots to obtain a historical record of their rotation velocity in the plane of the sky ($|V_T|$) and the amplitude of their velocity gradients along the line of sight ($A_{\rm grad}^{v_z}$) measured from the host. For the majority… ▽ More We present a kinematic study of a thousand of dwarf satellites of MW/M31-like hosts from the IllustrisTNG50 simulation. Internal kinematics were derived for all the snapshots to obtain a historical record of their rotation velocity in the plane of the sky ($|V_T|$) and the amplitude of their velocity gradients along the line of sight ($A_{\rm grad}^{v_z}$) measured from the host. For the majority of the satellites we initially detected rotation in the plane of the sky (65%) or velocity gradients (80%), and this was progressively reduced to 45% and 68% at $z = 0$ respectively. We find that the evolution of the rotation in the plane of the sky and the velocity gradients differs according to type of dwarfs, which could be explained in terms of their different masses and orbital histories. We observe that interaction with the host has an impact on the evolution of the internal kinematics of the satellites. The rotation signal of the satellites is progressively reduced during pericentric passages, the first pericentre being especially disruptive for the initial kinematics. We observe temporary increases in $A_{\rm grad}^{v_z}$ during pericentric passage caused by tidal interaction with the host, $A_{\rm grad}^{v_z}$ increasing as the satellites approach their pericentre and dropping as they move away. In summary, we conclude that the presence of detectable rotation in dwarf satellites is not uncommon, and that the evolution of their internal kinematics is clearly affected by their interaction with the host. △ Less

Submitted 30 October, 2023; v1 submitted 25 July, 2023; originally announced July 2023.

Comments: 12 pages, 9 figures, 2 tables. Accepted for publication in MNRAS

Merging galaxy clusters in IllustrisTNG

Authors: Ewa L. Lokas

Abstract: Mergers between galaxy clusters are an important stage in the formation of the large-scale structure of the Universe. Some of the mergers show a spectacular bow shock that formed as a result of recent passage of a smaller cluster through a bigger one, the classic example of this being the so-called bullet cluster. In this paper, I describe ten examples of interacting clusters identified among 200… ▽ More Mergers between galaxy clusters are an important stage in the formation of the large-scale structure of the Universe. Some of the mergers show a spectacular bow shock that formed as a result of recent passage of a smaller cluster through a bigger one, the classic example of this being the so-called bullet cluster. In this paper, I describe ten examples of interacting clusters identified among 200 of the most massive objects, with total masses above $1.4 \times 10^{14}$ M$_{\odot}$, from the IllustrisTNG300 simulation by searching for prominent bow shocks in their temperature maps. Despite different mass ratios of the two merging clusters, the events are remarkably similar in many respects. In all cases, the companion cluster passed close to the main one only once, between 0.9 and 0.3 Gyr ago, with the pericenter distance of 100-530 kpc and a velocity of up to 3400 km s$^{-1}$. The subcluster, typically an order of magnitude smaller in mass than the main cluster before the interaction, loses most of its dark matter and gas in the process. The displacement between the collisionless part of the remnant and the bow shock is such that the remnant typically lags behind the shock or coincides with it, with a single exception of the merger occurring with the largest velocity. Usually about 1% of the gas cells in the merging clusters are shocked, and the median Mach numbers of these gas cells are around two. Due to the relatively small size of the simulation box, no close analog of the bullet cluster was found, but I identified one case that is similar in terms of mass, velocity, and displacement. The presented cases bear more resemblance to less extreme observed interacting clusters such as A520 and Coma. △ Less

Submitted 26 April, 2023; originally announced April 2023.

Comments: 8 pages, 8 figures, accepted for publication in A&A

Journal ref: A&A 673, A131 (2023)

Formation of a barred galaxy in a major merger: The role of AGN feedback

Authors: Ewa L. Lokas

Abstract: Among the many processes involved in galaxy evolution, those of bar formation, quenching, and feedback from an active galactic nucleus (AGN) seem to be connected, however, the nature of these relations remains unclear. In this work, we aim to elucidate them by studying the formation of a barred galaxy in a major merger of two disks in the IllustrisTNG simulations. This merger involves a coalescenc… ▽ More Among the many processes involved in galaxy evolution, those of bar formation, quenching, and feedback from an active galactic nucleus (AGN) seem to be connected, however, the nature of these relations remains unclear. In this work, we aim to elucidate them by studying the formation of a barred galaxy in a major merger of two disks in the IllustrisTNG simulations. This merger involves a coalescence of two supermassive black holes and a sudden switch to the kinetic mode of AGN feedback implemented in the simulations, which leads to the removal of the gas from the inner part of the galaxy, followed by quenching of star formation and the formation of the bar. This causal relation between AGN feedback and bar formation explains a number of correlations observed in the data, such as the higher frequency of bars among red spirals and the presence of central gas holes in barred galaxies. In such a picture, the bars do not feed the black holes, so their presence does not increase the AGN strength, and they do not cause quenching. However, bars do form in regions characterized by a low gas fraction resulting from AGN feedback. This scenario is probably applicable to many barred galaxies, not only those formed in major mergers. △ Less

Submitted 24 November, 2022; originally announced November 2022.

Comments: 5 pages, 4 figures, accepted for publication in A&A Letters

Journal ref: A&A 668, L3 (2022)

The origin and properties of red spirals: Insights from cosmological simulations

Authors: Ewa L. Lokas

Abstract: A significant fraction of spiral galaxies are red, gas-poor, and have low star formation rates (SFRs). We study these unusual objects using the IllustrisTNG100 simulation. Among 1912 well-resolved disk galaxies selected from the last simulation output, we identify 377 red objects and describe their properties and origins using a few representative examples. The simulated red spirals turn out to be… ▽ More A significant fraction of spiral galaxies are red, gas-poor, and have low star formation rates (SFRs). We study these unusual objects using the IllustrisTNG100 simulation. Among 1912 well-resolved disk galaxies selected from the last simulation output, we identify 377 red objects and describe their properties and origins using a few representative examples. The simulated red spirals turn out to be typically very gas-poor, have very low SFRs, are more metal-rich, and have larger stellar masses than the remaining disks. Only about 13% of red spirals suffered strong mass loss and thus could have resulted from environmental quenching by ram pressure and tidal stripping of the gas, or similar processes. The majority of red disks were probably quenched by feedback from the active galactic nucleus (AGN). This conclusion is supported by the higher black hole masses and lower accretion rates of red disks, as well as the larger total AGN feedback energies injected into the surrounding gas in the kinetic feedback mode implemented in the IllustrisTNG simulations. The timescales of the gas loss correlate with the black hole growth for the AGN-quenched galaxies and with the dark-matter loss for the environmentally quenched ones. The red spirals are more likely to possess bars, and their bars are stronger than in the remaining disks, which is probably the effect of gas loss rather than the reason for quenching. △ Less

Submitted 23 September, 2022; originally announced September 2022.

Comments: 11 pages, 10 figures, accepted for publication in A&A

Journal ref: A&A 667, A27 (2022)

Lopsided galactic disks in IllustrisTNG

Authors: Ewa L. Lokas

Abstract: 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… ▽ More 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. △ Less

Submitted 4 April, 2022; originally announced April 2022.

Comments: 9 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 662, A53 (2022)

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

Authors: Klaudia Kowalczyk, Ewa L. Lokas

Abstract: 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 meth… ▽ More 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. △ Less

Submitted 1 January, 2022; originally announced January 2022.

Comments: 12 pages, 15 figures, accepted for publication in A&A

Journal ref: A&A 659, A119 (2022)

Lopsided galactic bars

Authors: Ewa L. Lokas

Abstract: Most of the observed and simulated galactic bars are symmetric in the face-on view. However, there are indeed cases of bars that are off-center with respect to the disk or have an asymmetric shape. The only well-known example showing both these features is the Large Magellanic Cloud. We report on the identification of several lopsided galactic bars in the Illustris TNG100 simulation found among a… ▽ More Most of the observed and simulated galactic bars are symmetric in the face-on view. However, there are indeed cases of bars that are off-center with respect to the disk or have an asymmetric shape. The only well-known example showing both these features is the Large Magellanic Cloud. We report on the identification of several lopsided galactic bars in the Illustris TNG100 simulation found among a sample of elongated bar-like galaxies studied in the past. The bars show a clear asymmetry in the face-on view, which is in the shape of a footprint. We measured the evolution of the different parameters of the bars' shape and asymmetry as a function of time and find that the asymmetry is preserved for a few Gyr. It can grow together with a bar or appear later, after bar formation. We considered two scenarios leading to the formation of lopsided bars using controlled simulations. In the first, a Milky Way-like galaxy interacts with a massive companion placed on a radial orbit in the plane of the disk and perpendicular to the orientation of the bar at the time of the first passage. In the second, the galaxy initially has an off-center disk and the growth of the bar and its asymmetry is more similar to the one found in IllustrisTNG galaxies, where it is also preceded by the presence of an asymmetric disk. It is thus possible that lopsided bars are born in lopsided disks, although in some cases, the time difference between the occurrence of the asymmetry in the two components is quite large. △ Less

Submitted 9 September, 2021; originally announced September 2021.

Comments: 12 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 655, A97 (2021)

Bar-like galaxies in IllustrisTNG

Authors: Ewa L. Lokas

Abstract: We study a sample of bar-like galaxies in the Illustris TNG100 simulation, in which almost the whole stellar component is in the form of a prolate spheroid. The sample is different from the late-type barred galaxies studied before. In addition to the requirement of a high enough stellar mass and resolution, the 277 galaxies were selected based on the single condition of a low enough ratio of the i… ▽ More We study a sample of bar-like galaxies in the Illustris TNG100 simulation, in which almost the whole stellar component is in the form of a prolate spheroid. The sample is different from the late-type barred galaxies studied before. In addition to the requirement of a high enough stellar mass and resolution, the 277 galaxies were selected based on the single condition of a low enough ratio of the intermediate to long axis of the stellar component. We followed the mass and shape evolution of the galaxies as well as their interactions with other objects and divided them into three classes based on the origin of the bar and the subsequent history. In galaxies of class A (comprising 28% of the sample), the bar was induced by an interaction with a larger object, most often a cluster or group central galaxy, and the galaxies were heavily stripped of dark matter and gas. In classes B and C (27% and 45% of the sample, respectively) the bars were induced by a merger or a passing satellite, or they were formed by disk instability. Class B galaxies were then partially stripped of mass, while those of class C evolved without strong interactions, thus retaining their dark matter and gas in the outskirts. We illustrate the properties of the different classes with three representative examples of individual galaxies. In spite of the different evolutionary histories, the bars are remarkably similar in strength, length, and formation times. The gas fraction in the baryonic component within two stellar half-mass radii at the time of bar formation is always below 0.4 and usually very low, which confirms in the cosmological context the validity of this threshold, which has previously been identified in controlled simulations. Observational counterparts of these objects can be found among early-type fast rotators, S0 galaxies, or red spirals with bars. △ Less

Submitted 5 February, 2021; originally announced February 2021.

Comments: 10 pages, 8 figures, 1 table, accepted for publication in A&A

Journal ref: A&A 647, A143 (2021)

Revealing the Structure and Internal Rotation of the Sagittarius Dwarf Spheroidal Galaxy with Gaia and Machine Learning

Authors: Andrés del Pino, Mark A. Fardal, Roeland P. van der Marel, Ewa L. Łokas, Cecilia Mateu, Sangmo Tony Sohn

Abstract: We present a detailed study of the internal structure and kinematics of the core of the Sagittarius dwarf spheroidal galaxy (Sgr). Using machine-learning techniques, we have combined the information provided by 3300 RR Lyrae stars, more than 2000 spectroscopically observed stars, and the Gaia second data release to derive the full phase space, i.e. 3D positions and kinematics, of more than… ▽ More We present a detailed study of the internal structure and kinematics of the core of the Sagittarius dwarf spheroidal galaxy (Sgr). Using machine-learning techniques, we have combined the information provided by 3300 RR Lyrae stars, more than 2000 spectroscopically observed stars, and the Gaia second data release to derive the full phase space, i.e. 3D positions and kinematics, of more than $1.2\times10^5$ member stars in the core of the galaxy. Our results show that Sgr has a bar structure $\sim 2.5$ kpc long, and that tidal tails emerge from its tips to form what it is known as the Sgr stream. The main body of the galaxy, strongly sheared by tidal forces, is a triaxial (almost prolate) ellipsoid with its longest principal axis of inertia inclined $43^\circ\pm6^\circ$ with respect to the plane of the sky and axis ratios of 1:0.67:0.60. Its external regions are expanding mainly along its longest principal axis, yet the galaxy conserves an inner core of about $500\times330\times300$ pc that shows no net expansion and is rotating at $v_{\rm rot} = 4.13 \pm 0.16$ ${\rm{ km \ s^{-1}}}$. The internal angular momentum of Sgr forms an angle $θ= 18^\circ\pm6^\circ$ with respect to its orbital angular momentum, meaning that the galaxy is in an inclined prograde orbit around the Milky Way. We compared our results with predictions from $N$-body models with spherical, pressure-supported progenitors and a model whose progenitor is a flattened rotating disk. Only the rotating model, based on preexisting simulations aimed at reproducing the line-of-sight velocity gradients observed in Sgr, was able to reproduce the observed properties in the core of the galaxy. △ Less

Submitted 15 February, 2021; v1 submitted 4 November, 2020; originally announced November 2020.

Comments: 37 pages, 26 figures, 2 tables. Accepted for publication in the Astrophysical Journal

Galaxies with kinematically distinct cores in Illustris

Authors: Ivana Ebrová, Ewa L. Łokas, Jiří Eliášek

Abstract: The growing amount of integral-field spectroscopic data creates an increased demand for understanding kinematic peculiarities that carry valuable information about the evolution of the host galaxies. For kinematically distinct cores (KDCs), a number of formation mechanisms have been proposed, but it is still unclear which of them commonly occur in the Universe. We aim to address the KDC formation… ▽ More The growing amount of integral-field spectroscopic data creates an increased demand for understanding kinematic peculiarities that carry valuable information about the evolution of the host galaxies. For kinematically distinct cores (KDCs), a number of formation mechanisms have been proposed, but it is still unclear which of them commonly occur in the Universe. We aim to address the KDC formation in the cosmological context. We used the publicly available data of the large-scale hydrodynamic cosmological simulation Illustris. We identify 134 KDCs, study their properties, and follow their evolution back in time. Examples of four galaxies hosting KDCs are presented and described in detail. The masses of the KDC hosts follow the general distribution of the Illustris galaxies, with a possible slight preference towards massive galaxies. KDCs can be long-lived features, with their formation epochs roughly uniformly distributed in look-back times 0-11.4 Gyr, and they can survive even major or multiple subsequent mergers. There is no single channel of KDC formation, but mergers seem to be the formation mechanism for about 60% of KDCs with a significant preference for major mergers and with the percentage being higher among massive hosts. Other KDCs formed during a pericentric passage or flyby of another galaxy, by precession of a previously formed rapidly rotating core, or without an obvious external cause. The mean mass-weighted stellar age inside the KDC radius is either about the same as the look-back time of the KDC formation or older. Although the radii of our KDCs are on average larger than observed, we find that younger stellar ages are typically associated with smaller KDCs. A significant fraction of KDC hosts possess stellar shells formed during mergers that led to KDCs within the last 5 Gyr, or double peaks in their velocity dispersion maps. △ Less

Submitted 21 December, 2020; v1 submitted 5 October, 2020; originally announced October 2020.

Comments: 18 pages, 16 figures, 1 table; accepted for publication in A&A

Journal ref: A&A 647, A103 (2021)

An interesting case of the formation and evolution of a barred galaxy in the cosmological context

Authors: Ewa L. Lokas

Abstract: Elongated, bar-like galaxies without a significant disk component, with little rotation support and no gas, often form as a result of tidal interactions with a galaxy cluster, as was recently demonstrated using the IllustrisTNG-100 simulation. Galaxies that exhibit similar properties are, however, also found to be infalling into the cluster for the first time. We use the same simulation to study i… ▽ More Elongated, bar-like galaxies without a significant disk component, with little rotation support and no gas, often form as a result of tidal interactions with a galaxy cluster, as was recently demonstrated using the IllustrisTNG-100 simulation. Galaxies that exhibit similar properties are, however, also found to be infalling into the cluster for the first time. We use the same simulation to study in detail the history of such a galaxy over cosmic time in order to determine its origin. The bar appears to be triggered at t=6.8 Gyr by the combined effect of the last significant merger with a subhalo and the first passage of another dwarf satellite, both ten times less massive than the galaxy. The satellites deposit all their gas in the galaxy, contributing to its third and last star-formation episode, which perturbs the disk and may also contribute to the formation of the bar. The galaxy then starts to lose its gas and dark matter due to its passage near a group of more massive galaxies. The strongest interaction involves a galaxy 22 times more massive, leaving the barred galaxy with no gas and half of its maximum dark matter mass. During this time, the bar grows steadily, seemingly unaffected by the interactions, although they may have aided its growth by stripping the gas. The studied galaxy, together with two other similar objects briefly discussed in this letter, suggest the existence of a new class of early-type barred galaxies and thereby demonstrate the importance of interactions in galaxy formation and evolution. △ Less

Submitted 28 September, 2020; originally announced September 2020.

Comments: 5 pages, 4 figures, accepted for publication in A&A Letters

Journal ref: A&A 642, L12 (2020)

Tidally induced warps of spiral galaxies in IllustrisTNG

Authors: Marcin Semczuk, Ewa L. Lokas, Elena D'Onghia, E. Athanassoula, Victor P. Debattista, Lars Hernquist

Abstract: 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 m… ▽ More 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. △ Less

Submitted 17 February, 2020; originally announced February 2020.

Comments: 13 pages, 11 figures, submitted to MNRAS

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

Authors: Ewa L. Lokas

Abstract: 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 peric… ▽ More 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. △ Less

Submitted 25 June, 2020; v1 submitted 3 February, 2020; originally announced February 2020.

Comments: 15 pages, 14 figures, final version published in A&A

Journal ref: A&A 638, A133 (2020)

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

Authors: Ewa L. Lokas

Abstract: 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 v… ▽ More 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. △ Less

Submitted 20 January, 2020; originally announced January 2020.

Comments: 11 pages, 12 figures, accepted for publication in A&A

Journal ref: A&A 634, A122 (2020)

Dark matter distribution and dynamics of dwarf spheroidal galaxies

Authors: Ewa L. Lokas

Abstract: 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 a… ▽ More 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. △ Less

Submitted 23 October, 2019; originally announced October 2019.

Comments: 8 pages, 3 figures, invited review for the IAU Symposium 353 "Galactic Dynamics in the Era of Large Surveys"

Disc galaxies formed from major mergers in Illustris

Authors: Nicolas Peschken, Ewa L. Łokas, E. Athanassoula

Abstract: We show how wet major mergers can create disc galaxies in a cosmological context, using the Illustris simulation. We select a sample of 38 disc galaxies having experienced a major merger in their history with no subsequent significant minor merger, and study how they transform into discs after the merger. In agreement with what was previously found in controlled simulations of such mergers, we fin… ▽ More We show how wet major mergers can create disc galaxies in a cosmological context, using the Illustris simulation. We select a sample of 38 disc galaxies having experienced a major merger in their history with no subsequent significant minor merger, and study how they transform into discs after the merger. In agreement with what was previously found in controlled simulations of such mergers, we find that their disc is built gradually from young stars formed after the merger in the disc region, while the old stars born before the merger form an ellipsoidal component. Focusing on one fiducial case from our sample, we show how the gas was initially dispersed in the halo region right after the merger, but is then accreted onto a disc to form stars, and builds the disc component. We then select a sample of major mergers creating elliptical galaxies, to show that those cases correspond mainly to dry mergers, where the lack of star formation prevents the formation of a disc in the remnant galaxy. The amount of gas in the remnant galaxy after the merger is therefore essential to determine the final outcome of a major merger. △ Less

Submitted 3 March, 2020; v1 submitted 3 September, 2019; originally announced September 2019.

Comments: 14 pages, 14 figures, 1 table, accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 493, Issue 1, March 2020, Pages 1375-1387

Anatomy of a buckling galactic bar

Authors: Ewa L. Lokas

Abstract: 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 ove… ▽ More 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 ν= 4 Ω$. 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. △ Less

Submitted 3 August, 2019; v1 submitted 10 June, 2019; originally announced June 2019.

Comments: 12 pages, 14 figures, minor corrections

Journal ref: A&A 629, A52 (2019)

Buckling instability in tidally induced galactic bars

Authors: Ewa L. Lokas

Abstract: 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… ▽ More 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. △ Less

Submitted 20 February, 2019; v1 submitted 19 February, 2019; originally announced February 2019.

Comments: 9 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 624, A37 (2019)

The major merger origin of the Andromeda II kinematics

Authors: Ivana Ebrová, Ewa L. Łokas, Sylvain Fouquet, Andrés del Pino

Abstract: Prolate rotation (i.e. rotation around the long axis) has been reported for two Local Group dwarf galaxies: Andromeda II, a dwarf spheroidal satellite of M31, and Phoenix, a transition type dwarf galaxy. The prolate rotation may be an exceptional indicator of a past major merger between dwarf galaxies. We showed that this type of rotation cannot be obtained in the tidal stirring scenario, in which… ▽ More Prolate rotation (i.e. rotation around the long axis) has been reported for two Local Group dwarf galaxies: Andromeda II, a dwarf spheroidal satellite of M31, and Phoenix, a transition type dwarf galaxy. The prolate rotation may be an exceptional indicator of a past major merger between dwarf galaxies. We showed that this type of rotation cannot be obtained in the tidal stirring scenario, in which the satellite is transformed from disky to spheroidal by tidal forces of the host galaxy. However, we successfully reproduced the observed Andromeda II kinematics in controlled, self-consistent simulations of mergers between equal-mass disky dwarf galaxies on a radial or close-to-radial orbit. In simulations including gas dynamics, star formation and ram pressure stripping, we are able to reproduce more of the observed properties of Andromeda II: the unusual rotation, the bimodal star formation history and the spatial distribution of the two stellar populations, as well as the lack of gas. We support this scenario by demonstrating the merger origin of prolate rotation in the cosmological context for sufficiently resolved galaxies in the Illustris large-scale cosmological hydrodynamical simulation. △ Less

Submitted 12 September, 2018; originally announced September 2018.

Comments: 4 pages, 3 figures, to appear in proceedings of IAU Symposium 344 "Dwarf Galaxies: From the Deep Universe to the Present"

Journal ref: Proc. IAU 14 (2018) 62-65

Schwarzschild dynamical model of the Fornax dwarf spheroidal galaxy

Authors: Klaudia Kowalczyk, Andrés del Pino, Ewa L. Lokas, Monica Valluri

Abstract: 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 $Υ$ varying with radius and found that Fornax is with high probability embedded in a massive and extended d… ▽ More 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 $Υ$ 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$σ$ level, but still require a high value of $Υ\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 $β$ 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$σ$ confidence level. △ Less

Submitted 13 November, 2018; v1 submitted 20 July, 2018; originally announced July 2018.

Comments: 11 pages, 9 figures, 3 tables, accepted for publication in MNRAS

Tidally induced bars in gas-rich dwarf galaxies orbiting the Milky Way

Authors: Grzegorz Gajda, Ewa L. Łokas, E. Athanassoula

Abstract: Bars in galaxies may form not only through instability, but also as a result of an interaction with another galaxy. In particular, they may appear in disky dwarf galaxies interacting with Milky Way-like galaxies. Here we report the results of N-body/SPH simulations of such dwarfs orbiting in the static potential of a larger galaxy. We used several models of the dwarf galaxy, all of the same mass,… ▽ More Bars in galaxies may form not only through instability, but also as a result of an interaction with another galaxy. In particular, they may appear in disky dwarf galaxies interacting with Milky Way-like galaxies. Here we report the results of N-body/SPH simulations of such dwarfs orbiting in the static potential of a larger galaxy. We used several models of the dwarf galaxy, all of the same mass, but covering a large range of gas fractions: 0, 30 and 70%. We also tested the impact of subgrid star formation processes. In all cases bars of similar length formed in the stellar disk of the dwarfs at the first pericenter passage. However, unexpectedly, the gaseous component remained approximately axisymmetric and unaffected by the bar potential. The bar properties did not change significantly between two consecutive pericenters. The impact of the later encounters with the host depends strongly on the exact orientation of the bar at the pericenter. When the bar is spun up by the tidal force torque, it is also shortened. Conversely, if it is slowed down, it gets longer. In the models with a low gas fraction the bars were more pronounced and survived until the end of the simulations, while in the dwarfs with a high gas fraction the bars were destroyed after the second or third pericenter passage. In terms of the ratio of the corotation radius to the bar length, the bars are slow, and remain so independently of the encounters with the host. △ Less

Submitted 23 October, 2018; v1 submitted 2 July, 2018; originally announced July 2018.

Comments: 19 pages, 18 figures, accepted to ApJ

Journal ref: ApJ (2018), 868, 100

Tidally induced bars in Illustris galaxies

Authors: Nicolas Peschken, Ewa L. Łokas

Abstract: We study barred galaxies selected from the Illustris cosmological simulation, focusing on tidally induced bars formed from flyby interactions. To guarantee high enough resolution we focus on high mass disc galaxies ($M_*>8.3 \times 10^{10} M_{\odot}$). We find that the fraction of barred galaxies among those (21% at redshift $z=0$) is lower in Illustris than observed in the local Universe, and the… ▽ More We study barred galaxies selected from the Illustris cosmological simulation, focusing on tidally induced bars formed from flyby interactions. To guarantee high enough resolution we focus on high mass disc galaxies ($M_*>8.3 \times 10^{10} M_{\odot}$). We find that the fraction of barred galaxies among those (21% at redshift $z=0$) is lower in Illustris than observed in the local Universe, and the fraction grows slightly with redshift. The bar fraction also increases with the stellar mass and decreases with the amount of gas in the disc. Only very few bars at redshift $z=0$ are formed in secular evolution ($\sim$ 7%) and most of them are triggered by external perturbers in mergers or flybys. Many of these bars disappear over time, mostly during secular evolution, which leads to a lower fraction of bars at redshift $z=0$. We then focus on the effect of flyby interactions on the disc and look at tidally induced bars created by a flyby, or pre-existing bars influenced by the passage of a perturber. In the latter case, the interaction can enhance or weaken the bar. During the interaction, the change in the bar strength occurs right after the pericentre passage. The resulting tidally induced bars tend to be stronger than the overall bar sample in Illustris. The preferred scenario to create or enhance a bar seems to be with a strong interaction involving a perturber on a prograde orbit. Furthermore, the strength of the created bar grows with the strength of the interaction. △ Less

Submitted 3 December, 2018; v1 submitted 17 April, 2018; originally announced April 2018.

Comments: 16 pages, 14 figures, accepted for publication in MNRAS

Tidally induced morphology of M33 in hydrodynamical simulations of its recent interaction with M31

Authors: Marcin Semczuk, Ewa L. Lokas, Jean-Baptiste Salomon, E. Athanassoula, Elena D'Onghia

Abstract: 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 interactio… ▽ More 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. △ Less

Submitted 18 July, 2018; v1 submitted 12 April, 2018; originally announced April 2018.

Comments: 22 pages, 17 figures, accepted for publication in ApJ

Dark halo structure in the Carina dwarf spheroidal galaxy: joint analysis of multiple stellar components

Authors: Kohei Hayashi, Michele Fabrizio, Ewa L. Łokas, Giuseppe Bono, Matteo Monelli, Massimo Dall'Ora, Peter B. Stetson

Abstract: Photometric and spectroscopic observations of the Carina dSph revealed that this galaxy contains two dominant stellar populations of different age and kinematics. The co-existence of multiple populations provides new constraints on the dark halo structure of the galaxy, because different populations should be in equilibrium in the same dark matter potential well. We develop non-spherical dynamical… ▽ More Photometric and spectroscopic observations of the Carina dSph revealed that this galaxy contains two dominant stellar populations of different age and kinematics. The co-existence of multiple populations provides new constraints on the dark halo structure of the galaxy, because different populations should be in equilibrium in the same dark matter potential well. We develop non-spherical dynamical models including such multiple stellar components and attempt to constrain the properties of the non-spherical dark halo of Carina. We find that Carina probably has a larger and denser dark halo than found in previous works and a less cuspy inner dark matter density profile, even though the uncertainties of dark halo parameters are still large due to small volume of data sample. Using our fitting results, we evaluate astrophysical factors for dark matter annihilation and decay and find that Carina should be one of the most promising detectable targets among classical dSph galaxies. We also calculate stellar velocity anisotropy profiles for both stellar populations and find that they are both radially anisotropic in the inner regions, while in the outer regions the older population becomes more tangentially biased than the intermediate one. This is consistent with the anisotropy predicted from tidal effects on the dynamical structure of a satellite galaxy and thereby can be considered as kinematic evidence for the tidal evolution of Carina. △ Less

Submitted 28 August, 2018; v1 submitted 5 April, 2018; originally announced April 2018.

Comments: 14 pages, 7 figures and 3 tables. Accepted for publication in MNRAS

Formation of tidally induced bars in galactic flybys: prograde versus retrograde encounters

Authors: Ewa L. Lokas

Abstract: 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 morphol… ▽ More 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. △ Less

Submitted 26 March, 2018; originally announced March 2018.

Comments: 12 pages, 10 figures, accepted for publication in ApJ

The effect of non-sphericity on mass and anisotropy measurements in dSph galaxies with Schwarzschild method

Authors: Klaudia Kowalczyk, Ewa L. Lokas, Monica Valluri

Abstract: 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 mer… ▽ More 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. △ Less

Submitted 16 February, 2018; v1 submitted 30 August, 2017; originally announced August 2017.

Comments: 14 pages, 12 figures, accepted for publication in MNRAS

Galaxies with prolate rotation in Illustris

Authors: Ivana Ebrova, Ewa L. Lokas

Abstract: 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 i… ▽ More 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. △ Less

Submitted 8 December, 2017; v1 submitted 10 August, 2017; originally announced August 2017.

Comments: 21 pages, 16 figures, published in ApJ; particle number corrected

Journal ref: ApJ, 850, 144 (2017)

Tidally induced bars in dwarf galaxies on different orbits around a Milky Way-like host

Authors: Grzegorz Gajda, Ewa L. Lokas, E. Athanassoula

Abstract: 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… ▽ More 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. △ Less

Submitted 22 May, 2017; v1 submitted 8 March, 2017; originally announced March 2017.

Comments: 17 pages, 14 figures, accepted to ApJ

Recovering the mass profile and orbit anisotropy of mock dwarf galaxies with Schwarzschild modelling

Authors: Klaudia Kowalczyk, Ewa L. Lokas, Monica Valluri

Abstract: 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 t… ▽ More 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. △ Less

Submitted 19 June, 2017; v1 submitted 20 February, 2017; originally announced February 2017.

Comments: 12 pages, 11 figures, accepted for publication in MNRAS

The structure of Andromeda II dwarf spheroidal galaxy

Authors: Andrés del Pino, Ewa L. Lokas, Sebastian L. Hidalgo, Sylvain Fouquet

Abstract: 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) foll… ▽ More 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$. △ Less

Submitted 15 May, 2017; v1 submitted 25 November, 2016; originally announced November 2016.

Comments: 17 pages, 17 figures, 6 tables. Revised version accepted for publication in MNRAS

Tidal origin of spiral arms in galaxies orbiting a cluster

Authors: Marcin Semczuk, Ewa L. Lokas, Andres del Pino

Abstract: 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 clu… ▽ More 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. △ Less

Submitted 3 November, 2016; v1 submitted 11 August, 2016; originally announced August 2016.

Comments: 15 pages, 15 figures, accepted for publication in ApJ

Damping of the Milky Way bar by manifold-driven spirals

Authors: Ewa L. Lokas

Abstract: 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… ▽ More 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. △ Less

Submitted 4 October, 2016; v1 submitted 28 July, 2016; originally announced July 2016.

Comments: 6 pages, 5 figures, revised version accepted for publication in ApJ Letters

Journal ref: ApJ, 830, L20 (2016)

The Carina Project. X. On the kinematics of old and intermediate-age stellar populations

Authors: M. Fabrizio, G. Bono, M. Nonino, E. L. Lokas, I. Ferraro, G. Iannicola, R. Buonanno, S. Cassisi, G. Coppola, M. Dall'Ora, R. Gilmozzi, M. Marconi, M. Monelli, M. Romaniello, P. B. Stetson, F. Thévenin, A. R. Walker

Abstract: 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… ▽ More 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/σ 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. △ Less

Submitted 11 July, 2016; originally announced July 2016.

Comments: 21 pages, 16 figures, 4 tables. Accepted for publication in ApJ

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

Authors: Sylvain Fouquet, Ewa L. Lokas, Andres del Pino, Ivana Ebrova

Abstract: 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… ▽ More 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. △ Less

Submitted 29 September, 2016; v1 submitted 29 June, 2016; originally announced June 2016.

Comments: 15 pages, 14 figures, 2 tables, accepted for publication in MNRAS

The orbital structure of a tidally induced bar

Authors: Grzegorz Gajda, Ewa L. Lokas, E. Athanassoula

Abstract: 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 stella… ▽ More 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. △ Less

Submitted 17 August, 2016; v1 submitted 1 June, 2016; originally announced June 2016.

Comments: 15 pages, 15 figures, accepted for publication in ApJ

Journal ref: ApJ, 830, 108 (2016)

Rotating stellar populations in the Fornax dSph galaxy

Authors: Andrés del Pino, Antonio Aparicio, Sebastian L. Hidalgo, Ewa L. Lokas

Abstract: 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 t… ▽ More 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. △ Less

Submitted 18 November, 2016; v1 submitted 30 May, 2016; originally announced May 2016.

Comments: 17 pages, 16 figures. Accepted for publication in MNRAS

Tidally induced bars of galaxies in clusters

Authors: Ewa L. Lokas, Ivana Ebrova, Andres del Pino, Agnieszka Sybilska, E. Athanassoula, Marcin Semczuk, Grzegorz Gajda, Sylvain Fouquet

Abstract: 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 is… ▽ More 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. △ Less

Submitted 30 May, 2016; v1 submitted 27 January, 2016; originally announced January 2016.

Comments: 13 pages, 13 figures, accepted for publication in ApJ

Journal ref: ApJ, 826, 227 (2016)

Orbit anisotropy of dark matter haloes with Schwarzschild modelling

Authors: Klaudia Kowalczyk, Ewa L. Lokas, Monica Valluri

Abstract: 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 disti… ▽ More 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. △ Less

Submitted 22 December, 2015; originally announced December 2015.

Comments: 3 pages, 2 figures, to be published in the Proceedings of the XXXVII Meeting of the Polish Astronomical Society

Simulations of tidally induced spiral arms

Authors: Marcin Semczuk, Ewa L. Lokas

Abstract: 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… ▽ More 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. △ Less

Submitted 18 December, 2015; originally announced December 2015.

Comments: 3 pages, 2 figures, to be published in the Proceedings of the XXXVII Meeting of the Polish Astronomical Society

Spectral analysis of stellar orbits in a tidally induced bar

Authors: Grzegorz Gajda, Ewa L. Lokas, E. Athanassoula

Abstract: 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 elonga… ▽ More 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. △ Less

Submitted 13 November, 2015; originally announced November 2015.

Comments: 3 pages, 2 figures, to be published in the Proceedings of the XXXVII Meeting of the Polish Astronomical Society

On the tidal radius of satellites on prograde and retrograde orbits

Authors: Grzegorz Gajda, Ewa L. Lokas

Abstract: 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 th… ▽ More 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. △ Less

Submitted 15 January, 2016; v1 submitted 13 August, 2015; originally announced August 2015.

Comments: 7 pages, 6 figures, significantly expanded, accepted for publication in ApJ

Journal ref: ApJ, 819, 20 (2016)

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

Authors: Ivana Ebrova, Ewa L. Lokas

Abstract: 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 merge… ▽ More 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. △ Less

Submitted 9 October, 2015; v1 submitted 20 May, 2015; originally announced May 2015.

Comments: 14 pages, 17 figures, two new sections added (on comparison with tidally stirred dwarfs and on comparison with And II data), accepted for publication in ApJ

Journal ref: ApJ, 813, 10 (2015)

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

Authors: Ewa L. Lokas, Marcin Semczuk, Grzegorz Gajda, Elena D'Onghia

Abstract: 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 differ… ▽ More 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. △ Less

Submitted 5 August, 2015; v1 submitted 5 May, 2015; originally announced May 2015.

Comments: 10 pages, 10 figures, minor revision, accepted for publication in ApJ

Journal ref: 2015, ApJ, 810, 100

Estimating masses of dwarf spheroidal galaxies

Authors: Klaudia Kowalczyk, Ewa L. Lokas

Abstract: 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 initial… ▽ More 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. △ Less

Submitted 22 August, 2014; originally announced August 2014.

Comments: 6 pages, 3 figures, to be published in the Proceedings of the XXXVI Annual Meeting of the Polish Astronomical Society

Tidal evolution of disky dwarf galaxies: prograde versus retrograde orbits

Authors: Ewa L. Lokas, Marcin Semczuk

Abstract: 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 respe… ▽ More 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. △ Less

Submitted 21 August, 2014; originally announced August 2014.

Comments: 6 pages, 3 figures, to be published in the Proceedings of the XXXVI Annual Meeting of the Polish Astronomical Society

Radial orbit instability in dwarf dark matter haloes

Authors: Grzegorz Gajda, Ewa L. Lokas, Radoslaw Wojtak

Abstract: 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… ▽ More 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. △ Less

Submitted 20 November, 2014; v1 submitted 24 July, 2014; originally announced July 2014.

Comments: 14 pages, 10 figures, accepted by MNRAS

Journal ref: MNRAS, 447, 97 (2015)

Andromeda II as a merger remnant

Authors: Ewa L. Lokas, Ivana Ebrova, Andres del Pino, Marcin Semczuk

Abstract: 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 m… ▽ More 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. △ Less

Submitted 30 July, 2014; v1 submitted 21 May, 2014; originally announced May 2014.

Comments: 5 pages, 6 figures, accepted for publication in MNRAS Letters

Journal ref: 2014, MNRAS, 445, L6