Search

We discuss parallel (additive) and sequential (multiplicative) variants of overlapping Schwarz methods for the Helmholtz equation in $\mathbb{R}^d$, with large real wavenumber and smooth variable wave speed. The radiation condition is approximated by a Cartesian perfectly-matched layer (PML). The domain-decomposition subdomains are overlapping hyperrectangles with Cartesian PMLs at their boundaries. In a recent paper ({\tt arXiv:2404.02156}), the current authors proved (for both variants) that, after a specified number of iterations -- depending on the behaviour of the geometric-optic rays -- the error is smooth and smaller than any negative power of the wavenumber $k$. For the parallel method, the specified number of iterations is less than the maximum number of subdomains, counted with their multiplicity, that a geometric-optic ray can intersect. The theory, which is given at the continuous level and makes essential use of semi-classical analysis, assumes that the overlaps of the subdomains and the widths of the PMLs are all independent of the wavenumber. In this paper we extend the results of {\tt arXiv:2404.02156} by experimentally studying the behaviour of the methods in the practically important case when both the overlap and the PML width decrease as the wavenumber increases. We find that (at least for constant wavespeed), the methods remain robust to increasing $k$, even for miminal overlap, when the PML is one wavelength wide., Comment: 8 pages, 2 figures

We prove sharp wavenumber-explicit error bounds for first- or second-type-N\'ed\'elec-element (a.k.a. edge-element) conforming discretisations, of arbitrary (fixed) order, of the variable-coefficient time-harmonic Maxwell equations posed in a bounded domain with perfect electric conductor (PEC) boundary conditions. The PDE coefficients are allowed to be piecewise regular and complex-valued; this set-up therefore includes scattering from a PEC obstacle and/or variable real-valued coefficients, with the radiation condition approximated by a perfectly matched layer (PML). In the analysis of the $h$-version of the finite-element method, with fixed polynomial degree $p$, applied to the time-harmonic Maxwell equations, the $\textit{asymptotic regime}$ is when the meshwidth, $h$, is small enough (in a wavenumber-dependent way) that the Galerkin solution is quasioptimal independently of the wavenumber, while the $\textit{preasymptotic regime}$ is the complement of the asymptotic regime. The results of this paper are the first preasymptotic error bounds for the time-harmonic Maxwell equations using first-type N\'ed\'elec elements or higher-than-lowest-order second-type N\'ed\'elec elements. Furthermore, they are the first wavenumber-explicit results, even in the asymptotic regime, for Maxwell scattering problems with a non-empty scatterer.

We study overlapping Schwarz methods for the Helmholtz equation posed in any dimension with large, real wavenumber and smooth variable wave speed. The radiation condition is approximated by a Cartesian perfectly-matched layer (PML). The domain-decomposition subdomains are overlapping hyperrectangles with Cartesian PMLs at their boundaries. The overlaps of the subdomains and the widths of the PMLs are all taken to be independent of the wavenumber. For both parallel (i.e., additive) and sequential (i.e., multiplicative) methods, we show that after a specified number of iterations -- depending on the behaviour of the geometric-optic rays -- the error is smooth and smaller than any negative power of the wavenumber. For the parallel method, the specified number of iterations is less than the maximum number of subdomains, counted with their multiplicity, that a geometric-optic ray can intersect. These results, which are illustrated by numerical experiments, are the first wavenumber-explicit results about convergence of overlapping Schwarz methods for the Helmholtz equation, and the first wavenumber-explicit results about convergence of any domain-decomposition method for the Helmholtz equation with a non-trivial scatterer (here a variable wave speed).

We show that for the Lindblad evolution defined using (at most) quadratically growing classical Hamiltonians and (at most) linearly growing classical jump functions (quantized into jump operators assumed to satisfy certain ellipticity conditions and modeling interaction with a larger system), the evolution of a quantum observable remains close to the classical Fokker--Planck evolution in the Hilbert--Schmidt norm for times vastly exceeding the Ehrenfest time (the limit of such agreement with no jump operators). The time scale is the same as in the recent papers by Hern\'andez--Ranard--Riedel but the statement and methods are different. The appendix presents numerical experiments illustrating the classical/quantum correspondence in Lindblad evolution and comparing it to the mathematical results., Comment: Main article by Jeffrey Galkowski and Maciej Zworski with an appendix by Zhen Huang and Maciej Zworski -- new appendix with numerical experiments and a new section providing estimates for the Hilbert--Schmidt norm under Lindblad evolution

