We revisit collisionless major and minor mergers of spheroidal galaxies in the context of the size evolution of elliptical galaxies. The simulations are performed as a series of mergers with mass-ratios of 1:1 and 1:10 for models representing pure bulges as well as bulges embedded in dark matter halos. For major and minor mergers, respectively, we identify and analyze two different processes, violent relaxation and stripping, leading to size evolution and a change of the dark matter fraction within the observable effective radius. Violent relaxation - which is the dominant mixing process for major mergers but less important for minor mergers - scatters relatively more dark matter particles than bulge particles to small. Stripping in minor mergers assembles stellar satellite particles at large radii in halo dominated regions of the massive host. This strongly increases the size of the bulge into regions with higher dark matter fractions leaving the inner host structure almost unchanged. A factor of two mass increase by minor mergers increases the dark matter fraction by 20 per cent. We present analytic corrections to simple one-component virial estimates for the evolution of the gravitational radii. If such a two-component system grows by minor mergers alone its size growth, $r_{\mathrm{e}} \propto M^\alpha$, reaches values of $\alpha \approx 2.4$, significantly exceeding the simple theoretical limit of $\alpha = 2$. For major mergers the sizes grow with $\alpha \lesssim 1$. Our results indicate that minor mergers of galaxies embedded in massive dark matter halos provide a potential mechanism for explaining the rapid size growth and the build-up of massive elliptical systems predicting significant dark matter fractions and radially biased velocity dispersions at large radii (abbreviated)

We use oblate axisymmetric dynamical models including dark halos to determine the orbital structure of intermediate mass to massive Coma early-type galaxies. We find a large variety of orbital compositions. Averaged over all sample galaxies the unordered stellar kinetic energy in the azimuthal and the radial direction are of the same order, but they can differ by up to 40 percent in individual systems. In contrast, both for rotating and non-rotating galaxies the vertical kinetic energy is on average smaller than in the other two directions. This implies that even most of the rotating ellipticals are flattened by an anisotropy in the stellar velocity dispersions. Using three-integral axisymmetric toy models we show that flattening by stellar anisotropy maximises the entropy for a given density distribution. Collisionless disk merger remnants are radially anisotropic. The apparent lack of strong radial anisotropy in observed early-type galaxies implies that they may not have formed from mergers of disks unless the influence of dissipational processes was significant.

We use two-dimensional kinematic maps of simulated binary disc mergers to investigate the $\lambda_R$-parameter, which is a luminosity weighted measure of projected angular momentum per unit mass. This parameter was introduced to subdivide the SAURON sample of early type galaxies in so called fast $\lambda_R > 0.1$ and slow rotators $\lambda_R < 0.1$. Tests on merger remnants reveal that $\lambda_R$ is a robust indicator of the true angular momentum content in elliptical galaxies. We find the same range of $\lambda_R$ values in our merger remnants as in the SAURON galaxies. The merger mass ratio is decisive in creating a slow or a fast rotator in a single binary merger, the former being created mostly in an equal mass merger. Slow rotators have a $\lambda_R$ which does not vary with projection. The confusion rate with face-on fast rotators is very small. Merger with low gas fractions form slow rotators with smaller ellipticities and are in much better agreement with the SAURON slow rotators. Remergers of merger remnants are slow rotators but tend to have too high ellipticities. Fast rotators maintain the angular momentum content from the progenitor disc galaxy if merger mass ratio is high. Some SAURON galaxies have values of $\lambda_R$ as high as our progenitor disc galaxies.

We examine the origin of photometrical properties of N-body merger remnants from the perpective of their orbital content. We show that disc mergers alone are unlikely to form a boxy-discy dichotomy in isophotal shape, because of the survival of disc-like orbits even in violent mergers. The shape of the different orbit families can vary strongly, depending on how violent the merger was. Minor axis tubes can become boxy and box orbits can be round. However, for edge-on projections isophotal shape, orbital content and line-of-sight velocity are well connected.

