Search

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.

In this article, we study propagation of defect measures for Schrödinger operators, $-h^2Δ_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.

We prove explicit asymptotics for the location of semiclassical scattering resonances in the setting of $h$-dependent delta-function potentials on $\mathbb{R}$. In the cases of two or three delta poles, we are able to show that resonances occur along specific lines of the form $\Im z \sim -\gamma h \log(1/h).$ More generally, we use the method of Newton polygons to show that resonances near the real axis may only occur along a finite collection of such lines, and we bound the possible number of values of the parameter $\gamma.$ We present numerical evidence of the existence of more and more possible values of $\gamma$ for larger numbers of delta poles., Comment: 16 pages, 3 figures, comments are welcome!

A static rotating cosmic string metric is singular along a timelike line and fails to be globally hyperbolic; these features make it difficult to solve the wave equation by conventional energy methods. Working on a single angular mode at a time, we use microlocal methods to construct forward parametrices for wave and Klein--Gordon equations on such backgrounds., 15 pages; improvements to exposition in response to referee comments

A crucial role in the theory of uncertainty quantification (UQ) of PDEs is played by the regularity of the solution with respect to the stochastic parameters; indeed, a key property one seeks to establish is that the solution is holomorphic with respect to (the complex extensions of) the parameters. In the context of UQ for the high-frequency Helmholtz equation, a natural question is therefore: how does this parametric holomorphy depend on the wavenumber $k$? The recent paper [Ganesh, Kuo, Sloan 2021] showed for a particular nontrapping variable-coefficient Helmholtz problem with affine dependence of the coefficients on the stochastic parameters that the solution operator can be analytically continued a distance $\sim k^{-1}$ into the complex plane. In this paper, we generalise the result in [Ganesh, Kuo, Sloan 2021] about $k$-explicit parametric holomorphy to a much wider class of Helmholtz problems with arbitrary (holomorphic) dependence on the stochastic parameters; we show that in all cases the region of parametric holomorphy decreases with $k$, and show how the rate of decrease with $k$ is dictated by whether the unperturbed Helmholtz problem is trapping or nontrapping. We then give examples of both trapping and nontrapping problems where these bounds on the rate of decrease with $k$ of the region of parametric holomorphy are sharp, with the trapping examples coming from the recent results of [Galkowski, Marchand, Spence 2021]. An immediate implication of these results is that the $k$-dependent restrictions imposed on the randomness in the analysis of quasi-Monte Carlo (QMC) methods in [Ganesh, Kuo, Sloan 2021] arise from a genuine feature of the Helmholtz equation with $k$ large (and not, for example, a suboptimal bound).

A subject of recent interest in inverse problems is whether a corner must diffract fixed frequency waves. We generalize this question somewhat and study cones $[0,\infty)\times Y$ which do not diffract high frequency waves. We prove that if $Y$ is analytic and does not diffract waves at high frequency then every geodesic on $Y$ is closed with period $2\pi$. Moreover, we show that if $\dim Y=2$, then $Y$ is isometric to either the sphere of radius 1 or its $\mathbb{Z}^2$ quotient, $\mathbb{R}\mathbb{P}^2$., Comment: Updated in response to referee comments. 13 pages

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$)., arXiv admin note: text overlap with arXiv:2102.13081

A rotating cosmic string spacetime has a singularity along a timelike curve corresponding to a one-dimensional source of angular momentum. Such spacetimes are not globally hyperbolic: they admit closed timelike curves near the string. This presents challenges to studying the existence of solutions to the wave equation via conventional energy methods. In this work, we show that semi-global forward solutions to the wave equation do nonetheless exist, but only in a microlocal sense. The main ingredient in this existence theorem is a propagation of singularities theorem that relates energy entering the string to energy leaving the string. The propagation theorem is localized in the fibers of a certain fibration of the blown-up string, but global in time, which means that energy entering the string at one time may emerge previously., Comment: 58 pgs; Lemma 9.3 has been updated

We obtain estimates on the rate of decay of a solution to the wave equation on a stationary spacetime that tends to Minkowski space at a rate $O(\lvert x \rvert^{-\kappa}),$ $\kappa \in (1,\infty) \backslash \mathbb{N}.$ Given suitably smooth and decaying initial data, we show a wave locally enjoys the decay rate $O(t^{-\kappa-2+\epsilon})$., Comment: 36 pages; updated with clarifying exposition and fixed typos

