The physical causes of multiple peaks obversed in thermal desorption spectra are investigated using chemical rate equations and a two binding site model. An exclusion term is introduced which adds the repulsive adatom-adatom hard sphere interactions into the rate equations. Finally, using a fundamental statistical mechanical approach it is demonstrated how general adatom-adatom interactions can be incorporated into the expression for the surface delay time.