A library of hybrid N-body + test particle simulations for dynamical modeling of galaxy interactions. These are intended for use within the Identikit environment.
This repository consists of bulge+disk+NFW halo galaxy models to facilitate exploration of encounters with various mass ratios. The "full mass" and "half mass" models are the same as those provided by Barnes & Hibbard (2009) and the makefiles are essentially identical to those provided with Identikit, but separated into different files.
The "quarter mass" and "eighth mass" models have been created by George Privon. These models were computed by further reducing the scale lengths by a factor sqrt(2) for each halving of the mass.
The initial conditions files given below provide sample encounters for parabolic orbits and a variety of pericenter seperations.
Initial Conditions and Run Files
The initial conditions files and run files follow a naming convention of
IdkZYXrun, where the numbers replacing ZYX correspond to:
Z= 1 indicating parabolic orbits.
Yindicating the sum of the fractional masses of the interacting galaxies.
2indicates an equal-mass encounter (
5indicates a 1:4 mass ratio (
Xindicates the pericenter separation at first pass, for the idealized Keplerian orbit. This corresponds to multiples of
Idk12X contains initial conditions for equal mass galaxy encounters on parabolic orbits.
The file has a range of first pericenter encounters (1/16 through 16/16, in simulation units).
Idk12Xrun file contains routines to evolve the initial conditions for a suitable time range and save the snapshots for comparison with data.
halo_...Halo mass profiles, accounting for the effects of gravitational softening.
bulge_...Hernquist bulge mass profiles, accounting for the effects of gravitational softening.
Users of these models should cite Barnes & Hibbard (2009).