Galactic spheroids can form as a result of galaxy interactions and mergers of disks. Detailed analyses of the photometric properties, the intrinsic orbital structure, the line-of-sight velocity distributions and the kinemetry of simulated merger remnants, which depend critically on the geometry and the gas content of the interacting progenitors, indicate that low and intermediate mass rotating ellipticals can form from mergers of disks. The masses and metallicities of all massive ellipticals and the kinematics of some massive non-rotating ellipticals cannot be explained by binary mergers. Thus these galaxies might have formed in a different way.

We model mock observations of collisionless N-body disc-disc mergers with an axisymmetric orbit superposition program that has been used to model Coma ellipticals. The remnants sample representatively the shapes of disc-disc mergers including prolate, triaxial and oblate objects. The aim is to better understand how the assumption of axial symmetry affects reconstructed masses and stellar motions of systems which are intrinsically not axisymmetric, whether it leads to a bias and how such a potential bias can be recognised in models of real galaxies. The mass recovery at the half-light radius depends on viewing-angle and intrinsic shape: edge-on views allow to reconstruct total masses with an accuracy between 20% (triaxial/prolate remnants) and 3% (oblate remnant). Masses of highly flattened, face-on systems are underestimated by up to 50%. Deviations in local mass densities can be larger where remnants are strongly triaxial or prolate. Luminous M/L are sensitive to box orbits in the remnants. Box orbits cause the central kinematics to vary with viewing-angle. Reconstructed luminous M/L and central masses follow this variation. Luminous M/L are always underestimated (up to a factor of 2.5). Respective dark halos in the models can be overestimated by about the same amount, depending again on viewing angle. Reconstructed velocity anisotropies depend on viewing angle and on the orbital composition of the remnant. We construct N-body realisations of the Schwarzschild models to discuss chaotic orbits and the virial equilibrium in our models. Apparently flattened, rotating ellipticals of intermediate mass are likely close to both, axial symmetry and edge-on orientation. Our results imply that Schwarzschild models allow a reconstruction of their masses and stellar anisotropies with high accuracy. (abridged)

It has been shown that the dissipative gas infall during galaxy formation has the capability to modify the shape of dark matter halos. In this paper we perform the first detailed analysis of particle orbits in a cosmological dark matter halos to understand {\it how} and {\it why} baryons alter its shape. We perform a series of numerical experiments where we grow a baryonic core inside a live dark matter halo extracted from a cosmological simulation. We follow the orbits of more than 300 particles with more that 50000 timesteps. Our results clearly show that the dissipational component is responsible for repeatedly deflecting orbits which visit often the center of the system. Moreover the gravitational potential time dependence associated with the growth of the baryonic mass, shifts the frequencies of the orbits, making them extremely chaotic. This randomization makes the orbits explore a large phase space. When this effect takes place for a significant number of orbits it will be manifested in the density distribution as an approach to a rounder configuration. As a consequence, the influence of the central mass on the shape of the phase space decreases with increasing distance from the center. We discuss the importance of future analysis of controlled experiments (i.e. using analytic potentials instead of live DM halos) to better decipher the dynamics of this phenomenon.

We present a two-dimensional kinematic analysis for a sample of simulated binary disc merger remnants with mass ratios 1:1 and 3:1. For the progenitor discs we used pure stellar models as well as models with 10% of their mass in gas. A multitude of phenomena also observed in real galaxies are found in the simulations. These include misaligned rotation, embedded discs, gas rings, counter-rotating cores and kinematic misaligned discs. Using the 2D maps we illustrate projection effects and the change in properties of a merger remnant when gas is included in the merger. We find that kinematic peculiar subsystems are preferably formed in equal mass mergers. Equal-mass collisionless remnants can show almost no rotation, regular rotation or strong kinematic misalignment. The inclusion of gas makes the remnants appear more round(1:1) and axisymmetric(3:1). Counter-Rotating Cores (CRCs) are almost exclusively formed in equal-mass mergers with a dissipational component. 3:1 remnants show a much more regular structure. We quantify these properties by applying the kinemetric methods recently developed by Krajnovi\'c et al. This work will help to understand observations of elliptical galaxies with 2D field spectrographs, like SAURON.

