Dynamics of barred galaxies: effects of disk height

TL;DR

This study investigates how disk thickness influences bar formation and evolution in disk galaxies within dark matter halos, highlighting the importance of numerical accuracy and code consistency in simulations.

Contribution

It demonstrates that disk thickness significantly affects bar properties and emphasizes the necessity of high-resolution, accurate numerical methods in N-body simulations.

Findings

Thick disks produce slow, long bars inconsistent with observations.

Thin disks generate faster, shorter bars with realistic rotation ratios.

Numerical effects can significantly impact simulation results, but consistent high-resolution methods yield similar outcomes.

Abstract

We study dynamics of bars in models of disk galaxies embeded in realistic dark matter halos. We find that disk thickness plays an important, if not dominant, role in the evolution and structure of the bars. We also make extensive numerical tests of different N-body codes used to study bar dynamics. Models with thick disks typically used in this type of modeling (height-to-length ratio hz/Rd=0.2) produce slowly rotating, and very long, bars. In contrast, more realistic thin disks with the same parameters as in our Galaxy (hz/Rd= 0.1) produce bars with normal length Rbar approx R_d, which rotate quickly with the ratio of the corotation radius to the bar radius 1.2-1.4 compatible with observations. Bars in these models do not show a tendency to slow down, and may lose as little as 2-3 percent of their angular momentum due to dynamical friction with the dark matter over cosmological time. We attribute the differences between the models to a combined effect of high phase-space density and smaller Jeans mass in the thin disk models, which result in the formation of a dense central bulge. Special attention is paid to numerical effects such as the accuracy of orbital integration, force and mass resolution. Using three N-body codes -- Gadget, ART, and Pkdgrav -- we find that numerical effects are very important and, if not carefully treated, may produce incorrect and misleading results. Once the simulations are performed with sufficiently small time-steps and with adequate force and mass resolution, all the codes produce nearly the same results: we do not find any systematic deviations between the results obtained with TREE codes (Gadget and Pkdgrav) and with the Adaptive-Mesh-Refinement (ART) code.