We present a detailed study of photoluminescence transients for neutral donor bound excitons ͑DBEs͒ in GaN, notably the O N donor DBE at 3.4714 eV and the Si Ga DBE at 3.4723 eV. The studied samples are thick strain free nominally undoped bulk GaN samples, with a spectroscopic linewidth Ͻ0.5 meV at 2 K. The photoluminescence ͑PL͒ decay curves for these no-phonon ͑NP͒ lines are strongly nonexponential, and do not allow a proper assessment of the characteristic BE decay time. The decay of the LO-phonon replicas as well as the so-called two-electron transitions ͑TETs͒ at lower energies show a nicely exponential behavior, and allow extraction of DBE decay times of about 1.1 ns for the Si DBE and 1.8 ns for the O DBE, respectively. The initial nonexponential decay behavior of the NP lines has been studied in both the common front surface excitation-detection mode and with detection in transmission through the sample. This initial decay is explained as related to scattering processes in the near surface region, involving the DBEs and free excitons ͑FEs͒. Light scattering processes may also contribute to this complex decay shape. The DBE-LO-phonon decay does not discriminate between the O and Si DBEs because of spectral overlap involving different LO modes. The TET decays at 2 K are very different for transitions related to the DBE ground state and DBE excited states ͑going to p-like donor final states͒, for T Ͼ 10 K thermalization between the DBE ground state and DBE excited states produces a common decay time. Thermalization between free and bound excitons appears to occur above about 20 K, when the DBE decay follows the FE decay. A simple two-level modeling of exciton capture and recombination for the PL decay curves of the FE and the DBEs, as commonly used in the literature, is shown to be generally inadequate. A broad PL background in the TET spectral region is suggested to be related to a radiative Auger process, where the DBEs recombine while leaving the donors ionized.