S. A. Rodionov

Bar formation and evolution in disc galaxies with gas and a triaxial halo: morphology, bar strength and halo properties

We follow the formation and evolution of bars in N-body simulations of disc galaxies with gas and/or a triaxial halo. We find that both the relative gas fraction and the halo shape play a major role in the formation and evolution of the bar. In gas-rich simulations, the disc stays near-axisymmetric much longer than in gas-poor ones, and, when the bar starts growing, it does so at a much slower rate. Due to these two effects combined, large-scale bars form much later in gas-rich than in gas-poor discs. This can explain the observation that bars are in place earlier in massive red disc galaxies than in blue spirals. We also find that the morphological characteristics in the bar region are strongly influenced by the gas fraction. In particular, the bar at the end of the simulation is much weaker in gas-rich cases. In no case did we witness bar destruction. Halo triaxiality has a dual influence on bar strength. In the very early stages of the simulation it induces bar formation to start earlier. On the other hand, during the later, secular evolution phase, triaxial haloes lead to considerably less increase of the bar strength than spherical ones. The shape of the halo evolves considerably with time. The inner halo parts may become more elongated, or more spherical, depending on the bar strength. The main body of initially triaxial haloes evolves towards sphericity, but in initially strongly triaxial cases it stops well short of becoming spherical. Part of the angular momentum absorbed by the halo generates considerable rotation of the halo particles that stay located relatively near the disc for long periods of time. Another part generates halo bulk rotation, which, contrary to that of the bar, increases with time but stays small.

An iterative method for constructing equilibrium phase models of stellar systems

We present a new method for constructing equilibrium phase models for stellar systems, which we call the iterative method. It relies on constrained, or guided evolution, so that the equilibrium solution has a number of desired parameters and/or constraints. This method is very powerful, to a large extent due to its simplicity. It can be used for mass distributions with an arbitrary geometry and a large variety of kinematical constraints. We present several examples illustrating it. Applications of this method include the creation of initial conditions for N-body simulations and the modelling of galaxies from their photometric and kinematic observations.

Estimating the dark halo mass from the relative thickness of stellar disks

We analyze the relationship between the mass of a spherical component and the minimum possible thickness of stable stellar disks. This relationship for real galaxies allows the lower limit on the dark halo mass to be estimated (the thinner the stable stellar disk is, the more massive the dark halo must be). In our analysis, we use both theoretical relations and numerical N-body simulations of the dynamical evolution of thin disks in the presence of spherical components with different density profiles and different masses. We conclude that the theoretical relationship between the thickness of disk galaxies and the mass of their spherical components is a lower envelope for the model data points. We recommend using this theoretical relationship to estimate the lower limit for the dark halo mass in galaxies. The estimate obtained turns out to be weak. Even for the thinnest galaxies, the dark halo mass within four exponential disk scale lengths must be more than one stellar disk mass.

Optimal choice of the softening length and time step in N-body simulations

A criterion for the choice of optimal softening length

Bending instability of stellar disks: The stabilizing effect of a compact bulge

\epsilon

Mechanisms of the vertical secular heating of a stellar disk

and time-step

An iterative method for the construction of equilibrium N-body models for stellar disks

dt

Dynamical models of the elliptical galaxy NGC 4494

for

Anisotropic models of dark halos

N

Bending instability in galactic discs: advocacy of the linear theory

-body simulations of a collisionless stellar system is analyzed. Plummer and Hernquist spheres are used as models to follow how changes in various parameters of an initially equilibrium stable model depend on