The regularity of integration kernels forces decay rates of singular values of associated integral operators. This is well-known for symmetric operators with kernels defined on (a,b)×(a,b), where (a,b) is an interval. Over time, many authors have studied this case in detail [9,10,4,15,12]. The case of spheres has also been resolved [3]. A few authors have examined the higher dimensional case or the case of manifolds [1,6]. Typically, these authors have provided decay estimates of singular values in lp norms, p≥1 or in case of faster decay due to regularity, lp quasi-norms, 0<p<1. With that approach, it is straightforward to show that their estimates are optimal using periodic kernels obtained from Fourier series. Our new approach for deriving decay estimates of these singular values uses Weyl's asymptotic formula for Neumann eigenvalues [21] that we combine to an appropriately defined inverse Laplacian. We obtain decay estimates in the form nα where for the n -th singular value where α depends on dimension and on the Sobolev regularity of the kernel. Since we are interested in optimal estimates in case of regular kernels, instead of writing an upper bound by a constant times nα, we use the singular values of the kernel obtained by differentiation. While in [1,6]lp estimates were proven to be optimal by simply considering periodic Fourier series, these Fourier series do not provide sharp results for our estimates. Instead, series of Neumann eigenfunctions for the Laplacian that are specific to the domain of interest are used to prove that our decay estimates are optimal. Finally, we cover the case of real analytic kernels where we are also able to derive optimal estimates.