Motivated by the study of flat bands in models of twisted bilayer graphene (TBG), we give abstract conditions which guarantee the existence of a discrete set of parameters for which periodic Hamiltonians exhibit flat bands. As an application, we show that a scalar operator derived from the chiral model of TBG has flat bands for a discrete set of parameters.

For $h$-FEM discretisations of the Helmholtz equation with wavenumber $k$, we obtain $k$-explicit analogues of the classic local FEM error bounds of [Nitsche, Schatz 1974], [Wahlbin 1991], [Demlow, Guzm\'an, Schatz 2011], showing that these bounds hold with constants independent of $k$, provided one works in Sobolev norms weighted with $k$ in the natural way. We prove two main results: (i) a bound on the local $H^1$ error by the best approximation error plus the $L^2$ error, both on a slightly larger set, and (ii) the bound in (i) but now with the $L^2$ error replaced by the error in a negative Sobolev norm. The result (i) is valid for shape-regular triangulations, and is the $k$-explicit analogue of the main result of [Demlow, Guzm\'an, Schatz, 2011]. The result (ii) is valid when the mesh is locally quasi-uniform on the scale of the wavelength (i.e., on the scale of $k^{-1}$) and is the $k$-explicit analogue of the results of [Nitsche, Schatz 1974], [Wahlbin 1991]. Since our Sobolev spaces are weighted with $k$ in the natural way, the result (ii) indicates that the Helmholtz FEM solution is locally quasi-optimal modulo low frequencies (i.e., frequencies $\lesssim k$). Numerical experiments confirm this property, and also highlight interesting propagation phenomena in the Helmholtz FEM error.

In this article, we study propagation of defect measures for Schr\"odinger operators, $-h^2\Delta_g+V$, on a Riemannian manifold $(M,g)$ of dimension $n$ with $V$ having conormal singularities along a hypersurface $Y$ in the sense that derivatives along vector fields tangent to $Y$ preserve the regularity of $V$. We show that the standard propagation theorem holds for bicharacteristics travelling transversally to the surface $Y$ whenever the potential is absolutely continuous. Furthermore, even when bicharacteristics are tangent to $Y$ at exactly first order, as long as the potential has an absolutely continuous first derivative, standard propagation continues to hold.

In the analysis of the $h$-version of the finite-element method (FEM), with fixed polynomial degree $p$, applied to the Helmholtz equation with wavenumber $k\gg 1$, the $\textit{asymptotic regime}$ is when $(hk)^p C_{\rm sol}$ is sufficiently small and the sequence of Galerkin solutions are quasioptimal; here $C_{\rm sol}$ is the $L^2 \to L^2$ norm of the Helmholtz solution operator, with $C_{\rm sol} \sim k$ for nontrapping problems. In the $\textit{preasymptotic regime}$, one expects that if $(hk)^{2p}C_{\rm sol}$ is sufficiently small, then (for physical data) the relative error of the Galerkin solution is controllably small. In this paper, we prove the natural error bounds in the preasymptotic regime for the variable-coefficient Helmholtz equation in the exterior of a Dirichlet, or Neumann, or penetrable obstacle (or combinations of these) and with the radiation condition $\textit{either}$ realised exactly using the Dirichlet-to-Neumann map on the boundary of a ball $\textit{or}$ approximated either by a radial perfectly-matched layer (PML) or an impedance boundary condition. Previously, such bounds for $p>1$ were only available for Dirichlet obstacles with the radiation condition approximated by an impedance boundary condition. Our result is obtained via a novel generalisation of the "elliptic-projection" argument (the argument used to obtain the result for $p=1$) which can be applied to a wide variety of abstract Helmholtz-type problems.

Let $(M,g)$ be a Zoll manifold, i.e., a smooth, compact, Riemannian manifold without boundary all of whose geodesics are closed with a minimal common period $T$. The positive definite Laplace-Beltrami operator has eigenvalues $\{\lambda_j^2\}_j$ which cluster around $\nu^2_\ell$ for some sequence $\nu_\ell\to \infty$. This article is concerned with the number of $\lambda_j$ in a window of fixed size $\mathrm{w}$ around $\nu_\ell$, denoted by $\mathbf{N}(\nu_\ell,\mathrm{w}):=\#\{j\,:\, \lambda_j\in[\nu_\ell-\mathrm{w},\nu_\ell+\mathrm{w}]\}.$ When the set of trajectories with period smaller than $T$ has zero measure, there is $c_{n}>0$, depending only on $n=\operatorname{dim} M$, such that $$ \mathbf{N}(\nu_\ell,\mathrm{w}) =c_n\operatorname{vol}_g(M)\nu_{\ell}^{n-1}+o(\nu_{\ell}^{n-1}), $$ as $\ell \to \infty$. However, for a general Zoll manifold this may not be the case. We show that, nevertheless, there is $N>0$, independent of $\ell$, such that $$ \sum_{j=0}^{N-1}\mathbf{N}(\nu_{\ell+j},\mathrm{w})= c_nN\operatorname{vol}_g(M)\nu_{\ell}^{n-1}+o(\nu_{\ell}^{n-1}), $$ as $\ell \to \infty$. In addition to asymptotics for the counting function, we study the kernel of the spectral projector for the Laplacian, $\Pi_{\ell,\mathrm{w}}(x,y)$ onto the spectrum in ${\bigcup_{j=0}^{N-1}[\nu_{\ell+j}-\mathrm{w},\nu_{\ell+j}+\mathrm{w}]}$. We show that for $x$ and $y$ in a shrinking neighborhood of a point with few loops of length smaller than $T$, $\Pi_{\ell,\mathrm{w}}(x,y)$ and its derivatives have the same asymptotics as those on the round sphere and flat torus., Comment: The main results in the current version are novel. Theorem 1 in this version is new. Theorem 1 in the previous version was misstated; specifically missing the assumption that the Zoll manifold be SC_T. Theorem 2 in the current posting replaces this assumption with a weaker one

We prove sharp lower bounds on the error incurred when approximating any oscillating function using piecewise polynomial spaces., Comment: Updated introductory material and included proof of optimality for the low-frequency estimates in full generality

The scattering phase, defined as $ \log \det S ( \lambda ) / 2\pi i $ where $ S ( \lambda ) $ is the (unitary) scattering matrix, is the analogue of the counting function for eigenvalues when dealing with exterior domains and is closely related to Krein's spectral shift function. We revisit classical results on asymptotics of the scattering phase and point out that it is never monotone in the case of strong trapping of waves. Perhaps more importantly, we provide the first numerical calculations of scattering phases for non-radial scatterers. They show that the asymptotic Weyl law is accurate even at low frequencies and reveal effects of trapping such as lack of monotonicity. This is achieved by using the recent high level multiphysics finite element software FreeFEM.

In this article, we consider the asymptotic behaviour of the spectral function of Schr\"odinger operators on the real line. Let $H: L^2(\mathbb{R})\to L^2(\mathbb{R})$ have the form $$ H:=-\frac{d^2}{dx^2}+V, $$ where $V$ is a formally self-adjoint first order differential operator with smooth coefficients, bounded with all derivatives. We show that the kernel of the spectral projector, $\mathbb{1}_{(-\infty,\rho^2]}(H)$, has a complete asymptotic expansion in powers of $\rho$. This settles the 1-dimensional case of a conjecture made by the last two authors., Comment: expanded exposition and simplified various steps proofs

We consider approximation of the variable-coefficient Helmholtz equation in the exterior of a Dirichlet obstacle using perfectly-matched-layer (PML) truncation; it is well known that this approximation is exponentially accurate in the PML width and the scaling angle, and the approximation was recently proved to be exponentially accurate in the wavenumber $k$ in [Galkowski, Lafontaine, Spence, 2021]. We show that the $hp$-FEM applied to this problem does not suffer from the pollution effect, in that there exist $C_1,C_2>0$ such that if $hk/p\leq C_1$ and $p \geq C_2 \log k$ then the Galerkin solutions are quasioptimal (with constant independent of $k$), under the following two conditions (i) the solution operator of the original Helmholtz problem is polynomially bounded in $k$ (which occurs for "most" $k$ by [Lafontaine, Spence, Wunsch, 2021]), and (ii) either there is no obstacle and the coefficients are smooth or the obstacle is analytic and the coefficients are analytic in a neighbourhood of the obstacle and smooth elsewhere. This $hp$-FEM result is obtained via a decomposition of the PML solution into "high-" and "low-frequency" components, analogous to the decomposition for the original Helmholtz solution recently proved in [Galkowski, Lafontaine, Spence, Wunsch, 2022]. The decomposition is obtained using tools from semiclassical analysis (i.e., the PDE techniques specifically designed for studying Helmholtz problems with large $k$).

Let $M$ be a smooth compact manifold of dimension $d$ without boundary. We introduce the concept of predominance for Riemannian metrics on $M$, a notion analogous to full Lebesgue measure which, in particular, implies density. We show that for a predominant metric, the number of closed geodesics of length smaller than $T$ has a stretched exponential upper bound in $T$. In addition, we study remainders in the Weyl law for predominant metrics. The Weyl law states that the number of Laplace-Beltrami eigenvalues smaller than $\lambda^2$ is asymptotic to $C\lambda^d$ with an $O(\lambda^{d-1})$ error. We show that, for a predominant metric, the estimate on the error can by improved by a power of $\log \lambda$. After an application of recent results of the authors in the case of the Weyl law, these estimates follow from a study of the non-degeneracy properties of nearly closed orbits for predominant sets of metrics.

In $d$ dimensions, accurately approximating an arbitrary function oscillating with frequency $\lesssim k$ requires $\sim k^d$ degrees of freedom. A numerical method for solving the Helmholtz equation (with wavenumber $k$ and in $d$ dimensions) suffers from the pollution effect if, as $k\to\infty$, the total number of degrees of freedom needed to maintain accuracy grows faster than this natural threshold (i.e., faster than $k^d$ for domain-based formulations, such as finite element methods, and $k^{d-1}$ for boundary-based formulations, such as boundary element methods). It is well known that the $h$-version of the finite element method (FEM) (where accuracy is increased by decreasing the meshwidth $h$ and keeping the polynomial degree $p$ fixed) suffers from the pollution effect, and research over the last $\sim$ 30 years has resulted in a near-complete rigorous understanding of how quickly the number of degrees of freedom must grow with $k$ to maintain accuracy. In contrast to the $h$-FEM, at least empirically, the $h$-version of the boundary element method (BEM) does $\textit{not}$ suffer from the pollution effect (recall that in the boundary element method the scattering problem is reformulated as an integral equation on the boundary of the scatterer, with this integral equation then solved numerically using a finite-element-type approximation space). However, the current best results in the literature on how quickly the number of degrees of freedom for the $h$-BEM must grow with $k$ to maintain accuracy fall short of proving this. In this paper, we prove that the $h$-version of the Galerkin method applied to the standard second-kind boundary integral equations for solving the Helmholtz exterior Dirichlet problem does not suffer from the pollution effect when the obstacle is nontrapping (i.e., does not trap geometric-optic rays).

We employ separation of variables to prove weighted resolvent estimates for the semiclassical Schr\"odinger operator $-h^2 \Delta + V(|x|) - E$ in dimension $n \ge 2$, where $h, \, E > 0$, and $V: [0, \infty) \to \mathbb{R}$ is $L^\infty$ and compactly supported. The weighted resolvent norm grows no faster than $\exp(Ch^{-1})$, while an exterior weighted norm grows $\sim h^{-1}$. We introduce a new method based on the Mellin transform to handle the two-dimensional case.

Let $(\Omega,g)$ be a compact, analytic Riemannian manifold with analytic boundary $\partial \Omega = M.$ We give $L^2$-lower bounds for Steklov eigenfunctions and their restrictions to interior hypersurfaces $H \subset \Omega^{\circ}$ in a geometrically defined neighborhood of $M$. Our results are optimal in the entire geometric neighborhood and complement the results on eigenfunction upper bounds in the author's previous work., Comment: 16 pages, 2 figures

We study a commonly-used second-kind boundary-integral equation for solving the Helmholtz exterior Neumann problem at high frequency, where, writing $\Gamma$ for the boundary of the obstacle, the relevant integral operators map $L^2(\Gamma)$ to itself. We prove new frequency-explicit bounds on the norms of both the integral operator and its inverse. The bounds on the norm are valid for piecewise-smooth $\Gamma$ and are sharp up to factors of $\log k$ (where $k$ is the wavenumber), and the bounds on the norm of the inverse are valid for smooth $\Gamma$ and are observed to be sharp at least when $\Gamma$ is smooth with strictly-positive curvature. Together, these results give bounds on the condition number of the operator on $L^2(\Gamma)$; this is the first time $L^2(\Gamma)$ condition-number bounds have been proved for this operator for obstacles other than balls.

We consider a wide variety of scattering problems including scattering by Dirichlet, Neumann, and penetrable obstacles. We consider a radial perfectly-matched layer (PML) and show that for any PML width and a steep-enough scaling angle, the PML solution is exponentially close, both in frequency and the tangent of the scaling angle, to the true scattering solution. Moreover, for a fixed scaling angle and large enough PML width, the PML solution is exponentially close to the true scattering solution in both frequency and the PML width. In fact, the exponential bound holds with rate of decay $c(w\tan\theta -C) k$ where $w$ is the PML width and $\theta$ is the scaling angle. More generally, the results of the paper hold in the framework of black-box scattering under the assumption of an exponential bound on the norm of the cutoff resolvent, thus including problems with strong trapping. These are the first results on the exponential accuracy of PML at high-frequency with non-trivial scatterers., Comment: 41 pages, 6 figures

Over the last ten years, results from [Melenk-Sauter, 2010], [Melenk-Sauter, 2011], [Esterhazy-Melenk, 2012], and [Melenk-Parsania-Sauter, 2013] decomposing high-frequency Helmholtz solutions into "low"- and "high"-frequency components have had a large impact in the numerical analysis of the Helmholtz equation. These results have been proved for the constant-coefficient Helmholtz equation in either the exterior of a Dirichlet obstacle or an interior domain with an impedance boundary condition. Using the Helffer-Sj\"ostrand functional calculus, this paper proves analogous decompositions for scattering problems fitting into the black-box scattering framework of Sj\"ostrand-Zworski, thus covering Helmholtz problems with variable coefficients, impenetrable obstacles, and penetrable obstacles all at once. These results allow us to prove new frequency-explicit convergence results for (i) the $hp$-finite-element method applied to the variable coefficient Helmholtz equation in the exterior of a Dirichlet obstacle, when the obstacle and coefficients are analytic, and (ii) the $h$-finite-element method applied to the Helmholtz penetrable-obstacle transmission problem. In particular, the result in (i) shows that the $hp$-FEM applied to this problem does not suffer from the pollution effect.

We consider GMRES applied to discretisations of the high-frequency Helmholtz equation with strong trapping; recall that in this situation the problem is exponentially ill-conditioned through an increasing sequence of frequencies. Under certain assumptions about the distribution of the eigenvalues, we prove upper bounds on how the number of GMRES iterations grows with the frequency. Our main focus is on boundary-integral-equation formulations of the exterior Dirichlet and Neumann obstacle problems in 2- and 3-d; for these problems, we investigate numerically the sharpness (in terms of dependence on frequency) of both our bounds and various quantities entering our bounds. This paper is therefore the first comprehensive study of the frequency-dependence of the number of GMRES iterations for Helmholtz boundary-integral equations under trapping.

We consider approximating the solution of the Helmholtz exterior Dirichlet problem for a nontrapping obstacle, with boundary data coming from plane-wave incidence, by the solution of the corresponding boundary value problem where the exterior domain is truncated and a local absorbing boundary condition coming from a Pad\'e approximation (of arbitrary order) of the Dirichlet-to-Neumann map is imposed on the artificial boundary (recall that the simplest such boundary condition is the impedance boundary condition). We prove upper- and lower-bounds on the relative error incurred by this approximation, both in the whole domain and in a fixed neighbourhood of the obstacle (i.e. away from the artificial boundary). Our bounds are valid for arbitrarily-high frequency, with the artificial boundary fixed, and show that the relative error is bounded away from zero, independent of the frequency, and regardless of the geometry of the artificial boundary., Comment: 86 pages

For the Helmholtz equation posed in the exterior of a Dirichlet obstacle, we prove that if there exists a family of quasimodes (as is the case when the exterior of the obstacle has stable trapped rays), then there exist near-zero eigenvalues of the standard variational formulation of the exterior Dirichlet problem (recall that this formulation involves truncating the exterior domain and applying the exterior Dirichlet-to-Neumann map on the truncation boundary). Our motivation for proving this result is that a) the finite-element method for computing approximations to solutions of the Helmholtz equation is based on the standard variational formulation, and b) the location of eigenvalues, and especially near-zero ones, plays a key role in understanding how iterative solvers such as the generalised minimum residual method (GMRES) behave when used to solve linear systems, in particular those arising from the finite-element method. The result proved in this paper is thus the first step towards rigorously understanding how GMRES behaves when applied to discretisations of high-frequency Helmholtz problems under strong trapping (the subject of the companion paper [Marchand, Galkowski, Spence, Spence, 2021]).

