The spectral and temporal dynamics of a multimode dye laser has been studied theoretically and experimentally. The analytical model includes quantum fluctuations as well as four-wave mixing due to population pulsations, stimulated Brillouin scattering, and Rayleigh scattering both in a standing-wave linear laser, and in a unidirectional ring laser. The nonlinearity found most important in the multimode dye laser is four-wave mixing due to pulsations of the population of the upper laser level. Numerical simulations show features that characterize this particular type of mode coupling: broadening of the emission spectrum, oscillations of the light flux in individual laser modes, suppression of certain beat notes. Observations of these features confirm population pulsations dominating the laser dynamics. Four-wave mixing due to population pulsations tends to arrange the phases of the laser modes such as to minimize the pulsations and to limit its own strength.