A Polynomial Rate of Convergence for the Dirichlet Problem on Orthodiagonal Maps
David Pechersky
TL;DR
The paper proves a polynomial rate of convergence for the Dirichlet problem on orthodiagonal maps to its continuum counterpart as the mesh $\varepsilon\to0$, for Hölder boundary data. The core method combines mollification of discrete harmonic functions to produce almost-harmonic functions, Beurling-type estimates to control boundary-influenced behavior, and Green-function representations to compare discrete and continuous solutions. A principal contribution is the explicit, though not sharp, rate bound $|h^{\bullet}(v)-h(v)| \le C (\varepsilon/\mathrm{diam}(\Omega))^{\lambda(\alpha,\beta)}$ with an exponent $\lambda(\alpha,\beta)$ defined via a minimax expression, plus a robust mesoscopic Lipschitz regularity result for discrete harmonic functions on OD maps. The results extend prior work by Gurel-Gurevich–Jerison–Nachmias and Bou-Rabee–Gwynne, and provide a framework for quantitative convergence in discrete complex analysis on 2D lattice discretizations, with implications for critical lattice models. Overall, the paper advances the quantitative understanding of discretization convergence in planar domains and introduces tools that may improve rates via a bootstrap of mesoscopic regularity.
Abstract
We extend recent work of Gurel-Gurevich--Jerison--Nachmias (2020) and Bou-Rabee--Gwynne (2024) by showing that as the mesh of our lattice tends to $0$, we have a polynomial rate of convergence for the Dirichlet problem on orthodiagonal maps with Hölder boundary data to its continuous counterpart. The key idea is that the convolution of a discrete harmonic function on an orthodiagonal map with a smooth mollifier has small Laplacian and so is ``almost harmonic." This also allows us to show that discrete harmonic functions on orthodiagonal maps are Lipschitz in the bulk on a mesoscopic scale.