In this article we consider asymptotics for the spectral function of Schr\"odinger operators on the real line. Let $P:L^2(\mathbb{R})\to L^2(\mathbb{R})$ have the form $$ P:=-\tfrac{d^2}{dx^2}+W, $$ where $W$ is a self-adjoint first order differential operator with certain modified almost periodic structure. We show that the kernel of the spectral projector, $\mathbb{1}_{(-\infty,\lambda^2]}(P)$ has a full asymptotic expansion in powers of $\lambda$. In particular, our class of potentials $W$ is stable under perturbation by formally self-adjoint first order differential operators with smooth, compactly supported coefficients. Moreover, it includes certain potentials with dense pure point spectrum. The proof combines the gauge transform methods of Parnovski-Shterenberg and Sobolev with Melrose's scattering calculus., Comment: Small typographical corrections in section 5

We obtain new quantitative estimates on Weyl Law remainders under dynamical assumptions on the geodesic flow. On a smooth compact Riemannian manifold $(M,g)$ of dimension $n$, let $\Pi_\lambda$ denote the kernel of the spectral projector for the Laplacian, $\mathbb{1}_{[0,\lambda^2]}(-\Delta_g)$. Assuming only that the set of near periodic geodesics over $W\subset M$ has small measure, we prove that as $\lambda \to \infty$ $$ \int_{W} \Pi_\lambda(x,x)dx=(2\pi)^{-n}\text{vol}_{\mathbb{R}^n}(B)\text{vol}_g(W)\,\lambda^n+O\Big(\frac{\lambda^{n-1}}{\log \lambda }\Big),$$ where $B$ is the unit ball. One consequence of this result is that the improved remainder holds on all product manifolds, in particular giving improved estimates for the eigenvalue counting function in the product setup. Our results also include logarithmic gains on asymptotics for the off-diagonal spectral projector $\Pi_\lambda(x,y)$ under the assumption that the set of geodesics that pass near both $x$ and $y$ has small measure, and quantitative improvements for Kuznecov sums under non-looping type assumptions. The key technique used in our study of the spectral projector is that of geodesic beams., Comment: Updated references

