The first reported two-dimensional diffusion-ordered spectroscopy (DOSY) experiments were recorded at low field (LF) on a benchtop NMR spectrometer using the BPP-STE-LED (bipolar pulse pair-stimulated echo sequence with a longitudinal eddy current delay) pulse sequence which limits phase anomalies and baseline discrepancies. A LF DOSY map was first obtained from a solution of a model pharmaceutical formulation containing a macromolecule and an active pharmaceutical ingredient. It revealed a clear separation between the components of the mixture and gave apparent diffusion coefficients (ADC) values consistent with those measured from the reference high field experiment. LF DOSY was then applied to a real esomeprazole medicine and several gradient sampling schemes (linear, exponential and semi-gaussian (SG)) were compared. With a pulsed field gradient range of 4–70%, the most reliable results were given by the SG ramp. The resulting LF DOSY map obtained after 2.84 h of acquisition confirmed that the diffusion dimension is of prime interest to facilitate the assignment of overcrowded LF spectra although relevant ADC values could not be obtained in part of the spectrum with highly overlapped signals.