Ten weeks of daily imaging of the giant elliptical galaxy M87 with the Hubble Space Telescope (HST) has yielded 41 nova light curves of unprecedented quality for extragalactic cataclysmic variables. We have recently used these light curves to demonstrate that the observational scatter in the so-called Maximum-Magnitude Rate of Decline (MMRD) relation for classical novae is so large as to render the nova-MMRD useless as a standard candle. Here we demonstrate that a modified Buscombe - de Vaucouleurs hypothesis, namely that novae with decline times t2 > 10 days converge to nearly the same absolute magnitude about two weeks after maximum light in a giant elliptical galaxy, is supported by our M87 nova data. For 13 novae with daily-sampled light curves, well determined times of maximum light in both the F606W and F814W filters, and decline times $t2 > 10 days we find that M87 novae display M(606W,15) = -6.37 +/- 0.46 and M(814W,15) = -6.11 +/- 0.43. If very fast novae with decline times t2 < 10 days are excluded, the distances to novae in elliptical galaxies with stellar binary populations similar to those of M87 should be determinable with 1 sigma accuracies of +/-20% with the above calibrations.

In this paper, the third in the series, we continue our study of combinatorics in chaotic Newtonian dynamics. We study the chaotic four-body problem in Newtonian gravity assuming finite-sized particles, and we focus on interactions that produce direct collisions between any two stars. Our long-term goal is to construct an equation that gives the probability of a given collision event occurring over the course of the interaction, as a function of the total encounter energy and angular momentum as well as the numbers and properties of the particles. In previous papers, we varied the number of interacting particles and the distribution of particle radii, for all equal mass particles. Here, we focus on the effects of different combinations of particle masses. We develop an analytic formalism for calculating the time-scales for different collision scenarios to occur. Our analytic time-scales reproduce the simulated time-scales when gravitational focusing is included. We present a method for calculating the relative rates for different types of collisions to occur, assuming two different limits for the particle orbits; radial and tangential. These limits yield relative collision probabilities that bracket the probabilities we obtain directly from numerical scattering experiments, and are designed to reveal important information about the (time-averaged) trajectories of the particles as a function of the interaction parameters. Finally, we present a Collision Rate Diagram (CRD), which directly compares the predictions of our analytic rates to the simulations and quantifies the quality of the agreement. The CRD will facilitate refining our analytic collision rates in future work, as we expand in to the remaining parameter space.