We give an elementary proof of weighted resolvent estimates for the semiclassical Schr\"odinger operator $-h^2 \Delta + V(x) - E$ in dimension $n \neq 2$, where $h, \, E > 0$. The potential is real-valued, $V$ and $\partial_r V$ exhibit long range decay at infinity, and may grow like a sufficiently small negative power of $r$ as $r \to 0$. The resolvent norm grows exponentially in $h^{-1}$, but near infinity it grows linearly. When $V$ is compactly supported, we obtain linear growth if the resolvent is multiplied by weights supported outside a ball of radius $CE^{-1/2}$ for some $C > 0$. This $E$-dependence is sharp and answers a question of Datchev and Jin., Comment: Update includes new section 5 and estimates when the potential may blow-up near 0

Gajic--Warnick have recently proposed a definition of scattering resonances based on Gevrey-2 regularity at infinity and introduced a new class of potentials for which resonances can be defined. We show that standard methods based on complex scaling apply to a slightly larger class of potentials and show existence of resonances in larger angles., Comment: 1 figure, 12 pages

We consider Keldysh-type operators, $ P = x_1 D_{x_1}^2 + a (x) D_{x_1} + Q (x, D_{x'} ) $, $ x = ( x_1, x') $ with analytic coefficients, and with $ Q ( x, D_{x'} ) $ second order, principally real and elliptic in $ D_{x'} $ for $ x $ near zero. We show that if $ P u =f $, $ u \in C^\infty $, and $ f $ is analytic in a neighbourhood of $ 0 $ then $ u $ is analytic in a neighbourhood of $ 0 $. This is a consequence of a microlocal result valid for operators of any order with Lagrangian radial sets. Our result proves a generalized version of a conjecture made by the second author and Lebeau and has applications to scattering theory., Comment: 24 pages, 1 figure

This work concerns $L^p$ norms of high energy Laplace eigenfunctions, $(-\Delta_g-\lambda^2)\phi_\lambda=0$, $\|\phi_\lambda\|_{L^2}=1$. In 1988, Sogge gave optimal estimates on the growth of $\|\phi_\lambda\|_{L^p}$ for a general compact Riemannian manifold. The goal of this article is to give general dynamical conditions guaranteeing quantitative improvements in $L^p$ estimates for $p>p_c$, where $p_c$ is the critical exponent. We also apply previous results of the authors to obtain quantitative improvements in concrete geometric settings including all product manifolds. These are the first results improving estimates for the $L^p$ growth of eigenfunctions that only require dynamical assumptions. In contrast with previous improvements, our assumptions are local in the sense that they depend only on the geodesics passing through a shrinking neighborhood of a given set in $M$. Moreover, the article gives a structure theorem for eigenfunctions which saturate the quantitatively improved $L^p$ bound. Modulo an error, the theorem describes these eigenfunctions as finite sums of quasimodes which, roughly, approximate zonal harmonics on the sphere scaled by $1/\sqrt{\log \lambda}$., Comment: 51 pages, 2 figures, added Theorem 2 describing the structure of near extremal eigenfunctions has been added

In this note, we prove weighted resolvent estimates for the semiclassical Schr\"odinger operator $-h^2 \Delta + V(x) : L^2(\mathbb{R}^n) \to L^2(\mathbb{R}^n)$, $n \neq 2$. The potential $V$ is real-valued, and assumed to either decay at infinity or to obey a radial $\alpha$-H\"older continuity condition, $0\leq \alpha \leq 1$, with sufficient decay of the local radial $C^\alpha$ norm toward infinity. Note, however, that in the H\"older case, the potential need \emph{not} decay. If the dimension $n \ge 3$, the resolvent bound is of the form $\exp \left(C h^{-1 - \frac{1 - \alpha}{3 + \alpha}} [(1-\alpha) \log(h^{-1})+c]\right)$, while for $n = 1$ it is of the form $\exp(Ch^{-1})$. A new type of weight and phase function construction allows us to reduce the necessary decay even in the pure $L^\infty$ case., Comment: 15 pages

We prove a uniform lower bound on Cauchy data on an arbitrary curve on a negatively curved surface using the Dyatlov-Jin(-Nonnenmacher) observability estimate on the global surface. In the process, we prove some further results about defect measures of restrictions of eigenfunctions to a hypersurface.

In this expository article we relate the presentation of weighted estimates in the book of Martinez to the Bergman kernel approach of Sj\"ostrand. It is meant as an introduction to the Helffer--Sj\"ostrand theory (designed for the study of quantum resonances) in the simplest setting and to its adaptations to compact manifolds.

Motivated by the work of Colin de Verdi\`ere and Saint-Raymond on spectral theory 0th order pseudodifferential operators on tori we consider viscosity limits in which 0th order operators $ P $ are replaced by $ P + i \nu \Delta $, $ \nu > 0 $. By adapting the Helffer--Sj\"ostrand theory of scattering resonances we show that in a complex neighbourhood of the continuous spectrum eigenvalues of $ P + i \nu \Delta $ have limits as viscosity $ \nu $ goes to 0. In the simplified setting of tori this justifies claims made in the physics literature., Comment: Proposition 6.3 added and some small issues corrected

This article concerns the asymptotic geometric character of the nodal set of the eigenfunctions of the Steklov eigenvalue problem $$ -\Delta \phi_{\sigma_j}=0,\quad\text{ on }\Omega,\qquad\qquad \partial_\nu \phi_{\sigma_j}=\sigma_j \phi_{\sigma_j}\quad \text{ on }\partial\Omega $$ in two-dimensional domains $\Omega$. In particular, this paper presents a dense family $\mathcal{A}$ of simply-connected two-dimensional domains with analytic boundaries such that, for each $\Omega\in \mathcal{A}$, the nodal set of the eigenfunction $\phi_{\sigma_j}$ "is $not$ dense at scale $\sigma_j^{-1}$". This result addresses a question put forth under "Open Problem 10" in Girouard and Polterovich, J. Spectr. Theory, 321-359 (2017). In fact, the results in the present paper establish that, for domains $\Omega\in \mathcal{A}$, the nodal sets of the eigenfunctions $\phi_{\sigma_j}$ associated with the eigenvalue $\sigma_j$ have starkly different character than anticipated: they are not dense at any shrinking scale. More precisely, for each $\Omega\in \mathcal{A}$ there is a value $r_1>0$ such that for each $j$ there is $x_j\in \Omega$ such that $\phi_{\sigma_j}$ does not vanish on the ball of radius $r_1$ around $x_j$., Comment: 21 Pages, 5 figures

In this article we develop new techniques for studying concentration of Laplace eigenfunctions $\phi_\lambda$ as their frequency, $\lambda$, grows. The method consists of controlling $\phi_\lambda(x)$ by decomposing $\phi_\lambda$ into a superposition of geodesic beams that run through the point $x$. Each beam is localized in phase-space on a tube centered around a geodesic whose radius shrinks slightly slower than $\lambda^{-\frac{1}{2}}$. We control $\phi_\lambda(x)$ by the $L^2$-mass of $\phi_\lambda$ on each geodesic tube and derive a purely dynamical statement through which $\phi_\lambda(x)$ can be studied. In particular, we obtain estimates on $\phi_\lambda(x)$ by decomposing the set of geodesic tubes into those that are non self-looping for time $T$ and those that are. This approach allows for quantitative improvements, in terms of $T$, on the available bounds for $L^\infty$ norms, $L^p$ norms, pointwise Weyl laws, and averages over submanifolds., Comment: 61 pages, 2 figures. Improved exposition and includes new explanatory material in the introduction as well as an examples section (1.5) and a full section on comparison with previous work (1.6). Appendices A.1 (Index of notation) and B were also added

We study the resolvent for nontrapping obstacles on manifolds with Euclidean ends. It is well known that for such manifolds, the outgoing resolvent satisfies $\|\chi R(k) \chi\|_{L^2\to L^2}\leq C{k}^{-1}$ for ${k}>1$, but the constant $C$ has been little studied. We show that, for high frequencies, the constant is bounded above by $2/\pi$ times the length of the longest generalized bicharacteristic of $|\xi|_g^2-1$ remaining in the support of $\chi.$ We show that this estimate is optimal in the case of manifolds without boundary. We then explore the implications of this result for the numerical analysis of the Helmholtz equation., Comment: 40 pages

Let $(M,g)$ be a compact Riemannian manifold and $P_1:=-h^2\Delta_g+V(x)-E_1$ so that $dp_1\neq 0$ on $p_1=0$. We assume that $P_1$ is quantum completely integrable in the sense that there exist functionally independent pseuodifferential operators $P_2,\dots P_n$ with $[P_i,P_j]=0$, $i,j=1,\dots ,n$. We study the pointwise bounds for the joint eigenfunctions, $u_h$ of the system $\{P_i\}_{i=1}^n$ with $P_1u_h=E_1u_h+o(1)$. We first give polynomial improvements over the standard H\"ormander bounds for typical points in $M$. In two and three dimensions, these estimates agree with the Hardy exponent $h^{-\frac{1-n}{4}}$ and in higher dimensions we obtain a gain of $h^{\frac{1}{2}}$ over the H\"ormander bound. In our second main result, under a real-analyticity assumption on the QCI system, we give exponential decay estimates for joint eigenfunctions at points outside the projection of invariant Lagrangian tori; that is at points $x\in M$ in the "microlocally forbidden" region $p_1^{-1}(E_1)\cap \dots \cap p_n^{-1}(E_n)\cap T^*_xM=\emptyset.$ These bounds are sharp locally near the projection of the invariant tori., Comment: 30 pages, 1 figure