Rate equations for sympathetic cooling of trapped bosons or fermions

We derive two different sets of rate equations for sympathetic cooling of harmonically trapped bosons or fermions. The rate equations are obtained from a master equation derived earlier by Lewenstein et al. [Phys. Rev. A 51, 4617 (1995)] by means of decoherence and ergodicity arguments. We show analytically that the thermal equilibrium state is a stationary solution of our rate equation. We present analytical results for the rate coefficients which are needed to solve the rate equations, and we give approximate formulas that permit their computation in practice. We solve the two sets of rate equations numerically and compare the results. The cooling times obtained in both approaches agree very well. The equilibration rates show fair agreement