BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Jeffrey Galkowski

Jeffrey Galkowski contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.APmath-phmath.MPmath.SPmath.NANumerical Analysismath.DGmath.DS

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

13 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