We present a large set of merger simulations of early-type disc galaxies with mass ratios of 1:1 and 3:1 and 10% of the total disc mass in gas. In contrast to the collisionless case equal-mass mergers with gas do not result in very boxy remnants which is caused by the suppression of box orbits and the change of the projected shape of minor-axis tube orbits in the more axisymmetric remnants. The isophotal shape of 3:1 remnants and the global kinematic properties of 1:1 and 3:1 remnants are only weakly affected by the presence of gas. 1:1 remnants are slowly rotating whereas 3:1 remnants are fast rotating and discy. The shape of the stellar LOSVD is strongly influenced by gas. The LOSVDs of collisionless remnants have broad leading wings while their gaseous counterparts show steep leading wings, more consistent with observations of elliptical galaxies. We show that this change is also caused by the suppressed populating of box orbits and it is amplified by the formation of extended gas discs in the merger remnants. If elliptical galaxies have formed from mergers our results indicate that massive, slowly rotating boxy elliptical galaxies can not have formed from dissipative mergers of discs. Pure stellar (dry) mergers are the more likely candidates. On the other hand lower mass, fast rotating and discy ellipticals can have formed from dissipative (wet) mergers of early-type discs. So far, only unequal-mass disc mergers with gas can successfully explain their observed substructure. This is consistent with the revised morphological classification scheme of increasing importance of gas dissipation when moving from boxy ellipticals to discy ellipticals and then to spiral galaxies, proposed by Kormendy & Bender (abbreviated).

We present a detailed investigation of the relation between the orbital content of merger remnants and observable properties of elliptical and S0 galaxies. Our analysis is based on the statistical sample of collisionless mergers of disk galaxies with different mass ratios and orbital parameters, published by Naab & Burkert. We use the spectral method by Carpintero & Aguilar to determine the orbital content of every remnant and correlate it with its intrinsic shape, and its projected kinematic and photometric properties. We discuss the influence of the bulge component and varying pericenter distances. The two most abundant orbit classes are the minor axis tubes and the box orbits. Their ratio seems to determine the basic properties of a remnant. On average, the fraction of minor axis tubes increases by a factor of two from a merger mass ratio of 1:1 to 4:1, whereas the fraction of box orbits decreases by 10%. At a given mass the central velocity dispersion of a remnant scales with the ratio of minor axis tubes to box orbits. Interestingly, the division line between rotational supported systems and pressure supported objects, $(v_{maj}/\sigma_0)^*=0.7$, turns out to coincide with a box to minor axis tube ratio of unity. The observed $h_3$-$v/\sigma$ anti-correlation for ellipticals can not be reproduced by collisionless merger remnants. We propose that this can only be reconciled by an additional physical process that significantly reduces the box orbit content. Remnants which are dominated by minor axis tube orbits have predominantly disky projections. Boxy remnants have always a box to minor axis tube ratio larger than one. This study will enable to identify observed ellipticals that could have formed, in the collisionless limit, by gas-poor disk mergers.

We use high resolution numerical simulations to investigate the adiabatic contraction of dark matter halos with a Hernquist density profile. We test the response of the halos to the growth of additional axisymmetric disk potentials with various central concentrations and the spherically symmetric potential of a softened point mass. Adding the potentials on timescales that are long compared to the dynamical time scale of the halo, the contracted halos have density profiles that are in excellent agreement with analytical predictions based on the conservation of the adiabatic invariant $M(r)r$. This is surprising as this quantity is strictly conserved only for particles on circular orbits and in spherically symmetric potentials. If the same potentials are added on timescales that are short compared to the dynamical timescale, the result depends strongly on the adopted potential. The adiabatic approximation still works for disk potentials. It does, however, fail for the central potential.