Extensions and applications of the iterative method

TL;DR

This paper develops an iterative algorithm to construct equilibrium initial conditions for realistic disk galaxy simulations, including gaseous components and triaxial halos, enabling more accurate N-body galaxy modeling.

Contribution

It introduces a new iterative method for creating equilibrium models of disk galaxies with complex components, advancing simulation accuracy.

Findings

Successfully constructed equilibrium models with gaseous disks and triaxial halos.

Models remain close to equilibrium after construction.

Method applicable to various galaxy configurations.

Abstract

Aims: We aim to develop an algorithm for constructing equilibrium initial conditions for simulations of disk galaxies with a triaxial halo and/or a gaseous component. This will pave the way for N-body simulations of realistic disk galaxies. Methods: We use the iterative method, which we presented in a previous article. The idea of this method is very simple. It relies on constrained evolution. Results: We develop an algorithm for constructing equilibrium models of disk galaxies including a gaseous disk and a triaxial or axisymmetric halo. We discuss two test models. The first model consists of a spherical halo, a stellar disk, and an isothermal gaseous disk. The second model consists of a triaxial halo, a stellar disk, and a star-forming gaseous disk. We demonstrate that both test models are very close to equilibrium, as we had intended.