Stellar evolution in N-body simulations of disk galaxies. I

The Kalnajs (1972, 1976) Omega models of global mass and velocity distributions are employed in the present two-dimensional N-body simulation, which allows for a spectrum of particle masses, stellar explosions, explosion remnant interactions with an interstellar medium, and the creation of new stars from the gas. Two sequences of runs using the Omega values of 0.8 and 0.9 examine the separate and combined effects of particle mass distribution, the gravitational influence of an interstellar gas distribution on the N-body particles, and stellar evolution, allowing for stellar explosions and star formation from the gas. It is found that both Omega values' nonequilibrium results dramatically change when evolution is allowed to occur. These results call for more realistic coupled N-body and evolution simulations in order to improve the understanding of disk galaxy evolution.