The intermittency of pressure and pressure gradient in stationary isotropic turbulence at low to moderate Reynolds numbers is studied by direct numerical simulation (DNS) and theoretically. The energy spectra scale in Kolmogorov units as required by the universal-equilibrium hypothesis, but the pressure spectra do not. It is found that the variances of the pressure and pressure gradient are larger than those computed using the Gaussian approximation for the fourth-order moments of velocity, and that the variance of the pressure gradient, normalized by Kolmogorov units, increases roughly as [Rscr ]1/2λ, where [Rscr ]λ is the Taylor microscale Reynolds number. A theoretical explanation of the Reynolds number dependence is presented which assumes that the small-scale pressure field is driven by coherent small-scale vorticity–strain domains. The variance of the pressure gradient given by the model is the product of the variance of ui,juj,i, the source term of the Poisson equation for pressure, and the square of an effective length of the small-scale coherent vorticity–strain structures. This length can be expressed in terms of the Taylor and Kolmogorov microscales, and the ratio between them gives the observed Reynolds number dependence. Formal asymptotic matching of the spectral scaling observed at small scales in the DNS with the classical scaling at large scales suggests that at high Reynolds numbers the pressure spectrum in these forced flows consists of three scaling ranges which are joined by two inertial ranges, the classical k−7/3 range and a k−5/3 range at smaller scale. It is not possible, within the classical Kolmogorov theory, to determine the length scale at which the inertial range transition occurs because information beyond the energy dissipation rate is required.