The influence of numerical parameters on tidally triggered bar formation

Abstract

The joint influence of numerical parameters such as the number of particles N, the gravitational softening length $\epsilon$ and the time-step $\Delta t$ is investigated in the context of galaxy simulations. For isolated galaxy models we have performed a convergence study and estimated the numerical parameters ranges for which the relaxed models do not deviate significantly from its initial configuration. By fixing N, we calculate the range of the mean interparticle separation $\lambda (r)$ along the disc radius. We have found that in the simulations with N=1310720 particles $\lambda$ varies by a factor of 6, and the corresponding final Toomre's parameters Q change by only about 5 per cent. By decreasing N, the $\lambda$ and Q ranges broaden. Large $\epsilon$ and small N cause an earlier bar formation. For a given set of parameters the disc heating is smaller with the Plummer softening than with the spline softening. For galaxy collision models numerical simulations indicate that the properties of the formed bars strongly depend upon the selection of N and $\epsilon$. Large values of the gravitational softening parameter and a small number of particles results in the rapid formation of a well defined, slowly rotating bar. On the other hand, small values of $\epsilon$ produce a small, rapidly rotating disc with tightly wound spiral arms, and subsequently a weak bar emerges. We have found that by increasing N, the bar properties converge and the effect of the softening parameter diminishes. Finally, in some cases short spiral arms are observed at the ends of the bar that change periodically from trailing to leading and vice-versa - the wiggle.

17 pages, 13 figures, 3 tables. Accepted for publication in A&A. A high resolution version of the paper is found at http://www.astro.inin.mx/ruslan/tidal_bars/gabbasov.pdf

Cited by in corpus (11)

• Star formation efficiency in galaxy interactions and mergers: a statistical study
• The rapid onset of stellar bars in the baryon-dominated centers of disk galaxies
• Tidally induced bars in gas-rich dwarf galaxies orbiting the Milky Way
• On the bar pattern speed determination of NGC 3367
• Cosmological simulations using a static scalar-tensor theory
• The effect of softening on dynamical simulations of galaxies
• Probing Yukawian gravitational potential by numerical simulations. I. Changing N-body codes
• An interaction scenario of the galaxy pair NGC 3893/96 (KPG 302). A single passage?
• Power spectrum of large-scale structure cosmological models in the framework of scalar-tensor theories
• A thick-disk galaxy model and simulations of equal-mass galaxy pair collisions
• Numerical Study of Bar Suppression in Galaxy Models Due to Disc Heating