It is well known that when the geometry and/or coefficients allow stable trapped rays, the outgoing solution operator of the Helmholtz equation (a.k.a. the resolvent of the Laplacian) grows exponentially through a sequence of real frequencies tending to infinity. In this paper we show that, even in the presence of the strongest-possible trapping, if a set of frequencies of arbitrarily small measure is excluded, the Helmholtz solution operator grows at most polynomially as the frequency tends to infinity. One significant application of this result is in the convergence analysis of several numerical methods for solving the Helmholtz equation at high frequency that are based on a polynomial-growth assumption on the solution operator (e.g. $hp$-finite elements, $hp$-boundary elements, certain multiscale methods). The result of this paper shows that this assumption holds, even in the presence of the strongest-possible trapping, for most frequencies.

We establish propagation of singularities for the semiclassical Schr\"odinger equation, where the potential is conormal to a hypersurface. We show that semiclassical wavefront set propagates along generalized broken bicharacteristics, hence reflection of singularities may occur along trajectories reaching the hypersurface transversely. The reflected wavefront set is weaker, however, by a power of $h$ that depends on the regularity of the potential. We also show that for sufficiently regular potentials, wavefront set may not stick to the hypersurface, but rather detaches from it at points of tangency to travel along ordinary bicharacteristics., Comment: 99 pages; exposition revised in response to referee comments

In this survey, we review some applications and extensions of the author's results with Richard Melrose on propagation of singularities for solutions to the wave equation on manifolds with conical singularities. These results mainly concern: the local decay of energy on noncompact manifolds with diffractive trapped orbits (joint work with Dean Baskin); singularities of the wave trace created by diffractive closed geodesics (joint work with G. Austin Ford); and the distribution of scattering resonances associated to such closed geodesics (joint work with Luc Hillairet)., Comment: 15 pages; contribution to Seminaire Laurent Schwartz proceedings

Let $(X,g)$ be a compact manifold with conic singularities. Taking $\Delta_g$ to be the Friedrichs extension of the Laplace-Beltrami operator, we examine the singularities of the trace of the half-wave group $e^{- i t \sqrt{ \smash[b]{\Delta_g}}}$ arising from strictly diffractive closed geodesics. Under a generic nonconjugacy assumption, we compute the principal amplitude of these singularities in terms of invariants associated to the geodesic and data from the cone point. This generalizes the classical theorem of Duistermaat-Guillemin on smooth manifolds and a theorem of Hillairet on flat surfaces with cone points., Comment: 46 pages; 3 figures. New version reflects correction to trace formula (factor of primitive length rather than length), as well as some corrections to discussions of Jacobi fields in section 1.3

We consider a non-trapping $n$-dimensional Lorentzian manifold endowed with an end structure modeled on the radial compactification of Minkowski space. We find a full asymptotic expansion for tempered forward solutions of the wave equation in all asymptotic regimes. The rates of decay seen in the asymptotic expansion are related to the resonances of a natural asymptotically hyperbolic problem on the "northern cap" of the compactification. For small perturbations of Minkowski space that fit into our framework, we show a rate of decay that improves on the Klainerman--Sobolev estimates., Comment: 67 pages, 2 figures; version 2 was a substantial revision with the main result (Theorem 1.1) improved; version 3 fixes a minor error in the iterative argument establishing the asymptotic expansion; version 4 includes a few cosmetic changes and adds a few references

Using a new local smoothing estimate of the first and third authors, we prove local-in-time Strichartz and smoothing estimates without a loss exterior to a large class of polygonal obstacles with arbitrary boundary conditions and global-in-time Strichartz estimates without a loss exterior to a large class of polygonal obstacles with Dirichlet boundary conditions. In addition, we prove a global-in-time local smoothing estimate in exterior wedge domains with Dirichlet boundary conditions and discuss some nonlinear applications., Comment: 15 pages. This paper has been substantially revised; the results have been improved to include global-in-time estimates in the Dirichlet case. Further revised to include more related scattering references thanks to a careful reading by Fabrice Planchon

We consider a family of surfaces of revolution, each with a single periodic geodesic which is degenerately unstable. We prove a local smoothing estimate for solutions to the linear Schr��dinger equation with a loss that depends on the degeneracy, and we construct explicit examples to show our estimate is saturated on a weak semiclassical time scale. As a byproduct of our proof, we obtain a cutoff resolvent estimate with a sharp polynomial loss., 26 pages, incorporated referee's comments

We consider the fundamental solution to the wave equation on a manifold with corners of arbitrary codimension. If the initial pole of the solution is appropriately situated, we show that the singularities which are diffracted by the corners (i.e., loosely speaking, are not propagated along limits of transversely reflected rays) are smoother than the main singularities of the solution. More generally, we show that subject to a hypothesis of nonfocusing, diffracted wavefronts of any solution to the wave equation are smoother than the incident singularities. These results extend our previous work on edge manifolds to a situation where the fibers of the boundary fibration, obtained here by blowup of the corner in question, are themselves manifolds with corners., Comment: Revision adding significant exposition, correcting many small errors; index of notation added

