Dynamic light scattering (DLS) from ultra-low concentration suspensions (the average number of particles in the scattering volume is less than ~100) gives rise to autocorrelation functions (ACFs) containing a non-Gaussian term due to particle number fluctuations. This term is difficult to characterize and account for and makes recovery of particle size distribution (PSD) information unreliable. We show that an initial analysis of the intensity ACF to determine parameters describing the amplitude and relaxation rate of the non-Gaussian term and then using these parameters to create a better theoretical non-Gaussian ACF model allows a more accurate recovery of the PSD. The modified model is consistent with the measured ACF data, and a reconstructed kernel matrix matching the measured data is obtained. When compared with the usual kernel function reconstruction (KFR) method, the proposed method gives significantly improved PSD recovery accuracy with experimental data. Furthermore, the PSDs obtained have no obvious differences to those obtained from measurements at normal particle concentrations. • Fitting the measured ACF data with a non-Gaussian ACF model of lower concentration. • Reducing focused beam waist value to obtain non-Gaussian model at lower concentration. • Providing a solution for nanoparticle size measurement at ultra-low concentration.