On learning higher-order cumulants in diffusion models