Molecular dynamics simulations were used to study relaxation of a vibrationally excited C6F6* molecule in a N 2 bath. Ab initio calculations were performed to develop N 2-N2 and N2-C6F6 intermolecular potentials for the simulations. Energy transfer from "hot" C6F6 is studied versus the bath density (pressure) and number of bath molecules. For the large bath limit, there is no heating of the bath. As C6F6* is relaxed, the average energy of C6F6* is determined versus time, i.e., 〈E(t)〉, and for each bath density 〈E(t)〉 is energy dependent and cannot be fit by a single exponential. In the long-time limit C6F6 is fully equilibrated with the bath. For a large bath and low pressures, the simulations are in the fixed temperature, independent collision regime and the simulation results may be compared with gas phase experiments of collisional energy transfer. The derivative d[〈E(t)〉]/ dt divided by the collision frequency ω of the N2 bath gives the average energy transferred from C6F6* per collision 〈ΔEc〉, which is in excellent agreement with experiment. For the ∼100-300 ps simulations reported here, energy transfer from C6F6* is to N2 rotation and translation in accord with the equipartition model, with no energy transfer to N2 vibration. The energy transfer dynamics from C6F 6* is not statistically sensitive to fine details of the N 2-C6F6 intermolecular potential. Tests, with simulation ensembles of different sizes, show that a relatively modest ensemble of only 24 trajectories gives statistically meaningful results.