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.

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.

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.

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

Using N-body simulations we study the 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.

Stellar shells observed in many giant elliptical and lenticular as well as a few spiral and dwarf galaxies presumably result from radial minor mergers of galaxies. We show that the line-of-sight velocity distribution of the shells has a quadruple-peaked shape. We found simple analytical expressions that connect the positions of the four peaks of the line profile with the mass distribution of the galaxy, namely, the circular velocity at the given shell radius and the propagation velocity of the shell. The analytical expressions were applied to a test-particle simulation of a radial minor merger, and the potential of the simulated host galaxy was successfully recovered. Shell kinematics can thus become an independent tool to determine the content and distribution of dark matter in shell galaxies up to ~100 kpc from the center of the host galaxy. Moreover we investigate the dynamical friction and gradual disruption of the cannibalized galaxy during the shell formation in the framework of a simulation with test particles. The coupling of both effects can considerably redistribute positions and luminosities of shells. Neglecting them can lead to significant errors in attempts to date the merger in observed shell galaxies.

Context. The elliptical galaxy NGC 3923 is surrounded by numerous stellar shells that are concentric arcs centered on the galactic core. They are very likely a result of a minor merger and they consist of stars in nearly radial orbits. For a given potential, the shell radii at a given time after the merger can be calculated and compared to observations. The Modified Newtonian Dynamics (MOND) is a theory that aims to solve the missing mass problem by modifying the laws of classical dynamics in the limit of small accelerations. Hernquist & Quinn(1987) claimed that the shell distribution of NGC 3923 contradicted MOND, but Milgrom(1988) found several substantial insufficiencies in their work. Aims. We test whether the observed shell distribution in NGC 3923 is consistent with MOND using the current observational knowledge of the shell number and positions and of the host galaxy surface brightness profile, which supersede the data available in the 1980s when the last (and negative) tests of MOND viability were performed on NGC 3923. Methods. Using the 3.6 um bandpass image of NGC 3923 from the Spitzer space telescope we construct the mass profile of the galaxy. The evolution of shell radii in MOND is then computed using analytical formulae. We use 27 currently observed shells and allow for their multi-generation formation, unlike the Hernquist & Quinn one-generation model that used the 18 shells known at the time. Results. Our model reproduces the observed shell radii with a maximum deviation of 5% for 25 out of 27 known shells while keeping a reasonable formation scenario. A multi-generation nature of the shell system, resulting from successive passages of the surviving core of the tidally disrupted dwarf galaxy, is one of key ingredients of our scenario supported by the extreme shell radial range. The 25 reproduced shells are interpreted as belonging to three generations.

Stellar shells observed in many giant elliptical and lenticular as well as a few spiral and dwarf galaxies, presumably result from galaxy mergers. Line-of-sight velocity distributions of the shells could, in principle, if measured with a sufficiently high S/N, constitute one of methods to constrain the gravitational potential of the host galaxy. Merrifield & Kuijken (1998) predicted a double-peaked line profile for stationary shells resulting from a nearly radial minor merger. In this paper, we aim at extending their analysis to a more realistic case of expanding shells, inherent to the merging process, whereas we assume the same type of merger and the same orbital geometry. We use analytical approach as well as test particle simulations to predict the line-of-sight velocity profile across the shell structure. Simulated line profiles are convolved with spectral PSFs to estimate the peak detectability. The resulting line-of-sight velocity distributions are more complex than previously predicted due to non-zero phase velocity of the shells. In principle, each of the Merrifield & Kuijken (1998) peaks splits into two, giving a quadruple-peaked line profile, which allows more precise determination of the potential of the host galaxy and, moreover, contains additional information. We find simple analytical expressions that connect the positions of the four peaks of the line profile and the mass distribution of the galaxy, namely the circular velocity at the given shell radius and the propagation velocity of the shell. The analytical expressions were applied to a test-particle simulation of a radial minor merger and the potential of the simulated host galaxy was successfully recovered. The shell kinematics can thus become an independent tool to determine the content and distribution of the dark matter in shell galaxies, up to ~100 kpc from the center of the host galaxy.

The model of a radial minor merger proposed by Quinn (1984), which successfully reproduces the observed regular shell systems in shell galaxies, is ideal for a test-particle simulation. We compare such a simulation with a self-consistent one. They agree very well in positions of the first generation of shells but potentially important effects -- dynamical friction and gradual decay of the dwarf galaxy -- are not present in the test-particle model, therefore we look for a proper way to include them.

We present an improved study of the expected shape of the line-of-sight velocity distribution in shell galaxies. We found a simple analytical expression connecting prominent and in principle observable characteristics of the line profile and mass-distribution of the galaxy. The prediction was compared with the results from a test-particle simulation of a radial merger.

Using N-body simulations of shell galaxies created in nearly radial minor mergers, we investigate the error of collision dating, resulting from the neglect of dynamical friction and of gradual disruption of the cannibalized dwarf.

With the goal to refine modelling of shell galaxies and the use of shells to probe the merger history, we develop a new method for implementing dynamical friction in test-particle simulations of radial minor mergers. The friction is combined with a gradual decay of the dwarf galaxy. The coupling of both effects can considerably redistribute positions and luminosities of shells; neglecting them can lead to significant errors in attempts to date the merger.

In the context of exploring mass distributions of dark matter haloes in giant ellipticals, we extend the analysis carried out Merrifield and Kuijken (1998) for stellar line profiles of shells created in nearly radial mergers of galaxies. We show that line-of-sight velocity distributions are more complex than previously predicted. We simulate shell formation and analyze the detectability of spectroscopic signatures of shells after convolution with spectral PSFs.