Lecture notes from 2008 CMI/ETH Summer School on Evolution Equations. These notes are an informal introduction to the applications of microlocal methods in the study of linear evolution equations and spectral theory. Calculi of pseudodifferential operators and Fourier integral operators are discussed and axiomatized, but not constructed: the focus is on how to apply these tools., Comment: Corrects an error in statement of Theorem 9.2 kindly pointed out by Amir Vig

We obtain the Strichartz inequalities $$ \| u \|_{L^q_t L^r_x([0,1] \times M)} \leq C \| u(0) \|_{L^2(M)}$$ for any smooth $n$-dimensional Riemannian manifold $M$ which is asymptotically conic at infinity (with either short-range or long-range metric perturbation) and non-trapping, where $u$ is a solution to the Schr��dinger equation $iu_t + {1/2} ��_M u = 0$, and $2 < q, r \leq \infty$ are admissible Strichartz exponents ($\frac{2}{q} + \frac{n}{r} = \frac{n}{2}$). This corresponds with the estimates available for Euclidean space (except for the endpoint $(q,r) = (2, \frac{2n}{n-2})$ when $n > 2$). These estimates imply existence theorems for semi-linear Schr��dinger equations on $M$, by adapting arguments from Cazenave and Weissler \cite{cwI} and Kato \cite{kato}. This result improves on our previous result in \cite{HTW}, which was an $L^4_{t,x}$ Strichartz estimate in three dimensions. It is closely related to the results of Staffilani-Tataru, Burq, Tataru, and Robbiano-Zuily, who consider the case of asymptotically flat manifolds., 50 pages, 2 figures

We obtain an $L^4$ space-time Strichartz inequality for any smooth three-dimensional Riemannian manifold $(M,g)$ which is asymptotically conic at infinity and non-trapping, where $u$ is a solution to the Schr\"odinger equation $iu_t + {1/2} \Delta_M u = 0$. The exponent $H^{1/4}(M)$ is sharp, by scaling considerations. In particular our result covers asymptotically flat non-trapping manifolds. Our argument is based on the interaction Morawetz inequality introduced by Colliander et al., interpreted here as a positive commutator inequality for the tensor product $U(t,z',z'') := u(t,z') u(t,z'')$ of the solution with itself. We also use smoothing estimates for Schr\"odinger solutions including a new one proved here with weight $r^{-1}$ at infinity and with the gradient term involving only one angular derivative., Comment: 43 pages, 1 figure; minor corrections to earlier version

We discuss the form of the propagator $U(t)$ for the time-dependent Schr��dinger equation on an asyptotically Euclidean, or, more generally, asymptotically conic, manifold with no trapped geodesics. In the asymptotically Euclidean case, if $��\in \mathcal{C}_0^\infty$, and with $\mathcal{F}$ denoting Fourier transform, $\mathcal{F}\circ e^{-ir^2/2t} U(t) ��$ is a Fourier integral operator for $t\neq 0.$ The canonical relation of this operator is a ``sojourn relation'' associated to the long-time geodesic flow. This description of the propagator follows from its more precise characterization as a ``scattering fibered Legendrian,'' given by the authors in a previous paper and sketched here. A corollary is a propagation of singularities theorem that permits a complete description of the wavefront set of a solution to the Schr��dinger equation, restricted to any fixed nonzero time, in terms of the oscillatory behavior of its initial data. We discuss two examples which illustrate some extremes of this propagation behavior.

Let g be a scattering metric on a compact manifold X with boundary, i.e., a smooth metric giving the interior of X the structure of a complete Riemannian manifold with asymptotically conic ends. An example is any compactly supported perturbation of the standard metric on Euclidean space. Consider the operator $H = \half \Delta + V$, where $\Delta$ is the positive Laplacian with respect to g and V is a smooth real-valued function on X vanishing to second order at the boundary. Assuming that g is non-trapping, we construct a global parametrix for the kernel of the Schroedinger propagator $U(t) = e^{-itH}$ and use this to show that the kernel of U(t) is, up to an explicit quadratic oscillatory factor, a class of `Legendre distributions' on $X \times X^{\circ} \times \halfline$ previously considered by Hassell-Vasy. When the metric is trapping, then the parametrix construction goes through microlocally in the non-trapping part of the phase space. We apply this result to obtain a microlocal characterization of the singularities of $U(t) f$, for any tempered distribution $f$ and any fixed $t \neq 0$, in terms of the oscillation of f near the boundary of X. If the metric is non-trapping, then we obtain a complete characterization; more generally we need to assume that f is microsupported in the nontrapping part of the phase space. This generalizes results of Craig-Kappeler-Strauss and Wunsch., Comment: 29 pages, no figures