A new approximate method is developed to allow S(&alpha, β) thermal collision contributions to next-event estimators in Monte Carlo calculations. The new technique is applicable to the more general problem of contributions to next-event estimators from any discrete probability distribution. The method is incorporated into the MCNP Monte Carlo radiation transport code.