Abstract Magnetic particle imaging (MPI) is a tracer-based tomographic imaging technique utilized in applications such as lung perfusion imaging, cancer diagnosis, stem cell tracking, etc. The goal of translating MPI to clinical use has prompted studies on further improving the spatial-temporal resolutions of MPI through various methods, including image reconstruction algorithm, scanning trajectory design, magnetic field profile design, and tracer design. Iron oxide magnetic nanoparticles (MNPs) are favored for MPI and magnetic resonance imaging (MRI) over other materials due to their high biocompatibility, low cost, and ease of preparation and surface modification. For core-shell MNPs, the tracers’ magnetic core size and non-magnetic coating layer characteristics can significantly affect MPI signals through dynamic magnetization relaxations. Most works to date have assumed an ensemble of MNP tracers with identical sizes, ignoring that artificially synthesized MNPs typically follow a log-normal size distribution, which can deviate theoretical results from experimental data. In this work, we first characterize the size distributions of four commercially available iron oxide MNP products and then model the collective magnetic responses of these MNPs for MPI applications. For an ensemble of MNP tracers with size standard deviations of σ, we applied a stochastic Langevin model to study the effect of size distribution on MPI imaging performance. Under an alternating magnetic field (AMF), i.e., the excitation field in MPI, we collected the time domain dynamic magnetizations (M-t curves), magnetization-field hysteresis loops (M-H curves), point-spread functions (PSFs, dM/dH-H curves), and higher harmonics from these MNP tracers. The intrinsic MPI spatial resolution, which is related to the full width at half maximum (FWHM) of the PSF profile, along with the higher harmonics, serve as metrics to provide insights into how the size distribution of MNP tracers affects MPI performance.