Unraveling the Complexity of Dwarf Galaxy Dynamics: A Study of Binary Orbital Motions

We investigate the impact of binary orbital motions on the dynamical modeling of dwarf galaxies with intrinsic line-of-sight velocity dispersions ( sigma(vr) ) of 1-9 km s(-1). Using dwarf galaxies from the auriga level-2 and level-3 simulations, we apply the Jeans Anisotropic Multi-Gaussian Expansion modeling to tracer stars before and after including binaries to recover the dynamical masses. The recovered total masses within the half-mass radius of tracers, M(< r(half)), are always inflated due to binary motions, with greater inflations occurring for smaller sigma vr . However, many dwarf galaxies experience central density deflated due to binary motions, with little dependence on sigma vr . This is due to the negative radial gradients in the velocity dispersion profiles, with the fractional inflation in sigma vr due to binaries more significant in outskirts. An extreme binary fraction of 70% can lead to central density deflation of up to 10%-20% at 3 km s(-1) < sigma vr < 8 km s-1, with M( < r(half)) inflated by 4% at 9 km s(-1) and up to 15% at 3 km s(-1). A lower binary fraction of 36% leads to similar deflations, with the inflations decreasing to approximately 10% at 3 km s(-1) and becoming statistically insignificant. The choice of binary orbit distribution models does not result in significant differences, and observational errors tend to slightly weaken the deflations in the recovered central density. Two observations separated by 1 yr to exclude binaries lead to almost zero inflations/deflations for a binary fraction of 36% over 3 km s(-1) < sigma(vr)<9 km s(-1). For sigma vr similar to 1 km s(-1) to 3 km s(-1), a binary fraction of 70% (36%) still results in 60% (30%) to 10% (1%) of inflations in M( < r(half)), even with two-epoch observation.