Pólya's conjecture up to $ε$-loss and quantitative estimates for the remainder of Weyl's law

Let $Ω\subset\mathbb{R}^n$ be a bounded Lipschitz domain. For any $ε\in (0,1)$ we show that for any Dirichlet eigenvalue $λ_k(Ω)>Λ(ε,Ω)$, it holds \begin{align*} k&\le (1+ε)\frac{|Ω|ω(n)}{(2π)^n}λ_k(Ω)^{n/2}, \end{align*} where $Λ(ε,Ω)$ is given explicitly. This reduces the $ε$-loss version of Pólya's conjecture to a computational problem. This estimate is based on a uniform estimate on the remainder of the Weyl law on Lipschitz domains, which appears to be the first quantitative estimate for the remainder of Weyl's law since Weyl's seminal work in the year 1911. We also provide in all dimensions $n\ge 2$, a class of domains that may even have rather irregular shapes or boundaries but satisfy Pólya's conjecture.