Estimating the masses of the spherical and disk components of galaxies via numerical simulations

A series of numerical N-body simulations is performed in order to dynamically model the properties of four galaxies (NGC 5603, NGC 3198, NGC 891, and NGC 1566) with known rotation curves, radial disk scales L, and velocity dispersions of old disk stars at various galactocentric distances r. Each model includes a three-dimensional collisionless disk and rigid spherical components, whose relative mass μ was treated as a free parameter that differed from simulation to simulation. The observed disk stellar velocity dispersions were assumed to be equal to or (in the general case) greater than the corresponding line-of-sight projections of the simulated values for the adopted μ after the initially unstable disk is heated and arrives at a steady state. A comparison of the simulated and observed rotational velocities and velocity dispersions provides evidence for “light” disks with μ≥2 in the disk (r<4L).