Let $R\left(\lambda \right)$ be the error term in Weyl’s law for a 3-dimensional Riemannian Heisenberg manifold. We prove that ${\int }_1^T{\left|R\left(t\right)\right|}^{2}dt=c{T}^{\frac{5}{2}}+O_{\delta }\left(T^{\frac{9}{4}+\delta }\right)$, where $c$ is a specific nonzero constant and $\delta$ is an arbitrary small positive number. This is consistent with the conjecture of Petridis and Toth stating that $R\left(t\right)=O_{\delta }\left(t^{\frac{3}{4}+\delta }\right)$.The idea of the proof is to use the Poisson summation formula to write the error term in a form which can be estimated by the method of the stationary phase. The similar result will be also proven in the $2n+1$-dimensional case.

Consider the family uₜ = Δu + G(u), t > 0, $x\in {\Omega }_{\epsilon }$, ${\partial }_{{\nu }_{\epsilon }}u=0$, t > 0, $x\in \partial {\Omega }_{\epsilon }$, $\left(E_{\epsilon }\right)$ of semilinear Neumann boundary value problems, where, for ε > 0 small, the set ${\Omega }_{\epsilon }$ is a thin domain in ${ℝ}^l$, possibly with holes, which collapses, as ε → 0⁺, onto a (curved) k-dimensional submanifold of ${ℝ}^l$. If G is dissipative, then equation $\left(E_{\epsilon }\right)$ has a global attractor ${}_{\epsilon }$. We identify a “limit” equation for the family $\left(E_{\epsilon }\right)$, prove convergence of trajectories and establish an upper semicontinuity result for the family ${}_{\epsilon }$ as ε → 0⁺.