BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jeffrey Galkowski

Jeffrey Galkowski appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.APmath-phmath.MPmath.SPmath.NANumerical Analysismath.DGmath.DS

Catalog footprint

What is connected

17works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2026arXiv

Sharp error bounds for edge-element discretisations of the high-frequency Maxwell equations

We prove sharp wavenumber-explicit error bounds for first- or second-family-Nédélec-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-family Nédélec elements or higher-than-lowest-order second-family Nédélec elements. Furthermore, they are the first wavenumber-explicit results, even in the asymptotic regime, for Maxwell scattering problems with a non-empty scatterer.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Decompositions of high-frequency Helmholtz solutions via functional calculus, and application to the finite element method

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östrand functional calculus, this paper proves analogous decompositions for scattering problems fitting into the black-box scattering framework of Sjöstrand-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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Does the Helmholtz boundary element method suffer from the pollution effect?

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).

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Logarithmic improvements in the Weyl law and exponential bounds on the number of closed geodesics are predominant

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 $λ^2$ is asymptotic to $Cλ^d$ with an $O(λ^{d-1})$ error. We show that, for a predominant metric, the estimate on the error can by improved by a power of $\log λ$. 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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Semiclassical resolvent bounds for long range Lipschitz potentials

We give an elementary proof of weighted resolvent estimates for the semiclassical Schrödinger operator $-h^2 Δ+ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Weyl remainders: an application of geodesic beams

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 $Π_λ$ denote the kernel of the spectral projector for the Laplacian, $\mathbb{1}_{[0,λ^2]}(-Δ_g)$. Assuming only that the set of near periodic geodesics over $W\subset M$ has small measure, we prove that as $λ\to \infty$ $$ \int_{W} Π_λ(x,x)dx=(2π)^{-n}\text{vol}_{\mathbb{R}^n}(B)\text{vol}_g(W)\,λ^n+O\Big(\frac{λ^{n-1}}{\log λ}\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 $Π_λ(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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Complete asymptotic expansions of the spectral function for symbolic perturbations of almost periodic Schrödinger operators in dimension one

In this article we consider asymptotics for the spectral function of Schrödinger 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,λ^2]}(P)$ has a full asymptotic expansion in powers of $λ$. 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Eigenfunction concentration via geodesic beams

In this article we develop new techniques for studying concentration of Laplace eigenfunctions $ϕ_λ$ as their frequency, $λ$, grows. The method consists of controlling $ϕ_λ(x)$ by decomposing $ϕ_λ$ 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 $λ^{-\frac{1}{2}}$. We control $ϕ_λ(x)$ by the $L^2$-mass of $ϕ_λ$ on each geodesic tube and derive a purely dynamical statement through which $ϕ_λ(x)$ can be studied. In particular, we obtain estimates on $ϕ_λ(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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Semiclassical resolvent bounds for weakly decaying potentials

In this note, we prove weighted resolvent estimates for the semiclassical Schrödinger operator $-h^2 Δ+ 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 $α$-Hölder continuity condition, $0\leq α\leq 1$, with sufficient decay of the local radial $C^α$ norm toward infinity. Note, however, that in the Hölder 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 - α}{3 + α}} [(1-α) \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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Viscosity limits for 0th order pseudodifferential operators

Motivated by the work of Colin de Verdière 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 νΔ$, $ ν> 0 $. By adapting the Helffer--Sjöstrand theory of scattering resonances we show that in a complex neighbourhood of the continuous spectrum eigenvalues of $ P + i νΔ$ have limits as viscosity $ ν$ goes to 0. In the simplified setting of tori this justifies claims made in the physics literature.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Distribution of Resonances in Scattering by Thin Barriers

We study high energy resonances for the operators $-Δ+δ_{\partialΩ}\otimes V$ and $-Δ+δ_{\partialΩ}'\otimes V\partial_ν$ where $Ω$ is strictly convex with smooth boundary, $V:L^2(\partialΩ)\to L^2(\partialΩ)$ may depend on frequency, and $δ_{\partialΩ}$ is the surface measure on $\partialΩ$. These operators are model Hamiltonians for quantum corrals and leaky quantum graphs. We give a quantum version of the Sabine Law from the study of acoustics for both the $δ$ and $δ'$ interactions. It characterizes the decay rates (imaginary parts of resonances) in terms of the system's ray dynamics. In particular, the decay rates are controlled by the average reflectivity and chord length of the barrier. For the $δ$ interaction we show that generically there are infinitely many resonances arbitrarily close to the resonance free region found by our theorem. In the case of the $δ'$ interaction, the quantum Sabine law gives the existence of a resonance free region that converges to the real axis at a fixed polynomial rate and is optimal in the case of the unit disk in the plane. As far as the author is aware, this is the only class of examples that is known to have resonances converging to the real axis at a fixed polynomial rate but no faster. The proof of our theorem requires several new technical tools. We adapt intersecting Lagrangian distributions to the semiclassical setting and give a description of the kernel of the free resolvent as such a distribution. We also construct a semiclassical version of the Melrose--Taylor parametrix for complex energies. We use these constructions to give a complete microlocal description of boundary layer operators and to prove sharp high energy estimates on the boundary layer operators in the case that $\partialΩ$ is smooth and strictly convex.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quantum Ergodicity for a Class of Mixed Systems

We examine high energy eigenfunctions for the Dirichlet Laplacian on domains where the billiard flow exhibits mixed dynamical behavior. (More generally, we consider semiclassical Schrodinger operators with mixed assumptions on the Hamiltonian flow.) Specificially, we assume that the billiard flow has an invariant ergodic component, U, and study defect measures, mu, of positivie density subsequences of eigenfunctions (and, more generally, of almost orthogonal quasimodes). We show that any defect measure associated to such a subsequence, when restricted to U, satisfies mu = c mu_L where mu_L is the Liouville measure. This proves part of a conjecture of Percival.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sharp norm estimates of layer potentials and operators at high frequency

In this paper, we investigate single and double layer potentials mapping boundary data to interior functions of a domain at high frequency $λ^2\to\infty$. For single layer potentials, we find that the $L^{2}(\partialΩ)\to{}L^{2}(Ω)$ norms decay in $λ$. The rate of decay depends on the curvature of $\partialΩ$: The norm is $λ^{-3/4}$ in general domains and $λ^{-5/6}$ if the boundary $\partialΩ$ is curved. The double layer potential, however, displays uniform $L^{2}(\partialΩ)\to{}L^{2}(Ω)$ bounds independent of curvature. By various examples, we show that all our estimates on layer potentials are sharp. The appendix by Galkowski gives bounds $L^{2}(\partialΩ)\to{}L^{2}(\partialΩ)$ for the single and double layer operators at high frequency that are sharp modulo $\log λ$. In this case, both the single and double layer operator bounds depend upon the curvature of the boundary.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nonlinear Instability in a Semiclassical Problem

We consider a nonlinear evolution problem with an asymptotic parameter and construct examples in which the linearized operator has spectrum uniformly bounded away from Re z >= 0 (that is, the problem is spectrally stable), yet the nonlinear evolution blows up in short times for arbitrarily small initial data. We interpret the results in terms of semiclassical pseudospectrum of the linearized operator: despite having the spectrum in Re z < -c < 0, the resolvent of the linearized operator grows very quickly in parts of the region Re z > 0. We also illustrate the results numerically.

0 citations0 reviewsNo code linkedOpen access