Comparable upper and lower bounds for boundary values of Neumann eigenfunctions and tight inclusion of eigenvalues

For smooth bounded domains in ℝⁿ, we prove upper and lower L² bounds on the boundary data of Neumann eigenfunctions, and we prove quasiorthogonality of this boundary data in a spectral window. The bounds are tight in the sense that both are independent of the eigenvalues; this is achieved by working with an appropriate norm for boundary functions, which includes a spectral weight, that is, a function of the boundary Laplacian. This spectral weight is chosen to cancel concentration at the boundary that can happen for whispering gallery-type eigenfunctions. These bounds are closely related to wave equation estimates due to Tataru. Using this, we bound the distance from an arbitrary Helmholtz parameterE > 0 to the nearest Neumann eigenvalue in terms of boundary normal derivative data of a trial function u solving the Helmholtz equation (Δ - E)u = 0. This inclusion bound improves over previously known bounds by a factor of E⁵⁼⁶, analogously to a recently improved inclusion bound in the Dirichlet case due to the first two authors. Finally, we apply our theory to present an improved numerical implementation of the method of particular solutions for computation of Neumann eigenpairs on smooth planar domains. We show that the new inclusion bound improves the relative accuracy in a computed Neumann eigenvalue (around the 42000th) from nine to fourteen digits, with negligible extra numerical effort.