\section{Introduction}\label{sec:intro}

\subsection{Boundary integral equations for Laplace's equation}
If an explicit expression for the fundamental solution of a linear PDE is known, then boundary value problems (BVPs) for that PDE can be converted to integral equations on the boundary of the domain.
The main advantage of this procedure is that the dimension of the problem is reduced; indeed, the problem is converted from one on a $d$-dimensional domain to one on a $(d-1)$-dimensional domain. Futhermore, if the original domain is the exterior of a bounded obstacle, then the problem is reduced from one on a $d$-dimensional \emph{infinite} domain, to one on a $(d-1)$-dimensional \emph{finite} domain.

This reduction to the boundary has both theoretical and practical benefits:~on the theoretical side, C.~Neumann famously used boundary integral equations (BIEs) to prove existence of the solution of the Dirichlet problem for Laplace's equation in convex domains in \cite{Ne:77} (see, e.g., the account in \cite[Chapter 1]{Mc:00}), and BIEs have a long history of use in the harmonic analysis literature to prove wellposedness of BVPs on rough domains (see, e.g., \cite{CoMcMe:82}, \cite{Ve:84}, \cite{Ca:85}, \cite[\S2.1]{Ke:94}, \cite{MitreaTaylor:99}, \cite[Chapter 15]{MeCo:00}, \cite[Chapter 4]{Ta:00}, \cite{MiMiTa:01}).
On the more practical side, numerical methods based on Galerkin, collocation, or numerical quadrature discretisation of BIEs, coupled with fast matrix-vector multiply and compression algorithms, and iterative solvers such as GMRES, provide spectacularly effective computational tools for solving a range of linear boundary value problems, for example in potential theory, elasticity, and acoustic and electromagnetic wave scattering (see, e.g., \cite{Rokhlin83,At:97, LaSc:99,BrKu:01,ContETAL02, ChSoCuVeHa:04,BoSa:04,XiTaWe:08,GrGuMaRo09,SaSc:11, ChDaLo:17}).

Let $\Phi(\mathbf{x},\mathbf{y})$ be the fundamental solution for Laplace's equation:
\begin{equation}\label{eq:fund}
\Phi(\mathbf{x},\mathbf{y}):=
\displaystyle{\frac{1}{2\pi} \log \left(\frac{a}{|\bx-\by|}\right),}  \quad d= 2,
\qquad :=
\dfrac{1}{(d-2)C_d|\mathbf{x}-\mathbf{y}|^{d-2}},  \quad d\geq 3,
\end{equation}
where $C_d$ is the surface area of the unit sphere $S^{d-1}\subset \mathbb{R}^d$ and  $a >0$.
With $\Gamma$ the boundary of a bounded Lipschitz domain,
the boundary integral operators (BIOs) $S$, $D$, $D'$, and $H$,
the \emph{single-layer}, \emph{double-layer}, \emph{adjoint double-layer}, and \emph{hypersingular} operators, respectively,
 are defined for $\phi\in{L^2(\Gamma)}$, $\psi\inH^1(\Gamma)$, and $\mathbf{x}\in\Gamma$ by
\begin{equation}\label{eq:bio1}
S_k \phi(\mathbf{x}) = \int_\Gamma \Phi_k}%{\rm MH}(\mathbf{x},\mathbf{y}) \phi(\mathbf{y})\, {\rm d} s(\mathbf{y}), \quad
D \phi(\mathbf{x}) = \int_\Gamma \pdiff{\Phi(\mathbf{x},\mathbf{y})}{n(\mathbf{y})} \phi(\mathbf{y})\, {\rm d} s(\mathbf{y}), \quad
\end{equation}
and
\begin{equation}\label{eq:bio2}
D' \phi(\mathbf{x}) = \int_\Gamma \pdiff{\Phi(\mathbf{x},\mathbf{y})}{n(\mathbf{x})} \phi(\mathbf{y})\, {\rm d} s(\mathbf{y}), \quad
H \psi(\mathbf{x}) = \pdiff{}{n(\mathbf{x})}\int_\Gamma \pdiff{\Phi_k}%{\rm MH}(\mathbf{x},\mathbf{y})}{n(\mathbf{y})} \psi(\mathbf{y})\, {\rm d} s(\mathbf{y}).
\end{equation}
When $\Gamma$ is Lipschitz, the integrals in $D$ and $D'$ are defined as Cauchy principal values, in general only for almost all $\mathbf{x}\in \Gamma$ with respect to the surface measure ${\rm d} s$.
The definition of $H$ on spaces larger than $H^1(\Gamma)$ is complicated (it must be understood either as a finite-part integral, or as the non-tangential limit of a potential; see \cite[Chapter 7]{Mc:00}, \cite[Page 113]{ChGrLaSp:12} respectively), but these details are not essential to the present paper. The standard mapping properties of $S, D, D'$, and $H$ on Sobolev spaces on $\Gamma$ are recalled in Appendix \ref{app:B} (see \eqref{eq:map}).

The BIE operators involved in the standard first- and second-kind BIEs for the Dirichlet and Neumann problems for  Laplace's equation are shown in Table \ref{tab:bies}; although we do not explicitly consider the Neumann problem in this paper, we use the information in this table in what follows.
For the details of the right-hand sides and unknowns for the integral equations corresponding to the operators in Table \ref{tab:bies}, see, e.g., \cite[\S3.4]{SaSc:11}, \cite[Chapter 7]{Mc:00}, \cite[Chapter 7]{St:08}, \cite[\S2.5]{ChGrLaSp:12}.
Recall that the adjective ``direct" in the table refers to equations where the unknown is either the Dirichlet or Neumann trace of the solution to the corresponding BVP, and the adjective ``indirect" refers to equations where the unknown does not have immediate physical relevance.

\begin{table}
\begin{center}
\begin{tabular}{|c|c|c|c|c|
\hlin
&Interior Dirichlet   & Interior Neumann & Exterior Dirichlet  & Exterior Neumann  \\
& problem & problem & problem & problem \\
\hline
 Direct  & $S$ & $H$ &  $S$ & $H$\\
 & $\displaystyle{ \frac{1}{2}I- D^\prime}$   & $\displaystyle{ \frac{1}{2}I+ D}$  &  $\displaystyle{ \frac{1}{2}I+ D'}$ & $\displaystyle{ \frac{1}{2}I- D}$\\
 \hline
 Indirect  & $S$ & $H$ &  $S$ & $H$\\
& $\displaystyle{ \frac{1}{2}I- D}$ & $\displaystyle{ \frac{1}{2}I+ D'}$& $\displaystyle{ \frac{1}{2}I+ D}$ & $\displaystyle{ \frac{1}{2}I- D'}$\\
\hlin
     \end{tabular}
     \end{center}
     \caption{
The integral operators involved in the standard boundary-integral-equation formulations of the interior and exterior Dirichlet and Neumann problems for Laplace's equation.\label{tab:bies}}
\end{table}

Following \cite[Pages 9 and 10]{SaSc:11},
we call BIEs \emph{first kind} where the unknown function only appears under the integral, and \emph{second kind} where the unknown function appears outside the integrand as well as inside; by this definition, the BIEs in the first and third row of Table \ref{tab:bies} are first kind, and in the second and fourth row second kind.
An alternative definition of second kind BIEs is that, in addition to the unknown function appearing outside the integrand as well as inside, the BIO is Fredholm of index zero (i.e., the Fredholm alternative applies to the BIE); see, e.g., \cite[\S1.1.4]{At:97}.
Every BIE that we describe in the paper as second-kind is second-kind in both senses above.

\subsection{The Galerkin method} \label{sec:Gal}

We focus on solving Laplace BIEs with the \emph{Galerkin method}. The Galerkin method  applied to the equation $A\phi = f$, where $\phi, f \in {\cal H}$, $A:{\cal H}\rightarrow {\cal H}$ is a continuous (i.e.~bounded) linear operator, and ${\cal H}$ is a complex\footnote{It is convenient, since we deal with non-self-adjoint operators and talk at some points about spectra and numerical ranges, to assume throughout that all Hilbert spaces and function spaces are complex. Of course results for the corresponding real case are easily deduced, if needed, from the complex function space case.}
Hilbert space, is:~given a sequence $({\cal H}_N)_{N=1}^\infty$ of finite-dimensional subspaces of ${\cal H}$ with $\dim({\cal H}_N)\rightarrow \infty$ as $N\rightarrow \infty$,
\begin{equation}\label{eq:G}
\text{ find } \phi_N \in{\cal H}_N \text{ such that } \big(A\phi_N,\psi_N)_{\cal H}=\big(f,\psi_N\big)_{\cal H} \quad \text{ for all } \psi_N\in {\cal H}_N.
\end{equation}
We say that the \emph{Galerkin method converges for the sequence $({\cal H}_N)_{N=1}^\infty$} if, for every $f\in {\cal H}$, the Galerkin equations \eqref{eq:G} have a unique solution for all sufficiently large $N$ and $\phi_N\to A^{-1}f$ as $N\to\infty$.
We say that \emph{$({\cal H}_N)_{N=1}^\infty$ is asymptotically dense in ${\cal H}$} if, for every $\phi\in {\cal H}$,
\begin{equation}\label{eq:space_converge}
\min_{\psi_N\in {\cal H}_N}\N{\phi-\psi_N}_{{\cal H}} \to 0 \quad \mbox{as} \quad N\to \infty.
\end{equation}

A necessary condition for the convergence of the Galerkin method is that $({\cal H}_N)_{N=1}^\infty$ is asymptotically dense in ${\cal H}$.
Indeed, a standard necessary and sufficient condition  (e.g., \cite[Chapter II, Theorem 2.1]{GoFe:74}) for convergence of the Galerkin method is that  $({\cal H}_N)_{N=1}^\infty$ is asymptotically dense and that, for some $N_0\in \mathbb{N}$ and $C_{\mathrm{dis}}>0$,
\begin{equation} \label{eq:d_infsup}
 \frac{\|{\cal P}_N A\psi_N\|_{{\cal H}}}{\|\psi_N\|_{\cal H}} \geq C_{\mathrm{dis}} \quad \mbox{for all non-zero }\psi_N\in {\cal H}_N \mbox{ and } N\geq N_0,
\end{equation}
where ${\cal P}_N$ is orthogonal projection of ${\cal H}$ onto ${\cal H}_N$. Importantly, if \eqref{eq:d_infsup} holds, then (\cite[Chapter II, Equation (2.5)]{GoFe:74} or see \cite[Theorem 4.2.1 and Remark 4.2.5]{SaSc:11})
\begin{equation} \label{eq:quasioptimal2}
\|\phi- \phi_N\|_{{\cal H}} \leq \left(1+ \frac{\|A\|_{{\cal H}\to{\cal H}}}{C_{\mathrm{dis}}}\right)\min_{\psi_N\in {\cal H}_N}\|\phi-\psi_N\|_{\cal H}, \quad \mbox{for } N\geq N_0,
\end{equation}
where $\phi=A^{-1}f$ and $\phi_N$ is the unique solution of the Galerkin equations \eqref{eq:G}.
We note that \eqref{eq:quasioptimal2} is known as a \emph{quasioptimal} error estimate.

We now recap the main abstract theorem on convergence of the Galerkin method; this theorem uses the definition that
an operator $A:{\cal H} \rightarrow {\cal H}$ is \emph{coercive}\footnote{In the literature, the property \eqref{eq:coer} (and its analogue for operators $A: {\cal H}\rightarrow {\cal H}'$, where ${\cal H}'$ is the dual of ${\cal H}$) is sometimes called ``$\cH$-ellipticity" (as in, e.g., \cite[Page 39]{SaSc:11}, \cite[\S3.2]{St:08}, and \cite[Definition 5.2.2]{HsWe:08}) or ``strict coercivity" (e.g., \cite[Definition 13.22]{Kr:89}),  with ``coercivity" then used to mean \emph{either} that $A$ is the sum of a coercive operator and a compact operator (as in, e.g., \cite[\S3.6]{St:08} and \cite[\S5.2]{HsWe:08}) \emph{or} that $A$ satisfies a G\aa rding inequality (as in \cite[Definition 2.1.54]{SaSc:11}).} if there exists $C_{\rm coer}>0$ such that
\begin{equation}\label{eq:coer}
\big|(A\psi,\psi)_{{\cal H}}\big|\geq C_{\rm coer} \N{\psi}^2_{{\cal H}} \quad \text{ for all } \psi \in {\cal H}.
\end{equation}

\begin{theorem}[The main abstract theorem on convergence of the Galerkin method.]\label{thm:Galerkin}

\

\begin{enumerate}
\item[(a)]
If $A$ is invertible then there exists a sequence $({\cal H}_N)_{N=1}^\infty$ for which the Galerkin method converges.
\item[(b)]
If $A$ is invertible then the following are equivalent:
\begin{enumerate}
\item[(i)]
The Galerkin method converges for every asymptotically-dense sequence $({\cal H}_N)_{N=1}^\infty$ in ${\cal H}$.
\item[(ii)] $A=A_0+K$ where $A_0$ is coercive and $K$ is compact.
\end{enumerate}
\item[(c)]
 If $A$ is coercive (i.e.~\eqref{eq:coer} holds) then, for every sequence $({\cal H}_N)_{N=1}^\infty$ and every $N\in \mathbb{N}$, the Galerkin equations \eqref{eq:G} have a unique solution $\phi_N$ and, where $\phi = A^{-1}f$,
\begin{equation}\label{eq:quasioptimal}
\big\|\phi-\phi_N\big\|_{\cal H} \leq \frac{\|A\|_{{\cal H}\rightarrow{\cal H}}}{C_{\rm coer}} \,  \min_{\psi \in {\cal H}_N}\big\|\phi-\psi\big\|_{\cal H},
\end{equation}
(so that $\phi_N\to \phi$ as $N\to\infty$ if $({\cal H}_N)_{N=1}^\infty$ is asymptotically dense in ${\cal H}$).
\end{enumerate}
\end{theorem}

\begin{proof}[References for the proof]
Part (a) was first proved in \cite[Theorem 1]{Ma:74}; see also \cite[Chapter II, Theorem 4.1]{GoFe:74}.
Part (b) was first proved in \cite[Theorem 2]{Ma:74}, with this result building on results in \cite{Va:65}; see also \cite[Chapter II, Lemma 5.1 and Theorem 5.1]{GoFe:74}.
Part (c) is C\'ea's Lemma, first proved in \cite{Ce:64}.
\end{proof}

\subsection{The rationale for using second-kind BIEs posed in ${L^2(\Gamma)}$}\label{sec:rationale}

The BIOs in Table \ref{tab:bies} are coercive in the trace spaces $H^{\pm 1/2}(\Gamma)$ (or certain subspaces of these) for Lipschitz $\Gamma$, thus insuring convergence of the associated Galerkin methods by Part (c) of Theorem \ref{thm:Galerkin}; this coercivity theory was established for first-kind equations by N\'ed\'elec and Planchard \cite{NePl:73}, Le Roux \cite{LeRo:74}, \cite{LeRo:77}, and Hsiao and Wendland \cite{HsWe:77},
and for second-kind equations by Steinbach and Wendland \cite{StWe:01}. These arguments involve transferring boundedness/coercivity properties of the PDE solution operator to the associated boundary integral operators
via the trace map and layer potentials; the generality of these arguments is why coercivity holds with $\Gamma$ only assumed to be Lipschitz, and Costabel \cite{Co:07} highlighted how these ideas can be traced back to the work of Gauss and Poincar\'e.

Despite convergence of the associated Galerkin methods, using the first-kind formulations in the trace spaces has the disadvantage that
the condition numbers of the Galerkin matrices grow as the discretisation is refined; e.g., for the $h$-version of the Galerkin method (where convergence is obtained by decreasing the mesh-width $h$ and keeping the polynomial degree fixed), the condition numbers grow like $h^{-1}$; see, e.g. \cite[\S4.5]{SaSc:11}.
Furthermore, using the second-kind formulations in the trace spaces has the disadvantage that the inner products on $H^{\pm 1/2}(\Gamma)$ are non-local and non-trivial to evaluate; even if the basis functions $\phi_N$ and $\psi_N$ in \eqref{eq:G} have supports only on a subset of $\Gamma$, $(A\phi_N,\psi_N)_{{\cal H}}$ is an integral over all of $\Gamma$, and the calculation of the Galerkin matrix in this case is impractical.

For the second-kind BIEs, an attractive alternative to working in the trace spaces is to work in ${L^2(\Gamma)}$.
When $\Gamma$ is $C^1$, $D$ and $D'$ are compact in ${L^2(\Gamma)}$ by the results of Fabes, Jodeit, and Rivi\`ere \cite[Theorems 1.2 and 1.9]{FaJoRi:78} and thus each of the second-kind BIOs $\frac{1}{2} I \pm D$ and $\frac{1}{2} I \pm D'$ is the sum of
a coercive operator and a compact operator, and convergence of the associated Galerkin methods in ${L^2(\Gamma)}$ is ensured by Part (b) of Theorem \ref{thm:Galerkin}. Since the ${L^2(\Gamma)}$ norm is local, $(A\phi_N,\psi_N)_{{\cal H}}$ is an integral over the support of $\psi_N$, and the Galerkin matrix is much more easily computable. Furthermore, when $D$ and $D'$ are compact, the condition numbers of the Galerkin matrices of  $\frac{1}{2} I \pm D$ and $\frac{1}{2} I \pm  D'$ are independent of the discretisation (without preconditioning); see \cite[\S3.6.3]{At:97}, \cite[\S4.5.5]{Ha:95}.

\subsection{Convergence of the Galerkin method in $L^2(\Gamma)$ for the standard second-kind integral equations
 on polyhedral and Lipschitz domains.
}\label{sec:openbook}

The disadvantage of using second-kind BIEs in ${L^2(\Gamma)}$ is that convergence of the Galerkin method is harder to establish when
$\Gamma$ is only Lipschitz, or Lipschitz polyhedral. Indeed, in these cases $D$ and $D'$ are not compact; e.g., when $\Gamma$ has a corner or edge their spectra are not discrete; see, e.g. \cite[\S8.1.3]{At:97}.
When $\Gamma$ is only Lipschitz, $D$ and $D'$  are bounded on $L^2(\Gamma)$ by  the results on boundedness of the Cauchy integral on Lipschitz $\Gamma$ of Coifman, McIntosh, and Meyer \cite{CoMcMe:82} (following earlier work by Calder\'on \cite{Ca:77} on boundedness for $\Gamma$ with small Lipschitz character).
Verchota \cite{Ve:84}
showed that the operators
$\frac{1}{2} I \pm D$ and $\frac{1}{2} I \pm D'$
are Fredholm of index zero on $L^2(\Gamma)$; when $\Gamma$ is connected, $\frac{1}{2} I - D$ and $\frac{1}{2} I- D'$ are invertible on ${L^2(\Gamma)}$ and $\frac{1}{2} I + D$ and $\frac{1}{2} I+ D'$ invertible on $L^2_0(\Gamma)$, the set of $\phi\in{L^2(\Gamma)}$ with mean value zero; see \cite[Theorems 3.1 and 3.3(i)]{Ve:84}.
\footnote{The invertibility of $\frac{1}{2} I-D'$ on ${L^2(\Gamma)}$
implies that the bilinear form of the associated least-squares formulation
\begin{equation*}
a(\phi,\psi)= \left(\left(\frac{1}{2} I - D'\right)\phi, \left(\frac{1}{2} I - D'\right)\psi\right)_{L^2(\Gamma)}
\end{equation*}
is coercive. This formulation, however, suffers from the same disadvantages as the variational formulation of $\frac{1}{2} I- D'$ in ${H^{-1/2}(\Gamma)}$, including that computing the entries of the Galerkin matrix requires computing integrals over all of $\Gamma$, even when the basis functions have support on (small) subsets of $\Gamma$.
}

A long-standing open question has been
\begin{quotation}
\noindent Can $\frac{1}{2} I \pm D$ and $\frac{1}{2} I \pm  D'$ be written as the sum of a coercive operator and a compact operator in the space $L^2(\Gamma)$ when $\Gamma$ is only assumed to be Lipschitz?
\end{quotation}
By Part (b) of Theorem \ref{thm:Galerkin}, this question is equivalent to the question: does the Galerkin method
applied to  $\frac{1}{2} I \pm D$ and $\frac{1}{2} I \pm  D'$ in ${L^2(\Gamma)}$ converge
for every asymptotically-dense sequence of subspaces
 when $\Gamma$ is only assumed to be Lipschitz?

Until recently, this question was answered only in the following two cases, both in the affirmative: (i) $\Gamma$ is a 2d curvilinear polygon with each side $C^{1,\alpha}$ for some $0<\alpha<1$ and with each corner angle in the range $(0,2\pi)$. (ii) $\Gamma$ is Lipschitz, with sufficiently small Lipschitz character. Regarding (i): this result was announced by Shelepov in \cite{Sh:69}, with details of the proof given in \cite{Sh:91}, and with the analogous result for polygons following from the result of Chandler \cite[\S3]{Ch:84}; see, e.g. \cite[Lemma 1.5]{BoMaNiPa:90}.
Regarding (ii): Wendland \cite[\S4.2]{We:09} recognised that the results of I.~Mitrea \cite[Lemma 1, Page 392]{Mi:99} about the essential spectral radius could be adapted to prove this result, with this result proved in full in \cite[Corollary 3.5]{ChSp:22}; for more discussion on both (i) and (ii), see \cite[\S1]{ChSp:22}.

The recent paper \cite{ChSp:22} finally settled the question above negatively by giving examples of 2-d Lipschitz domains and 3-d star-shaped Lipschitz polyhedra for which $\frac{1}{2} I \pm D$ and $\frac{1}{2} I \pm  D'$ \emph{cannot} be written as the sum of a coercive operator and a compact operator in the space $L^2(\Gamma)$.
The 3-d star-shaped Lipschitz polyhedra are defined in \cite[Definition 5.7]{ChSp:22}, and called the \emph{open-book polyhedra}; see Figure \ref{fig:book} for an example, where we use the notation that $\Omega_{\theta,n}$ is the open-book polyhedron with $n$ pages and opening angle $\theta$.
Given $\epsilon>0$ there exists $\theta_0\in (0,\pi]$ such that the essential numerical range of $D$ in ${L^2(\Gamma)}$ contains the interval $[-\sqrt{n}/2 + \epsilon, \sqrt{n}/2-\epsilon]$ \cite[Theorem 1.3]{ChSp:22}. By the definition of the essential numerical range (see, e.g., \cite[Equation 2.3]{ChSp:22}), this result implies that if $\theta$ is sufficiently small and $n\geq 2$, then $\frac{1}{2} I \pm D$ and $\frac{1}{2} I \pm  D'$ \emph{cannot} be written as the sum of a coercive operator and a compact operator in the space $L^2(\Gamma)$ when $\Gamma = \partial \Omega_{\theta,n}$.

\begin{figure}
\includegraphics[width=.5\textwidth]{NewBook1}\includegraphics[width=.5\textwidth]{NewBook2}
\caption{\label{fig:book} Views from above and below of the open-book polyhedron $\Omega_{\theta,n}$ of \cite[Definition 5.7]{ChSp:22}, with $n=4$ pages and opening angle $\theta = \pi/2$}
\end{figure}

Nevertheless, Part (b) of Theorem \ref{thm:Galerkin} only shows that the Galerkin method applied to these domains does not converge for every asymptotically dense sequence $({\cal H}_N)_{N=1}^\infty \subset {L^2(\Gamma)}$, leaving opening the possibility that all Galerkin methods used in practice (based on boundary element method discretisation \cite{St:08, SaSc:11}) are in fact convergent. However, the following result from \cite{ChSp:22} clarifies that this is not the case.

\begin{theorem}\textbf{\emph{(\cite[Theorem 1.4]{ChSp:22})}}\label{thm:GalGenNew} Suppose that $A$ is invertible but $A$ cannot be written in the form $A=A_0+K$, where $A_0$ is coercive and $K$ is compact, and that $({\cal H}^*_N)_{N=1}^\infty$ is a sequence of finite-dimensional subspaces of ${\cal H}$, with ${\cal H}^*_1\subset {\cal H}^*_2 \subset ...$,
for which the Galerkin method converges.
Then there exists a sequence $({\cal H}_N)_{N=1}^\infty$ of finite-dimensional subspaces of ${\cal H}$, with ${\cal H}_1\subset {\cal H}_2 \subset ...$, such that:
\begin{enumerate}
\item[(a)] the Galerkin method does not converge
for the sequence $({\cal H}_N)_{N=1}^\infty$; and
\item[(b)] for each $N\in \mathbb{N}$,
\begin{equation}\label{eq:embed}
{\cal H}^*_N \subset {\cal H}_N \subset {\cal H}^*_{M_N},\quad\text{ for some } M_N\in \mathbb{N}.
\end{equation}
\end{enumerate}
\end{theorem}

We can apply this result when $({\cal H}^*_N)_{N=1}^\infty$ is a sequence of boundary element subspaces that is asymptotically dense in ${L^2(\Gamma)}$, in which case $({\cal H}_N)_{N=1}^\infty$, satisfying \eqref{eq:embed}, is also a sequence of boundary element subspaces (since ${\cal H}_N\subset {\cal H}_{M_N}^*$) and is also asymptotically dense in ${L^2(\Gamma)}$ (since ${\cal H}^*_N\subset {\cal H}_N$).

In summary, the results of \cite{ChSp:22} show that there exist Lipschitz and polyhedral boundaries $\Gamma$ for which there are Galerkin methods for solving BIEs involving $\frac{1}{2} I \pm D$ and $\frac{1}{2} I \pm  D'$
that do \emph{not} converge, with these methods based on asymptotically-dense sequences $({\cal H}_N)_{N=1}^\infty \subset{L^2(\Gamma)}$ of boundary element subspaces.

\subsection{Motivation for the present paper and summary of the main results}

Given the negative results of  \cite{ChSp:22} about convergence of the Galerkin method for the standard second-kind formulations, a natural question is therefore
\begin{quotation}
\noindent Do there exist second-kind BIE formulations in $L^2(\Gamma)$ of Laplace's equation where, with $\Gamma$ only assumed to be Lipschitz, the operators are continuous, invertible, and can be written as the sum of a coercive operator and a compact operator?
\end{quotation}
In this paper we answer this question in the affirmative for the Laplace interior and exterior Dirichlet problems. We present new BIE formulations that are continuous and in fact \emph{coercive} (i.e., not just the sum of a coercive and a compact operator) in the space $L^2(\Gamma)$, with $\Gamma$ only assumed to be Lipschitz; thus
convergence of the Galerkin method in ${L^2(\Gamma)}$
for every asymptotically-dense sequence $({\cal H}_N)_{N=1}^\infty$, plus the explicit error estimate \eqref{eq:quasioptimal},
is ensured by Part (c) of Theorem \ref{thm:Galerkin}. Furthermore, the strong property of coercivity allows us to prove that, if the Galerkin matrices are preconditioned by a specified diagonal matrix, then
the number of GMRES iterations required to solve the associated linear systems to a prescribed accuracy does \emph{not} increase as the discretisation is refined and $N$ increases.

\paragraph{Outline of the paper.} \S\ref{sec:notation} defines more precisely the Laplace BVPs we consider. \S\ref{sec:recap} recaps results about a non-standard layer potential introduced in \cite{Ca:85} and its non-tangential limits. \S\ref{sec:main_results} states the main results. \S\ref{sec:idea} discusses the ideas behind the main results, and the links to other work in the literature. \S\ref{sec:proofs} proves the main results, except the parts of the proofs that are related to the wellposedness and regularity of the Laplace oblique Robin problem, with these given in \S\ref{sec:oblique}. \S\ref{sec:Helmholtz} presents results for the Helmholtz exterior Dirichlet problem (with these results corollaries of the Laplace results in \S\ref{sec:main_results}).

\subsection{Notation and statement of the BVPs}\label{sec:notation}

Let $\Omega_-\subset \mathbb{R}^d$, $d\geq 2$, be a bounded (not necessarily connected) Lipschitz open set, and let $\Oe:= \mathbb{R}^d\setminus \overline{\Omega_-}$ and $\Gamma :=\partial \Omega_-$.
Let $\mathbf{n}$ be the outward-pointing unit normal vector to $\Omega_-$ (so $\mathbf{n}$ points out of $\Omega_-$ and into $\Oe$).
For $v \in H^1(\Omega_-)$ let $\gamma^- v$ denote its Dirichlet trace. For $v \in H^1(\Omega_-,\Delta):= \{ w : w \in H^1(\Omega_-), \Delta w \in L^2(\Omega_-)\}$ let $\partial_n^- v$ denote its Neumann trace; recall that, if $v \in H^2(\Omega_-)$ then $\partial_n^- v = \mathbf{n} \cdot \gamma^- \nabla v$.
Similarly, for $v\in H^1_{{\rm loc}}(\Oe):=\{ w: \Oe\rightarrow \mathbb{R} : \chi w \in H^1(\Oe) \text{ for all } \chi \in C^\infty_{\rm comp}(\mathbb{R}^d)\}$,
let $\gamma^+ v$ denote its Dirichlet trace. For
$v\in H^1_{{\rm loc}}(\Oe, \Delta):=\{ w: \Oe\to \mathbb{R} : \chi w \in H^1(\Oe), \chi \Delta w \in L^2(\Oe) \text{ for all } \chi \in C^\infty_{\rm comp}(\mathbb{R}^d)\}$, let $\partial_n^+ v$ denote its Neumann trace.

\begin{definition}[Laplace interior Dirichlet problem (IDP)]\label{def:idp}
Given $g_D\in H^{1/2}(\Gamma)$, we say that $u\in H^1(\Omega_-)$ satisfies the {\em interior Dirichlet problem (IDP)} if $\Delta u =0$ in $\Omega_-$ and $\gamma^- u = g_D$ on $\Gamma$.
\end{definition}

\begin{definition}[Laplace exterior Dirichlet problem (EDP)]\label{def:edp}
With $\Omega_-$ and $\Oe$ as above, assume further that $\Oe$ is connected. Given $g_D\in H^{1/2}(\Gamma)$, we say that $u\in H^1_{\rm{loc}}(\Oe)$ satisfies the {\em exterior Dirichlet problem (EDP)} if $\Delta u =0$ in $\Oe$, $\gamma^+ u = g_D$ on $\Gamma$, and $u(\mathbf{x})= O(1)$ when $d=2$ and $u(\mathbf{x}) = o(|\mathbf{x}|^{3-d})$ when $d\geq 3$ as $|\mathbf{x}| \rightarrow \infty$ (uniformly in all directions $\mathbf{x}/|\mathbf{x}|$).
\end{definition}

\noindent We make three remarks.

(i) Recall that, by elliptic regularity (see, e.g., \cite[Theorem 4.16]{Mc:00}),
the solution of the IDP and EDP are $C^\infty$ in $\Omega_-$ and $\Oe$ respectively. Therefore, the pointwise conditions at infinity imposed in the EDP
make sense.

(ii) For the IDP and EDP, uniqueness of the solution is shown in, e.g., \cite[Corollary 8.3]{Mc:00} and \cite[Theorems 8.9 and 8.10]{Mc:00} respectively\footnote{
\cite[Theorem 8.10]{Mc:00} proves uniqueness under the condition $u(\mathbf{x}) = O(|\mathbf{x}|^{2-d})$ as $|\mathbf{x}| \rightarrow \infty$, but when $d\geq 3$ the proof still goes through when $u(\mathbf{x}) = o(|\mathbf{x}|^{3-d})$ as $|\mathbf{x}| \rightarrow \infty$.
}.
Existence then follows from Fredholm theory and, e.g., \cite[Theorems 7.5, 7.6, and 7.15]{Mc:00}.

(iii) The Neumann traces of the solutions of both the IDP and EDP are in ${H^{-1/2}(\Gamma)}$; see, e.g., \cite[Lemma 4.3]{Mc:00}. Later, we consider both these BVPs when the Dirichlet data is in $H^1(\Gamma)$. The regularity result of Ne\v{c}as \cite[\S5.1.2]{Ne:67} (restated as Theorem \ref{thm:Necas} below) then implies that $\partial_n^- u$ and $\partial_n^+ u$ (in Definitions \ref{def:idp} and \ref{def:edp} respectively) are both in ${L^2(\Gamma)}$, as opposed to just in ${H^{-1/2}(\Gamma)}$.

\

The IDP and EDP can equivalently be formulated in terms of non-tangential limits, with these alternative formulations standard in the harmonic-analysis literature (see, e.g., \cite[Corollary 3.2]{Ve:84}, \cite[\S3]{Ca:85}, \cite[Theorem 2]{MeCo:00}, \cite[Proposition 5.1]{Ta:00}). We state these alternative formulations, and recall their equivalence, so that we can easily use results from the harmonic-analysis literature (summarised in Appendix \ref{app:HA} below).

Given $\mathbf{x}\in \Gamma$, let $\Theta^{\pm}(\mathbf{x})$ be the \emph{non-tangential approach set to $\mathbf{x}$ from $\Omega^\pm$} defined by
\begin{equation}\label{eq:ntapproach2}
\Theta^{\pm}(\mathbf{x}):= \Big\{ \mathbf{y} \in \Omega^{\pm} \, : \, |\mathbf{x}-\mathbf{y}|\leq \min\big\{c, C \mathrm{dist}(\mathbf{y},\Gamma)\big\}\Big\},
\end{equation}
for some $c>0$ and some $C>1$ sufficiently large depending on the Lipschitz character of $\Omega^\pm$.\footnote{E.g., $C> M+1$ is large enough when $M$ is the Lipschitz character.}
If $u\in C(\Omega^{\pm})$, its non-tangential maximal function $u^*: \Gamma \rightarrow [0,\infty]$ is defined by
\begin{equation}\label{eq:ntmax}
u^*(\mathbf{x}) := \sup_{\mathbf{y} \in \Theta^{\pm}(\mathbf{x})}|u(\mathbf{y})|, \quad \mathbf{x}\in\Gamma.
\end{equation}
Define the non-tangential limit
\begin{equation}\label{eq:ntlim}
\widetilde{\gamma}^\pm u(\mathbf{x}) := \lim_{\mathbf{y}\rightarrow \mathbf{x},\, \mathbf{y} \in \Theta^\pm(\mathbf{x})}u(\mathbf{y}).
\end{equation}
If $u\in C^2(\Omega^\pm)$, $\Delta u=0$, and $u^* \in {L^2(\Gamma)}$, then $\widetilde{\gamma}^\pm u(x)$ is well-defined for almost all $x\in \Gamma$ and $\widetilde{\gamma}^\pm u \in {L^2(\Gamma)}$ by \cite[Corollary 5.5]{JeKe:95} (restated as Part (ii) of Theorem \ref{thm:JeKe:95} below). Furthermore, if $u\in H^s_{\rm loc}(\Omega^\pm)$ with $s>1/2$, then $\widetilde{\gamma}^\pm u =\gamma^\pm u$ by \cite[Lemma A.9]{ChGrLaSp:12} (restated as Lemma \ref{lem:A9} below).

\begin{definition}[Laplace IDP formulated via non-tangential limits]\label{def:idp2}
Given $g_D\in L^2(\Gamma)$, we say that $u\in C^2(\Omega_-)$ with $u^* \in L^2(\Gamma)$ satisfies the IDP if $\Delta u =0$ in $\Omega_-$ and $\widetilde{\gamma}^- u = g_D$ on $\Gamma$.
\end{definition}

\begin{definition}[Laplace EDP formulated via non-tangential limits]\label{def:edp2}
With $\Omega_-$ and $\Oe$ as above, assume further that $\Oe$ is connected.
Given $g_D\in L^2(\Gamma)$, we say that $u\in C^2(\Oe)$ with $u^* \in L^2(\Gamma)$ satisfies the EDP if $\Delta u =0$ in $\Oe$, $\widetilde{\gamma}^+ u = g_D$ on $\Gamma$, and $u(\mathbf{x})= O(1)$ when $d=2$ and $u(\mathbf{x}) = o(|\mathbf{x}|^{3-d})$ when $d\geq 3$ as $|\mathbf{x}| \rightarrow \infty$ (uniformly in all directions $\mathbf{x}/|\mathbf{x}|$).
\end{definition}

Existence and uniqueness of the solutions of these formulations of the IDP and EDP go back to the work of Dahlberg \cite{Da:79}, and are given explicitly in, e.g., \cite[Corollary 3.2 and Lemma 3.7]{Ve:84}, \cite[\S3]{Ca:85}.
The following equivalence result is proved in Appendix \ref{app:BVPequiv}.

\begin{theorem}[Equivalence of the different formulations of the IDP and EDP]\label{thm:BVPequiv}
If $g_D\in H^{1/2}(\Gamma)$, then the solution of the IDP in the sense of Definition \ref{def:idp} is the solution of the IDP in the sense of Definition \ref{def:idp2}, and vice versa.

Similarly, if $g_D\in H^{1/2}(\Gamma)$, then the solution of the EDP in the sense of Definition \ref{def:edp} is the solution of the EDP in the sense of Definition \ref{def:edp2}, and vice versa.
\end{theorem}

\subsection{Recap of results about layer potentials and their non-tangential limits}\label{sec:recap}

Recall that the \emph{surface gradient} operator on $\Gamma$ is the unique operator such that, when $v \in C^1(\overline{\Omega_-})$, $\nabla v = \mathbf{n} (\mathbf{n}\cdot \nabla v) + \nabla_{\Gamma} (\gamma^- v)$ on $\Gamma$ (and similarly for $v \in C^1(\overline{\Oe})$); for an explicit expression for $\nabla_{\Gamma}$ in terms of a parametrisation of $\Gamma$, see, e.g., \cite[Equation A.14]{ChGrLaSp:12}.

The following results all rely on the harmonic-analysis results in \cite{CoMcMe:82} and \cite{Ve:84} (see also the accounts in \cite[Chapter 15]{MeCo:00}, \cite[Chapter 4]{Ta:00}, \cite[Chapter 2, Section 2]{Ke:94}).
Define
\begin{equation*}
\nabla_{\Gamma} S \phi(\mathbf{x}):= \int_\Gamma \nabla_{\Gamma,\mathbf{x}}\Phi(\mathbf{x},\mathbf{y})\phi(\mathbf{y}){\rm d} s(\mathbf{y}) \quad\text{ for a.e. } \mathbf{x}\in \Gamma,
\end{equation*}
where the integral is understood in the principal-value sense. By \cite[Theorem 1.6]{Ve:84}, $\nabla_{\Gamma} S:{L^2(\Gamma)}\rightarrow{L^2(\Gamma)}$, with this mapping continuous, and $(\nabla_{\Gamma} S)\phi= \nabla_{\Gamma} (S\phi)$.
The following potential was introduced in \cite[\S2]{Ca:85}; given  $\bZ\in (L^\infty(\Gamma))^d$ that is real-valued (which we assume throughout), let
\begin{equation}\label{eq:Caldp}
\mathcal{K}_{\vfb} \phi(\mathbf{x}):= \int_\Gamma \bZ(\mathbf{y}) \cdot \nabla_\mathbf{y} \Phi(\mathbf{x},\mathbf{y}) \phi(\mathbf{y})\, {\rm d} s(\mathbf{y}) \quad
\quad\text{ for }\mathbf{x} \in \mathbb{R}^d\setminus \Gamma,
\end{equation}
and let
\begin{equation}\label{eq:Caldb}
K_{\vfb} \phi(\mathbf{x}):= \int_\Gamma  \bZ(\mathbf{y})\cdot\nabla_\mathbf{y}\Phi(\mathbf{x},\mathbf{y})\phi(\mathbf{y}) \, {\rm d} s(\mathbf{y}) \quad\quad\text{ for a.e. }\mathbf{x} \in \Gamma,
\end{equation}
where the integral in \eqref{eq:Caldb} is understood in the principal-value sense.
The results of \cite{CoMcMe:82} and \cite{Ve:84} imply that $K_{\vfb}:{L^2(\Gamma)}\rightarrow {L^2(\Gamma)}$ and, for $\phi\in{L^2(\Gamma)}$,
$\mathcal{K}_{\vfb} \phi \in C^2(\Omega^\pm)$, $\mathcal{K}_{\vfb} \phi$ satisfies Laplace's equation, and
$(\mathcal{K}_{\vfb}\phi)^*\in {L^2(\Gamma)}$ with
\begin{equation}\label{eq:jumpCald}
\widetilde{\gamma}^{\pm}
\mathcal{K}_{\vfb} \phi(\mathbf{x}) = \pm \frac{1}{2} \big( \bZ(\mathbf{x})\cdot \mathbf{n}(\mathbf{x})\big) \phi(\mathbf{x}) + K_{\vfb} \phi (\mathbf{x}) \quad\quad\text{ for a.e. } \mathbf{x}\in \Gamma.
\end{equation}

Observe that (i) when $\bZ= \mathbf{n}$, $\mathcal{K}_{\vfb} = \mathcal{D}$, $K_{\vfb}= D$, and \eqref{eq:jumpCald} is the usual jump relation for the double-layer potential, and (ii) we can rewrite $K_{\vfb}$ as
\begin{equation}\label{eq:Caldb2}
K_{\vfb} \phi(\mathbf{x}) = \int_\Gamma \left( \bZ(\mathbf{y})\cdot \mathbf{n}(\mathbf{y}) \pdiff{\Phi(\mathbf{x},\mathbf{y})}{n(\mathbf{y})} \phi(\mathbf{y}) + \bZ(\mathbf{y}) \cdot \nabla_{\Gamma} \Phi(\mathbf{x},\mathbf{y}) \phi(\mathbf{y})\right) {\rm d} s(\mathbf{y}).
\end{equation}
In a similar way to how the $L^2$ adjoint of $D$ is $D'$ (see, e.g., \cite[Chapter 15, text around Equation 4.10]{MeCo:00}), the $L^2$ adjoint of $K_{\vfb}$ is
\begin{equation}\label{eq:Caldb'}
K_{\vfb}' \phi(\mathbf{x}) := \big(\bZ(\mathbf{x}) \cdot \mathbf{n}(\mathbf{x}) \big)D' \phi(\mathbf{x}) + \bZ(\mathbf{x}) \cdot \nabla_{\Gamma} S \phi(\mathbf{x}).
\end{equation}
The significance of the operator $K_{\vfb}'$ is that it appears in the inner product of $\bZ$ with the non-tangential limit of $\nabla\mathcal{S}$, where ${\cal S}$ is the single-layer potential defined for $\phi\in {L^2(\Gamma)}$ by
\begin{equation}\label{eq:SLP}
{\cal S}\phi(\mathbf{x}):= \int_\Gamma \Phi(\mathbf{x},\mathbf{y})\phi(\mathbf{y})\, {\rm d} s(\mathbf{y}) \quad\quad\text{ for } \mathbf{x} \in\mathbb{R}^d\setminus \Gamma.
\end{equation}
Indeed, by \cite[Theorems 1.6 and 1.11]{Ve:84} (see also \cite[Theorem 5]{MeCo:00}, \cite[Equation 2.30]{ChGrLaSp:12}), for almost every $\mathbf{x}\in \Gamma$,
\begin{equation}\label{eq:gradS}
\widetilde{\gamma}^\pm
\nabla {\cal S}\phi(\mathbf{x}) = \mathbf{n}(\mathbf{x}) \left( \mp \frac{1}{2} I + D'\right) \phi(\mathbf{x}) + \nabla_{\Gamma}(S\phi)(\mathbf{x}),
\end{equation}
so that
\begin{equation}\label{eq:gradS2}
 \bZ(\mathbf{x})\cdot\widetilde{\gamma}^\pm\nabla {\cal S}\phi(\mathbf{x}) = \left( \mp\frac{1}{2}(\mathbf{Z}(\mathbf{x})\cdot\mathbf{n}(\mathbf{x}))I + K_{\vfb}'\right)\phi(\mathbf{x}).
\end{equation}

\section{Statement of the main results}\label{sec:main_results}

\subsection{New boundary integral equations for the Laplace interior and exterior Dirichlet problems on general Lipschitz domains for $d\geq 3$}\label{sec:dgeq3}

We focus on the case $d\geq 3$, since the question of whether or not there exist BIE formulations of the Laplace IDP and EDP that are coercive, or coercive up to a compact perturbation, on Lipschitz domains is more pressing when $d=3$ than $d=2$ (because of
the existing convergence theory for $\pm\frac{1}{2} I + D$ and $\pm\frac{1}{2} I+ D'$ on curvilinear polygons \cite{Sh:69, Ch:84, Sh:91}
 but negative results for these operators for certain 3-d star-shaped polyhedra \cite{ChSp:22} recapped in \S\ref{sec:openbook}).
Results for $d=2$ are given in \S\ref{sec:d=2}.

\subsubsection{The interior Dirichlet problem}

Given $\bZ\in (L^\infty(\Gamma))^d$ and $\alpha \in \mathbb{R}$,
define the integral operators $A'_{I,\bZ,\alpha}$,
$A_{I,\bZ,\alpha}$,
and $B_{I,\bZ,\alpha}$ by
\begin{align}
\label{eq:opLIZ}
&A'_{I,\bZ,\alpha}:= \frac{1}{2} (\bZ\cdot \mathbf{n})I - K_{\vfb}' + \alpha S, \qquad
A_{I,\bZ,\alpha}:= \frac{1}{2} (\bZ\cdot \mathbf{n})I - K_{\vfb} + \alpha S, \\
&B_{I,\bZ,\alpha}:=-(\bZ\cdot \mathbf{n}) H  - \bZ\cdot\nabla_{\Gamma} \left( \frac{1}{2} I + D\right) + \alpha\left(\frac{1}{2} I + D\right),
\end{align}
with the subscript $I$ standing for ``interior", and the $'$ superscript indicating that $A'_{I,\bZ,\alpha}$ is the $L^2$ adjoint of $A_{I,\bZ,\alpha}$.

\begin{theorem}[New integral equations for Laplace IDP with $d\geq 3$]\label{thm:Laplace_int}

\

(i) \textbf{\emph{Direct formulation.}} Let $u$ be the solution of the Laplace IDP of Definition \ref{def:idp} with $d\geq 3$ and additionally $g_D\in H^1(\Gamma)$.
Then $\partial_n^- u$ satisfies
\begin{equation}\label{eq:new_int}
A'_{I,\bZ,\alpha} \partial_n^- u =B_{I,\bZ,\alpha} g_D.
\end{equation}

(ii) \textbf{\emph{Indirect formulation.}} Given $g_D\in {L^2(\Gamma)}$, if $\phi\in{L^2(\Gamma)}$ satisfies
\begin{equation}\label{eq:new_int_indirect}
A_{I,\bZ,\alpha} \phi =- g_D,
\end{equation}
then
\begin{equation}\label{eq:indirectansatz1}
u:= (\mathcal{K}_{\vfb} - \alpha {\cal S})\phi
\end{equation}
is the solution of the Laplace IDP of Definition \ref{def:idp2}.

(iii) \textbf{\emph{Continuity.}}
$A'_{I,\bZ,\alpha}: {L^2(\Gamma)}\rightarrow{L^2(\Gamma)}$,
$A_{I,\bZ,\alpha}: {L^2(\Gamma)}\rightarrow{L^2(\Gamma)}$, and $B_{I,\bZ,\alpha}:H^1(\Gamma)\rightarrow{L^2(\Gamma)}$, and these mappings are continuous.

 (iv) \textbf{\emph{Coercivity up to compact perturbation.}} If $\mathbf{Z}\in (C(\Gamma))^d$ and there exists $c>0$ such that
\begin{equation}\label{eq:c}
\bZ(\mathbf{x}) \cdot \mathbf{n}(\mathbf{x}) \geq c \quad \text{ for almost every } \mathbf{x}\in \Gamma,
\end{equation}
 then both $A'_{I,\bZ,\alpha}$ and $A_{I,\bZ,\alpha}$ are the sum of a coercive operator and a compact operator on ${L^2(\Gamma)}$.

 (v) \textbf{\emph{Invertibility for all $\alpha>0$.}} If $\alpha>0$, $\mathbf{Z}\in (C^{0,\beta}(\Gamma))^d$ for some $0<\beta<1$, and there exists $c>0$ such that
\eqref{eq:c} holds,
then both $A'_{I,\bZ,\alpha}:{L^2(\Gamma)}\to{L^2(\Gamma)}$ and $A_{I,\bZ,\alpha}:{L^2(\Gamma)}\to{L^2(\Gamma)}$ are invertible.

 (vi) \textbf{\emph{Coercivity for sufficiently large $\alpha$.}} If $\mathbf{Z}$ satisfies \eqref{eq:c}, $\mathbf{Z} \in (C^{0,1}(\Gamma))^d$ with Lipschitz constant $L_{\bZ}$, and
 \begin{equation}\label{eq:alpha_bound_boundary}
 2\alpha \geq 3d L_{\bZ},
 \end{equation}
  then both
$A'_{I,\bZ,\alpha}$ and $A_{I,\bZ,\alpha}$ are coercive on ${L^2(\Gamma)}$ with coercivity constant $c/2$ (with $c$ defined by \eqref{eq:c});
indeed,
\begin{equation}\label{eq:BIEcoercive}
\big(A'_{I,\bZ,\alpha} \psi,\psi\big)_{{L^2(\Gamma)}} \geq \frac{c}{2} \N{\psi}^2_{{L^2(\Gamma)}} \quad\text{ for all real-valued } \psi \in {L^2(\Gamma)},
\end{equation}
and similarly for $A_{I,\bZ,\alpha}$.
 \end{theorem}

Recall that if $A$ is real and $(A\psi,\psi)\geq C_{\rm coer} \|\psi\|^2_{{L^2(\Gamma)}}$ for all real-valued $\psi\in {L^2(\Gamma)}$, then $\Re(A\phi,\phi)\geq C_{\rm coer} \|\phi\|^2_{{L^2(\Gamma)}}$ for all complex-valued $\phi\in{L^2(\Gamma)}$; thus \eqref{eq:BIEcoercive} implies that $A'_{I,\bZ,\alpha}$ and $A_{I,\bZ,\alpha}$ are coercive on complex-valued ${L^2(\Gamma)}$.

For any bounded Lipschitz  open set $\Omega_-$ there exists $\mathbf{Z} \in (C^{0,1}(\Gamma))^d$ such that \eqref{eq:c} holds; for completeness, and because a formula for $\mathbf{Z}$ is needed for implementation, we include this result and a concrete, constructive proof in \S\ref{sec:Z} below.
The combination of this result and Parts (iii) and (vi) of Theorem \ref{thm:Laplace_int} imply that, for any bounded Lipschitz open set, there exists a BIE formulation of the Laplace IDP that is continuous and coercive in ${L^2(\Gamma)}$.

 \subsubsection{The exterior Dirichlet problem}

Given  $\bZ\in (L^\infty(\Gamma))^d$ and $\alpha \in \mathbb{R}$,
define the integral operators $A'_{E,\bZ,\alpha}$,
$A_{E,\bZ,\alpha}$,
and $B_{E,\bZ,\alpha}$ by
\begin{align}
\label{eq:opLEZ}
&A'_{E,\bZ,\alpha}:= \frac{1}{2} (\bZ\cdot \mathbf{n})I + K_{\vfb}' + \alpha S, \qquad
A_{E,\bZ,\alpha}:= \frac{1}{2} (\bZ\cdot \mathbf{n})I + K_{\vfb} + \alpha S, \\
&B_{E,\bZ,\alpha}:=(\bZ\cdot \mathbf{n}) H  + \bZ\cdot\nabla_{\Gamma} \left( -\frac{1}{2} I + D\right) + \alpha\left(-\frac{1}{2} I + D\right),
\end{align}
with the subscript $E$ standing for ``exterior".

\begin{theorem}[New integral equations for Laplace EDP with $d\geq 3$]\label{thm:Laplace_ext}

\

(i) \textbf{\emph{Direct formulation.}} Let $u$ be the solution of the Laplace EDP of Definition \ref{def:edp} with $d=3$ and additionally $g_D\in H^1(\Gamma)$. Then $\partial_n^+ u$ satisfies
\begin{equation}\label{eq:new_ext}
A'_{E,\bZ,\alpha}\partial_n^+ u =B_{E,\bZ,\alpha} g_D.
\end{equation}

(ii) \textbf{\emph{Indirect formulation.}} Given $g_D\in {L^2(\Gamma)}$, if $\phi\in {L^2(\Gamma)}$ satisifes
\begin{equation}\label{eq:new_ext_indirect}
A_{E,\bZ,\alpha} \phi = g_D,
\end{equation}
then
\begin{equation}\label{eq:indirect_ext_u}
u:= (\mathcal{K}_{\vfb} + \alpha {\cal S})\phi
\end{equation}
is the solution of the Laplace EDP of Definition \ref{def:edp2}.

  (iii) \textbf{\emph{Continuity.}}
 $A'_{E,\bZ,\alpha}: {L^2(\Gamma)}\rightarrow{L^2(\Gamma)}$,
$A_{E,\bZ,\alpha}: {L^2(\Gamma)}\rightarrow{L^2(\Gamma)}$,
and $B_{E,\bZ,\alpha}:H^1(\Gamma)\rightarrow{L^2(\Gamma)}$, and these mappings are continuous.

  (iv) \textbf{\emph{Coercivity up to compact perturbation.}} If $\mathbf{Z}\in (C(\Gamma))^d$ and there exists $c>0$ such that \eqref{eq:c} holds,
 then both $A'_{E,\bZ,\alpha}$ and $A_{E,\bZ,\alpha}$ are the sum of a coercive operator and a compact operator on ${L^2(\Gamma)}$.

 (v) \textbf{\emph{Invertibility for all $\alpha>0$}.} If $\alpha>0$, $\mathbf{Z}\in (C^{0,\beta}(\Gamma))^d$ for some $0<\beta<1$, and there exists $c>0$ such that
\eqref{eq:c} holds, then
 both $A'_{E,\bZ,\alpha}:{L^2(\Gamma)}\to{L^2(\Gamma)}$ and $A_{E,\bZ,\alpha}:{L^2(\Gamma)}\to{L^2(\Gamma)}$ are invertible.

 (vi) \textbf{\emph{Coercivity for sufficiently large $\alpha$.}} If $\mathbf{Z}\in (C^{0,1}(\Gamma))^d$ with Lipschitz constant $L_{\bZ}$ and \eqref{eq:alpha_bound_boundary} holds, then both
$A'_{E,\bZ,\alpha}$ and $A_{E,\bZ,\alpha}$ are coercive on ${L^2(\Gamma)}$ with coercivity constant $c/2$ (with $c$ defined by \eqref{eq:c}), in that \eqref{eq:BIEcoercive} holds with $A'_{I,\bZ,\alpha}$ replaced by either $A'_{E,\bZ,\alpha}$ or $A_{E,\bZ,\alpha}$.
 \end{theorem}

Similar to the case of the IDP, the result in \S\ref{sec:Z} and Parts (iii) and (vi) of Theorem \ref{thm:Laplace_ext} imply that, for any bounded Lipschitz open set $\Omega_-$ such that $\Oe$ is connected, there exists a BIE formulation of the Laplace EDP that is continuous and coercive in ${L^2(\Gamma)}$.

\subsubsection{The new formulations of the IDP and EDP for $d\geq 3$ on domains that are star-shaped with respect to a ball}\label{sec:Laplace_star}

When $\Omega_-$ is star-shaped with respect to a ball, the coercivity results in Theorems \ref{thm:Laplace_int} and \ref{thm:Laplace_ext} take a particularly simple form.

\begin{definition}
(i) $D$ is \emph{star-shaped with respect to the point $\mathbf{x}_0$} if, whenever $\mathbf{x} \in D$, the segment $[\mathbf{x}_0,\mathbf{x}]\subset D$.

\noindent (ii) $D$ is \emph{star-shaped with respect to the ball $B_{\kappa}(\mathbf{x}_0)$} if it is star-shaped with respect to every point in $B_{\kappa}(\mathbf{x}_0)$.

\end{definition}

\begin{lemma}\textbf{\emph{(\cite[Lemma 5.4.1]{Mo:11})}}\label{lem:star} If $D$ is Lipschitz with outward unit normal vector $\boldsymbol{\nu}$, then
$D$ is star-shaped with respect to $B_{\kappa}(\mathbf{x}_0)$, for some $\kappa>0$, if and only if $(\mathbf{x}-\mathbf{x}_0) \cdot \boldsymbol{\nu}(\mathbf{x}) \geq {\kappa}$ for all  $\mathbf{x} \in \partial D$ for which $\boldsymbol{\nu}(\mathbf{x})$ is defined.
\end{lemma}

\noindent From now on, if a domain $D$ is star-shaped with respect to $\mathbf{x}_0$, we assume (without loss of generality) that $\mathbf{x}_0=\mathbf{0}$.

\begin{theorem}[Coercivity for star-shaped domains]\label{thm:star}
Let $\Omega_- \subset \mathbb{R}^d$, $d\geq 3$, be a bounded Lipschitz domain that is star-shaped with respect to a ball of radius $\kappa$, i.e.
\begin{equation}\label{eq:ssg}
\kappa:= \mathop{{\rm ess} \inf}_{\mathbf{x}\in\Gamma}( \mathbf{x} \cdot \mathbf{n}(\mathbf{x})),
\end{equation}
Then
\begin{equation*}
A^I_{\bx} \,\,\text{ and }\,\,A_{I,\bx,\alpha}, \,\,\text{ with } \alpha\geq -(d-2)/2,
\end{equation*}
and
\begin{equation*}
\quadA^E_{\bx}\,\, \text{ and }\,\, A_{E,\bx,\alpha},\,\,\text{ with } \alpha\geq (d-2)/2,
\end{equation*}
are all coercive on ${L^2(\Gamma)}$ with coercivity constant $\kappa/2$, in that \eqref{eq:BIEcoercive} holds with $c$ replaced by $\kappa$, and $A'_{I,\bZ,\alpha}$ replaced by any one of $A^I_{\bx}, A_{I,\bx,\alpha}, A^E_{\bx}$, or $A_{E,\bx,\alpha}$.
\end{theorem}

\subsection{Convergence and conditioning of the associated Galerkin methods}\label{sec:Galerkin_imp}

We now show how Theorems \ref{thm:Laplace_int} and \ref{thm:Laplace_ext} imply that (i) the associated Galerkin methods converge (see \S\ref{sec:convergence}), and (ii)
the associated Galerkin matrices are provably well-conditioned
as the discretisation is refined, without the need for operator preconditioning (see \S\ref{sec:conditioning}).
We focus on the case $d\geq 3$ and the new BIE formulations for the IDP and EDP (appearing in Theorems \ref{thm:Laplace_int} and \ref{thm:Laplace_ext}), but analogous results hold for the BIEs for star-shaped domains in \S\ref{sec:Laplace_star}
and also for the new BIEs for $d=2$ in \S\ref{sec:d=2} below.

\subsubsection{Convergence of the Galerkin method for the new formulations}\label{sec:convergence}

\begin{corollary}[Convergence of the Galerkin method]\label{cor:Galerkin}
Let $({\cal H}_N)_{N=1}^\infty$ denote any sequence of finite-dimensional subsets of ${\cal H} := L^2(\Gamma)$ that is asymptotically dense in $L^2(\Gamma)$ in the sense defined in \S\ref{sec:Gal}.

(a) If $\mathbf{Z}\in (C^{0,\beta}(\Gamma))^d$ for some $0<\beta<1$, and there exists $c>0$ such that
\eqref{eq:c} holds then, for all $\alpha >0$, the Galerkin method \eqref{eq:G} applied to any one of the BIEs \eqref{eq:new_int}, \eqref{eq:new_int_indirect}, \eqref{eq:new_ext} or \eqref{eq:new_ext_indirect} converges (in the sense defined in \S\ref{sec:intro}).

(b) If, additionally, $\mathbf{Z}$ is Lipschitz and $\alpha$ satisfies \eqref{eq:alpha_bound_boundary}, then, additionally, the Galerkin solution exists for every finite-dimensional subspace ${\cal H}_N\subset{L^2(\Gamma)}$ and satisfies the quasioptimal error estimate \eqref{eq:quasioptimal}, with constant $2\|A^\dag_{\bZ,\alpha}\|_{{\LtG\rightarrow \LtG}}/c$, where $A^\dag_{\bZ,\alpha} :=A_{I,\bZ,\alpha}$ for the BIEs \eqref{eq:new_int} and \eqref{eq:new_int_indirect}, $A^\dag_{\bZ,\alpha} :=A_{E,\bZ,\alpha}$ for the BIEs \eqref{eq:new_ext} and \eqref{eq:new_ext_indirect}.
\end{corollary}

Since the proof is so short, we include it here.

\

\begin{proof}[Proof of Corollary \ref{cor:Galerkin}]
(a) This follows from Parts (iv) and (v) of Theorem \ref{thm:Laplace_int}/Theorem \ref{thm:Laplace_ext}
and Part (b) of Theorem \ref{thm:Galerkin}.
(b) This follows from Part (vi) of Theorem \ref{thm:Laplace_int}/Theorem \ref{thm:Laplace_ext} and Part (c) of Theorem \ref{thm:Galerkin}.
\end{proof}

\

We highlight that Corollary \ref{cor:Galerkin} is the first time convergence of the Galerkin method for a BIE posed in ${L^2(\Gamma)}$ used to solve a boundary-value problem for Laplace's equation has been proved with the only assumption on $\Gamma$ that it is Lipschitz; the same is true if $\Gamma$ is assumed to be Lipschitz polyhedral in 3-d.

\begin{remark}[Bounding the best approximation and Galerkin errors] \label{rem:BestApprox}
For 3-d Lipschitz polyhedra  the smoothness of the solution, in particular its singularities at corners and edges, is well understood (see, e.g., \cite{PeSt:90}) for the direct formulations \eqref{eq:new_int} and \eqref{eq:new_ext}, where the solution of the integral equation is $\phi=\partial_n^\pm u$. Moreover, it is well understood how to design effective $h$- and $hp$-boundary element approximation spaces ${\cal H}_N$
based on graded, anisotropic meshes so as to obtain optimal best approximation error estimates (see, e.g., \cite{PeSt:90,El:92a,El:95,LuNi:95,MaSt:97}), indeed exponential convergence of $\min_{\psi \in {\cal H}_N}\big\|\phi-\psi\big\|_{{L^2(\Gamma)}}$ as a function of $M_N:=\dim({\cal H}_N)$ if the Dirichlet data $g_D$ is the restriction to $\Gamma$ of an analytic function (see \cite[Theorem 3.1]{MaSt:97}). Further,  by Part (a) of Corollary \ref{cor:Galerkin} and the quasioptimality \eqref{eq:quasioptimal2}, the same rates of convergence follow for the Galerkin error $\|\phi-\phi_N\|_{L^2(\Gamma)}$ as long as $\alpha>0$.
\end{remark}

\subsubsection{Solution of the Galerkin linear systems of the new formulations}\label{sec:conditioning}

Let ${\cal H}_N = {\rm span} \{\psi^{N}_1,\ldots, \psi^{N}_{M_N}\}$, with $M_N=\dim({\cal H}_N)$ and $\{\psi^{N}_1,\ldots, \psi^{N}_{M_N}\}$ a basis for ${\cal H}_N$.
The Galerkin method \eqref{eq:G} applied to \eqref{eq:new_int_indirect}
is then equivalent to the linear system
\begin{equation}\label{eq:system0}
{\mymatrix{A}} \mathbf{x} = \mathbf{b}
\end{equation}
where
\begin{equation}\label{eq:matrix}
({\mymatrix{A}})_{ij} := \big(A^\dag_{\bZ,\alpha} \psi^{N}_j,\psi^{N}_i\big)_{L^2(\Gamma)} \quad\text{ and }\quad b_i := -(g_D,\psi^{N}_i)_{L^2(\Gamma)}, \quad i,j=1,\ldots, M_N,
\end{equation}
with $A^\dag_{\bZ,\alpha}:= A_{I,\bZ,\alpha}$, and with the Galerkin solution $\phi_N$ given by $\phi_N= \sum_{j=1}^{M_N} x^{N}_j \psi^{N}_j$, where $\mathbf{x} = (x_1^N,\ldots,x_{M_N}^N)^T$. The Galerkin method applied to \eqref{eq:new_int}, \eqref{eq:new_ext_indirect}, or \eqref{eq:new_ext}, respectively, is also equivalent to \eqref{eq:system0} with $A^\dag_{\bZ,\alpha}:=A'_{I,\bZ,\alpha}$, $A_{E,\bZ,\alpha}$, or $A'_{E,\bZ,\alpha}$  in \eqref{eq:matrix} and with correspondingly different definitions of the right-hand side components $b_i$.

In each case, whether $A^\dag_{\bZ,\alpha}=A_{I,\bZ,\alpha}$, $A'_{I,\bZ,\alpha}$, $A_{E,\bZ,\alpha}$, or $A'_{E,\bZ,\alpha}$, the matrix ${\mymatrix{A}}$ defined in \eqref{eq:matrix} is non-symmetric, and a popular method for solving such non-symmetric linear systems is
the \emph{generalised minimum residual method (GMRES)} \cite{SaSc:86}, which we now briefly recall.
Consider the abstract  linear system
${\mymatrix{C}} \mathbf{x} = \mathbf{d}$
in $\mathbb{C}^{M_N}$, where ${\mymatrix{C}} \in \mathbb{C}^{M_N\times M_N}$ is invertible. Let $\mathbf{x}_{0}$ be an  initial guess for $\mathbf{x}$, and define the corresponding initial residual $\mathbf{r}^0 :={\mymatrix{C}} \mathbf{x}^0-\mathbf{d}$ and the corresponding standard Krylov spaces by
\begin{equation*}
{\cal K}^m({\mymatrix{C}}, \mathbf{r}_0) := \mathrm{span}\big\{{\mymatrix{C}}^j \mathbf{r}_0 : j = 0, \ldots, m-1\big\}.
\end{equation*}
For $m \geq 1$, define the $m$th GMRES iterate $\mathbf{x}_m$  to be the unique element of ${\cal K}^m$ such that its residual $\mathbf{r}_m:= {\mymatrix{C}} \mathbf{x}_m-\mathbf{d}$ satisfies  the   minimal residual  property
\begin{equation*}
 \Vert \mathbf{r}_m \Vert_2 =  \min_{\mathbf{y} \in {\cal K}_m({\mymatrix{C}}, \mathbf{r}^0)} \Vert   {{\mymatrix{C}}} {\mathbf{y}} -{\mathbf{d}}\Vert_2.
\end{equation*}

The main result of this subsection (Theorem \ref{thm:conditioning} below) is a result about the convergence of GMRES applied to \eqref{eq:system0} preconditioned by diagonal matrices. This result is proved under the following assumption in which (and subsequently) for every $v_N\in {\cal H}_N$ we denote by $\mathbf{v}\in \mathbb{C}^{M_N}$ the unique vector $\mathbf{v}=(v^N_1,\ldots, v^N_{M_N})^T$ such that $v_N = \sum_{j=1}^{M_N} v^N_j \psi^{N}_j$.

\begin{assumption}\label{ass:D}
$({\cal H}_N)_{N=1}^\infty$, and the associated bases $(\{\psi^N_1,\ldots,\psi^N_{M_N}\})_{N=1}^\infty$, are such that there exists a sequence of diagonal matrices $({\mymatrix{D}}_N)_{N=1}^\infty$ and $C_1, C_2>0$, independent of $N$, such that
\begin{equation}\label{eq:norm_equiv}
C_1 \big\|{\mymatrix{D}}_N^{1/2} \mathbf{w}\big\|_2 \leq \N{w_N}_{{L^2(\Gamma)}} \leq C_2 \big\|{\mymatrix{D}}_N^{1/2} \mathbf{w}\big\|_2 \quad\text{ for all } w_N \in {\cal H}_N.
\end{equation}
\end{assumption}

\begin{remark}[Relation of $C_1$ and $C_2$ in \eqref{eq:norm_equiv} to the mass matrix]
Let ${\mymatrix{M}}_N$ be the mass matrix defined by
\begin{equation}\label{eq:massmatrix}
({\mymatrix{M}}_N)_{ij} := \big(\psi_j^{N},\psi^{N}_i\big)_{L^2(\Gamma)}.
\end{equation}
Since $({\mymatrix{M}}_N\mathbf{w},\mathbf{w})_2 = \|w_N\|^2_{L^2(\Gamma)}$,
and thus $\|{\mymatrix{M}}_N^{1/2}\mathbf{w}\|_2 = \|w_N\|_{L^2(\Gamma)}$,
 \eqref{eq:norm_equiv} implies that
\begin{equation*}
C_1 \| \mathbf{v}\|_2 \leq \big\| {\mymatrix{M}}^{1/2}_N {\mymatrix{D}}^{-1/2}_N \mathbf{v}\big\|_{{L^2(\Gamma)}} \leq C_2 \|\mathbf{v}\|_2 \quad\text{ for all } v_N \in {\cal H}_N;
\end{equation*}
i.e., ${\mymatrix{D}}_N^{1/2}$ can be considered as a right-preconditioner for ${\mymatrix{M}}_N^{1/2}$, removing the $N$-dependence of the norms of ${\mymatrix{M}}_N^{1/2}$ and ${\mymatrix{M}}_N^{-1/2}$.
\end{remark}

\begin{remark}[When does Assumption \ref{ass:D} hold?]
If $(\psi^{N}_j)_{j=1}^{M_N}$ is an orthogonal basis of ${\cal H}_N$, then Assumption \ref{ass:D} is satisfied with ${\mymatrix{D}}_N= {\mymatrix{M}}_N$ and $C_1=C_2=1$; therefore, if $(\psi^{N}_j)_{j=1}^{M_N}$ is an orthonormal basis of ${\cal H}_N$, then Assumption \ref{ass:D} is satisfied with ${\mymatrix{D}}_N$ equal the identity matrix.

Lemma \ref{lem:meshes} below shows that Assumption \ref{ass:D} is satisfied (and specifies the matrices ${\mymatrix{D}}_N$) when $({\cal H}_N)_{N=1}^\infty$ are piecewise-polynomial subspaces allowing discontinuities across elements, under mild constraints on the sequence of meshes; in particular
Lemma \ref{lem:meshes} covers nodal basis functions on highly anisotropic meshes, such as the meshes highlighted in Remark \ref{rem:BestApprox}. We highlight that the assumption that discontinuities are allowed is made so that we can assume in the proof that each basis function is supported on only one element, but we expect the result to hold more generally. In particular, if $d=3$ and the sequence of meshes is regular, shape-regular, and quasi-uniform (in the sense of \cite[Definitions 4.1.4, 4.1.12, and 4.1.13]{SaSc:11}, respectively)
on a polyhedral or piecewise curved domain (in the sense of \cite[Assumptions 4.3.17 and 4.3.18]{SaSc:11}, respectively),
then Assumption \ref{ass:D} holds for a general nodal basis with ${\mymatrix{D}}_N = h^{d-1}{\mymatrix{I}}_N$ by \cite[Theorem 4.4.7]{SaSc:11}.
\end{remark}

Let $\mathbf{y}_m$ be the $m$th iterate when the linear system
\begin{equation}\label{eq:pclinear}
({\mymatrix{D}}_N^{-1/2}{\mymatrix{A}} {\mymatrix{D}}_N^{-1/2})\mathbf{y} = {\mymatrix{D}}_N^{-1/2}\mathbf{b}
\end{equation}
is solved using GMRES with zero initial guess. (Since ${\mymatrix{D}}_N$ is diagonal, the cost of calculating the action of ${\mymatrix{D}}_N^{-1/2}{\mymatrix{A}} {\mymatrix{D}}_N^{-1/2}$ is dominated by the cost of calculating the action of ${\mymatrix{A}}$.)
Let
\begin{equation}\label{eq:phiNm}
\phi_N^m :=\sum_{j=1}^{M_N} ({\mymatrix{D}}_N^{-1/2}\mathbf{y}_m)_j \psi^{N}_j,
\end{equation}
and observe that, by \eqref{eq:pclinear}, \eqref{eq:system0}, and \eqref{eq:matrix}, the Galerkin solution $\phi_N$ is given by
\begin{equation}\label{eq:phiN}
\phi_N=\sum_{j=1}^{M_N} ({\mymatrix{D}}_N^{-1/2}\mathbf{y})_j \psi^{N}_j.
\end{equation}

\begin{theorem}[Convergence of GMRES applied to the linear system \eqref{eq:pclinear}]\label{thm:conditioning}
Assume that $\mathbf{Z}$ is Lipschitz, there exists $c>0$ such that
\eqref{eq:c} holds, $\alpha$ satisfies \eqref{eq:alpha_bound_boundary}, and Assumption \ref{ass:D} holds.
With $C_1$ and $C_2$ as in  \eqref{eq:norm_equiv}, and where $\|A^\dag_{\bZ,\alpha}\|$ denotes $\|A^\dag_{\bZ,\alpha}\|_{L^2(\Gamma)\to L^2(\Gamma)}$,
 let $\beta\in [0,\pi/2)$ be defined such that
\begin{equation}\label{eq:cosbeta}
\cos\beta= \frac{c}{2\|A^\dag_{\bZ,\alpha}\|}\left(\frac{C_1}{C_2}\right)\quad\text{ and let } \quad \gamma_\beta:= 2 \sin\left(\frac{\beta}{4-2\beta/\pi}\right)
\end{equation}
(observe that $\cos \beta$ is indeed $\leq 1$ since, by definition, $C_1\leq C_2$ and $c/2 \leq \|A^\dag_{\bZ,\alpha}\|$).
Given $\varepsilon>0$, if
\begin{equation}\label{eq:iteration_bound}
m\geq  \left( \log\left(\frac{1}{\gamma_\beta}\right)\right)^{-1}
\left[
\log \left( \frac{24 \|A^\dag_{\bZ,\alpha}\|}{c}\left(\frac{C_2}{C_1}\right)^3\right)+ \log\left(\frac{1}{\varepsilon}\right) \right],
\end{equation}
then
\begin{equation}\label{eq:full_discrete}
\frac{
\N{\phi-\phi_N^m}_{{L^2(\Gamma)}}}
{\N{\phi}_{{L^2(\Gamma)}}}
\leq (1+\varepsilon)\frac{2\|A^\dag_{\bZ,\alpha}\|}{c}\left( \min_{\psi\in{\cal H}_N}\frac{\N{\phi-\psi}_{{L^2(\Gamma)}}}{\N{\phi}_{{L^2(\Gamma)}}}\right)+ \varepsilon
\end{equation}
(compare to \eqref{eq:quasioptimal}).
\end{theorem}

The key point about Theorem \ref{thm:conditioning} is that \emph{both} the bound on the number of iterations \eqref{eq:iteration_bound} \emph{and} the terms on the right-hand side of
\eqref{eq:full_discrete} other than the best-approximation error are independent of the dimension $M_N$. Therefore, the number of iterations required to solve systems involving ${\mymatrix{D}}_N^{-1/2}{\mymatrix{A}} {\mymatrix{D}}_N^{-1/2}$  to a prescribed accuracy
does \emph{not} increase as the discretisation is refined and
$M_N$ increases.
The same property holds when the conjugate-gradient method is applied to
sequences of $M_N\times M_N$ symmetric, positive-definite matrices whose condition number is bounded independently of $M_N$.

\begin{remark}[Bounds on the condition number]
Recall that, in general, a bound on the condition number
for a nonnormal matrix cannot be used to rigorously prove results about the convergence of GMRES applied to that matrix; see, e.g., \cite[Page 165]{LiTi:04}, \cite[Page 3]{Em:99}. We have no reason to expect that ${\mymatrix{A}}$ is normal, so to prove Theorem \ref{thm:conditioning} we crucially use the coercivity of $A_{I,\bZ,\alpha}$.

Nevertheless, since there is a long history of studying the condition numbers of second-kind integral equations posed on ${L^2(\Gamma)}$, we record that in the course of proving Theorem \ref{thm:conditioning} we prove that, where $C:= \|A^\dag_{\bZ,\alpha}\|_{L^2(\Gamma)\to L^2(\Gamma)}$,
\begin{equation*}
\cond\big( {\mymatrix{D}}_N^{-1/2}{\mymatrix{A}} {\mymatrix{D}}_N^{-1/2}\big) \leq \frac{2 C}{c}\left(\frac{C_2}{C_1}\right)^2,
\end{equation*}
where $\cond({\mymatrix{B}}):= \|{\mymatrix{B}}\|_2 \|{\mymatrix{B}}^{-1}\|_2$ (see \eqref{eq:GMRES3} and \eqref{eq:GMRES3a} below); i.e.,
$\cond( {\mymatrix{D}}_N^{-1/2}{\mymatrix{A}} {\mymatrix{D}}_N^{-1/2})$ is bounded independently of the dimension $M_N$.
Furthermore, by the arguments in \cite[\S III]{BeChGrLaLi:11}, \cite[Equation B.8]{MaGaSpSp:21}, and \cite[Equation 3.6.166]{At:97},
\begin{equation}\label{eq:conditionnumber}
\cond({\mymatrix{A}}) \leq\frac{2C}{c}\,  \cond({\mymatrix{M}}_N);
\end{equation}
recall that for a piecewise polynomial boundary element approximation space (in dimensions $d=2$ or $3$) using nodal basis functions on a quasiuniform mesh  (with these terms defined in \S\ref{sec:conditioning2a}), $\cond({\mymatrix{M}}_N)$ is independent of the dimension $M_N$ (see, e.g., the proof of Part (i) of Lemma \ref{lem:meshes} or \cite[Remark 4.5.3]{SaSc:11}).
\end{remark}

\begin{remark}[Calculating the entries of the Galerkin matrices for the new BIEs]
Calculating the entries of the Galerkin matrices for the new BIEs requires evaluating integrals involving \emph{only} the operators
$S$, $D$, and $D'$.
Indeed, for the direct BIE \eqref{eq:new_int}, the expression \eqref{eq:Caldb'} shows that constructing the Galerkin matrix requires evaluating integrals involving the operators above, and evaluating integrals of the form
\begin{equation}\label{eq:ibps}
\int_\Gamma \bZ(\mathbf{x})\cdot \nabla_{\Gamma} \big(S \psi_j(\mathbf{x})\big) \psi_i (\mathbf{x})\, {\rm d} s(\mathbf{x})
\end{equation}
where $\psi_i, \psi_j \in {\cal H}_N$.
It is shown in \cite[\S4.3]{SpChGrSm:11} that, using integration by parts, the integral \eqref{eq:ibps} can be evaluated in terms of integrals involving derivatives of $\psi_i$ and values (but not derivatives) of $S\psi_j$.
Constructing the Galerkin matrix of the indirect BIE \eqref{eq:new_int_indirect} requires evaluating the integral
\begin{equation}\label{eq:ibps2}
\int_\Gamma \big(K_{\vfb} \psi_j(\mathbf{x})\big) \psi_i (\mathbf{x})\, {\rm d} s(\mathbf{x}),
\end{equation}
where $K_{\vfb}$ is defined by \eqref{eq:Caldb}.
Using the expansion $\mathbf{Z} = \sum_{i=1}^d Z_i \mathbf{e}_i$ in the definition of $K_{\vfb}$, we have
\begin{equation*}
K_{\vfb} \phi(\mathbf{x})= \sum_{i=1}^d \mathbf{e}_i \cdot\int_\Gamma  \nabla_\mathbf{y}\Phi(\mathbf{x},\mathbf{y})Z_i(\mathbf{y})\phi(\mathbf{y}) \, {\rm d} s(\mathbf{y}).
\end{equation*}
Given $\mathbf{x}\in \Gamma$, $\mathbf{e}_i = (\mathbf{e}_i\cdot \mathbf{n}(\mathbf{x})) \mathbf{n}(\mathbf{x}) + \mathbf{e}_T(\mathbf{x})$, where $\mathbf{e}_T(\mathbf{x})$ is tangent to $\Gamma$ at $\mathbf{x}$.
Thus,
\begin{equation*}
K_{\vfb} \phi(\mathbf{x}) = \sum_{i=1}^d\Big( (\mathbf{e}_i\cdot \mathbf{n}(\mathbf{x})) D'\big(Z_i \phi\big)(\mathbf{x}) + \mathbf{e}_i\cdot\nabla_\Gamma S \big(Z_i \phi\big)(\mathbf{x})\Big);
\end{equation*}
the integral \eqref{eq:ibps2} can therefore be evaluated in terms of integrals only involving $D'$ and (by the discussion above regarding \eqref{eq:ibps}) $S$.
\end{remark}

\subsection{New boundary integral equations for the Laplace interior and exterior Dirichlet problems on general 2-d Lipschitz domains}\label{sec:d=2}

The biggest difference in going from $d\geq 3$ to $d=2$ is that the single-layer potential is no longer $o(1)$ at infinity, and is only $O(1)$ for a restricted class of densities; see \eqref{eq:Sasym}, \eqref{eq:Sasym2} below.
In this section, we first outline what parts of the $d\geq 3$ results in \S\ref{sec:dgeq3} immediately carry over to $d=2$.
We then present modifications of the integral equations in Theorem \ref{thm:Laplace_int} and Theorem \ref{thm:Laplace_ext} that are coercive for general 2-d Lipschitz domains when $\alpha$ is sufficiently large.

Inspecting the proof of Theorem \ref{thm:Laplace_int} in \S\ref{sec:proofs}, we see that
Parts (i), (ii), (iii), and (iv) hold when $d=2$ (i.e.,~everything apart from invertibility (v) and coercivity for sufficiently large $\alpha$ (vi)).

Similarly, inspecting the proof of Theorem \ref{thm:Laplace_ext} in \S\ref{sec:proofs}, we see that
Parts (i), (iii), (iv), and (v) hold when $d=2$ (i.e.,~everything apart from the indirect formulation (ii) and coercivity for sufficiently large $\alpha$ (vi)), although, firstly,
$\alpha u_\infty$ must be added to the right-hand side of the BIE \eqref{eq:new_ext}, where $u_\infty$ is the limit of $u$ at infinity\footnote{This is because Green's integral representation for the solution of the Laplace EDP with $d=2$ takes the form $u(\mathbf{x}) = -{\cal S}\partial_n^+ u (\mathbf{x}) + {\cal D} \gamma^+ u (\mathbf{x}) + u_\infty$ for $\mathbf{x} \in \Oe$.} and, secondly, Part (v) holds when $d=2$ provided the constant $a$ in the fundamental solution \eqref{eq:fund} is not equal to the capacity of $\Gamma$, $\mathrm{Cap}_\Gamma$ (defined in, e.g., \cite[Page 263]{Mc:00}), which holds, in particular, if $a>\mathop{{\rm diam}}(\Gamma)$.

Let
\begin{equation}\label{eq:P}
P_\Gamma \phi (\mathbf{x}):= \frac{1}{|\Gamma|} \int_\Gamma \phi (\mathbf{y})\, {\rm d} s(\mathbf{y})= \frac{1}{|\Gamma|} \big(\phi,1\big)_{L^2(\Gamma)} \quad \text{ for }\mathbf{x}\in\Gamma;
\end{equation}
i.e., $P_\Gamma \phi$ is the mean value of $\phi$. Observe that $P_\Gamma^2=P_\Gamma$ and $P_\Gamma'=P_\Gamma$. Let $Q_\Gamma:=I-P_\Gamma$.

We give two theorems: the first for general 2-d Lipschitz domains, the second for 2-d star-shaped Lipschitz domains. Recall that we are assuming throughout that $\bZ$ is real-valued.

\begin{theorem}[New integral equations for Laplace IDP and EDP in 2-d]\label{thm:2d} Suppose that $\bZ\in (L^\infty(\Gamma))^2$ and $\alpha,\beta\in \mathbb{R}$.

(i) \textbf{\emph{IDP direct formulation.}} Let $u$ be the solution of the Laplace IDP of Definition \ref{def:idp} with $d=2$ and $g_D\in H^1(\Gamma)$.
Let
\begin{equation}\label{eq:TIprime}
T'_{I, \bZ,\alpha,\beta}:= Q_\Gamma A'_{I,\bZ,\alpha} Q_\Gamma + \beta P_\Gamma.
\end{equation}
Then $\partial_n^- u$ satisfies
\begin{equation}\label{eq:2d_IDP_direct}
T'_{I,\bZ,\alpha,\beta}(\partial_n^- u) =Q_\Gamma B_{I,\bZ,\alpha} \,g_D.
\end{equation}

(ii) \textbf{\emph{IDP indirect formulation.}}
Let
\begin{equation*}
T_{I, \bZ,\alpha,\beta}:= Q_\Gamma A_{I,\bZ,\alpha} Q_\Gamma + \beta P_\Gamma.
\end{equation*}
Given $g_D\in {L^2(\Gamma)}$, if $\phi$ satisifes
\begin{equation}\label{eq:2d_IDP_indirect}
T_{I,\bZ,\alpha,\beta}\phi =- g_D,
\end{equation}
then, if $d=2$,
\begin{equation}\label{eq:u2d_int}
u:= (\mathcal{K}_{\vfb}-\alpha{\cal S}) Q_\Gamma\phi + P_\Gamma A_{I,\bZ,\alpha} Q_\Gamma \phi - \beta P_\Gamma\phi
\end{equation}
is the solution of the Laplace IDP of Definition \ref{def:idp2}.

(iii) \textbf{\emph{EDP direct formulation.}} Let $u$ be the solution of the Laplace EDP of Definition \ref{def:edp} with $d=2$ and $g_D\in H^1(\Gamma)$.
Let
\begin{equation}\label{eq:TEprime}
T'_{E, \bZ,\alpha,\beta}:= Q_\Gamma A'_{E,\bZ,\alpha} Q_\Gamma + \beta P_\Gamma.
\end{equation}
Then $\partial_n^+ u$ satisfies
\begin{equation}\label{eq:2d_EDP_direct}
T'_{E,\bZ,\alpha,\beta}(\partial_n^+ u) =Q_\Gamma B_{E,\bZ,\alpha} \,g_D.
\end{equation}

(iv) \textbf{\emph{EDP indirect formulation.}}
Let
\begin{equation*}
T_{E, \bZ,\alpha,\beta}:= Q_\Gamma A_{E,\bZ,\alpha} Q_\Gamma + \beta P_\Gamma.
\end{equation*}
Given $g_D\in {L^2(\Gamma)}$, if $\phi$ satisifes
\begin{equation}\label{eq:2d_EDP_indirect}
T_{E,\bZ,\alpha,\beta}\phi =g_D,
\end{equation}
then, if $d=2$,
\begin{equation}\label{eq:u2d_ext}
u:= (\mathcal{K}_{\vfb}+\alpha{\cal S}) Q_\Gamma\phi - P_\Gamma A_{I,\bZ,\alpha} Q_\Gamma \phi + \beta P_\Gamma\phi
\end{equation}
is the solution of the Laplace EDP of Definition \ref{def:edp2}.

(v) \textbf{\emph{Coercivity.}} If $\bZ\in (C^{0,1}(\Gamma))^2$ satisfies \eqref{eq:c}, $\alpha$ satisfies \eqref{eq:alpha_bound_boundary}, and $\beta=c/2$, then
$T'_{I,\bZ,\alpha,\beta}$, $T_{I,\bZ,\alpha,\beta}$, $T'_{E,\bZ,\alpha,\beta}$, and $T_{E,\bZ,\alpha,\beta}$ are all coercive on ${L^2(\Gamma)}$ with coercivity constant $c/2$.
\end{theorem}

\begin{theorem}[New integral equations for 2-d star-shaped domains]\label{thm:2dstar}

\

(i) \textbf{\emph{IDP direct formulation.}} Let $u$ be the solution of the Laplace IDP of Definition \ref{def:idp} with $d=2$ and $g_D\in H^1(\Gamma)$.
Then $\partial_n^- u$ satisfies
\begin{equation}\label{eq:2dstar1}
\left(A'_{I,\mathbf{x},0} - \frac{|\Gamma|}{4\pi}P_\Gamma\right) \partial_n^- u =B_{I,\mathbf{x},0} \,g_D.
\end{equation}

(ii) \textbf{\emph{IDP indirect formulation.}}
Given $g_D\in {L^2(\Gamma)}$, if $\phi$ satisifes
\begin{equation}\label{eq:2dstar2}
\left(A_{I,\mathbf{x},0} - \frac{|\Gamma|}{4\pi}P_\Gamma\right)
\phi =- g_D,
\end{equation}
then, if $d=2$,
\begin{equation*}
u:= \mathcal{K}_{\vfb} \phi + \frac{|\Gamma|}{4\pi}P_\Gamma\phi
\end{equation*}
is the solution of the Laplace IDP of Definition \ref{def:idp2}.

(iiii) \textbf{\emph{EDP direct formulation.}} Let $u$ be the solution of the Laplace EDP of Definition \ref{def:idp} with $d=2$ and $g_D\in H^1(\Gamma)$.
Then $\partial_n^+ u$ satisfies
\begin{equation}\label{eq:2dstar3}
\left(A'_{E,\mathbf{x},0} + \frac{|\Gamma|}{4\pi}P_\Gamma\right) \partial_n^+ u =B_{E,\mathbf{x},0} \,g_D.
\end{equation}

(iv) \textbf{\emph{EDP indirect formulation.}}
Given $g_D\in {L^2(\Gamma)}$, if $\phi$ satisifes
\begin{equation}\label{eq:2dstar4}
\left(A_{E,\mathbf{x},0} + \frac{|\Gamma|}{4\pi}P_\Gamma\right)
\phi = g_D,
\end{equation}
then, if $d=2$,
\begin{equation*}
u:= \mathcal{K}_{\vfb} \phi + \frac{|\Gamma|}{4\pi}P_\Gamma\phi
\end{equation*}
is the solution of the Laplace EDP of Definition \ref{def:edp2}.

(v) \textbf{\emph{Coercivity.}} If $\Omega_-$  is star-shaped with respect to a ball of radius $\kappa$ (i.e.~\eqref{eq:ssg} holds),
then each of the integral operators on the left-hand sides of \eqref{eq:2dstar1}, \eqref{eq:2dstar2}, \eqref{eq:2dstar3}, and \eqref{eq:2dstar4} is coercive on ${L^2(\Gamma)}$ with coercivity constant $\kappa/2$.
\end{theorem}

\section{Discussion of the ideas behind the new BIEs and links to previous work}\label{sec:idea}

\subsection{How the BIEs arise}
The indirect BIE \eqref{eq:new_int_indirect} for the IDP arises from imposing the boundary condition on the ansatz $u= (\mathcal{K}_{\vfb} -\alpha {\cal S})\phi$ via taking the nontangential limit.
Similarly, the indirect BIE \eqref{eq:new_ext_indirect} for the EDP arises from the ansatz $u= (\mathcal{K}_{\vfb} +\alpha {\cal S})\phi$.
For the indirect BIEs for $d=2$ in Theorem \ref{thm:2d}, the idea is the same, except now a) the density in the ansatz is not a general $L^2(\Gamma)$ function (so that ${\cal S}\phi$ has the correct behaviour at infinity), and b) extra terms are added to the ansatz to ensure that the resulting BIE is still coercive on ${L^2(\Gamma)}$.

For the direct BIE \eqref{eq:new_int} for the IDP, recall that $u={\cal S} \partial_n^- u - {\cal D} \gamma^- u$ by Green's integral representation. The direct BIE \eqref{eq:new_int} then arises from considering
\begin{equation*}
-\bZ \cdot \widetilde{\gamma}^- (\nabla u ) + \alpha \gamma^- u.
\end{equation*}
Similarly, the direct BIE \eqref{eq:new_ext} for the EDP arises from considering
\begin{equation*}
\bZ \cdot \widetilde{\gamma}^+ (\nabla u ) + \alpha \gamma^+ u,
\end{equation*}
with $u = -{\cal S}\partial_n^+ u + {\cal D} \gamma^+ u$. Alternatively, since (informally) $\bZ \cdot \nabla = (\bZ \cdot \mathbf{n}) \partial_n + \bZ \cdot \nabla_{\Gamma}$, the direct BIE \eqref{eq:new_int} can be obtained by adding (i) $(\bZ\cdot \mathbf{n})$ multiplied by the standard direct second-kind BIE
\begin{equation}\label{eq:inteq2}
\left(\frac{1}{2} I - D'\right)\partial_n^- u = -H\gamma^- u,
\end{equation}
(ii) $-\bZ\cdot \nabla_{\Gamma}$ applied to the standard direct first-kind BIE
\begin{equation}\label{eq:inteq1}
S\partial_n^- u= \left(\frac{1}{2} I + D\right) \gamma^- u,
\end{equation}
and (iii) $\alpha$ multiplied by \eqref{eq:inteq1}. Similar considerations hold for the direct BIE \eqref{eq:new_ext}, and the 2-d direct BIEs of Theorems \ref{thm:2d} and \ref{thm:2dstar}, where, additionally, one uses that $P_\Gamma (\partial^{\pm}_n u)=0$ (see Lemma \ref{lem:Neumann_trace}).

\subsection{The other BVPs solved by the new BIEs}
The BIEs introduced in \S\ref{sec:main_results} to solve the Dirichlet problem can be used to solve other Laplace BVPs. Although the focus of this paper is on solving the Dirichlet problem, we highlight this fact here since
these other BVPs affect the properties of the new BIEs.

For example, the BIO $A'_{E,\bZ,\alpha}$ used to solve the EDP in Theorem \ref{thm:Laplace_ext} can also be used to solve the Laplace interior oblique Robin problem, i.e., the problem of finding $u$ in $\Omega_-$ satisfying $\Delta u=0$ and
\begin{equation}\label{eq:obliqueLintro}
(\mathbf{Z}\cdot\mathbf{n}) \partial_n^- u + \mathbf{Z} \cdot \nabla_\Gamma( \gamma^- u) + \alpha \, \gamma^- u =  g \quad\text{ on } \Gamma;
\end{equation}
see Definition \ref{prob:ioipL} and Theorem \ref{thm:star_in_obliqueLE} below. Similarly, the BIO $A'_{I,\bZ,\alpha}$ used to solve the IDP in Theorem \ref{thm:Laplace_int} can also be used to solve the Laplace exterior oblique Robin problem; see Definition \ref{prob:eoipL} and Theorem \ref{thm:star_in_obliqueL}  below.
This relationship means that the injectivity results implicit in Part (v) of Theorems \ref{thm:Laplace_int} and \ref{thm:Laplace_ext} are obtained by proving uniqueness of these oblique Robin problems; see \S\ref{subsec:oblique3}.

\subsection{The use of similar BIEs by Calder\'on \cite{Ca:85} and Medkov\'a \cite{Me:18}}
Calder\'on \cite{Ca:85} used indirect versions of the BIEs in Theorems \ref{thm:Laplace_int} and \ref{thm:Laplace_ext} with $\alpha=0$ to prove wellposedness results about the Dirichlet problem and the oblique derivative problem (i.e., \eqref{eq:obliqueLintro} with $\alpha=0$) with data in $L^p(\Gamma)$.
Indeed, \cite{Ca:85} posed the ansatz $u= \mathcal{K}_{\vfb} \phi$ for the IDP, which gives the BIE $A_{I, \bZ, 0} \phi =-g_D$ \cite[Page 39]{Ca:85}, and posed the ansatz $u={\cal S} \phi$ for the oblique derivative problem, which gives the BIE $A'_{E,\bZ,0} \phi = g$ \cite[Page 45]{Ca:85}.
Furthermore, Medkov\'a \cite[\S5.23]{Me:18}
posed the ansatz $u = {\cal S}\phi$ for the interior oblique Robin problem, resulting in $A'_{I,\bZ,\alpha}\phi = -g$.

In both \cite{Ca:85} and \cite{Me:18}, the BIOs are proved to be Fredholm of index zero on $L^2(\Gamma)$; see \cite[Page 39]{Ca:85} (where the result is proved to hold on a slightly wider range of $L^p(\Gamma)$ spaces) and \cite[Proposition 5.23.2]{Me:18}.

\subsection{The main new properties of the BIEs of this paper:~coercivity for appropriate $\alpha$}\label{sec:Rellich_recap}
Building on the work of Calder\'on and Medkova, we show that the BIOs are not only Fredholm of index zero on $L^2(\Gamma)$, but \emph{invertible} for general Lipschitz domains as soon as $\alpha>0$, and, crucially, \emph{coercive} if $\alpha$ is chosen appropriately (so also \emph{coercive plus compact} for all $\alpha>0$). For star-shaped domains this coercivity can be proved using a simple modification of Calder\'on's proof that the BIOs are Fredholm of index zero (see Lemmas \ref{lem:Calderon} and \ref{lem:Calderon2} below).
For general domains this coercivity (for appropriate $\alpha$) is proved using Rellich-type identities (with this method also giving an alternative proof of coercivity for star-shaped domains). Recall that identities arising from multiplying $\Delta u$ by a derivative of $u$ are associated with the name Rellich, due to Rellich's introduction of the multiplier $\mathbf{x} \cdot \nabla u$ for the Helmholtz equation in \cite{Re:40}; these identities have been well-used in the study of the Laplace, Helmholtz, and other elliptic equations, see, e.g., the overviews in \cite[Pages 111 and 112]{Ke:94}, \cite[\S5.3]{ChGrLaSp:12}, \cite[\S1.4]{MoSp:14}. Verchota \cite{Ve:84} famously used Rellich identities to prove invertibility of $\frac{1}{2} I - D$ and $\frac{1}{2} I - D'$ on ${L^2(\Gamma)}$ (see Remark \ref{rem:Verchota} below) and  Medkov\'a  \cite[\S5.23]{Me:18} also used Rellich identities to prove that $A'_{I, \bZ, \alpha}$ is invertible for sufficiently large $\alpha$ \cite[Lemma 5.23.1, Prop.~5.23.2, Theorem 5.23.4]{Me:18}.

Our coercivity results are proved using the identity arising from multiplying $\Delta u$ by $\bZ \cdot \nabla u + \alpha u$ (see Lemma \ref{lem:Rellich} below); our use of a multiplier that is a linear combination of $u$ and a derivative of $u$ is inspired by the use of such multipliers by Morawetz \cite{Mo:61, MoLu:68, Mo:75}, and the particular identity we use also appears as \cite[Equation 2.28]{KuMaVa:02}.
As recalled in \S\ref{sec:rationale}, the idea of proving coercivity of Laplace BIOs in the trace spaces goes back to N\'ed\'elec and Planchard \cite{NePl:73}, Le Roux \cite{LeRo:74}, Hsiao and Wendland \cite{HsWe:77}, and Steinbach and Wendland \cite{StWe:01}, with this method based on using Green's identity (i.e. multiplying $\Delta u$ by $u$).
The idea of proving coercivity of second-kind BIOs in ${L^2(\Gamma)}$ using Rellich-type identities was introduced in
\cite{SpChGrSm:11} for a particular  Helmholtz BIE on star-shaped domains and then further developed in  \cite{SpKaSm:15} for the standard second-kind Helmholtz BIE on smooth convex domains. The main differences between \cite{SpChGrSm:11, SpKaSm:15} and the present paper are that (i) \cite{SpChGrSm:11, SpKaSm:15} only consider direct BIEs for the exterior Helmholtz Dirichlet problem whereas the present paper considers direct and indirect BIEs for the interior and exterior Laplace Dirichlet problems and (ii) \cite{SpChGrSm:11, SpKaSm:15} only prove coercivity under geometric restrictions on $\Gamma$ (which is somewhat expected for the high-frequency Helmholtz equation; see \cite{BeSp:11}, \cite[\S6.3.2]{ChSpGiSm:20}), namely star-shapedness with respect to a ball for \cite{SpChGrSm:11} and strict convexity and a piecewise analytic $C^3$ boundary for  \cite{SpKaSm:15}, whereas the present paper proves coercivity of Laplace BIOs for general Lipschitz domains.

\subsection{Combined-potential ansatz for solutions of Laplace's equation}
A key difference between the indirect BIEs in the present paper and those in \cite{Ca:85} is that ours arise from the ansatz $u= (\mathcal{K}_{\vfb} -\alpha {\cal S})\phi$ for the solution of the Laplace IDP, whereas \cite{Ca:85} poses the ansatz $u = \mathcal{K}_{\vfb}\phi$. We saw in the discussion above that the presence of the parameter $\alpha$ -- i.e., the fact that we use a combined-potential ansatz -- is crucial for proving coercivity of our BIOs.

The combined-potential ansatz is also crucial to proving uniqueness for cases where coercivity does not hold. Indeed, using a linear combination of double- and single-layer potentials to find solutions of the Helmholtz equation is standard, and goes back to \cite{BrWe:65, Le:65, Pa:65}, with the motivation to ensure uniqueness at all wavenumbers (see \S\ref{sec:Helmholtz}).
Using such a combination for Laplace's equation is less common, but this was done by D.~Mitrea in \cite[Theorem 5.1]{MitreaD:97} and subsequently by Medkov\'a in \cite{Me:99}.
The rationale for this combined ansatz is similar, namely that the standard indirect second-kind equations (based on a double-layer-potential ansatz) have non-trivial null spaces for multiply connected domains (with these characterised in \cite{KulkarniEtal:05, MitreaD:97}) but the BIOs resulting from a combined double- and single-layer potential ansatz are invertible no matter the topology of $\Omega_-$; see \cite[Theorem 5.15.2]{Me:18} (for $d\geq 3$) and \cite[Theorem 5.15.3]{Me:18} (for $d=2$). The BIOs in \S\ref{sec:main_results} are also invertible (and even, for appropriate $\alpha$, coercive) no matter the topology of $\Omega_-$.

\section{Proofs of the main results}\label{sec:proofs}

In this section we prove all of the results in \S\ref{sec:main_results} \emph{apart from} the invertibility results in Part (v) of Theorem \ref{thm:Laplace_int}/\ref{thm:Laplace_ext}.
As discussed in \S\ref{sec:idea}, these invertibility results are equivalent to uniqueness of the Laplace interior and exterior oblique Robin problems, and these uniqueness results are proved in \S\ref{sec:oblique}.
Indeed, Part (v) of Theorem \ref{thm:Laplace_int} follows from Corollary \ref{cor:oblique:unique2}, and Part (v) of Theorem \ref{thm:Laplace_ext} follows from Corollary \ref{cor:oblique:unique}.

\subsection{Proofs of Parts (i), (ii), and (iii) of Theorems \ref{thm:Laplace_int} and \ref{thm:Laplace_ext}}\label{sec:4.1}

For Part (i) of Theorem \ref{thm:Laplace_int}, first recall that the standard direct BIEs for the IDP (corresponding to the top left of Table \ref{tab:bies}) are
\eqref{eq:inteq1} and \eqref{eq:inteq2}.
If $g_D\in H^1(\Gamma)$, then $\partial_n^- u\in{L^2(\Gamma)}$ (by Theorem \ref{thm:Necas}), and then the mapping properties \eqref{eq:map1} of $S$ and $D$  imply that both sides of \eqref{eq:inteq1} are in $H^1(\Gamma)$. Taking the surface gradient, $\nabla_{\Gamma}$, of \eqref{eq:inteq1} yields the (vector) integral equation in $({L^2(\Gamma)})^d$
\begin{equation}\label{eq:inteq3}
\nabla_{\Gamma} S}%\pdiff{S_k}{T}(\partial_n^- u)= \nabla_{\Gamma}\left(\frac{1}{2} I + D\right) \gamma^- u.
\end{equation}
Taking $(\bZ\cdot\mathbf{n})$ times the scalar equation \eqref{eq:inteq2}, minus $\bZ$ dot the vector equation \eqref{eq:inteq3}, plus $\alpha$ times \eqref{eq:inteq1} yields \eqref{eq:new_int}. The proof of Part (i) of Theorem \ref{thm:Laplace_ext} (i.e., that \eqref{eq:new_ext} holds) is very similar.

For Part (ii) of Theorem \ref{thm:Laplace_int}, first recall that $\mathcal{K}_{\vfb}\phi$ and ${\cal S}\phi$ are both in $C^2(\Omega_-)$ and satisfy Laplace's equation
(for $\mathcal{K}_{\vfb}$ this was recalled in \S\ref{sec:recap}). When $\phi \in {L^2(\Gamma)}$, ${\cal S}\phi \in H^{3/2}(\Omega_-)$ by \eqref{eq:LPmap} and then $({\cal S}\phi)^*\in {L^2(\Gamma)}$ by Part (iii) of Theorem \ref{thm:JeKe:95}.
As recalled in \S\ref{sec:recap}, $(\mathcal{K}_{\vfb}\phi)^* \in {L^2(\Gamma)}$ by \cite{Ve:84}, and thus $u$ defined by \eqref{eq:indirectansatz1} satisfies $u^*\in {L^2(\Gamma)}$.
To show that  $\phi$ satisfies the BIE \eqref{eq:new_int_indirect}, we take the non-tangential limit of \eqref{eq:indirectansatz1}, using \eqref{eq:jumpCald} and that, by Lemma \ref{lem:A9}, $\widetilde{\gamma}^- ({\cal S}\phi) = \gamma^-({\cal S}\phi)$, where $\gamma^-({\cal S}\phi)$ is given by the first jump relation in
\begin{equation}\label{eq:jump1}
\gamma^{\pm}{\cal S} = S, \qquad \partial_n^{\pm}{\cal S} = \mp \frac{1}{2} I + D'.
\end{equation}
(see, e.g., \cite[Page 219]{Mc:00} or \cite[Equation 2.41]{ChGrLaSp:12}).

Part (ii) of Theorem \ref{thm:Laplace_ext} follows in an analogous way, except that we now need to show that $u$ defined by \eqref{eq:indirect_ext_u}
satisfies $u(\mathbf{x})= o(|\mathbf{x}|^{3-d})$ when $d=3$ as $|\mathbf{x}|\rightarrow \infty$; these asymptotics follow from
the first bound in \eqref{eq:Sasym} and the bound
\begin{equation}\label{eq:Caldasym}
|\mathcal{K}_{\vfb}\phi(\mathbf{x})| = O(|\mathbf{x}|^{1-d}) \quad\text{ as } |\mathbf{x}|\rightarrow \infty,
\end{equation}
which is proved in a similar way to the bound on the double-layer potential in \cite[Equation 3.23]{SaSc:11}.

Part (iii) of both Theorems \ref{thm:Laplace_int} and \ref{thm:Laplace_ext} follows from the mapping properties \eqref{eq:map} with $s=1/2$, and from
the fact, recalled in \S\ref{sec:recap}, that $K_{\vfb}:{\LtG\rightarrow \LtG}$ and $\nabla_{\Gamma} S:{\LtG\rightarrow \LtG}$, with these mappings continuous.

\subsection{Proofs of Part (iv) of Theorems  \ref{thm:Laplace_int} and \ref{thm:Laplace_ext} (coercivity up to a compact perturbation).}

\begin{lemma}\label{lem:Calderon}
If $d\geq 2$, $\Gamma$ is Lipschitz and $\bZ \in (C(\Gamma))^d$ then
$K_{\vfb} + K_{\vfb}'$
is compact in ${L^2(\Gamma)}$. Thus
there exists a compact operator $ C:{L^2(\Gamma)}\rightarrow {L^2(\Gamma)}$ such that
\begin{equation*}
\big( K_\bZ \phi,\phi\big)_{{L^2(\Gamma)}}= \big(C \phi,\phi\big)_{{L^2(\Gamma)}} \quad\mbox{for all real-valued } \phi\in {L^2(\Gamma)}.
\end{equation*}
\end{lemma}

Part (iv) of both Theorems \ref{thm:Laplace_int} and \ref{thm:Laplace_ext} follow by combining Lemma \ref{lem:Calderon} with the assumption \eqref{eq:c} and the fact that $S$ is compact on ${L^2(\Gamma)}$ (via the mapping property in \eqref{eq:map1} with $s=1/2$).

\

\begin{proof}[Proof of Lemma \ref{lem:Calderon}]
Since $\Phi(\mathbf{x},\mathbf{y})$ is a function of $|\mathbf{x}-\mathbf{y}|$, $\nabla_\mathbf{x} \Phi(\mathbf{x},\mathbf{y})= - \nabla_\mathbf{y}\Phi(\mathbf{x},\mathbf{y})$; the definitions of $K_\bZ$ \eqref{eq:Caldb} and $K'_\bZ$ \eqref{eq:Caldb'} then imply that, for all $\phi\in {L^2(\Gamma)}$,
\begin{equation}\label{eq:K3}
\big( K_\bZ + K_{\bZ}'\big)\phi(\mathbf{x})= \int_\Gamma \big(\bZ(\mathbf{y})-\bZ(\mathbf{x})\big)\cdot\nabla_\mathbf{y} \Phi(\mathbf{x},\mathbf{y}) \phi(\mathbf{y})\, {\rm d} s(\mathbf{y}).
\end{equation}
If $\bZ \in (C^{0,\beta}(\Gamma))^d$ for $\beta>0$, then the kernel of the integral on the right-hand side of \eqref{eq:K3} is weakly singular, and thus the operator is compact on ${L^2(\Gamma)}$ by, e.g., the combination of \cite[Part 3 of the theorem on Page 49]{Pr:91} and
 the Riesz-Thorin interpolation theorem (see, e.g., \cite[Theorem 6.27]{Fo:84a}), where the latter is used to verify the hypothesis of the former. Therefore, the result of this lemma follows if we can show that if, for all $\beta>0$, $K_{\bZ} + K'_{\bZ}$ is compact for all $\bZ \in (C^{0,\beta}(\Gamma))^d$, then $K_{\bZ} + K'_{\bZ}$ is compact for all $\bZ \in (C(\Gamma))^d$.

Given $\bZ \in (C(\Gamma))^d$, there exist $\beta>0$ and $\bZ_\ell \in (C^{0,\beta}(\Gamma))^d$ for all $\ell\in \mathbb{N}$ such that $\|\bZ_\ell-\bZ\|_{L^\infty}\rightarrow 0$ as $\ell \rightarrow \infty$.
By \eqref{eq:Caldb'}, the operator $K_{\mathbf{Z}}'$ can be written $K_\bZ' = \bZ \cdot \mathbf{T}$, where $\mathbf{T}:L^2(\Gamma)\to (L^2(\Gamma))^d$ is bounded by the results of \cite{CoMcMe:82} and \cite{Ve:84} (as discussed in \S\ref{sec:recap}).
Let $K_{\bZ_\ell}' = \bZ_\ell\cdot \mathbf{T}$; then
\begin{equation*}
\N{K_{\bZ_\ell}' - K_{\bZ}'}_{L^2(\Gamma)\to {L^2(\Gamma)}}= \N{\bZ_\ell \cdot \mathbf{T} - \bZ\cdot \mathbf{T}}_{{L^2(\Gamma)}\to{L^2(\Gamma)}} \leq \N{\bZ_\ell -\bZ}_{(L^\infty(\Gamma))^d} \N{\mathbf{T}}_{{L^2(\Gamma)}\to{L^2(\Gamma)}}\to 0
\end{equation*}
as $\ell\to\infty$.
Therefore also $K_{\bZ_\ell} \to K_{\bZ}$, so that $K_{\bZ_\ell} + K_{\bZ_\ell}' \to K_{\bZ}+ K_{\bZ}'$. Since the space of compact operators is closed, $K_{\bZ}+ K_{\bZ}'$ is compact.
\end{proof}

\subsection{Proof of Theorem \ref{thm:star} (coercivity for $\Omega_-$ that are  star-shaped with respect  to a ball)}

Theorem \ref{thm:star} is an immediate consequence of combining (i) the following special case of Lemma \ref{lem:Calderon}, (ii) the definitions of $A'_{I,\bZ,\alpha}$ and $A_{I,\bZ,\alpha}$ in \eqref{eq:opLIZ} and
$A'_{E,\bZ,\alpha}$ and $A_{E,\bZ,\alpha}$ in \eqref{eq:opLEZ}, and (iii) the inequality $(S\phi,\phi)_{{L^2(\Gamma)}}\geq 0$ for all $\phi \in {L^2(\Gamma)}$. The inequality in (iii) is well-known, following from Green's identity, and is a special case of Lemma \ref{lem:key} below with $\widetilde{\bZ} ={\bf 0}$.

\begin{lemma}[Key lemma for coercivity for star-shaped $\Omega_-$]\label{lem:Calderon2}
Let $\Gamma$ be Lipschitz. If $d\geq 3$ then
\begin{equation}\label{eq:symmetry3d}
K_\mathbf{x} + K_\mathbf{x}'+ (d-2)S=0 \,\,\text{and thus}\,\, \big(K_\mathbf{x} \phi,\phi\big)_{{L^2(\Gamma)}} + \frac{d-2}{2}\big(S\phi,\phi\big)_{L^2(\Gamma)}=0\,\,\mbox{for all real-valued } \phi \in{L^2(\Gamma)}.
\end{equation}
If $d=2$ then
\begin{equation}\label{eq:symmetry2d}
K_\mathbf{x} + K_\mathbf{x}'+ \frac{|\Gamma|}{2\pi}P_\Gamma =0 \quad\text{and thus}
\,\,\left( \left( K_\mathbf{x} + \frac{|\Gamma|}{4\pi}P_\Gamma\right) \phi,\phi\right)_{{L^2(\Gamma)}}= 0 \,\, \mbox{for all real-valued } \phi \in {L^2(\Gamma)},
\end{equation}
where $P_\Gamma$ is defined by \eqref{eq:P}.
\end{lemma}

\begin{proof}[Proof of Lemma \ref{lem:Calderon2}]
By \eqref{eq:fund}, when $d\geq 3$, $(\mathbf{y}-\mathbf{x})\cdot \nabla_\mathbf{y} \Phi(\mathbf{x},\mathbf{y}) = -(d-2)\Phi(\mathbf{x},\mathbf{y})$, and when $d=2$, $(\mathbf{y}-\mathbf{x})\cdot \nabla_\mathbf{y} \Phi(\mathbf{x},\mathbf{y}) =-1/2\pi$. The results then follow from \eqref{eq:K3} with $\bZ(\mathbf{x})=\mathbf{x}$.
\end{proof}

\begin{remark}[Link with the work of Fabes, Sand, and Seo \cite{FaSaSe:92}]
The analogue of \eqref{eq:symmetry3d} when $\Gamma$ is the graph of a function (i.e., the boundary of a hypograph) appears in the first sentence after the first displayed equation on \cite[Page 133]{FaSaSe:92}.
Indeed, the analogue of the operator $K_{\vfb}'$ for the hypograph with $\bZ = {\bf e}_d$ (i.e., the unit vector pointing in the $x_d$ direction) arises in \cite{FaSaSe:92} when they apply the Rellich identity \eqref{eq:id1} below with $u= {\cal S} \phi$, as part of their proof that $\lambda I - D'$ is invertible on $L^2(\Gamma)$ for $\lambda \in \mathbb{R}$ with $|\lambda|\geq 1/2$.
\end{remark}

\subsection{Proof of Part (vi) of Theorems \ref{thm:Laplace_int} and \ref{thm:Laplace_ext} (coercivity for general $\Omega_-$)}

\begin{lemma}[Key lemma for coercivity for general $\Omega_-$]\label{lem:key}
Suppose that $\Omega_- \subset \mathbb{R}^d$ is Lipschitz, $\widetilde{\bZ} \in W^{1,\infty}(\mathbb{R}^d)$ with compact support, and $\alpha\in \mathbb{R}$ satisfies the lower bound
\begin{equation}\label{eq:alpha_bound}
2\alpha \geq 2\left(\sup_{\mathbf{x} \in \mathbb{R}^d} \big\|D\widetilde{\bZ}(\mathbf{x})
\big\|_2 \right)+ \big\|\nabla \cdot \widetilde{\bZ}\big\|_{L^\infty(\mathbb{R}^d)}
\end{equation}
(where $D\widetilde{\bZ}$ is the matrix with $(i,j)$th element $\partial_i \widetilde{Z}_j$ and $\|\cdot\|_2$ denotes the operator norm on $\mathbb{R}^d\times \mathbb{R}^d$ induced by the Euclidean norm on $\mathbb{R}^d$).
If $d\geq 3$ then
\begin{equation}\label{eq:key}
\pm
\big( K_{\widetilde{\bZ}} \phi,\phi)_{{L^2(\Gamma)}}
+ \alpha
(S \phi,\phi)_{{L^2(\Gamma)}} \geq 0
\end{equation}
for all real-valued $\phi\in{L^2(\Gamma)}$. If $d=2$, then \eqref{eq:key} holds for all real-valued $\phi\in {L^2(\Gamma)}$ with $P_\Gamma\phi =0$, where $P_\Gamma$ is defined by \eqref{eq:P}.
\end{lemma}

We first show how the coercivity results of Theorems \ref{thm:Laplace_int} and \ref{thm:Laplace_ext} are a consequence of Lemma \ref{lem:key} combined with the following lemma.

\begin{lemma}\label{lem:Lipschitz}
Given $\bZ\in (C^{0,1}(\Gamma))^d$ with non-zero Lipschitz constant, there exists a compactly supported $\widetilde{\bZ}_{\rm ext}\in (C^{0,1}(\mathbb{R}^d))^d$ with the same Lipschitz constant as $\bZ$ and such that $\widetilde{\bZ}_{\rm ext}|_\Gamma = \bZ$.
\end{lemma}

The proof of Lemma \ref{lem:Lipschitz} is given in  Appendix \ref{app:Lipschitz}. Note that, by the Kirszbraun theorem \cite{Ki:34}, \cite{Va:45},
$\bZ\in (C^{0,1}(\Gamma))^d$ can be extended to a function $\bZ_{\rm ext}\in (C^{0,1}(\mathbb{R}^d))^d$ with the same (non-zero) Lipschitz constant, so to prove Lemma \ref{lem:Lipschitz} we only need to show that there exists an extension with compact support.

\

\begin{proof}[Proof of Part (vi) of Theorems \ref{thm:Laplace_int} and \ref{thm:Laplace_ext} assuming Lemmas \ref{lem:key} and \ref{lem:Lipschitz}]
Given $\bZ$, by Lemma \ref{lem:Lipschitz} there exists a compactly-supported  Lipschitz extension of $\bZ$ to $\mathbb{R}^d$ with the same Lipschitz constant; call this $\widetilde{\bZ}$.
This $\widetilde{\bZ}$ then satisfies the assumptions of Lemma \ref{lem:key}, and
the inequality \eqref{eq:alpha_bound_boundary} then ensures that \eqref{eq:alpha_bound} holds (where we have used the inequality $\|A\|_2^2 \leq \sum_i \sum_j |(A)_{ij}|^2$ to show that $\sup_\mathbf{x} \|D\widetilde{\bZ}(\mathbf{x})\|_2\leq d L_{\bZ}$).
 Thus \eqref{eq:key} holds (with $K_{\widetilde{\bZ}}$ replaced by $K_{\bZ}$) and the coercivity results follow from the definitions of $A'_{I,\bZ,\alpha}$ and $A_{I,\bZ,\alpha}$ \eqref{eq:opLIZ} and
$A'_{E,\bZ,\alpha}$ and $A_{E,\bZ,\alpha}$ \eqref{eq:opLEZ} and the inequality \eqref{eq:c} on $\bZ\cdot\mathbf{n}$.
\end{proof}

\

The proof of Lemma \ref{lem:key} is based on the following identity. The relationship of this identity to other similar identities
in the literature  is discussed in \S\ref{sec:Rellich_recap}, and we note, in particular, that this identity appears as \cite[Equation 2.28]{KuMaVa:02}; for completeness we include the short proof.

\begin{lemma}[Rellich-type identity]
\label{lem:Rellich}
Let $v$ be a real-valued $C^2$ function on some open set $D\subset \mathbb{R}^d$, $d\geq 2$.
Let $\widetilde{\bZ}\in (C^1(D))^d$ and $\alpha\in C^1(D)$
and let both be real-valued.
Then, with the summation convention,
\begin{align}
2 \,{\cal Z} v \Delta v  =  \nabla \cdot \Big[ 2\, {\cal Z} v \nabla v  - |\nabla v|^2 \widetilde{\bZ}\Big] - \big(2 \alpha-\nabla \cdot \widetilde{\bZ}  \big)|\nabla v|^2  - 2  \partial_i \widetilde{Z}_j \partial_i v \partial_j v - 2 v\nabla \alpha\cdot \nabla v,
\label{eq:id1}
\end{align}
where
\begin{equation}\label{cZ}
{\cal Z} v := \big(\widetilde{\bZ}\cdot \nabla v + \alpha v\big).
\end{equation}
\end{lemma}

\begin{proof}
Splitting ${\cal Z} v$ into its component parts, we see that
the identity \eqref{eq:id1} is the sum of the identities
\begin{align}
2 \,\widetilde{\bZ}\cdot\nabla v \Delta v  =  \nabla \cdot \left[ 2\,  (\widetilde{\bZ}\cdot \nabla v)\,\nabla v  - |\nabla v|^2 \widetilde{\bZ}\right]  + \big(\nabla \cdot \widetilde{\bZ} \big)|\nabla v|^2  - 2 \partial_i \widetilde{Z}_j \partial_i v \partial_j v \label{A}
\end{align}
and
\begin{equation}\label{C}
2\alpha v \Delta v  = \nabla \cdot \left[ 2 \alpha v\nabla v \right] - 2\alpha |\nabla v|^2  - 2v\nabla\alpha \cdot \nabla v.
\end{equation}
To prove \eqref{C}, expand the divergence on the right-hand side.
The identity \eqref{A} is obtained by combining the identities
\begin{equation}\label{basic}
\widetilde{\bZ}\cdot \nabla v \Delta v = \nabla \cdot \left[(\widetilde{\bZ}\cdot \nabla v) \nabla v\right] - \partial_i \widetilde{Z}_j \partial_i v \partial_j v - \nabla v \cdot  (\widetilde{\bZ}\cdot \nabla) \nabla v
\end{equation}
and
\begin{equation}\label{Melenktrick}
2\nabla v \cdot \Big( \widetilde{\bZ}\cdot \nabla\Big) \nabla v= \nabla \cdot \left( |\nabla v|^2 \widetilde{\bZ} \right) - (\nabla \cdot \widetilde{\bZ}) |\nabla v|^2,
\end{equation}
which can both be proved by expanding the divergences on the right-hand sides.
\end{proof}

\

For the proof of Lemma \ref{lem:key}, we need the identity \eqref{eq:id1} integrated over a Lipschitz domain when $v$ is the single-layer potential. As a step towards this, we prove the following lemma.

\begin{lemma}[Integrated version of the identity]\label{lem:int} Let $D$ be a Lipschitz domain with outward-pointing unit normal vector $\boldsymbol{\nu}$.
Define
\begin{equation}\label{eq:V}
V(D):=\Big\{v: \,v\in H^1(D), \; \Delta v \in L^2(D),\; \gamma v \in H^1(\partial D),\; \partial_{\nu} v \in L^2(\partial D)
\Big\}.
\end{equation}
If $v\in V(D)$, $\widetilde{\bZ}\in (W^{1,\infty}(D))^d$ (i.e. $\widetilde{Z}_i$ and $\partial_i \widetilde{Z}_j \in L^\infty(D)$ for $i,j=1,\ldots,n $),  $\alpha \in W^{1,\infty}(D)$, and $v$, $\widetilde{\bZ}$, and $\alpha$ are all real-valued,
then
\begin{align}\nonumber
&\int_{\partial D}
\left[(\widetilde{\bZ}\cdot \boldsymbol{\nu}) \left( \left(\partial_{\nu} v\right)^2  - |\nabla_{\Gamma} v|^2 \right) + 2
\Big(\widetilde{\bZ}\cdot \nabla_{\Gamma} (\gamma v) + \alpha (\gamma v)\Big)
\partial_{\nu} v\right]{\rm d} s\\
&\qquad= \int_D \bigg(2 \,{\cal Z} v  \Delta v +2 \, \partial_i \widetilde{Z}_j \partial_i v \partial_j v + 2 \, v \nabla \alpha\cdot \nabla v
+\left(2 \alpha-\nabla \cdot \widetilde{\bZ}  \right)|\nabla v|^2\bigg) {\rm d} \mathbf{x}. \label{eq:int1}
\end{align}
\end{lemma}

\noindent Recall that, when  $D$ is Lipschitz, we can identify $W^{1,\infty}(D)$ with $C^{0,1}(\overline D)$ (see, e.g., \cite[\S4.2.3, Theorem 5]{EvGa:92}), and understand $\widetilde{\bZ}$ and $\alpha$ on $\partial D$ in \eqref{eq:int1} by restriction without needing a trace operator.

\

\begin{proof}[Proof of Lemma \ref{lem:int}]
We first assume that $\widetilde{\bZ}$ and $\alpha$ are as in the statement of the theorem, but $v\in {\cal D}(\overline{D}):=\{ U|_\Omega : U \in C^\infty(\mathbb{R}^d)\}$.
Recall that the divergence theorem $\int_D \nabla \cdot \textbf{F}\,{\rm d} \mathbf{x}
= \int_{{\partial D}} \textbf{F} \cdot \boldsymbol{\nu}\,{\rm d} s$
is valid when $\textbf{F} \in (C^1(\overline{D}))^d$ \cite[Theorem 3.34]{Mc:00}, and thus for $\mathbf{F}\in (H^1(D))^d$ by the density of $C^1(\overline{D})$ in $H^1(D)$ \cite[Theorem 3.29]{Mc:00} and the continuity of trace operator from $H^1(D)$ to $H^{1/2}({\partial D})$ \cite[Theorem 3.37]{Mc:00}.
Recall also that the product of an $H^1(D)$ function and a $W^{1,\infty}(D)$ function is in $H^1(D)$,
and the usual product rule for differentiation holds for such functions. Thus
$\mathbf{F}= 2 {\cal Z} v \nabla v  - |\nabla v|^2 \widetilde{\bZ}$ is in $(H^1(D))^d$
and then \eqref{eq:id1} implies that $\nabla\cdot\mathbf{F}$ is given by the integrand on the left-hand side of \eqref{eq:int1}.
Furthermore,
\begin{align*}
\gamma \mathbf{F}\cdot \boldsymbol{\nu} =(\widetilde{\bZ}\cdot \boldsymbol{\nu}) \left( \left(\pdiff{v}{\nu}\right)^2  - |\nabla_{\Gamma} v|^2 \right) + 2\,
\big(\widetilde{\bZ}\cdot \nabla_{\Gamma} v + \alpha v\big)
\pdiff{v}{\nu}
\end{align*}
on $\partial D$, where we have used the fact that $\nabla v = \boldsymbol{\nu} (\partial v/\partial \nu) + \nabla_{\Gamma} v$ on $\partial D$ for $v\in {\cal D}(\overline{D})$; the identity \eqref{eq:int1} then follows from the divergence theorem.

The result for $v\in V(D)$ then follows from (i) the density of $\cD(\overline{\Omega})}%C_{\text{comp}}^\infty(\overline{\Omega})$ in $V(D)$
\cite[Lemmas 2 and 3]{CoDa:98} and (ii) the fact that \eqref{eq:int1} is continuous in $v$ with respect to the topology of $V(D)$
\end{proof}

\

\begin{proof}[Proof of Lemma \ref{lem:key}]
As discussed in \S\ref{sec:idea}, our strategy is to
mimic the classic method of ``transferring" coercivity properties of the PDE formulation to the BIOs in the trace spaces, but with Green's identity
\begin{equation}\label{eq:Green_int}
-\int_D u \Delta u \, {\rm d} \mathbf{x} = \int_D |\nabla u|^2\,{\rm d} \mathbf{x} - \int_{\partial D} \gamma u\, \partial_n u \, {\rm d} s,
\end{equation}
replaced by the integrated version of the Rellich-type identity \eqref{eq:id1}. That is, we apply the integrated version of \eqref{eq:id1}, namely \eqref{eq:int1}, with $v$ replaced by $u={\cal S}\phi$ (with $\phi \in {L^2(\Gamma)}$),
and $D$ first equal to $\Omega_-$, and then equal to $\Oe \cap B_R$, where $R> \sup_{\mathbf{x}\in \Omega_-}|\mathbf{x}|$.
At this stage we let $\widetilde{\bZ}$ be a general real-valued $W^{1,\infty}(\mathbb{R}^3)$ vector field with compact support, and let $\alpha$ be an arbitrary real constant.
That \eqref{eq:int1} holds with $v$ replaced by $u={\cal S}\phi$, with $\phi$ real-valued, can be  justified by using the results of \cite{JeKe:95} and \cite[Appendix A]{ChGrLaSp:12} recapped in \S\ref{app:HA}.
Indeed,
when $\phi \in {L^2(\Gamma)}$, $u= {\cal S}\phi \in H^{3/2}(D)$ when $D= \Omega_-$ or $\Oe\cap B_R$ by the first mapping property in \eqref{eq:LPmap}; then $u\in V(D)$ by Corollary \ref{cor:CoDa} and \eqref{eq:int1} holds by Lemma \ref{lem:int}.
\footnote{A common alternative to justify that \eqref{eq:int1} holds with $v$ replaced by $u={\cal S}\phi$ is to
(a) approximate $\Omega_-$ by a sequence of subdomains (often they are assumed to be smoother than $\Omega_-$, but this is not necessary),
(b)  apply the identity \eqref{eq:int1} to each member of the sequence,
(c) pass to the limit using the facts that (i) the non-tangential limits of $u$ and $\nabla u$ exist, and (ii) the maximal functions $u^*$ defined by \eqref{eq:ntmax} and $(\nabla u)^*$ (defined analogously) are in ${L^2(\Gamma)}$;
for examples of this argument, see, e.g.,  \cite[Proof of Theorem 2.1]{Ve:84}, \cite[Proof of Theorem 1]{MeCo:00}, \cite[Proof of Proposition 4.2]{Ta:00}, and \cite[Proof of Lemma 4.5]{SpChGrSm:11}.
 }

We have therefore established that \eqref{eq:int1} holds when $D= \Omega_-$ or $\Oe\cap B_R$ and $u= {\cal S}\phi$ for $\phi \in{L^2(\Gamma)}$ that is real-valued.
That is, with the identity \eqref{eq:id1} written as $\nabla \cdot \mathbf{Q} =P$,
\begin{equation}\label{apofmoi}
\int_\Gamma \mathbf{Q}_-\cdot \mathbf{n} \,{\rm d} s = \int_{\Omega_-} P \,{\rm d} \mathbf{x}
\end{equation}
and
\begin{equation}\label{apofmoe}
-\int_\Gamma \mathbf{Q}_+\cdot \mathbf{n} \,{\rm d} s + \int_{{\Gamma_R}} Q_R \,{\rm d} s = \int_{\Oe \cap B_R} P \,{\rm d} \mathbf{x},
\end{equation}
where (remembering that $\Delta u =0$ and $\alpha$ is a constant)
\begin{align}\label{P}
P=  2 \, \partial_i \widetilde{Z}_j \partial_i u \partial_j u + \big(2 \alpha - \nabla \cdot \widetilde{\bZ} \big)|\nabla u|^2,
\end{align}
\begin{equation}\label{Qn}
\mathbf{Q}_{\pm}\cdot \mathbf{n} = (\widetilde{\bZ} \cdot \mathbf{n}) \left( \big(\partial_n^{\pm}u\big)^2 - |\nabla_{\Gamma}(\gamma^\pm u)|^2\right) + 2\, \Big(\widetilde{\bZ}\cdot \nabla_{\Gamma}(\gamma^\pm u) +\alpha \gamma^\pm u\Big)\partial_n^{\pm}u.
\end{equation}
If $R$ is chosen large enough so that $\mathrm{supp} \,\widetilde{\bZ} \subset B_R$, then
\begin{equation}\label{QR}
Q_R = \mathbf{Q}\cdot \widehat{\mathbf{x}}= 2\alpha \,u\pdiff{u}{r}\quad\text{ for } \mathbf{x}\in {\Gamma_R},
\end{equation}
 where we have used the fact that $u$ is $C^\infty$ in a neighbourhood of ${\Gamma_R}$ (either by elliptic regularity or directly by the definition of the single-layer potential \eqref{eq:SLP}) to justify writing $\partial u/\partial r$ in place of some appropriate trace.

Adding \eqref{apofmoi} and \eqref{apofmoe} yields
\begin{equation*}
\int_\Gamma(\mathbf{Q}_- - \mathbf{Q}_+) \cdot \mathbf{n} \,{\rm d} s + \int_{{\Gamma_R}} Q_R \,{\rm d} s = \int_{\Omega_-} P \, {\rm d} \mathbf{x} + \int_{\Oe \cap B_R} P \, {\rm d} \mathbf{x}.
\end{equation*}
Now if $d\geq 3$ and $\phi \in{L^2(\Gamma)}$, then
\begin{equation}\label{eq:Sasym}
|{\cal S}\phi(\mathbf{x})|= O(|\mathbf{x}|^{2-d})\quad\text{ and }\quad |\nabla{\cal S}\phi(\mathbf{x})| = O(|\mathbf{x}|^{1-d})
\end{equation}
as $|\mathbf{x}|\rightarrow \infty$, uniformly in all directions $\mathbf{x}/|\mathbf{x}|$.
If $d=2$ then
\begin{equation}\label{eq:Sasym2}
{\cal S}\phi(\mathbf{x})= \frac{1}{2\pi}\log \left(\frac{a}{|\mathbf{x}|}\right) \left(\phi,1\right)_{{L^2(\Gamma)}} + O(|\mathbf{x}|^{-1})\quad\text{ and }\quad \nabla{\cal S}\phi(\mathbf{x})  =
-\frac{1}{2\pi |\mathbf{x}|} \left(\phi,1\right)_{{L^2(\Gamma)}} + O(|\mathbf{x}|^{-2})
\end{equation}
as $|\mathbf{x}|\rightarrow \infty$, uniformly in all directions $\mathbf{x}/|\mathbf{x}|$;
 these asymptotics are proved for $d=2,3$ in, e.g., \cite[Lemma 6.21]{St:08} (see also \cite[Equations 3.22 and 3.23]{SaSc:11}  for $d=3$); the proof of \eqref{eq:Sasym} for $d\geq 4$ is analogous.
Recalling the definition of $P_\Gamma$ \eqref{eq:P} and the assumption that $P_\Gamma\phi=0$ when $d=2$, we see that,
by \eqref{QR},
$\int_{{\Gamma_R}}Q_R \,{\rm d} s=O(R^{2-d})$ for $d\geq 3$ and
$\int_{{\Gamma_R}}Q_R \,{\rm d} s=O(R^{-2})$ for $d=2$
as $R\rightarrow \infty$. Thus, in this limit,
\begin{equation}\label{eq:keygamma-1}
\int_\Gamma(\mathbf{Q}_- - \mathbf{Q}_+) \cdot \mathbf{n} \,{\rm d} s  = \int_{\Omega_-\cup \Oe} P \, {\rm d} \mathbf{x}.
\end{equation}
The expressions for $\mathbf{Q}_\pm \cdot \mathbf{n}$ \eqref{Qn} and the single-layer potential jump relations
\eqref{eq:jump1}
then imply that
\begin{equation}\label{eq:keygamma}
\int_\Gamma(\mathbf{Q}_- - \mathbf{Q}_+) \cdot \mathbf{n} \,{\rm d} s = 2\Big( (\widetilde{\bZ}\cdot \mathbf{n})D' + \widetilde{\bZ}\cdot \nabla_{\Gamma} S}%\pdiff{S_k}{T}+ \alpha S) \phi, \phi\Big)_{{L^2(\Gamma)}}.
\end{equation}
A key identity to help one see this is
\begin{eqnarray} \nonumber
\left(\partial_n^- u(\mathbf{x})\right)^2-\left(\partial_n^+ u(\mathbf{x})\right)^2 = 2\,\phi(\mathbf{x})\,\big(D^\prime\phi(\mathbf{x})\big) \quad \text{ for a.e. }\mathbf{x}\in\Gamma,\label{jump5}
\end{eqnarray}
which can be established using $a^2 -b^2 = (a-b)(a+b)$ and the jump relations \eqref{eq:jump1} for $\partial^{\pm}_n u$.

Combining \eqref{eq:keygamma-1}, \eqref{eq:keygamma}, and \eqref{eq:P}, we therefore have that
\begin{equation}\label{eq:1}
2 \Big( (\widetilde{\bZ}\cdot \mathbf{n})D' + \widetilde{\bZ}\cdot \nabla_{\Gamma} S}%\pdiff{S_k}{T}+ \alpha S) \phi, \phi\Big)_{{L^2(\Gamma)}} = \int_{\Omega_- \cup \Oe}\Big( 2 \partial_i \widetilde{Z}_j \partial_i u \partial_j u + \big(2 \alpha - \nabla \cdot \widetilde{\bZ} \big)|\nabla u|^2 \Big) {\rm d} \mathbf{x}.
\end{equation}
Using the Cauchy-Schwarz inequality and the definition of the matrix 2-norm for the term involving $2 \, \partial_i \widetilde{Z}_j \partial_i u \partial_j u= 2 \,\nabla u\cdot(D\widetilde{\bZ}\, \nabla u)$, and then standard results about integrals for both this term and the term involving $\nabla\cdot\widetilde{\bZ}$, we find that the right-hand side of \eqref{eq:1} is
\begin{equation*}
\geq \left( 2\alpha - \left( 2\sup_{\mathbf{x} \in \mathbb{R}^d} \big\|D\widetilde{\bZ}\big\|_2 + \big\|\nabla\cdot \widetilde{\bZ}\big\|_{L^\infty(\mathbb{R}^d)}\right)\right)\int_{\Omega_- \cup \Oe} |\nabla u|^2 \, {\rm d} \mathbf{x}.
\end{equation*}
Therefore, choosing $\alpha$ to satisfy the lower bound \eqref{eq:alpha_bound} establishes the lemma with the $+$ sign in \eqref{eq:key}. Multiplying \eqref{eq:1} by $-1$ and letting $\alpha \mapsto -\alpha$ we see again that if $\alpha$ satisfies \eqref{eq:alpha_bound} then this modified right-hand side is $\geq 0$, which establishes the lemma with the $-$ sign  in \eqref{eq:key}.
\end{proof}

\begin{remark}[Recovering the results of Lemma \ref{lem:Calderon2} for $d\geq 3$]
If $d\geq 3$ and $\widetilde{\bZ}=\mathbf{x}$, \eqref{eq:1} becomes
\begin{equation}\label{eq:1a}
2 \Big( (\mathbf{x}\cdot \mathbf{n})D' + \mathbf{x}\cdot \nabla_{\Gamma} S}%\pdiff{S_k}{T}+ \alpha S) \phi, \phi\Big)_{{L^2(\Gamma)}} = (2+ 2\alpha-d)\int_{\Omega_- \cup \Oe}|\nabla u|^2 \, {\rm d} \mathbf{x}.
\end{equation}
This is because, despite the additional terms in the analogue of \eqref{QR} coming from $\widetilde{\bZ}$ no longer having compact support,
it turns out that
$\int_{{\Gamma_R}}Q_R \,{\rm d} s
=O(R^{2-d})$ as $R\rightarrow \infty$ as before.
Letting $\alpha=(d-2)/2$ in \eqref{eq:1a} and recalling the definition \eqref{eq:Caldb'} of $K_{\vfb}'$, we obtain the second equality in \eqref{eq:symmetry3d}.
\end{remark}

\begin{remark}[Link with Verchota's proof of invertibility of $\frac{1}{2} I- D'$ on ${L^2(\Gamma)}$]\label{rem:Verchota}
Verchota's proof that $\frac{1}{2} I-D'$ is invertible on ${L^2(\Gamma)}$ when $\Gamma$ is Lipschitz in \cite[Theorem 3.1]{Ve:84} relies on the inequalities
\begin{equation}\label{eq:Verchota}
\N{\left(\frac{1}{2} I - D'\right)\phi}_{L^2(\Gamma)} \lesssim \N{\left(\frac{1}{2} I + D'\right)\phi}_{L^2(\Gamma)} \lesssim \N{\left(\frac{1}{2} I - D'\right)\phi}_{L^2(\Gamma)},
\end{equation}
which hold for all $\phi\in{L^2(\Gamma)}$ for $d\geq 3$ and for all $\phi \in {L^2(\Gamma)}$ with $P_\Gamma\phi=0$ for $d=2$, and where the omitted constants depend only on the Lipschitz character of $\Omega_-$. (Note that \cite[Theorem 2.1]{Ve:84} proves the slightly weaker result that
\begin{equation*}
\N{\left(\frac{1}{2} I \pm D_0'\right)\phi}_{L^2(\Gamma)} \lesssim \N{\left(\frac{1}{2} I \mp D_0'\right)\phi}_{L^2(\Gamma)} + \left|\int_\Gamma S_0\phi\, {\rm d} s\right|
\end{equation*}
but the final term on the right-hand side can be eliminated; see \cite[Chapter 15, Corollary 1, Page 273]{MeCo:00} when $\Gamma$ is the graph of a function and \cite[Corollary 2.20]{AmKa:07} for $\Omega_-$ bounded.)

The inequalities in \eqref{eq:Verchota} can be obtained by applying the following Dirichlet-to-Neumann and Neumann-to-Dirichlet map bounds with $u={\cal S} \phi$ and using the jump relations \eqref{eq:jump1}.

\noindent (i) If $u\in H^1(\Omega_-)$ is such that $\Delta u=0$ in $\Omega_-$, $\gamma^- u\in H^1(\Gamma)$, and $\partial_n^- u\in{L^2(\Gamma)}$, then
\begin{equation}\label{eq:V1}
\N{\nabla_{\Gamma}(\gamma^- u)}_{L^2(\Gamma)} \lesssim \N{\partial_n^- u}_{{L^2(\Gamma)}} \lesssim \N{\nabla_{\Gamma}(\gamma^- u)}_{L^2(\Gamma)}.
\end{equation}
(ii) If $u\in H^1_{\rm{loc}}(\Oe)$ is such that $\Delta u=0$ in $\Oe$, $\gamma^+ u\in H^1(\Gamma)$, $\partial_n^+ u\in{L^2(\Gamma)}$, and $u(\mathbf{x})= O(|\mathbf{x}|^{2-d})$ for $d\geq 3$ and $u(\mathbf{x})= O(|\mathbf{x}|^{-1})$ for $d=2$, then
\begin{equation}\label{eq:V2}
\N{\nabla_{\Gamma}(\gamma^+ u)}_{L^2(\Gamma)} \lesssim \N{\partial_n^+ u}_{{L^2(\Gamma)}} \lesssim \N{\nabla_{\Gamma}(\gamma^+ u)}_{L^2(\Gamma)}.
\end{equation}

\noindent The link with our proofs of coercivity of our new BIEs comes from the fact that the bounds \eqref{eq:V1} and \eqref{eq:V2} can be proved using the identity \eqref{eq:id1} with $\alpha=0$ and $\widetilde{\bZ}$ the vector field of Lemma \ref{lem:vf}; see, e.g., \cite[Corollary 2.20]{AmKa:07}.
\end{remark}

\subsection{Proofs of Theorems \ref{thm:2d} and \ref{thm:2dstar} (the 2d results)}

\begin{lemma}\label{lem:Neumann_trace}
If $u$ is the solution of the IDP then $P_\Gamma(\partial_n^- u)=0$. If $u$ is the solution of the EDP and $d=2$, then $P_\Gamma(\partial_n^+ u)=0$.
\end{lemma}

\begin{proof}
The result for the IDP follows from applying Green's second identity to $u$ and the constant function.
The result for the EDP when $d=2$ follows in a similar way, using the arguments in the proof of \cite[Theorem 8.9]{Mc:00} to deal with the integral at infinity. Alternatively, the result for the EDP when $d=2$ is proved in \cite[Proof of Theorem 6.10]{Kr:89}; see \cite[Equation 6.10]{Kr:89}.
\end{proof}

\

\begin{proof}[Proof of Theorem \ref{thm:2d}]
For Parts (i) and (iii), arguing exactly as in the proofs of Part (i) of Theorems \ref{thm:Laplace_int} and \ref{thm:Laplace_ext} gives
\begin{equation}\label{eq:England2}
A'_{I,\bZ,\alpha} \partial_n^- u =B_{I,\bZ,\alpha} g_D \quad\text{ and }\quad A'_{E,\bZ,\alpha}\partial_n^+ u =B_{E,\bZ,\alpha} g_D + \alpha u_\infty,
\end{equation}
where $u_\infty$ is the limit of the solution of the EDP at infinity
and we use Green's integral representation $u(\mathbf{x}) = -{\cal S}\partial_n^+ u (\mathbf{x}) + {\cal D} \gamma^+ u (\mathbf{x}) + u_\infty$ for $\mathbf{x} \in \Oe$ and $d=2$.
The BIEs \eqref{eq:2d_IDP_direct} and \eqref{eq:2d_EDP_direct} then follow by applying $Q_\Gamma = I-P_\Gamma$ to the equations in \eqref{eq:England2} and then using that $P_\Gamma \partial^{\pm}_n u=0$  by Lemma \ref{lem:Neumann_trace}, so that $\partial^{\pm}_n u = Q_\Gamma \partial^{\pm}_n u$.

For Part (ii), taking the non-tangential limit of $u$ defined by \eqref{eq:u2d_int} and using the jump relations \eqref{eq:jumpCald} and \eqref{eq:jump1}
(similar to the proof of Part (ii) of Theorem \ref{thm:Laplace_int}) and the fact that $Q_\Gamma= I-P_\Gamma$, we obtain that $\gamma_-u =g_D$ if the BIE \eqref{eq:2d_IDP_indirect} holds. Exactly as in the analogous proof for $d\geq 3$ in \S\ref{sec:4.1}, $\mathcal{K}_{\vfb}\psi$ and ${\cal S}\psi$
with $\psi\in{L^2(\Gamma)}$
 are in $C^2(\Omega_-)$, have non-tangential maximal functions in ${L^2(\Gamma)}$, and satisfy Laplace's equation; therefore $u$ defined by \eqref{eq:u2d_int} inherits these properties.

The proof of Part (iv) is very similar to the proof of Part (ii), except that now need to show that $u$ defined by \eqref{eq:u2d_ext}
satisfies $u(\mathbf{x}) = O(1)$ as $|\mathbf{x}|\rightarrow \infty$; these asymptotics follow from
the first bound in \eqref{eq:Sasym} (since $P_\Gamma Q_\Gamma \phi=0$) and the bound \eqref{eq:Caldasym}.

To see Part (v), arguing as in the proof of Part (vi) of Theorems \ref{thm:Laplace_int} and \ref{thm:Laplace_ext} below Lemma \ref{lem:Lipschitz}, but using \eqref{eq:key} with $d=2$, we see that $(A\psi,\psi)_{{L^2(\Gamma)}}\geq (c/2)\|\psi\|^2_{{L^2(\Gamma)}}$ for all real-valued $\psi\in L^2_0(\Gamma):= \{\phi\in {L^2(\Gamma)}:P_\Gamma\phi=0\}$ if $\alpha$ satisfies \eqref{eq:alpha_bound_boundary}, where $A$ denotes any of  $A_{I,\bZ,\alpha}$, $A'_{I,\bZ,\alpha}$, $A_{E,\bZ,\alpha}$, or $A'_{E,\bZ,\alpha}$.
 Part (v) then follows from the fact that if $(A\psi,\psi)_{{L^2(\Gamma)}}\geq (c/2)\|\psi\|^2_{{L^2(\Gamma)}}$ for all real-valued $\psi\in L^2_0(\Gamma)$ (so that $A$ is coercive on $L^2_0(\Gamma)$ with coercivity constant $c/2$),
then $Q_\Gamma A Q_\Gamma + cP_\Gamma/2$ is coercive on ${L^2(\Gamma)}$ with coercivity constant $c/2$. Indeed, since $Q_\Gamma'=Q_\Gamma$, $P_\Gamma^2=P_\Gamma$,  $P_\Gamma'=P_\Gamma$,
and $P_\GammaQ_\Gamma=0$,
it follows that, for all real-valued $\psi\in{L^2(\Gamma)}$, $Q_\Gamma\psi\in L^2_0(\Gamma)$ and
\begin{align*}
\left(\left( Q_\Gamma A Q_\Gamma + \frac{c}{2}P_\Gamma \right) \psi,\psi\right)_{{L^2(\Gamma)}}
&=\big( A Q_\Gamma \psi,Q_\Gamma\psi\big)_{{L^2(\Gamma)}}  +\frac{c}{2}\big(P_\Gamma^2 \psi,\psi\big)_{{L^2(\Gamma)}} \\
&\geq \frac{c}{2}\N{Q_\Gamma \psi}^2_{{L^2(\Gamma)}} + \frac{c}{2}\N{P_\Gamma\psi}^2_{{L^2(\Gamma)}}
= \frac{c}{2}\N{\psi}^2_{L^2(\Gamma)}.
\end{align*}
\end{proof}

\begin{proof}[Proof of Theorem \ref{thm:2dstar}] For Parts (i) and (iii), taking $\bZ=\mathbf{x}$ and $\alpha=0$ in \eqref{eq:England2} yields
\begin{equation*}
A'_{I,\mathbf{x},0} \partial_n^+ u =B_{I,\mathbf{x},0} \,g_D \quad\text{ and }\quad A'_{E,\mathbf{x},0} \partial_n^+ u =B_{E,\mathbf{x},0} \,g_D.
\end{equation*}
Since $P_\Gamma\partial^{\pm}_n u=0$ by Lemma \ref{lem:Neumann_trace}, the BIEs \eqref{eq:2dstar1} and \eqref{eq:2dstar3} follow.

The proofs of Parts (ii) and (iv) follow in the same way as the proofs of Parts (ii) and (iv) of Theorem \ref{thm:2d}, namely by taking non-tangential limits of $u$, using the jump relations \eqref{eq:jumpCald} and \eqref{eq:jump1}, and using the asymptotics \eqref{eq:Caldasym} for the exterior problem.

Part (v) follows immediately from using the second equation in \eqref{eq:symmetry2d}.
\end{proof}

\subsection{Proof of the results in \S\ref{sec:conditioning} (the conditioning results)}\label{sec:conditioning2}

\subsubsection{Proof of Theorem \ref{thm:conditioning}.}\label{sec:conditioning3}

Theorem \ref{thm:conditioning} is a special case of the following general theorem about GMRES applied to Galerkin linear systems of a continuous and coercive operator on a Hilbert space. We first establish some notation.

As in \S\ref{sec:Gal}, we consider the Galerkin method  applied to the equation $A\phi = f$, where $\phi, f \in {\cal H}$, $A:{\cal H}\rightarrow {\cal H}$ is a continuous (i.e.~bounded) linear operator, and ${\cal H}$ is a Hilbert space over $\mathbb{C}$. Let ${\cal H}_N \subset {\cal H}$ be such that ${\cal H}_N = {\rm span} \{\psi^{N}_1,\ldots, \psi^{N}_{M_N}\}$, with $M_N=\dim({\cal H}_N)$ and $\{\psi^{N}_1,\ldots, \psi^{N}_{M_N}\}$ a basis for ${\cal H}_N$.
The Galerkin matrix of $A$ is then defined by $({\mymatrix{A}})_{ij} :=(A\psi^{N}_j,\psi^{N}_i)_{{\cal H}}$, $i,j=1,\ldots, M_N$ (compare to \eqref{eq:matrix}).

The rest of the set up of \S\ref{sec:conditioning} then holds exactly as stated; i.e., we consider the equation ${\mymatrix{A}}\mathbf{x}= {\bf b}$, let $\mathbf{y}_m$ be the $m$th iterate when the linear system \eqref{eq:pclinear}
is solved using GMRES with zero initial guess, let $\mathbf{r}_m$ denote the corresponding residual, and let $\phi_N^{m}$ be defined by \eqref{eq:phiNm}, so that the Galerkin solution $\phi_N$ is given by \eqref{eq:phiN}.

\begin{theorem}\mythmname{Convergence of GMRES applied to the Galerkin linear system of a continuous and coercive operator}\label{thm:conditioning2}
Suppose that $A:{\cal H}\to{\cal H}$ is coercive (i.e., there exists $C_{\rm coer}>0$ such that \eqref{eq:coer} holds) and
Assumption \ref{ass:D} holds with $\|\cdot\|_{{L^2(\Gamma)}}$ in \eqref{eq:norm_equiv} replaced by $\|\cdot\|_{{\cal H}}$.
With $C_1$ and $C_2$ as in  \eqref{eq:norm_equiv}, let $\beta\in [0,\pi/2)$ be defined such that
\begin{equation}\label{eq:cosbeta_gen}
\cos\beta= \frac{C_{\rm coer}}{\|A\|_{{\cal H}\to{\cal H}}}\left(\frac{C_1}{C_2}\right)\quad\text{ and let } \quad \gamma_\beta:= 2 \sin\left(\frac{\beta}{4-2\beta/\pi}\right)
\end{equation}
Given $\varepsilon>0$, if
\begin{equation}\label{eq:iteration_bound_gen}
m\geq  \left( \log\left(\frac{1}{\gamma_\beta}\right)\right)^{-1}
\left[
\log \left( \frac{12 \|A\|_{{\cal H}\to{\cal H}}}{C_{\rm coer}}\left(\frac{C_2}{C_1}\right)^3\right)+ \log\left(\frac{1}{\varepsilon}\right) \right],
\end{equation}
then \begin{equation}\label{eq:full_discrete_gen}
\frac{
\N{\phi-\phi_N^m}_{{L^2(\Gamma)}}}
{\N{\phi}_{{L^2(\Gamma)}}}
\leq (1+\varepsilon)\frac{\|A^\dag_{\bZ,\alpha}\|_{{\cal H}\to{\cal H}}}{C_{\rm coer}}\left( \min_{\psi\in{\cal H}_N}\frac{\N{\phi-\psi}_{{L^2(\Gamma)}}}{\N{\phi}_{{L^2(\Gamma)}}}\right)+ \varepsilon.
\end{equation}
\end{theorem}

The first step in proving Theorem \ref{thm:conditioning2} is to establish the following relationship between the error $\N{\phi-\phi_N^m}_{{\cal H}}$, the GMRES relative residual $\N{\mathbf{r}_m}_2/\N{\mathbf{r}_0}_2$, and the Galerkin error $\N{\phi-\phi_N}_{{\cal H}}$.

\begin{lemma}\label{lem:conditioning1}
Suppose that $A:{\cal H}\to{\cal H}$ is coercive (i.e., there exists $C_{\rm coer}>0$ such that \eqref{eq:coer} holds) and
Assumption \ref{ass:D} holds with $\|\cdot\|_{{L^2(\Gamma)}}$ in \eqref{eq:norm_equiv} replaced by $\|\cdot\|_{{\cal H}}$.
If $C_1$ and $C_2$ are as in  \eqref{eq:norm_equiv} and $\phi_N^m$ is defined by \eqref{eq:phiNm}, then
\begin{equation}\label{eq:full_discrete1}
\frac{
\N{\phi-\phi_N^m}_{{\cal H}}}
{\N{\phi}_{{\cal H}}}
\leq
\left( 1 + \frac{\|A\|_{{\cal H}\to{\cal H}}}{C_{\rm coer}}\left(\frac{C_2}{C_1}\right)^3 \frac{\N{\mathbf{r}_m}_2}{\N{\mathbf{r}_0}_2}\right)
\frac{\N{\phi-\phi_N}_{{\cal H}}}{\N{\phi}_{{\cal H}}} +
\frac{\|A\|_{{\cal H}\to{\cal H}}}{C_{\rm coer}}\left(\frac{C_2}{C_1}\right)^3 \frac{\N{\mathbf{r}_m}_2}{\N{\mathbf{r}_0}_2}.
\end{equation}
\end{lemma}

\

The right-hand side of \eqref{eq:full_discrete1} contains the relative residual
$\N{\mathbf{r}_m}_2/\N{\mathbf{r}_0}_2$. The following bound, from \cite{BeGoTy:06}, gives sufficient conditions on $m$ for this relative residual
 to be controllably small; recall that this bound is an improvement of the so-called ``Elman estimate" from  \cite{El:82, EiElSc:83}.

\begin{theorem}\textbf{\emph{(Elman-type estimate for GMRES from \cite{BeGoTy:06})}} \label{thm:Elman}
Let ${\mymatrix{C}}$ be an $M_N\times M_N$ matrix with $0\notin W({\mymatrix{C}})$, where
\begin{equation*}
W({\mymatrix{C}}):= \big\{ \langle {\mymatrix{C}} \mathbf{v}, \mathbf{v}\rangle : \mathbf{v} \in \mathbb{C}^{M_N}, \|\mathbf{v}\|_2=1\big\}
\end{equation*}
is the \emph{field of values}, also called the \emph{numerical range}, of ${\mymatrix{C}}$.
Let $\beta\in [0,\pi/2)$ be such that
\begin{equation}\label{eq:cosbeta2}
\cos \beta \leq \frac{\mathrm{dist}\big(0, W({\mymatrix{C}})\big)}{\| {\mymatrix{C}}\|_2}
\end{equation}
(observe that $\cos\beta$ is indeed $\leq 1$ by the definition of $W({\mymatrix{C}})$)
and, given $\beta$, let
\begin{equation*}
 \gamma_\beta:= 2 \sin \left( \frac{\beta}{4-2\beta/\pi}\right).
\end{equation*}
Let  $\mathbf{r}_m$  be the $m$th GMRES residual, as defined in \S\ref{sec:conditioning}.
Then
\begin{equation}\label{eq:Elman2}
\frac{\|\mathbf{r}_m\|_{2}}{\|\mathbf{r}_0\|_{2}} \leq \left(2 + \frac{2}{\sqrt{3}}\right)\big(2+ \gamma_\beta\big) \,(\gamma_\beta)^m \leq 12 (\gamma_\beta)^m.
\end{equation}
\end{theorem}

\

\begin{proof}[Proof of Lemma \ref{lem:conditioning1}]
We first use continuity and coercivity of $A$ to obtain bounds on the norm of ${\mymatrix{D}}_N^{-1/2}{\mymatrix{A}} {\mymatrix{D}}_N^{-1/2}$ and its inverse.
By the definition \eqref{eq:matrix},
\begin{equation*}
\big({\mymatrix{A}} \mathbf{v},\mathbf{w}\big)_2 = \big(A v_N,w_N\big)_{{\cal H}} \quad\text{ for all } v_N, w_N \in {\cal H}_N.
\end{equation*}
Using this, along with the norm equivalence \eqref{eq:norm_equiv}, we find that, for all $\mathbf{v}, \mathbf{w}\in \mathbb{C}^{M_N}$,
\begin{equation*}
\big|\big({\mymatrix{A}} \mathbf{v},\mathbf{w}\big)_2\big| \leq \N{A}_{{\cal H}\to{\cal H}} \N{v_N}_{{\cal H}} \N{w_N}_{{\cal H}}
\leq  \N{A}_{{\cal H}\to{\cal H}}(C_2)^2  \big\|{\mymatrix{D}}_N^{1/2} \mathbf{v}\big\|_2 \big\|{\mymatrix{D}}_N^{1/2} \mathbf{w}\big\|_2
\end{equation*}
and
\begin{equation*}
\big|\big({\mymatrix{A}} \mathbf{v},\mathbf{v}\big)_2\big| \geq C_{\rm coer} \N{v_N}_{{\cal H}}^2 \geq C_{\rm coer}(C_1)^2 \big\|{\mymatrix{D}}_N^{1/2}\mathbf{v}\big\|_2^2.
\end{equation*}
Letting $\widetilde{\mathbf{v}}= {\mymatrix{D}}_N^{1/2}\mathbf{v}$ and $\widetilde{\mathbf{w}}= {\mymatrix{D}}_N^{1/2}\mathbf{w}$, we therefore have that, for all $\widetilde{\mathbf{v}}, \widetilde{\mathbf{w}}\in \mathbb{C}^{M_N}$,
\begin{equation}\label{eq:GMRES1}
\big| \big( {\mymatrix{D}}_N^{-1/2}{\mymatrix{A}}{\mymatrix{D}}_N^{-1/2}\widetilde{\mathbf{v}},\widetilde{\mathbf{w}}\big)_2 \big| \leq \N{A}_{{\cal H}\to{\cal H}}(C_2)^2  \N{\widetilde{\mathbf{v}}}_2 \N{\widetilde{\mathbf{w}}}_2
\end{equation}
and
\begin{equation}\label{eq:GMRES2}
\big| \big( {\mymatrix{D}}_N^{-1/2}{\mymatrix{A}}{\mymatrix{D}}_N^{-1/2}\widetilde{\mathbf{v}},\widetilde{\mathbf{v}}\big)_2 \big| \geq C_{\rm coer}  (C_1)^2 \N{\widetilde{\mathbf{v}}}_2^2.
\end{equation}
The inequalities \eqref{eq:GMRES1} and \eqref{eq:GMRES2} then imply that
\begin{equation}\label{eq:GMRES3}
\big\|{\mymatrix{D}}_N^{-1/2}{\mymatrix{A}} {\mymatrix{D}}_N^{-1/2}\big\|_2 \leq \N{A}_{{\cal H}\to{\cal H}} (C_2)^2 \quad\text{ and }\quad \mathrm{dist}\Big(0, W({\mymatrix{D}}_N^{-1/2}{\mymatrix{A}} {\mymatrix{D}}_N^{-1/2})\Big) \geq C_{\rm coer}(C_1)^2,
\end{equation}
with the second inequality and the Lax--Milgram theorem then implying that
\begin{equation}\label{eq:GMRES3a}
\big\|\big({\mymatrix{D}}_N^{-1/2}{\mymatrix{A}} {\mymatrix{D}}_N^{-1/2}\big)^{-1}\big\|_2 \leq \frac{1}{C_{\rm coer} (C_1)^2}.
\end{equation}

We now prove \eqref{eq:full_discrete1}. By the definitions of $\mathbf{y}_m$ (see \eqref{eq:pclinear}) and $\mathbf{r}_m$,
\begin{equation*}
\mathbf{r}_m ={\mymatrix{D}}_N^{-1/2}{\mymatrix{A}} {\mymatrix{D}}_N^{-1/2}(\mathbf{y}_m-\mathbf{y})
\end{equation*}
and (since $\mathbf{y}_0={\bf 0}$) $\mathbf{r}_0 = -{\mymatrix{D}}_N^{-1/2}\mathbf{b}= -{\mymatrix{D}}_N^{-1/2}{\mymatrix{A}} {\mymatrix{D}}_N^{-1/2}\mathbf{y}$. Therefore, by  \eqref{eq:GMRES3a} and the first bound in \eqref{eq:GMRES3},
\begin{align}\nonumber
\N{\mathbf{y}_m-\mathbf{y}}_2 \leq \big\|\big({\mymatrix{D}}_N^{-1/2}{\mymatrix{A}} {\mymatrix{D}}_N^{-1/2}\big)^{-1}\big\|_2 \N{\mathbf{r}_m}_2
&\leq \frac{1}{C_{\rm coer}(C_1)^2}
\left(\frac{\N{\mathbf{r}_m}_2}{\N{\mathbf{r}_0}_2}\right)
\N{\mathbf{r}_0}_2 \\ \nonumber
&\leq \frac{1}{C_{\rm coer}(C_1)^2} \left(\frac{\N{\mathbf{r}_m}_2}{\N{\mathbf{r}_0}_2}\right)\big\|{\mymatrix{D}}_N^{-1/2}{\mymatrix{A}} {\mymatrix{D}}_N^{-1/2}\big\|_2 \N{\mathbf{y}}_2\\
&\leq \frac{\N{A}_{{\cal H}\to{\cal H}}}{C_{\rm coer}}\left(\frac{C_2}{C_1}\right)^2 \left(\frac{\N{\mathbf{r}_m}_2}{\N{\mathbf{r}_0}_2}\right)\N{\mathbf{y}}_2.
\label{eq:GMRES4}
\end{align}
Next, the definition of $\phi_N^m$ \eqref{eq:phiNm}, the expression for $\phi_N$ \eqref{eq:phiN}, and the norm equivalence \eqref{eq:norm_equiv} imply that
\begin{equation*}
\frac{
\N{\phi_N^m-\phi_N}_{{\cal H}}
}{
\N{\phi_N}_{{\cal H}}
}
\leq
 \frac{C_2}{C_1}
\frac{
\big\|{\mymatrix{D}}_N^{1/2}({\mymatrix{D}}_N^{-1/2}(\mathbf{y}_m-\mathbf{y}))\big\|_2
}{
\big\|{\mymatrix{D}}_N^{1/2}{\mymatrix{D}}_N^{-1/2}\mathbf{y}\big\|_2
}
=
 \frac{C_2}{C_1}
\frac{
\big\|\mathbf{y}_m-\mathbf{y}\big\|_2
}{
\big\|\mathbf{y}\big\|_2
},
\end{equation*}
and then combining this with \eqref{eq:GMRES4} we obtain
\begin{equation*}
\frac{
\N{\phi_N^m-\phi_N}_{{\cal H}}
}{
\N{\phi_N}_{{\cal H}}
}
\leq \frac{\N{A}_{{\cal H}\to{\cal H}}}{C_{\rm coer}}\left(\frac{C_2}{C_1}\right)^3 \left(\frac{\N{\mathbf{r}_m}_2}{\N{\mathbf{r}_0}_2}\right).
\end{equation*}
Combining this last inequality with the triangle inequality, we obtain that
\begin{align*}
\N{\phi-\phi_N^m}_{{\cal H}} &\leq \N{\phi-\phi_N}_{{\cal H}} + \N{\phi_N-\phi_N^m}_{{\cal H}},\\
& \leq \N{\phi-\phi_N}_{{\cal H}}+  \frac{\N{A}_{{\cal H}\to{\cal H}}}{C_{\rm coer}}\left(\frac{C_2}{C_1}\right)^3 \left(\frac{\N{\mathbf{r}_m}_2}{\N{\mathbf{r}_0}_2}\right)\N{\phi_N}_{{\cal H}},
\end{align*}
and then the result \eqref{eq:full_discrete1} follows by another use of the triangle inequality.
\end{proof}

\

\begin{proof}[Proof of Theorem \ref{thm:conditioning2}]
By Part (c) of Theorem \ref{thm:Galerkin}, the Galerkin error $\N{\phi-\phi_N}_{{\cal H}}$ satisfies the quasioptimal error estimate \eqref{eq:quasioptimal}.
The definition of $\beta$ in \eqref{eq:cosbeta_gen} and the bounds on ${\mymatrix{D}}_N^{-1/2}{\mymatrix{A}} {\mymatrix{D}}_N^{-1/2}$ in \eqref{eq:GMRES3} imply that \eqref{eq:cosbeta2} is satisfied with ${\mymatrix{C}}= {\mymatrix{D}}_N^{-1/2}{\mymatrix{A}} {\mymatrix{D}}_N^{-1/2}$; note that here it is important that ${\cal H}$ is a Hilbert space over $\mathbb{C}$, so that continuity and coercivity of $A$ control $W({\mymatrix{A}})$ (which involves ${\mymatrix{A}}$ applied to vectors in $\mathbb{C}^{M_N}$).

Using both \eqref{eq:quasioptimal} and
the relative-residual bound \eqref{eq:Elman2}
in \eqref{eq:full_discrete1}, we obtain that
\begin{align*}
\frac{
\N{\phi-\phi_N^m}_{{\cal H}}}
{\N{\phi}_{{\cal H}}}
&\leq
\left( 1 + \frac{12\N{A}_{{\cal H}\to{\cal H}}}{C_{\rm coer}}\left(\frac{C_2}{C_1}\right)^3
(\gamma_\beta)^m\right)
\frac{\|A\|_{{\cal H}\to{\cal H}}}{C_{\rm coer}}\
\min_{\psi\in {\cal H}_N}\frac{\N{\phi-\psi}_{{\cal H}}}{\N{\phi}_{{\cal H}}}
\\
&\qquad\qquad
+
\frac{12\N{A}_{{\cal H}\to{\cal H}}}{C_{\rm coer}}\left(\frac{C_2}{C_1}\right)^3 (\gamma_\beta)^m.
\end{align*}
Given $\varepsilon>0$, if $m$ satisfies
\eqref{eq:iteration_bound_gen}, then
\begin{equation*}
\frac{12\N{A}_{{\cal H}\to{\cal H}}}{C_{\rm coer}}\left(\frac{C_2}{C_1}\right)^3 (\gamma_\beta)^m \leq \varepsilon
\end{equation*}
and thus the bound \eqref{eq:full_discrete_gen} holds.
\end{proof}

\subsubsection{Conditions under which Assumption \ref{ass:D} holds.}\label{sec:conditioning2a}

Our result about the convergence of GMRES applied to the Galerkin matrices of the new formulations, namely Theorem \ref{thm:conditioning}, is proved under Assumption \ref{ass:D}, which is an assumption about the sequence of finite-dimensional subspaces $({\cal H}_N)_{N=1}^\infty$ and their associated bases. Recall from \S\ref{sec:conditioning} that Assumption \ref{ass:D} holds, indeed with ${\mymatrix{D}}_N$ the identity matrix, for any sequence $({\cal H}_N)_{N=1}^\infty$ (and in any dimension $d\geq 2$) provided that the bases we choose are orthonormal. But many standard implementations of boundary element approximation methods use non-orthogonal bases, particularly bases of so-called nodal basis functions (e.g., \cite{AiMcTr:99,GrMc:06}, \cite[Page 216]{St:08}, \cite[Pages 205 and 280]{SaSc:11}.
We show as Lemma \ref{lem:meshes} below that Assumption \ref{ass:D} holds (moreover specifying the diagonal matrices ${\mymatrix{D}}_N$) under mild constraints on the sequence of meshes  when the approximation space allows discontinuities across elements. In particular, Lemma \ref{lem:meshes} holds when nodal basis functions are used, including for  sequences of highly anisotropic meshes.

To specify the conditions under which Assumption \ref{ass:D} holds, we recall the notion of a surface mesh on $\Gamma$, and aspects of the standard implementation of boundary element methods, including the notation of a reference element (for the moment, until we indicate otherwise, our results hold for any dimension $d\geq 2$). Following, e.g., \cite[Defn.~4.1.2]{SaSc:11}, we call $\mathcal{G}$ a {\em mesh} of $\Gamma$ if $\mathcal{G}$ is a set of finitely many disjoint, relatively open, topologically regular\footnote{By {\em topologically regular} we mean that the relative interior of the closure of $\tau$ is again $\tau$.} subsets of $\Gamma$ that cover $\Gamma$ in the sense that $\Gamma= \cup_{\tau\in \mathcal{G}} \overline{\tau}$, and are such that the relative boundary of each $\tau\in \mathcal{G}$ has zero surface measure. We call the elements of $\mathcal{G}$ the {\em (boundary) elements} of the mesh and, for $\tau\in \mathcal{G}$, set
$$
h_\tau:= \mathop{{\rm diam}}(\tau) \quad \mbox{and} \quad  s_\tau := |\tau|,
$$
where $|\tau|$ denotes the $(d-1)$-dimensional surface measure of $\tau$, and set $h:= \max_{\tau\in \mathcal{G}} h_\tau$.

We assume moreover that, for each $\tau\in \mathcal{G}$, there exists a mapping $\chi_\tau:\widehat{\tau}\to \tau$, for some $\widehat{\tau}\in \mathcal{R}$, the finite set of {\em reference elements}, that is bijective and at least bi-Lipschitz, so that $\chi_\tau^{-1}:\tau\to \widehat{\tau}$ is well-defined and also bi-Lipschitz.\footnote{It is standard (e.g., \cite{AiMcTr:99,GrMc:06}) to assume more smoothness for $\chi_\tau$, e.g.~that $\chi_\tau \in C^r(\overline{\widehat{\tau}})$ for some $r\in \mathbb{N}$, in which case also $\chi_\tau^{-1}\in C^r(\overline{\tau})$. In the important case when $d=3$, $\Gamma$ is a polyhedron, and $\mathcal{R}=\{\widehat{\tau}\}$, with $\widehat{\tau}$ the unit simplex (a triangle), it is usual (e.g., \cite[Chap.~10]{St:08}, \cite[Defn.~4.1.2]{SaSc:11}) for each $\tau$ to be a triangle and for $\chi_\tau$ to be affine.} Here, by a {\em reference element}, $\widehat{\tau}$, we mean, generically, some bounded, open, topologically regular subset of $\mathbb{R}^{d-1}$, but with the idea that, in practical implementations, $\widehat{\tau}$ is a polyhedron, usually the unit cube $\widehat{\tau} = (0,1)^{d-1}$ or the unit simplex $\widehat{\tau} = \{\widehat{\mathbf{x}}\in (0,1)^{d-1}:\widehat{\mathbf{x}}_1+\ldots + \widehat{\mathbf{x}}_{d-1} < 1\}$.
(In the case $d=2$ it is usual to take $\mathcal{R}=\{\hat\tau\}$ with $\widehat{\tau} = (0,1)$.) For each $\tau\in \mathcal{G}$ let $J_\tau\in (L^\infty(\Gamma))^{d\times (d-1)}$ denote the Jacobian of $\chi_\tau$. Importantly, for every $f\in L^1(\tau)$, where $\widehat{\tau} \in \mathcal{R}$ is the domain of $\chi_\tau$,
\begin{equation} \label{eq:cov}
\int_\tau f(\mathbf{x}) \,{\rm d} s(\mathbf{x}) = \int_{\widehat{\tau}} f(\chi_\tau(\widehat{\mathbf{x}}))g_\tau(\widehat{\mathbf{x}})\, {\rm d} \widehat{\mathbf{x}}, \quad \mbox{where} \quad g_\tau:= \big(\det(J_\tau^TJ_\tau)\big)^{1/2}\in L^\infty(\Gamma);
\end{equation}
in particular $s_\tau = \int_{\widehat{\tau}} g_\tau(\widehat{\mathbf{x}})\, {\rm d} \widehat{\mathbf{x}}$, so that
\begin{equation} \label{eq:gpm}
g^-_\tau |\widehat{\tau}| \leq s_\tau \leq g^+_\tau |\widehat{\tau}|, \quad \mbox{where} \quad g^+_\tau:= \mathop{{\rm ess} \sup}_{\widehat{\mathbf{x}}\in \widehat{\tau}} g_\tau(\widehat{\mathbf{x}}), \quad g^-_\tau:= \mathop{{\rm ess} \inf}_{\widehat{\mathbf{x}}\in \widehat{\tau}} g_\tau(\widehat{\mathbf{x}}).
\end{equation}

For $p\in \mathbb{N}_0$ and $\widehat \tau\in \mathcal{R}$ let $\mathbb{P}_p^{\widehat \tau}$ denote some finite-dimensional set of polynomials $\psi:\widehat{\tau}\to\mathbb{R}$ that
contains the polynomials of (total) degree $\leq p$.
When $\widehat{\tau}$ is a simplex one usually takes  $\mathbb{P}_p^{\widehat \tau}$ to be the set of polynomials of total degree $\leq p$; when $\widehat{\tau}$ is a cube one usually takes $\mathbb{P}_p^{\widehat \tau}$ to be the set of polynomials of coordinate degree $\leq p$; see, e.g., \cite[Page 1494, penultimate displayed equation]{GrMc:06}.
Following \cite[Defn.~4.1.17]{SaSc:11}, given a mesh $\mathcal{G}$ on $\Gamma$, define the {\em boundary element approximation space, $S_\mathcal{G}^p$,} of discontinuous piecewise polynomials of degree $\leq p$ on $\mathcal{G}$, by
\begin{equation} \label{eq:SpG}
S^p_\mathcal{G}:= \big\{\psi\in L^\infty(\Gamma):\psi|_\tau \circ \chi_\tau \in \mathbb{P}_p^{\widehat \tau}, \mbox{ for all } \tau\in \mathcal{G}, \mbox{ where }\widehat \tau \mbox{ is the domain of }\chi_\tau\big\}.
\end{equation}
Where $P_{\widehat \tau}=\dim(\mathbb{P}_p^{\widehat \tau})$ and $M=\dim(S^p_{\mathcal{G}})$, we equip $S^p_\mathcal{G}$ with a basis $\{\psi_1,\ldots, \psi_M\}$ constructed as follows.  For each $\widehat{\tau}\in \mathcal{R}$ choose a basis $\{\psi^{\widehat{\tau}}_1, \ldots,\psi^{\widehat{\tau}}_{P_{\widehat \tau}}\}$ for $\mathbb{P}_p^{\widehat \tau}$ (for example, a nodal basis as in \cite{AiMcTr:99,GrMc:06}). For each $\tau\in \mathcal{G}$, where $\widehat{\tau}$ is the domain of $\chi_\tau$, define $\psi^{\tau}_j\in L^\infty(\Gamma)$, for $j=1,\ldots,P_{\widehat \tau}$, by
\begin{equation}\label{eq:psij}
\psi^\tau_j(\mathbf{x}) :=
\begin{cases}
\psi^{\widehat{\tau}}_j\big(\chi_\tau^{-1}(\mathbf{x})\big), & \mathbf{x}\in \tau,\\
  0, & \mathbf{x}\in \Gamma\setminus \tau.
  \end{cases}
\end{equation}
Then set
\begin{equation} \label{eq:GlBa}
\big\{\psi_1,\ldots, \psi_M\big\} = \big\{\psi^\tau_j:\tau\in \mathcal{G}, j\in \{1,\ldots,P_{\widehat \tau}\}\big\},
\end{equation}
noting that (see, e.g., \cite[p.~1495]{GrMc:06}) $\{\psi_1,\ldots, \psi_M\}$ is a nodal basis if each $\{\psi^{\widehat{\tau}}_1, \ldots,\psi^{\widehat{\tau}}_{P_{\widehat \tau}}\}$, $\widehat{\tau}\in \mathcal{R}$, is a nodal basis.

Consider now the case that we keep $\mathcal{R}$, $p$, and the bases $\{\psi^{\widehat{\tau}}_1, \ldots,\psi^{\widehat{\tau}}_{P_{\widehat \tau}}\}$, $\widehat{\tau}\in \mathcal{R}$, fixed but use a sequence of meshes $\mathcal{G}_N$, $N\in \mathbb{N}$, with associated approximation spaces ${\cal H}_N:= S^p_\mathcal{G_N}$ that are such that $h_N:= \max_{\tau\in \mathcal{G}_N} h_\tau\to 0$ as $N\to\infty$, i.e.~we consider the \emph{$h$-version} of the boundary-element method. Lemma \ref{lem:meshes} below applies in this regime under the following assumption on the constants  $g^\pm_\tau$ defined by \eqref{eq:gpm}
(this assumption is the first half, Equation 3.5a, of \cite[Assumption 3.1]{GrMc:06}).

\begin{assumption} \label{ass:mild} There exists a constant $c_1\geq 1$ such that, for every $N\in \mathbb{N}$ and $\tau\in \mathcal{G}_N$,
\begin{equation} \label{eq:mild}
g^+_\tau \leq c_1 g^-_\tau;
\end{equation}
equivalently, there exists a constant $c_2\geq 1$ such that, for every $N\in \mathbb{N}$ and $\tau\in \mathcal{G}_N$,
\begin{equation} \label{eq:mild2}
c_2^{-1}s_\tau \leq g_\tau(\widehat{\mathbf{x}}) \leq c_2 s_\tau \quad \mbox{for almost all } \widehat{\mathbf{x}}\in \hat\tau.
\end{equation}
\end{assumption}
We make two remarks about Assumption \ref{ass:mild}.

(i) The claimed equivalence of \eqref{eq:mild} and \eqref{eq:mild2}
follows from \eqref{eq:gpm} (precisely, if \eqref{eq:mild} holds then \eqref{eq:mild2} holds with $c_2= c_1\max(|\widehat{\tau}|, |\widehat{\tau}|^{-1})$, and if \eqref{eq:mild2} holds then \eqref{eq:mild} holds with $c_1= c_2^2$).

(ii) Because $\chi_\tau$ is bi-Lipschitz, \eqref{eq:mild} holds for every $\tau\in \mathcal{G}_N$ for some $c_1\geq 1$ (not necessarily independent of $\tau$ and $N$). In particular \eqref{eq:mild} holds with $c_1 =1$  if  each $\chi_\tau$ is affine, so that Assumption \ref{ass:mild} holds in that case (see also the discussion below \cite[Assumption 3.1]{GrMc:06}).

\

For the following lemma, recall that the matrix ${\mymatrix{A}}$ is given by \eqref{eq:matrix} with $A^\dag_{\bZ,\alpha}$ equal to one of $A_{I,\bZ,\alpha}$, $A'_{I,\bZ,\alpha}$, $A_{E,\bZ,\alpha}$, or $A'_{E,\bZ,\alpha}$.

\begin{lemma}[Conditions under which Assumption \ref{ass:D} holds]\label{lem:meshes}
Suppose that, while keeping $\mathcal{R}$, $p$, and the bases $\{\psi^{\widehat{\tau}}_1, \ldots,\psi^{\widehat{\tau}}_{P_{\widehat \tau}}\}$, $\widehat{\tau}\in \mathcal{R}$, fixed, we use a sequence of meshes $\mathcal{G}_N$, $N\in \mathbb{N}$, with associated approximation spaces ${\cal H}_N:= S^p_\mathcal{G_N}$ and bases \eqref{eq:GlBa} that are such that $h_N:= \max_{\tau\in \mathcal{G}_N} h_\tau\to 0$ as $N\to\infty$ and Assumption \ref{ass:mild} holds. Then the following is true.

(i) Assumption \ref{ass:D} holds with $({\mymatrix{D}}_N)_{ii} := s_{i}$, $i=1,\ldots, M_N$, where $s_i:=s_\tau$ if $\psi_i$ is supported in $\tau$, with $C_1:=c_2^{-1/2}c_{\mathcal{R}}^{-1/2}$, $C_2 := c_2^{1/2}c_{\mathcal{R}}^{1/2}$, where $c_{\mathcal{R}}\geq 1$ depends only on the bases  $\{\psi^{\widehat{\tau}}_1, \ldots,\psi^{\widehat{\tau}}_{P_{\widehat{\tau}}}\}$, $\widehat{\tau}\in \mathcal{R}$.

(ii) If, in addition, $d\geq 3$, \eqref{eq:c} holds for some $c>0$, and $\alpha$ satisfies \eqref{eq:alpha_bound_boundary}, then Assumption \ref{ass:D} holds also with $({\mymatrix{D}}_N)_{ii} := |{\mymatrix{A}}_{ii}|$, for $i=1,\ldots, M_N$, with
\begin{equation} \label{eq:C1C2}
C_1 := C_-C^{-1/2} c_2^{-1}\quad\text{ and }\quad C_2 := C_+c^{-1/2}c_2,
\end{equation}
where $C:=\|A^\dag_{\bZ,\alpha}\|_{L^2(\Gamma)\to L^2(\Gamma)}$ and $C_\pm>0$ depend only on the bases  $\{\psi^{\widehat{\tau}}_1, \ldots,\psi^{\widehat{\tau}}_{P_{\widehat \tau}}\}$, $\widehat{\tau}\in \mathcal{R}$.
\end{lemma}

Part (ii) of Lemma \ref{lem:meshes} is proved for $d=3$ using the coercivity results of Part (vi) of both Theorems \ref{thm:Laplace_int} and \ref{thm:Laplace_ext}. An analogous result holds for $d=2$ using the coercivity results of Part (v) of Theorem \ref{thm:2d}, but we omit this for brevity.

\

\begin{proof}[Proof of Lemma \ref{lem:meshes}]
(i) Given $w_N = \sum_{\ell=1}^{M_N}w^N_{\ell} \psi_{\ell}^N \in {\cal H}_N$,
by \eqref{eq:GlBa},
\begin{equation}\label{eq:expand1}
w_N= \sum_{\tau\in \mathcal{G}_N}\sum_{j=1}^{P_{\widehat \tau}} w_{j}^\tau \psi_{j}^\tau
\end{equation}
for some coefficients $\{w_j^\tau\}_{j=1}^{P_{\widehat \tau}}$, and where $\widehat{\tau}$ is the domain of $\chi_\tau$. Thus, by \eqref{eq:cov}, for all $\tau \in \mathcal{G}_N$,
\begin{equation}\label{eq:L2wN}
\big\|w_N|_\tau\big\|^2_{L^2(\tau)} = \int_{\widehat{\tau}} |\widehat{\psi}(\widehat{\mathbf{x}})|^2 g_\tau(\widehat{\mathbf{x}})\, {\rm d} \widehat{\mathbf{x}},
\end{equation}
where
\begin{equation}\label{eq:psihat}
\widehat{\psi}(\widehat{\mathbf{x}}) := w_N\big(\chi_\tau(\widehat{\mathbf{x}})\big) = \sum_{j=1}^{P_{\widehat \tau}} w_j^\tau \psi^{\widehat{\tau}}_j(\widehat{\mathbf{x}}) \quad\text{ for }\widehat{\mathbf{x}}\in \widehat{\tau}
\end{equation}
(and we have used \eqref{eq:expand1} and \eqref{eq:psij}).
For every $\widehat{\tau}\in \mathcal{R}$, every $\widehat{\phi}\in \mathbb{P}_p^{\widehat \tau}$ can be written as $\widehat{\phi}= \sum_{j=1}^{P_{\widehat \tau}} a_j\psi^{\widehat{\tau}}_j$ for some unique vector $\mathbf{a}=(a_1,\ldots,a_{P_{\widehat \tau}})^T$. With $\mathbb{P}_p^{\widehat \tau}$ equipped with the norm $\|\cdot\|_{\widehat{\tau}}$ defined by
\begin{equation}\label{eq:normtau}
\|\widehat{\phi}\|^2_{\widehat{\tau}} = \sum_{j=1}^{P_{\widehat \tau}} a_j^2,
\end{equation}
since $\mathbb{P}_p^{\widehat \tau}$ is finite-dimensional,   there exists $c_{\widehat{\tau}}\geq 1$ such that
\begin{equation}\label{eq:ctau}
c_{\widehat{\tau}}^{-1}\big\|\widehat{\phi}\big\|^2_{\widehat{\tau}}\leq \big\|\widehat{\phi}\big\|^2_{L^2(\widehat{\tau})}\leq c_{\widehat{\tau}}\big\|\widehat{\phi}\big\|^2_{\widehat{\tau}} \quad\text{ for all } \widehat{\phi}\in \mathbb{P}_p^{\widehat \tau}.
\end{equation}
Therefore, using \eqref{eq:mild2} and \eqref{eq:ctau} (with $\widehat{\phi}=\widehat{\psi}$) in \eqref{eq:L2wN}, we have that, for all $\tau \in \mathcal{G}_N$,
\begin{equation}\label{eq:temp4}
c^{-1}_{\widehat{\tau}} c^{-1}_2s_\tau\big\|\widehat{\psi}\big\|^2_{\widehat{\tau}}\leq c_2^{-1}s_\tau\big\|\widehat{\psi}\big\|^2_{L^2(\widehat{\tau})} \leq \big\|w_N|_\tau\big\|^2_{L^2(\tau)} \leq c_2s_\tau\big\|\widehat{\psi}\big\|^2_{L^2(\widehat{\tau})} \leq c_2s_\tau c_{\widehat{\tau}}\big\|\widehat{\psi}\big\|^2_{\widehat{\tau}}.
\end{equation}
Furthermore, by \eqref{eq:psihat} and \eqref{eq:normtau},
\begin{equation}\label{eq:temp5}
\|\widehat{\psi}\|^2_{\widehat{\tau}} = \sum_{j=1}^{P_{\widehat \tau}} (w_j^\tau)^2.
\end{equation}
If $({\mymatrix{D}}_N)_{ii} := s_{i}$, $i=1,\ldots, M_N$, then, by \eqref{eq:expand1},
\begin{equation}\label{eq:temp6}
\big\|{\mymatrix{D}}_N^{1/2} \mathbf{w}\big\|^2_{2} = \sum_{\tau\in \mathcal{G}_N}s_\tau \sum_{j=1}^{P_{\widehat \tau}} (w_j^\tau)^2.
\end{equation}
Therefore, combining \eqref{eq:temp4}, \eqref{eq:temp5}, and \eqref{eq:temp6}, we see that, with the choice $({\mymatrix{D}}_N)_{ii} := s_{i}$, $i=1,\ldots, M_N$, Assumption \ref{ass:D} holds with $C_1=c^{-1/2}_2c_{\mathcal{R}}^{-1/2}$ and $C_2=c_2^{1/2}c^{1/2}_{\mathcal{R}}$, where $c_{\mathcal{R}}:= \max_{\widehat{\tau} \in \mathcal{R}} c_{\widehat{\tau}}$.

(ii)
By the coercivity and continuity of $A^\dag_{\bZ,\alpha}$ from Theorem \ref{thm:Laplace_int} or Theorem \ref{thm:Laplace_ext},
$$
\frac{c}{2}\|\psi^{N}_i\|^2_{L^2(\Gamma)}\leq \big|\big(A^\dag_{\bZ,\alpha} \psi^{N}_i,\psi^{N}_i\big)_{L^2(\Gamma)}\big| \leq C\|\psi^{N}_i\|^2_{L^2(\Gamma)};
$$
therefore, since $|{\mymatrix{A}}_{ii}|=|(A^\dag_{\bZ,\alpha} \psi^{N}_i,\psi^{N}_i)_{L^2(\Gamma)}|$,
\begin{equation}\label{eq:temp2}
\frac{c}{2}\|\psi^{N}_i\|^2_{L^2(\Gamma)}\leq |{\mymatrix{A}}_{ii}|\leq C\|\psi^{N}_i\|^2_{L^2(\Gamma)}.
\end{equation}
Further, where $\tau$ is the support of $\psi^N_i$ and $\widehat{\tau}$ is the domain of $\chi_\tau$, by \eqref{eq:cov},
$$
\|\psi^{N}_i\|^2_{L^2(\Gamma)} = \int_{\tau} |\psi^\tau_j(\mathbf{x})|^2 \, {\rm d} s(\mathbf{x}) = \int_{\widehat{\tau}} |\psi^{\widehat{\tau}}_j(\widehat{\mathbf{x}})|^2 g_\tau(\widehat{\mathbf{x}})\, {\rm d} \widehat{\mathbf{x}},
$$
for some $j\in \{1,...,P_{\widehat \tau}\}$, so that, by \eqref{eq:mild2},
\begin{equation}\label{eq:temp3}
c_- c^{-1}_2 s_\tau\leq \|\psi^{N}_i\|^2_{L^2(\Gamma)}\leq c_+ c_2 s_\tau,
\end{equation}
 where
$$
c_+ := \max_{\widehat{\tau} \in \mathcal{R}, j=1,\ldots,P_{\widehat \tau}} \|\psi^{\widehat{\tau}}_j\|_{L^2(\widehat{\tau})}^2 \quad\text{ and }\quad c_- := \min_{\widehat{\tau} \in \mathcal{R}, j=1,\ldots,P_{\widehat \tau}} \|\psi^{\widehat{\tau}}_j\|_{L^2(\widehat{\tau})}^2.
$$
Thus, combining \eqref{eq:temp2} and \eqref{eq:temp3}, we have that
$$
\frac{c}{2}c_- c^{-1}_2 s_\tau\leq |{\mymatrix{A}}_{ii}|\leq Cc_+ c_2 s_\tau \quad \mbox{ for } i=1,\ldots, M_N.
$$
Therefore, if $({\mymatrix{D}}_N)_{ii} := |{\mymatrix{A}}_{ii}|$, for $i=1,\ldots, M_N$, then, by \eqref{eq:expand1},
$$
 \frac{c}{2}c_-c_2^{-1}\sum_{\tau\in \mathcal{G}_N}s_\tau \sum_{j=1}^{P_{\widehat \tau}} (w_j^\tau)^2 \leq \big\|{\mymatrix{D}}_N^{1/2} \mathbf{w}\big\|^2_{2} \leq Cc_+c_2 \sum_{\tau\in \mathcal{G}_N}s_\tau \sum_{j=1}^{P_{\widehat \tau}} (w_j^\tau)^2.
$$
By \eqref{eq:temp5} and \eqref{eq:temp4}, Assumption \ref{ass:D} therefore holds with $C_1$ and $C_2$ given by \eqref{eq:C1C2}.
\end{proof}

\begin{remark}[The novelty of Lemma \ref{lem:meshes}]
Similar results to Lemma \ref{lem:meshes} are given in \cite{GrMc:06}, where, for a continuous, coercive, and symmetric sesquilinear form, scaling by the diagonal part of the Galerkin matrix (as in Part (ii) of Lemma \ref{lem:meshes}) is used to remove the ill-conditioning of the Galerkin matrix due to mesh degeneracy; see \cite[Equations 1.5-1.7]{GrMc:06}.

The advantage of the results of \cite{GrMc:06} compared to those of Lemma \ref{lem:meshes} is that \cite{GrMc:06} works in $H^s(\Gamma)$ for $|s|\leq 1$, whereas Lemma \ref{lem:meshes} only works in $L^2(\Gamma)$. However, Lemma \ref{lem:meshes} works with rather general meshes in arbitrary dimensions, subject only to Assumption \ref{ass:mild}, whereas \cite{GrMc:06} imposes the following conditions on the mesh: (i)
the mesh is regular in the sense of \cite[Definition 4.1.4]{SaSc:11}, see \cite[Page1495]{GrMc:06}, and (ii) the mesh satisfies
\cite[Assumptions 3.1 and 3.2]{GrMc:06}, with the latter requiring, e.g., that neighbouring mesh elements have comparable aspect ratios.
\end{remark}

\section{Wellposedness and regularity results for the Laplace interior and exterior oblique  Robin problems}\label{sec:oblique}

\subsection{Statement of the Laplace interior and exterior oblique Robin problems}

\begin{definition}[The Laplace interior oblique Robin problem (IORP)]\label{prob:ioipL}
With $\Omega_-$
as in \S\ref{sec:notation}, given $g\in L^2(\Gamma)$, $\mathbf{Z}\in (L^\infty(\Gamma))^{d}$, and $\alpha\in L^\infty(\Gamma)$, find $u\in H^1(\Omega_-)$ with $\gamma^- u\in H^1(\Gamma)$ and $\partial_n^- u\in {L^2(\Gamma)}$
such that
$\Delta u =0$ in $\Omega_-$ and
\begin{equation}\label{eq:obliqueL}
(\mathbf{Z}\cdot\mathbf{n}) \partial_n^- u + \mathbf{Z} \cdot \nabla_\Gamma( \gamma^- u) + \alpha \, \gamma^- u =  g \quad\text{ on } \Gamma.
\end{equation}
\end{definition}

\begin{definition}[The Laplace exterior oblique Robin problem (EORP)]\label{prob:eoipL}
With $\Oe$
as in \S\ref{sec:notation}, given $g\in L^2(\Gamma)$, $\mathbf{Z}\in (L^\infty(\Gamma))^{d}$, and $\alpha\in \L^\infty(\Gamma)$, find $u\in H^1_{\mathrm{loc}}(\Oe)$ with $\gamma^+ u \in H^1(\Gamma)$
and $\partial_n^+ u\in {L^2(\Gamma)}$
 such that
$\Delta u =0$ in $\Oe$,
\begin{equation}\label{eq:obliqueLext}
(\mathbf{Z}\cdot\mathbf{n}) \partial_n^+ u + \mathbf{Z} \cdot \nabla_\Gamma( \gamma^+ u) - \alpha \, \gamma^+ u =  g \quad\text{ on } \Gamma,
\end{equation}
and,  as $|\mathbf{x}| \rightarrow \infty$, $u(\mathbf{x})= O(1)$ when $d=2$ and $u(\mathbf{x}) = o(|\mathbf{x}|^{3-d})$ when $d\geq 3$ (uniformly in all directions $\mathbf{x}/|\mathbf{x}|$).
\end{definition}

A regularity result of Ne\v{c}as \cite{Ne:67} (stated as Theorem \ref{thm:Necas} below) implies that either of the requirements $\partial_n^- u\in {L^2(\Gamma)}$ and $\gamma^- u \in H^1(\Gamma)$ in Definition \ref{prob:ioipL} can be removed; similarly in Definition \ref{prob:eoipL}.

The IORP and EORP can also be formulated in terms of non-tangential maximal functions and non-tangential limits
(similar to the case of the Dirichlet problem discussed in \S\ref{sec:notation}).
We now give this alternative formulation for the IORP and prove that it is equivalent to Definition \ref{prob:ioipL}; this equivalence is necessary to use results from the harmonic-analysis literature on the standard Laplace oblique derivative problem (see Theorem \ref{thm:RegOblSt} below).
The alternative formulation for the EORP and proof of equivalence to Definition \ref{prob:eoipL} are completely analogous and are omitted.

\begin{definition}[The Laplace IORP via non-tangential limits]\label{prob:ioipL2}
With $\Omega_-$
 as in \S\ref{sec:notation}, given $g\in L^2(\Gamma)$, $\mathbf{Z}\in (L^\infty(\Gamma))^{d}$, and $\alpha\in L^\infty(\Gamma)$,
find $u\in C^2(\Omega_-)$ with $(\nabla u)^*\in (L^2(\Gamma))^d$ such that
$\Delta u =0$ in $\Omega_-$ and
\begin{equation}\label{eq:obliqueL2}
\mathbf{Z} \cdot \widetilde{\gamma}^- (\nabla u ) + \alpha \widetilde{\gamma}^-u = g
 \quad\text{ on } \Gamma,
\end{equation}
where $\widetilde{\gamma}^-$ is the non-tangential limit defined by \eqref{eq:ntlim}.
\end{definition}

\begin{theorem}[Equivalence of the different formulations of the IORP]\label{thm:BVPequiv2}
The formulations of the IORP in Definition \ref{prob:ioipL} and \ref{prob:ioipL2} are equivalent (i.e., if $u$ is a solution to the IORP in the sense of Definition \ref{prob:ioipL}, then it is a solution in the sense of Definition \ref{prob:ioipL2}, and vice versa).
\end{theorem}

\begin{proof}
If $u$ is a solution of the IORP in the sense of Definition \ref{prob:ioipL}, then $u\in C^\infty(\Omega_-)$ by elliptic regularity. Furthermore, $u\in H^{3/2}(\Omega_-)$ by Lemma \ref{lem:A10}, and then $(\nabla u)^* \in L^2(\Gamma)$ by Part (iii) of Theorem \ref{thm:JeKe:95}. By Lemma \ref{lem:A9}, $\widetilde{\gamma}^-u = \gamma^- u$, and, by Lemma \ref{lem:A10},
\begin{equation}\label{eq:moderna1}
\widetilde{\gamma}^-(\nabla u) = \mathbf{n} \partial_n^- u + \nabla_\Gamma(\gamma^- u) \quad\text{ almost everywhere on } \Gamma.
\end{equation}
Therefore
\begin{equation}\label{eq:moderna2}
\bZ\cdot \widetilde{\gamma}^- (\nabla u) + \alpha \widetilde{\gamma}^- u = (\bZ\cdot \mathbf{n}) \partial_n^- u + \bZ \cdot\nabla_\Gamma (\gamma^- u) + \alpha \gamma^- u,
\end{equation}
so that the boundary condition \eqref{eq:obliqueL2} is equivalent to \eqref{eq:obliqueL}; therefore, $u$ is a solution of the IORP in the sense of Definition \ref{prob:ioipL2}.

Conversely, if $u$ is the solution of the IORP in the sense of Definition \ref{prob:ioipL2}, then $u\in H^{3/2}(\Omega_-)$ by Part (iii) of Theorem \ref{thm:JeKe:95}. Then
Lemma \ref{lem:A10} implies that $\partial_n^- u \in {L^2(\Gamma)}$, $\gamma^- u\in H^1(\Gamma)$,
and \eqref{eq:moderna1} holds. Hence \eqref{eq:moderna2} holds and the boundary condition \eqref{eq:obliqueL} is equivalent to \eqref{eq:obliqueL2}; therefore, $u$ is a solution of the IORP in the sense of Definition \ref{prob:ioipL}.
\end{proof}

\subsection{Link between the IORP/EORP and the BIEs in Theorems \ref{thm:Laplace_int}, \ref{thm:Laplace_ext}}

\begin{theorem}[$A^\prime_{I,\mathbf{Z},\alpha}$ can be used to solve the EORP for $d\geq 3$] \label{thm:star_in_obliqueLE} If $d\geq 3$ then the single-layer potential $u={\cal S} \phi$ with density $\phi\in L^2(\Gamma)$ satisfies the exterior oblique Robin problem (Definition \ref{prob:eoipL})  if and only if
\begin{equation} \label{eq:BIE_eoipL}
A^\prime_{I,\mathbf{Z},\alpha}\phi = -g.
\end{equation}
Conversely, if $d\geq 3$ and $u$ satisfies the EORP, then $u={\cal S} \phi$ for some $\phi\in L^2(\Gamma)$ that satisfies \eqref{eq:BIE_eoipL}.
\end{theorem}

\begin{theorem}[$A^\prime_{E,\mathbf{Z},\alpha}$ can be used to solve the IORP for $d\geq 2$] \label{thm:star_in_obliqueL} The single-layer potential $u={\cal S} \phi$, with density $\phi\in L^2(\Gamma)$, satisfies the IORP (Definition \ref{prob:ioipL}) if and only if
\begin{equation} \label{eq:BIE_ioipL}
A^\prime_{E,\mathbf{Z},\alpha}\phi = g.
\end{equation}
Conversely, if $u$ satisfies the IORP, then, provided $a\neq \mathrm{Cap}_\Gamma$ when $d=2$, $u={\cal S} \phi$, where $\phi\in L^2(\Gamma)$ satisfies \eqref{eq:BIE_ioipL}.
\end{theorem}

\begin{proof}[Proof of Theorem \ref{thm:star_in_obliqueLE}]  If $d \geq 3$ and $u={\cal S} \phi$ with $\phi\in L^2(\Gamma)$, then by, e.g., \cite[Theorem 2.14]{ChGrLaSp:12} $u\in C^2(\Oe)$ and $\Delta u =0$ in $\Oe$, and, by \eqref{eq:Sasym}, $u(\mathbf{x})=O(|\mathbf{x}|^{2-d})$ as $|\mathbf{x}|\to\infty$, uniformly in $\mathbf{x}/|\mathbf{x}|$.
By, e.g., \cite[Theorem 2.14]{ChGrLaSp:12}, $u\in H^1_{\mathrm{loc}}(\Oe)$ and, by the jump relations \eqref{eq:jump1} and the definition of $K_{\vfb}'$ \eqref{eq:Caldb'}, \eqref{eq:obliqueLext} holds if and only if $\phi$ satisfies \eqref{eq:BIE_eoipL}. Conversely, if $u$ satisfies the EORP, then, by the invertibility of $S$ recalled in Lemma \ref{lem:SLP} below,
  $\phi:= S^{-1} \gamma^+ u  \in L^2(\Gamma)$. Defining $v:={\cal S} \phi$, $v$ satisfies the Laplace exterior Dirichlet problem with boundary data $\gamma^+ v=\gamma^+{\cal S} \phi = S\phi = \gamma^+ u$, so that $v=u$ by uniqueness for the EDP. As established in the first part of the proof, since $u$ satisfies the EORP, $\phi$ satisfies \eqref{eq:BIE_eoipL}.
\end{proof}

\

\begin{proof}[Proof of Theorem \ref{thm:star_in_obliqueL}]
This is very similar to the proof of Theorem \ref{thm:star_in_obliqueLE}, except that now we can also consider $d=2$, since (by definition) there are no conditions at infinity imposed on the solution of the IORP.
\end{proof}

\begin{theorem}\label{thm:inverse_formula_L}
Let $P_{\rm DtN}^\pm :H^1(\Gamma)\rightarrow {L^2(\Gamma)}$ denote the Dirichlet-to-Neumann maps for Laplace's equation in $\Omega^\pm$; i.e., the maps $g_D\mapsto \partial_n^\pm u$ for $u$ as in Definitions \ref{def:idp}/\ref{def:edp} respectively.
Let $P_{\rm ItD}^{-, \alpha, \bZ}: {L^2(\Gamma)}\rightarrow H^1(\Gamma)$ denote the map $g\mapsto \gamma^- u$ where $u$ is as in Definition \ref{prob:ioipL}.
Let $P_{\rm ItD}^{+, \alpha, \bZ}: {L^2(\Gamma)}\rightarrow H^1(\Gamma)$ denote the map $g\rightarrow \gamma^+ u$ where $u$ is as in Definition \ref{prob:eoipL}.
Then, as operators on ${L^2(\Gamma)}$,
\begin{equation}\label{eq:inverse1}
\big( A'_{E,\bZ,\alpha} \big)^{-1} = \frac{1}{\bZ\cdot \mathbf{n}}I - \left( P^+_{\rm DtN} + \frac{1}{ \bZ\cdot \mathbf{n}}\big( \alpha +\bZ\cdot\nabla_{\Gamma} \big) \right)
P^{-, \alpha, \bZ}_{\rm ItD}
\end{equation}
and
\begin{equation}\label{eq:inverse2}
\big( A'_{I,\bZ,\alpha} \big)^{-1} = \frac{1}{\bZ\cdot \mathbf{n}}I - \left( P^-_{\rm DtN} + \frac{1}{ \bZ\cdot \mathbf{n}}\big( -\alpha +\bZ\cdot\nabla_{\Gamma} \big) \right)
P^{+, \alpha, \bZ}_{\rm ItD}.
\end{equation}
\end{theorem}

\begin{proof}
We first prove \eqref{eq:inverse1}. Suppose $A'_{E,\bZ,\alpha} \phi =g$ with $\phi,g\in {L^2(\Gamma)}$ and let $u:= {\cal S}\phi$. Then $\gamma^+ u=\gamma^- u= P^{-,\alpha,\bZ}_{\rm ItD} g$ by the first jump relation in \eqref{eq:jump1} and Theorem \ref{thm:star_in_obliqueL}. By the second jump relation in \eqref{eq:jump1}, the definition of $P_{\rm DtN}^+$, and the boundary condition \eqref{eq:obliqueL},
\begin{align*}
\phi &= \partial_n^- u - \partial_n^+ u, \\
&= \frac{1}{\bZ\cdot\mathbf{n}} \big( g - \bZ \cdot \nabla_{\Gamma}(\gamma^- u) - \alpha \gamma^- u\big) - P_{\rm DtN}^+ \gamma^+ u
= \frac{1}{\bZ\cdot\mathbf{n}}g - \left( P_{\rm DtN}^+ + \frac{\alpha + \bZ\cdot\nabla_{\Gamma}}{\bZ \cdot\mathbf{n}}\right) P_{\rm ItD}^{-,\alpha,\bZ}g,
\end{align*}
which implies \eqref{eq:inverse1}.
The proof of \eqref{eq:inverse2} is then very similar, using Theorem \ref{thm:star_in_obliqueLE} instead of
 Theorem \ref{thm:star_in_obliqueL}.
\end{proof}

\subsection{Statement of the wellposedness results and implications for the BIEs in Theorems \ref{thm:Laplace_int} and \ref{thm:Laplace_ext}} \label{subsec:oblique3}

\begin{theorem}[Uniqueness for the IORP] \label{thm:oblique:unique} Suppose that, for some $\beta\in (0,1]$, $\mathbf{Z}\in (C^{0,\beta}(\Gamma))^d$ and $\alpha\in C^{0,\beta}(\Gamma)$ and that, for some constants $c, c_0>0$,
\begin{equation} \label{eq:lowerbounds}
\mathbf{Z}(\mathbf{x})\cdot \mathbf{n}(\mathbf{x}) \geq c  \quad \mbox{for almost every } \mathbf{x}\in \Gamma \quad \mbox{and} \quad \alpha(\mathbf{x}) \geq c_0 \quad\text{ for } \mathbf{x}\in \Gamma.
\end{equation}
Then the IORP has at most one solution.
\end{theorem}

\begin{corollary}[Existence for the IORP and invertibility of $A'_{E,\bZ,\alpha}$]\label{cor:oblique:unique}
If the assumptions of Theorem \ref{thm:oblique:unique} hold and $a\neq \mathrm{Cap}_\Gamma$ when $d=2$, then
 $A'_{E,\bZ,\alpha}$ is invertible and the IORP has exactly one solution.
\end{corollary}

\begin{theorem}[Uniqueness for the EORP] \label{thm:oblique:uniqueE} Suppose that, for some $\beta\in (0,1]$, $\mathbf{Z}\in (C^{0,\beta}(\Gamma))^d$ and $\alpha\in C^{0,\beta}(\Gamma)$, and that \eqref{eq:lowerbounds} holds, for some constants $c, c_0>0$.
Then the EORP has at most one solution.
\end{theorem}

\begin{corollary}[Existence for the EORP and invertibility of $A'_{I,\bZ,\alpha}$]\label{cor:oblique:unique2}
If the assumptions of Theorem \ref{thm:oblique:uniqueE} hold and $d\geq 3$, then the EORP has exactly one solution and $A'_{I,\bZ,\alpha}$ is invertible.
\end{corollary}

\subsection{Proofs of Theorems \ref{thm:oblique:unique} and \ref{thm:oblique:uniqueE}}\label{sec:obliqueproofs}

Recall that, for $1\leq p\leq\infty$, $H^{1,p}(\Gamma):=\{\phi\in L^p(\Gamma):\nabla_\Gamma \phi \in L^p(\Gamma)\}$ is a Banach space with the norm $\|\phi\|_{H^{1,p}(\Gamma)} := \|\phi\|_{L^p(\Gamma)} + \|\nabla_\Gamma \phi\|_{L^p(\Gamma)} $. Note that $H^1(\Gamma)=H^{1,2}(\Gamma)$, with equivalence of norms.

The following result is standard in the theory of potential theory on Lipschitz domains; see, e.g., \cite[Page 203]{Ve:07}.

\begin{lemma}\label{lem:obliqueness}
Suppose that $\mathbf{Z}\in (C(\Gamma))^d$ and the first of the bounds \eqref{eq:lowerbounds} holds for some $c>0$. Then, for each $\mathbf{x}\in \Gamma$ there exists $R>0$ and $F\in C^{0,1}(\mathbb{R}^{d-1})$ and a rotated coordinate system $0\tilde x_1...\tilde x_d$, with origin at $\mathbf{x}$ and with the $\tilde x_d$ axis pointing in the direction $\mathbf{Z}(\mathbf{x})$, such that, where $\tilde y^\prime := (\tilde y_1,...,\tilde y_{d-1})$,
$$
B_R(\mathbf{x})\cap \Oe = B_R(\mathbf{x})\cap \{\mathbf{y} = (\tilde y^\prime,\tilde y_d):\tilde y_d > F(\tilde y^\prime)\}, \quad B_R(\mathbf{x})\cap \Omega_- = B_R(\mathbf{x})\cap \{\mathbf{y} = (\tilde y^\prime,\tilde y_d):\tilde y_d < F(\tilde y^\prime)\}.
$$
\end{lemma}

The following key regularity estimate follows immediately from \cite{Pi:87,Ve:07}.

\begin{theorem}[Regularity for the interior oblique derivative problem.] \label{thm:RegOblSt}

\

Suppose that $\mathbf{Z}\in (C^{0,\beta}(\Gamma))^d$ for some $\beta\in (0,1]$, the first inequality in \eqref{eq:lowerbounds} holds for some constant $c>0$,
and $u$ satisfies the \emph{Laplace oblique derivative problem} (i.e., the IORP in the special case $\alpha=0$) with data $g$.

(i) If $g\in C^{0,\beta}(\Gamma)$, then
$u\in C^{1,\gamma}(\overline{\Omega_-})$ for some $\gamma\in (0,\beta]$ depending only on $\Omega_-$.

(ii) If $g\in L^p(\Gamma)$ with $2\leq p <\infty$, then $(\nabla u)^*\in (L^p(\Gamma))^d$.
\end{theorem}

\begin{proof}
(i) It is known from \cite[Section 4]{Ca:85}, \cite{Pi:87, KePi:88, Ve:07} that
if the first inequality in \eqref{eq:lowerbounds} holds, then the Laplace oblique derivative problem has a solution if and only if $g$ satisfies finitely-many linear conditions (i.e., conditions of the form $(g,\phi_j)_\Gamma=0$, $j=1,...,N$, for some $N\in \mathbb{N}$ and $\phi_1,...,\phi_N\in L^2(\Gamma)$).
 If $\mathbf{Z}\in (C^{0,\beta}(\Gamma))^d$ and $u$ is a solution for particular data $g\in C^{0,\beta}(\Gamma)$, the finitely-many linear conditions on $g$ are satisfied, and $u$ can be written as $u=u_P+u_H$ where $u_P$ is the particular solution studied in \cite{Pi:87}, which is shown in \cite[\S3]{Pi:87} to satisfy $u_P
 \in C^{1,\gamma}(\Omega_-)$ for some $\gamma\in (0,\beta]$ (dependent on $\Omega_-$),
 and $u_H$ is a solution of the homogeneous oblique derivative problem, which is shown in \cite[Corollary 2.7]{Ve:07} to be constant in $\Omega_-$.

(ii) This follows from arguing as in (i), but replacing the results of \cite{Pi:87} for H\"older continuous $g$ by those of \cite{Ca:85} for $g\in L^p(\Gamma)$ with $2-\epsilon<p<2+\epsilon$ (for some $\epsilon>0$ dependent on $\Omega_-$) and \cite{KePi:88} for $g\in L^p(\Gamma)$ with $p>2$ (note that while the results of \cite{Ca:85, KePi:88} only require that $\mathbf{Z}$ is continuous, \cite[Corollary 2.7]{Ve:07} requires $\mathbf{Z}$ to be H\"older continuous).
\end{proof}

\begin{theorem}[Regularity for the IORP] \label{thm:RegObl} Suppose that $\mathbf{Z}\in (C^{0,\beta}(\Gamma))^d$ and $\alpha\in C^{0,\beta}(\Gamma)$  for some $\beta\in (0,1]$, $\mathbf{Z}$ satisfies the first inequality in \eqref{eq:lowerbounds} for some $c>0$,
and $u$ satisfies the Laplace interior oblique Robin problem with data $g\in C^{0,\beta}(\Gamma)$. Then $u\in C^{1,\gamma}(\overline{\Omega_-})$ for some $\gamma\in (0,\beta]$.
\end{theorem}
\begin{proof}
Suppose that the conditions of the theorem are satisfied, in particular that $u$ satisfies the IORP with data $g\in C^{0,\beta}(\Gamma)$, for some $\beta\in (0,1]$. Suppose also that $2\leq p <\infty$ and that  $\gamma^- u \in L^p(\Gamma)$. Then since, clearly, $g\in L^p(\Gamma)$, $u$ is a solution of the Laplace oblique derivative problem with data in $L^{p}(\Gamma)$. Therefore, by Part (ii) of Theorem \ref{thm:RegOblSt}, $(\nabla u)^* \in (L^{p}(\Gamma))^d$ and thus $\gamma^- u \in H^{1,p}(\Gamma)$ by Corollary \ref{cor:Lp}. This implies, by the Sobolev embedding theorem \cite[Chapter V, Equations 6 and 4]{Adams}, that, if $p\geq d-1$, then
$\gamma^- u \in L^{q}(\Gamma)$ for all $2\leq q<\infty$, while, if $p<d-1$, then $\gamma^- u \in L^{q}(\Gamma)$ for $2\leq q \leq p_0$ where $1/p_0 = 1/p -1/(d-1)$. Since, certainly, $\gamma^- u \in {L^2(\Gamma)}$ (as $\gamma^- u\in H^1(\Gamma)$ by definition of the IORP), applying the above argument at most a finite number of times leads to the conclusion that $\gamma^- u \in H^{1,q}(\Gamma)$ for all $2\leq q<\infty$.
But this implies, by the Sobolev embedding theorem \cite[Chapter V, Equation 9]{Adams}, that $\gamma^- u \in C^{0,\beta'}(\Gamma)$ for all $0<\beta'<1$. Thus $u$ is a solution of the Laplace oblique derivative problem with data in $C^{0,\beta}(\Gamma)$, and thus the result that $u\in C^{1,\gamma}(\overline{\Omega_-})$ for some $\gamma\in (0,\beta]$ follows from Part (i) of Theorem \ref{thm:RegOblSt}.
\end{proof}

\begin{corollary}[Regularity for the EORP] \label{thm:RegOblE} Suppose that $\mathbf{Z}\in (C^{0,\beta}(\Gamma))^d$ and $\alpha\in C^{0,\beta}(\Gamma)$ for some $\beta\in (0,1]$, $\mathbf{Z}$ satisfies the first inequality in \eqref{eq:lowerbounds} for some $c>0$, and $u$ satisfies the Laplace exterior oblique Robin problem with data $g\in C^{0,\beta}(\Gamma)$. Then $u\in C^{1,\gamma}(\overline{\Oe\cap B_R})$ for all $R>0$ and for some $\gamma\in (0,\beta]$.
\end{corollary}

\begin{proof}
Since $\Omega_-$ is bounded, $\Gamma \subset B_r$ for some $r>0$. Suppose that $u$ satisfies the EORP and choose $R_2>R_1>R_0>r$ and $\chi\in C_{\rm comp}^\infty(\mathbb{R}^d)$ with $\chi(\mathbf{x})=1$ for $|\mathbf{x}|\leq R_0$ and $\chi(\mathbf{x})=0$ for $|\mathbf{x}|\geq R_1$. Let $v(\mathbf{x}):= \chi(\mathbf{x}) u(\mathbf{x})$ for $\mathbf{x}\in G:=\Oe\cap B_{R_2}$, so that, in particular, $v=u$ in $\Oe\cap B_{R_0}$.
The idea now is to create a solution of an IORP on $G$, and then use the interior regularity result of Theorem \ref{thm:RegObl}.
Since $u$ is harmonic in $\Oe$,
$$
\Delta v = F:= 2\nabla \chi\cdot \nabla u + u \Delta \chi \quad \mbox{in } G.
$$
Since $\chi \in C_{\rm comp}^\infty(\mathbb{R}^d)$ and $u\in C^{\infty}(\Oe)$, $F\in C_{\rm comp}^\infty(G)$. Therefore
$$
\widehat v(\mathbf{x}) := -\int_{G} \Phi(\mathbf{x},\mathbf{y})F(\mathbf{y})\, {\rm d} \mathbf{y}, \quad \text{ for } \mathbf{x} \in G,
$$
satisfies $\widehat v\in C^2(\overline{G})$ and $\Delta \widehat v = F$ in $G$. Let $w(\mathbf{x}):= v(\mathbf{x})-\widehat v(\mathbf{x})$ for $\mathbf{x}\in G$, and define $\widetilde \mathbf{Z}\in (C^{0,\beta}(\partial G))^d$ and $\widetilde \alpha\in C^{0,\beta}(\partial G)$ by  $\widetilde \mathbf{Z}:= - \mathbf{Z}$ on $\Gamma$, $\widetilde \mathbf{Z}(\mathbf{x}):= \mathbf{x}$, for $\mathbf{x}\in \partial B_{R_2}$, and $\widetilde \alpha := \alpha$ on $\Gamma$, $\widetilde \alpha := 0$ on $\partial B_{R_2}$. Then $w\in H^1(G)$ with trace $\gamma w\in H^1(\partial G)$, $\Delta w=0$ in $G$, and
\begin{equation*}
(\widetilde \mathbf{Z}\cdot\mathbf{n}) \partial_\mathbf{n} w + \widetilde \mathbf{Z} \cdot \nabla_\Gamma( \gamma w) + \widetilde \alpha \, \gamma w =  \widetilde g \text{ on } \partial G,
\end{equation*}
where $\mathbf{n}$ is the unit normal pointing out of $G$ and $\widetilde g\in C^{0,\beta}(\partial G)$ is defined by $\widetilde g := -(\widetilde \mathbf{Z}\cdot \nabla \widehat v +\widetilde \alpha \widehat v) - g$ on $\Gamma$, and by $\widetilde g := -(\widetilde \mathbf{Z}\cdot \nabla \widehat v +\widetilde \alpha \widehat v)$ on $\partial B_{R_2}$. Theorem \ref{thm:RegObl} implies that, for some $\gamma\in (0,\beta]$, $w\in C^{1,\gamma}(\overline{G})$, so that $v \in C^{1,\gamma}(\overline{G})$ and $u \in C^{1,\gamma}(\overline{\Oe\cap B_{R_0}})$. Since $u$ is harmonic in $\Oe$, $u \in C^{1,\gamma}(\overline{\Oe\cap B_{R}})$ for every $R>0$.
\end{proof}

\

We can now prove Theorems \ref{thm:oblique:unique} and \ref{thm:oblique:uniqueE} and Corollaries \ref{cor:oblique:unique} and \ref{cor:oblique:unique2}.

\

\begin{proof}[Proof of Theorem \ref{thm:oblique:unique}] Suppose that $u$ satisfies the IORP with $g=0$ and that, without loss of generality, $u$ is real-valued. To show that $u=0$ it is enough to show that $u\leq 0$ in $\Omega_-$, since this implies, by the same argument applied to $-u$, that also $u\geq 0$, and hence $u=0$. By Theorem \ref{thm:RegObl}, $u\in C^1(\overline{\Omega_-})$ (indeed $\nabla u$ is H\"older continuous). By the maximum principle, since $u\in C^2(\Omega_-)\cap C(\overline{\Omega_-})$ is harmonic in $\Omega_-$, the maximum value of $u$ in $\overline{\Omega_-}$ is attained at some point $\mathbf{x}_0\in\Gamma$. Since $u\in C^1(\overline{\Omega_-})$ it follows from \eqref{eq:obliqueL} with $g=0$ that
$$
\alpha(\mathbf{x}_0)u(\mathbf{x}_0) = -\mathbf{Z}(\mathbf{x}_0)\cdot \nabla u(\mathbf{x}_0)=-\lim_{h\to 0^+}\frac{u(\mathbf{x}_0)-u(\mathbf{x}_0-h\mathbf{Z}(\mathbf{x}_0))}{h}.
$$
Since $\mathbf{Z}$ is continuous and $\mathbf{Z}\cdot \mathbf{n}\geq c\geq 0$ almost everywhere on $\Gamma$,
$\mathbf{x}_0-h\mathbf{Z}(\mathbf{x}_0) \in \overline \Omega_-$ for all sufficiently small $h>0$ by Lemma \ref{lem:obliqueness},
so that $\mathbf{Z}(\mathbf{x}_0)\cdot \nabla u(\mathbf{x}_0)\geq 0$ since $\mathbf{x}_0$ is the global maximum. Since $\alpha(\mathbf{x}_0) >0$, it follows that $u(\mathbf{x}_0)\leq 0$, so that $u\leq 0$ in $\Omega_-$.
\end{proof}

\

\begin{proof}[Proof of Theorem \ref{thm:oblique:uniqueE}] Suppose that $u$ satisfies the EORP with $g=0$ and, without loss of generality, is real-valued. As in the proof of Theorem \ref{thm:oblique:unique}, it is enough to show that $u\leq 0$ in $\Oe$. We recall that, when $d=2$, the condition that
$u$ is bounded on $\Oe$ implies that, for some $u_\infty\in \mathbb{R}$,
$$
u(\mathbf{x}) = u_\infty + O(|\mathbf{x}|^{-1}) \quad \mbox{as } |\mathbf{x}|\to \infty,
$$
uniformly in $\mathbf{x}/|\mathbf{x}|$, and that
\begin{equation} \label{eq:average}
u_\infty = \frac{1}{2\pi R}\int_{\partial B_R} u\, {\rm d} s
\end{equation}
if $\Gamma\subset B_R$ \cite[Equation 6.11]{Kr:89}.

By Corollary \ref{thm:RegOblE}, $u\in C^1(\overline {\Oe})$. By the maximum principle, since $u\in C^2(\Oe)\cap C(\overline{\Oe})$ is harmonic in $\Oe$, the maximum value of $u$ in $\overline{\Oe}$ is attained on $\Gamma$ or, when $d=2$, $u(\mathbf{x})\leq u_\infty$ for $\mathbf{x}\in \Oe$.
If the maximum is attained on $\Gamma$, the result that $u\leq 0$ follows by arguing as in the proof of Theorem \ref{thm:oblique:unique}.
Therefore, it is sufficient to prove that the maximum is attained on $\Gamma$ when $d=2$. If $u(\mathbf{x})\leq u_\infty$ for $\mathbf{x}\in \Oe$, then \eqref{eq:average} implies that $u(\mathbf{x})= u_\infty$ for $|\mathbf{x}|\geq R$ if $\Gamma\subset B_R$, so that the maximum is attained in $\Oe$. The maximum principle (see, e.g., \cite[Theorem 6.8]{Kr:89}) then implies that $u$ is constant in $\Oe$, so that the maximum is also attained on $\Gamma$.
\end{proof}

\

The following proofs of Corollaries \ref{cor:oblique:unique} and \ref{cor:oblique:unique2} use the fact that, when $\alpha\in L^\infty(\Gamma)$, $\mathbf{Z}$ is continuous, and \eqref{eq:c} (i.e., the first lower bound in \eqref{eq:lowerbounds}) holds, then $A'_{I,\bZ,\alpha}$ and $A'_{E,\bZ,\alpha}$ are Fredholm of index zero by Parts (iii) and (iv) of Theorem \ref{thm:Laplace_int} and \ref{thm:Laplace_ext} respectively. Although these two theorems are for $d=3$, Parts (iii) and (iv) also hold when $d=2$ (as noted at the beginning of \S\ref{sec:d=2}).

\

\begin{proof}[Proof of Corollary \ref{cor:oblique:unique}]
If we can prove invertibility of $A'_{E,\bZ,\alpha}$, then existence of a solution to the IORP follows from Theorem \ref{thm:star_in_obliqueLE}. Since $A'_{E,\bZ,\alpha}$ is Fredholm of index zero on ${L^2(\Gamma)}$, by the Fredholm alternative (see, e.g., \cite[Theorem 2.27]{Mc:00}), to prove invertibility it is sufficient to prove injectivity.
Assume that  $A'_{E,\bZ,\alpha}\phi=0$ for $\phi\in L^2(\Gamma)$. By Theorem \ref{thm:star_in_obliqueL},
$u:= {\cal S}\phi$ satisfies the IORP, and by Theorem \ref{thm:oblique:unique} $u=0$ in $\Omega_-$.
Therefore $\gamma^- u=0$ and the first jump relation in \eqref{eq:jump1} implies that $S\phi =0$. Lemma \ref{lem:SLP} then implies that $\phi=0$ and the proof is complete.
\end{proof}

\

\begin{proof}[Proof of Corollary \ref{cor:oblique:unique2}]
This is very similar to that of Corollary \ref{cor:oblique:unique} except that now we only work in $d\geq 3$, since Theorem \ref{thm:star_in_obliqueLE} requires $d\geq 3$.
\end{proof}

\begin{remark}[The results of \cite{Li:87a}]
Although not directly used to prove the results in this section, the results of \cite{Li:87a} concern the Laplace IORP in Lipschitz domains with
H\"older continuous $\widetilde{\bZ}$ and $g$, and we comment here on their relevance to the results above.

The results of \cite{Li:87a} give an alternative route for obtaining uniqueness of the IORP (i.e., proving Theorem \ref{thm:oblique:unique}). Indeed, in the proof of Theorem \ref{thm:oblique:unique}, once we have established that $u\in C^1(\overline{\Omega_-})$ (by using Theorem \ref{thm:RegObl}), then uniqueness follows from \cite[Theorem 3.2]{Li:87a}.
The reason we argue as we do in the proof of Theorem \ref{thm:oblique:unique} is that this argument easily carries over to the proof of uniqueness for the EORP (Theorem \ref{thm:oblique:uniqueE}), whereas \cite[Theorem 3.2]{Li:87a} concerns only the IORP.

Furthermore, \cite[Theorem 3.2]{Li:87a} implies that, under the assumptions of Theorem \ref{thm:RegObl}, there exists $\beta_0<1$, depending only on the Lipschitz constant of $\Omega_-$, such that if $\beta< \beta_0$ then $u \in C^{1,\beta}(\overline{\Omega_-})$.
\end{remark}

\begin{remark}[Additional uniqueness results for the EORP with $\bZ=\mathbf{x}$]
The coercivity result of Theorem \ref{thm:star} allows us to extend the range of $\alpha$ for which the EORP is unique when $\bZ=\mathbf{x}$ and $d\geq 3$.
Indeed, Theorem \ref{thm:star} implies that $A'_{I,\bZ,\alpha}$ is injective when $\bZ=\mathbf{x}, \alpha(\mathbf{x})\geq -(d-2)/2$ for almost every $\mathbf{x}\in \Gamma$, and $d\geq 3$. Then, using Theorem \ref{thm:star_in_obliqueLE} and arguing as at the end of the proof of Corollary \ref{cor:oblique:unique}, we see that the solution of the EORP is unique under these conditions.
This result proves uniqueness for certain non-positive values of $\alpha$, which are not covered by Theorem \ref{thm:oblique:uniqueE}.
\end{remark}

\subsection{Link between the IORP/EORP and the BIEs in Theorem \ref{thm:2d}}

\begin{lemma}\label{lem:2doblique1}
Given $g\in {L^2(\Gamma)}$, if $\phi$ satisfies
\begin{equation}\label{eq:2doblique1}
T_{E,\bZ,\alpha,\beta}' \phi = g,
\end{equation}
(with $T_{E,\bZ,\alpha,\beta}'$ defined by \eqref{eq:TEprime})
and $d=2$, then
\begin{equation}\label{eq:2doblique2}
u = {\cal S} Q_\Gamma \phi + \frac{\beta}{\alpha} P_\Gamma \phi - \frac{1}{\alpha}P_\Gamma A'_{E,\bZ,\alpha} Q_\Gamma \phi
\end{equation}
satisfies the IORP.
\end{lemma}

\begin{lemma}\label{lem:2doblique2}
Given $g\in {L^2(\Gamma)}$, if $\phi$ satisfies
\begin{equation}\label{eq:2doblique1a}
T_{I,\bZ,\alpha,\beta}' \phi = -g,
\end{equation}
(with $T_{I,\bZ,\alpha,\beta}'$ defined by \eqref{eq:TIprime})
and $d=2$, then
\begin{equation}\label{eq:2doblique2a}
u = {\cal S} Q_\Gamma \phi + \frac{\beta}{\alpha} P_\Gamma \phi - \frac{1}{\alpha}P_\Gamma A'_{I,\bZ,\alpha} Q_\Gamma \phi
\end{equation}
satisfies the EORP.
\end{lemma}

\begin{proof}[Proofs of Lemmas \ref{lem:2doblique1}, \ref{lem:2doblique2}]
The fact that $u$ given by \eqref{eq:2doblique2}/\eqref{eq:2doblique2a} is $C^2$ and satisfies Laplace's equation follows from, e.g., \cite[Theorem 2.14]{ChGrLaSp:12}. The condition that $u= O(1)$ at infinity for the EORP follows from the asymptotics \eqref{eq:Sasym2}, the definition of $P_\Gamma$ \eqref{eq:P}, and that $Q_\Gamma:= I- P_\Gamma$.
The BIEs \eqref{eq:2doblique2}/\eqref{eq:2doblique2a} follow from the jump relations \eqref{eq:jump1} and the definitions of $A'_{I,\bZ,\alpha}$ \eqref{eq:opLIZ}, $A'_{E,\bZ,\alpha}$ \eqref{eq:opLEZ}, and $K_{\vfb}'$ \eqref{eq:Caldb'}.
\end{proof}

\begin{remark}[Link with the work of Medkova \cite{Me:18}]\label{rem:Medkova_oblique}
In \cite[Theorem 5.23.5]{Me:18}, the solution of the IORP is sought as \eqref{eq:2doblique2} without the final term on the right-hand side,
resulting in the BIE $(A'_{E,\bZ,\alpha} Q_\Gamma + \betaP_\Gamma)\phi =g$; this BIO is then proved to be invertible on ${L^2(\Gamma)}$ if $\beta=\alpha$ and $\alpha$ is sufficiently large \cite[Theorem 5.23.5]{Me:18}.
The advantage of including the final term on the right-hand side of \eqref{eq:2doblique2} is that, by Theorem \ref{thm:2d}, the resulting BIO $T_{E,\bZ,\alpha,\beta}' $ is not just invertible when $\alpha$ is sufficiently large, but also coercive by Part (v) of Theorem \ref{thm:2d}.
\end{remark}

\section{New formulations of the Helmholtz exterior Dirichlet problem}\label{sec:Helmholtz}

\subsection{Statement of the new formulations}\label{sec:Helmholtz_state}

As a corollary of Theorem \ref{thm:Laplace_ext}, we obtain results about BIE formulations of the Helmholtz exterior Dirichlet problem. We state the results in this subsection but defer the more substantial proofs to \S\ref{sec:Helmholtz_prove} (and see Remark \ref{rem:earlier}).

\begin{definition}[Helmholtz EDP]\label{def:Hedp}
With $\Omega_-$ and $\Oe$ as in \S\ref{sec:notation}, assume further that $\Oe$ is connected.
Given $g_D\in H^{1/2}(\Gamma)$, we say $u\in H^1_{\rm{loc}}(\Oe)$ satisfies the Helmholtz exterior Dirichlet problem (EDP) if $\Delta u +k^2u  =0$ in $\Oe$, $\gamma^+ u = g_D$ on $\Gamma$, and
 $u$ satisfies the Sommerfeld radiation condition
\begin{equation}\label{eq:src}
\pdiff{u}{r}(\mathbf{x}) - {\rm i} k u(\mathbf{x}) = o\left( \frac{1}{r^{(d-1)/2}}\right)
\end{equation}
as $r:= |\mathbf{x}|\rightarrow \infty$, uniformly in $\mathbf{x}/r$.
\end{definition}

Existence and uniqueness of the solution of the Helmholtz EDP for Lipschitz $\Omega_-$ is shown in, e.g., \cite[Theorem 2.10]{ChGrLaSp:12}.
Theorem \ref{thm:Hedp} below considers $g_D\in H^1(\Gamma)$; this holds, for example, for the \emph{sound-soft scattering problem} when $u$ corresponds to the scattered field and $g_D$ is the Dirichlet trace of the ($C^\infty$) incident field -- see, e.g., \cite[Proof of Theorem 2.12]{ChGrLaSp:12}.

\begin{definition}[Helmholtz EDP formulated via non-tangential limits]\label{def:Hedp2}
With $\Omega_-$ and $\Oe$ as  in \S\ref{sec:notation}, assume further that $\Oe$ is connected.
Given $g_D\in L^2(\Gamma)$, we say that $u\in C^2(\Oe)$ with $u^* \in L^2(\Gamma)$ satisfies the Helmholtz EDP if $\Delta u +k^2u=0$ in $\Oe$, $\widetilde{\gamma}^+ u = g_D$ on $\Gamma$, and $u$ satisfies the Sommerfeld radiation condition \eqref{eq:src}.
\end{definition}

The solution of the EDP in the sense of Definition \ref{def:Hedp} exists and is unique by \cite[Theorem 5.6, Part (ii)]{ToWe:93}.
The following equivalence result is proved in Appendix \ref{app:BVPequiv}.

\begin{theorem}
\label{thm:BVPequivH}
If $g_D\in H^{1/2}(\Gamma)$, then the solution of the Helmholtz EDP in the sense of Definition \ref{def:edp} is the solution of the Helmholtz EDP in the sense of Definition \ref{def:edp2}, and vice versa.
\end{theorem}

The fundamental solution for the Helmholtz equation is defined by
\begin{equation}\label{eq:fundH}
\Phi_k(\mathbf{x},\mathbf{y}):=
\frac{{\rm i}}{4} \left( \frac{k}{2\pi|\bx-\by|}\right)^{(d-2)/2} H^{(1)}_{(d-2)/2}\big(k|\bx-\by|\big);
\end{equation}
see, e.g., \cite[Equation 5.118]{Sta67}; the definition of the Hankel function $H^{(1)}_{(d-2)/2}$ implies that $\Phi_k(\mathbf{x},\mathbf{y})= ({\rm i}/4)H_0^{(1)}(k|\bx-\by|)$ for $d=2$ and $\exp({\rm i} k |\bx-\by|)/ (4\pi|\bx-\by|)$ for $d=3$.

The Helmholtz
BIOs $S_k, D_k, D'_k,$ $H_k$, $K_{\bZ,k}$, and $K_{\bZ,k}'$, are defined by \eqref{eq:bio1}, \eqref{eq:bio2}, \eqref{eq:Caldb}, and \eqref{eq:Caldb'}, respectively, with $\Phi$ replaced by $\Phi_k$.
Similarly, the potentials ${\cal S}_k$ and ${\cal K}_{\bZ,k}$ are defined by \eqref{eq:SLP} and \eqref{eq:Caldp}, respectively, with $\Phi$ replaced by $\Phi_k$.

Given $\bZ\in (L^\infty
(\Gamma))^d$ and $\eta \in L^\infty(\Gamma)$,
in analogy with \eqref{eq:opLEZ}, define the integral operators
 $A^\prime_{k,\eta,\bZ}$
 $A_{k,\eta,\bZ}$,
 and $B_{k,\eta,\bZ}$ by
\begin{align} \label{eq:AketaZpr_def}
&A^\prime_{k,\eta,\bZ} :=  \frac{1}{2} (\bZ\cdot \mathbf{n}) I + K_{\bZ,k}' - {\rm i} \eta S_k, \qquad
A_{k,\eta,\bZ} :=  \frac{1}{2} (\bZ\cdot \mathbf{n}) I + K_{\bZ,k} - {\rm i} \eta S_k, \\
&B_{k,\eta,\bZ} := (\bZ\cdot \mathbf{n}) H_k + \bZ\cdot\nabla_\Gamma \left(- \frac{1}{2}I + D_k\right) - {\rm i} \eta  \left(-{\frac{1}{2}}I+ D_k\right).
\label{eq:BketaZ_def}
\end{align}
Regarding notation, since we are only dealing with the exterior Helmholtz Dirichlet problem, we omit the subscript $E$ present in the analogous operators for Laplace's equation, but add a subscript $k$ to highlight the $k$ dependence. We use the notation $-{\rm i}\eta$ rather than $\alpha$ for consistency with standard notation for Helmholtz BIEs -- see \eqref{eq:opA} below -- and allow $\eta$ to be a function, rather than just a constant, to make a link to existing results -- see Theorem \ref{thm:Hoblique} below.

\begin{theorem}[New integral equations for the Helmholtz EDP]\label{thm:Hedp}

\

(i) \textbf{\emph{Direct formulation.}} Let $u$ be the solution of the  Helmholtz EDP of Definition \ref{def:Hedp} with additionally $g_D\in H^1(\Gamma)$. Then $\partial_n^+ u$ satisfies
\begin{equation}\label{eq:Hnew_ext}
A^\prime_{k,\eta,\bZ}\partial_n^+ u =B_{k,\eta,\bZ} g_D.
\end{equation}

(ii) \textbf{\emph{Indirect formulation.}} Given $g_D\in {L^2(\Gamma)}$, if $\phi$ satisifes
\begin{equation}\label{eq:Hnew_ext_indirect}
A_{k,\eta,\bZ}\phi = g_D,
\end{equation}
then
\begin{equation*}
u:= ({\cal K}_{\bZ,k} -{\rm i} \eta {\cal S}_k)\phi
\end{equation*}
is the solution of the Helmholtz EDP of Definition \ref{def:Hedp2}.

  (iii) \textbf{\emph{Continuity.}}
$A^\prime_{k,\eta,\bZ} : {L^2(\Gamma)}\rightarrow{L^2(\Gamma)}$,
$A_{k,\eta,\bZ}: {L^2(\Gamma)}\rightarrow{L^2(\Gamma)}$,
and $B_{k,\eta,\bZ}:H^1(\Gamma)\rightarrow{L^2(\Gamma)}$, and these mappings are continuous.

 (iv) \textbf{\emph{Coercivity up to compact perturbation for all $\eta$.}} If $\mathbf{Z}\in (C(\Gamma))^d$ and there exists $c>0$ such that \eqref{eq:c} holds,
 then, for all $\eta \in L^\infty(\Gamma)$, both $A^\prime_{k,\eta,\bZ}$ and $A_{k,\eta,\bZ}$ are the sum of a coercive operator and a compact operator on ${L^2(\Gamma)}$.
\end{theorem}

If $\bZ=\mathbf{n}$ then $A^\prime_{k,\eta,\bZ}=A'_{k,\eta}$, $A_{k,\eta,\bZ}=A_{k,\eta}$, and $B_{k,\eta,\bZ}=B_{k,\eta}$ where
\begin{equation}\label{eq:opA}
A'_{k,\eta} := \frac{1}{2} I + D'_k - {\rm i} \eta S_k,
\qquad A_{k,\eta} := \frac{1}{2} I + D_k - {\rm i} \eta S_k,
\quad\text{ and }\quad B_{k,\eta}:=H_k + {\rm i} \eta \left(\frac{1}{2} I - D_k \right)
\end{equation}
are the operators appearing in the standard ``combined field" or ``combined potential" BIEs for the Helmholtz EDP; see, e.g., \cite[Equations 2.68 and 2.69]{ChGrLaSp:12}.
Recall from \S\ref{sec:openbook} that there exist Lipschitz polyhedra such that each of $A'_{k,\eta}$ and $A_{k,\eta}$ is \emph{not} the sum of a coercive operator and a compact operator on ${L^2(\Gamma)}$.

By Part (b) of Theorem \ref{thm:Galerkin}, to prove that the Galerkin method converges when applied to either of the BIEs \eqref{eq:Hnew_ext} or \eqref{eq:Hnew_ext_indirect},
it is sufficient to show that  $A^\prime_{k,\eta,\mathbf{Z}}$ and $A_{k,\eta,\mathbf{Z}}$ are injective.
Analytic Fredholm theory and the fact that $A^\prime_{k,\eta,\mathbf{Z}}$ is a compact perturbation of $A'_{E,\bZ,\alpha}$ imply that $A^\prime_{k,\eta,\bZ}$ is invertible for all except at most a discrete set of wavenumbers.

\begin{lemma}\label{lem:analyticFredholm}
If $\bZ$ is Lipschitz, then
there exists a discrete set $D \subset \mathbb{C}$ (i.e.~each element of $D$ is isolated) such that
$A^\prime_{k,\eta,\bZ}$ is invertible on ${L^2(\Gamma)}$
  for all $k\in \mathbb{C}\setminus D$ ($d\geq 3$ and odd) or $k\in (\mathbb{C}\setminus \mathbb{R}^-)\setminus D$ ($d$ even), where $\mathbb{R}^-:= (-\infty,0]$.
\end{lemma}

Exactly as in the Laplace case, injectivity of $A^\prime_{k,\eta,\mathbf{Z}}$ and $A_{k,\eta,\mathbf{Z}}$ is equivalent to uniqueness of the solution of an interior oblique Robin problem.

\begin{definition}[The Helmholtz interior oblique impedance problem]\label{prob:ioip}
Let $\Omega_-$ be a bounded Lipschitz open set. Given $g\in L^2(\Gamma)$, $\bZ\in (L^\infty(\Gamma))^{d}$, and $\eta\in L^\infty(\Gamma)$, find $u\in H^1(\Omega_-)$ with $\gamma^- u \in H^1(\Gamma)$ and $\partial_n^- u\in {L^2(\Gamma)}$ such that
$\Delta u +k^ 2u =0$ in $\Omega_-$ and
\begin{equation}\label{eq:oblique}
(\bZ\cdot\mathbf{n}) \partial_n^- u + \bZ \cdot \nabla_\Gamma( \gamma^- u) -{\rm i} \eta \, \gamma^- u =  g \text{ on } \Gamma.
\end{equation}
\end{definition}

\begin{lemma}[$A^\prime_{k,\eta,\bZ}$ can be used to solve the interior oblique impedance problem] \label{thm:star_in_oblique} The single-layer potential $u={\cal S}_k \phi$, with density $\phi\in L^2(\Gamma)$, satisfies the interior oblique impedance problem \eqref{prob:ioip} if
\begin{equation} \label{eq:BIE_ioip1}
A^\prime_{k,\eta,\bZ}\phi = g.
\end{equation}
Conversely, if $u$ satisfies the interior oblique impedance problem, then $u={\cal S}_k \phi$, with $\phi\in L^2(\Gamma)$ satisfying \eqref{eq:BIE_ioip1}.
\end{lemma}

\begin{corollary}\label{cor:star_ioip_equiv}
 For $\eta\in L^\infty(\Gamma)$ and $\bZ\in (L^\infty(\Gamma))^{d}$, the operator $A^\prime_{k,\eta,\bZ}:L^2(\Gamma)\to L^2(\Gamma)$ is injective if and only if the interior oblique impedance problem with $g=0$ has only the trivial solution, and is surjective if and only if the interior oblique impedance problem has a solution for every $g\in L^2(\Gamma)$.
\end{corollary}

Whereas BVPs involving oblique derivatives have been well-studied for Laplace's equation (see the references in Theorem \ref{thm:RegOblSt} and also \cite{Li:13a}, \cite[\S\S5.23, 5.24, 6.19]{Me:18}), there do not appear to be any results in the literature on the unique solvability of the interior oblique impedance problem for the Helmholtz equation for all wavenumbers and general Lipschitz $\Omega_-$.
The following theorem collects three situations in which the solution of the interior oblique impedance problem is known to be unique, and hence, by the results above, the Galerkin method
applied to either of the BIEs \eqref{eq:Hnew_ext} or \eqref{eq:Hnew_ext_indirect}
is provably convergent
for every asymptotically-dense sequence of subspaces.

\begin{theorem}\label{thm:Hoblique}
The solution of the Helmholtz interior oblique impedance problem is unique in the following situations.

(i) $\Omega_-$ is Lipschitz and star-shaped with respect to a ball centred (without loss of generality) at the origin,
$k>0$,
$\bZ(\mathbf{x})= \mathbf{x}$, and
\begin{equation} \label{eq:special_eta}
\eta(\mathbf{x}) = k|\mathbf{x}| + {\rm i} \alpha \quad\text{ with }\quad \alpha \geq (d-1)/2.
\end{equation}

(ii) $\Omega_-$ is Lipschitz and is a  finite union of domains as in (i), $k>0$, and, on each connected part of $\Gamma$, where $\mathbf{x}_j$ is the point from which that connected component of $\Omega_-$ is star-shaped with respect to a ball, $\bZ(\mathbf{x})=\mathbf{x}-\mathbf{x}_j$ and
\begin{equation*}
\eta(\mathbf{x}) = k|\mathbf{x}-\mathbf{x}_j| + {\rm i} \alpha  \quad\text{ with }\quad \alpha \geq (d-1)/2.
\end{equation*}

(iii) $\Omega_-$ is $C^\infty$, $k$ is sufficiently large, $\eta= c_1 k$, where $c_1>0$ is independent of $k$, and $\bZ\in (C^\infty(\Gamma))^d$ is such that, for all $\mathbf{x} \in \Gamma$,
\begin{equation}\label{eq:GLS1}
\bZ(\mathbf{x})\cdot\mathbf{n}(\mathbf{x}) >0 \quad\text{ and } \quad |\bZ(\mathbf{x})\cdot \mathbf{t}(\mathbf{x}) - c_1 |>0
\end{equation}
for all unit tangent vectors $\mathbf{t}(\mathbf{x})$ at the point $\mathbf{x}$.
\end{theorem}

The following theorem is the Helmholtz analogue of Theorem \ref{thm:inverse_formula_L}.

\begin{theorem}\label{thm:inverse_formula_H}
Assume that the solution of the Helmholtz interior oblique impedance problem of Definition \ref{prob:ioip} exists and is unique, and let $P^{-, \eta, \bZ}_{\rm ItD}: L^2(\Gamma)\rightarrow H^1(\Gamma)$ denote the map that takes $g$ to $\gamma^- u$, where $u$ is as in Definition \ref{prob:ioip}.
Let $P^+_{\rm DtN}: H^1(\Gamma)\rightarrow {L^2(\Gamma)}$ denote the exterior Dirichlet-to-Neumman map for the Helmholtz exterior Dirichlet problem. Then, as an operator on ${L^2(\Gamma)}$,
\begin{equation}\label{eq:inverse_formula_H}
\big( A'_{k,\eta,\bZ}\big)^{-1} = \frac{1}{\bZ\cdot \mathbf{n}}I - \left( P^+_{\rm DtN} - \frac{1}{ \bZ\cdot \mathbf{n}}\big( {\rm i} \eta - \bZ\cdot\nabla_{\Gamma} \big) \right)
P^{-, \eta, \bZ}_{\rm ItD}.
\end{equation}
\end{theorem}

\begin{remark}[The star-combined operator]\label{rem:star}
When $\Omega_-$ is star-shaped with respect to a ball, one can choose $\bZ$ and $\eta$ so that $A^\prime_{k,\eta,\bZ}$ is coercive (not just coercive up to a compact perturbation). Indeed if
$\bZ(\mathbf{x})=\mathbf{x}$ and $\eta(\mathbf{x})$ is given by \eqref{eq:special_eta},
then \cite[Theorem 1.1]{SpChGrSm:11} showed that $A^\prime_{k,\eta,\bZ}$ is coercive for all $k>0$ with coercivity constant given by $\kappa$ in \eqref{eq:ssg}
\footnote{Strictly speaking, \cite[Theorem 1.1]{SpChGrSm:11} proves coercivity of $A^\prime_{k,\eta,\bZ}$ with $\eta$ as in \eqref{eq:special_eta} with $\alpha=(d-1)/2$. However, using the second inequality in \cite[Equation 2.8]{ChMo:08} in the proof of \cite[Lemma 2.4]{SpChGrSm:11}, one obtains the result for $\alpha \geq (d-1)/2$.
Furthermore, \cite{SpChGrSm:11} considers only $d=2,3$, but the proofs go through in exactly the same way for $d\geq 4$.
}. This operator was named the \emph{star-combined operator} and given the notation $\mathscr{A}_k$. Note that the coercivity of the star-combined operator implies the uniqueness result about the interior oblique impedance problem
in Part (i) of Theorem \eqref{thm:Hoblique}.
\end{remark}

\begin{remark} \label{rem:earlier}
Part (iv) of Theorem \ref{thm:Hedp} was stated (without proof) as \cite[Theorem 2.40]{ChGrLaSp:12}, with reference given to the present paper.
Lemma \ref{thm:star_in_oblique}, Corollary \ref{cor:star_ioip_equiv}, and Theorem \ref{thm:inverse_formula_H} were stated (with proof) as \cite[Theorem 2.38]{ChGrLaSp:12}, \cite[Corollary 2.39]{ChGrLaSp:12}, and \cite[Theorem 2.42]{ChGrLaSp:12} but are included here for completeness.
\end{remark}

\subsection{Discussion of the Helmholtz results and related literature}\label{sec:Helm_diss}

\paragraph{The ideas behind the results in \S\ref{sec:Helmholtz_state}.}
The proofs of Parts (i), (ii), and (iii) of Theorem \ref{thm:Hedp} are very similar to their Laplace analogues in Theorem \ref{thm:Laplace_ext}. The proof of Part (iv) (coercivity up to a compact perturbation on ${L^2(\Gamma)}$)
follows from Part (iv) of Theorem \ref{thm:Laplace_ext} (the analogous result for $A'_{E,\bZ,\alpha}$) since the difference $A'_{k,\eta,\bZ}- A'_{E,\bZ,\alpha}$ is compact (for any $\alpha\in \mathbb{R}$).
Indeed,
\begin{equation}\label{eq:Helm_diff}
A^\prime_{k,\eta,\bZ}- A'_{E,\bZ,\alpha}= (\bZ\cdot\mathbf{n})(D'_k-D') + \bZ \cdot \nabla_{\Gamma}(S_k-S) - {\rm i} \eta S_k - \alpha S.
\end{equation}
The bounds on $\Phi_k(\mathbf{x},\mathbf{y})- \Phi(\mathbf{x},\mathbf{y})$ in, e.g., \cite[Equation 2.25]{ChGrLaSp:12} imply that $D'_k-D'$ and $\nabla_{\Gamma}(S_k-S)$ map ${L^2(\Gamma)}$ continuously to $H^1(\Gamma)$, and are thus compact. Furthermore, the mapping property \eqref{eq:map2} holds with $S$ replaced by $S_k$ (again because of properties of $\Phi_k(\mathbf{x},\mathbf{y})- \Phi(\mathbf{x},\mathbf{y})$) and thus $S_k$ is compact on ${L^2(\Gamma)}$.

The link between $A'_{k,\eta,\bZ}$ and the Helmholtz interior oblique impedance problem expressed in Theorem \ref{thm:star_in_oblique} and Corollary \ref{cor:star_ioip_equiv} is proved in an analogous way to the link between $A'_{E,\bZ,\alpha}$ and the Laplace IORP in Theorem \ref{thm:star_in_obliqueL} and Corollary \ref{cor:oblique:unique}.

\paragraph{Discussion of our results in the context of related literature.}

The summary is that there does not yet exist a BIE posed in ${L^2(\Gamma)}$ for solving the Helmholtz EDP that, for all Lipschitz $\Omega_-$ and all $k>0$, is bounded, invertible, and the sum of a coercive operator and a compact operator. The standard BIOs $A'_{k,\eta}$ and $A_{k,\eta}$ \eqref{eq:opA} fail to be the sum of a coercive operator and a compact operator on the Lipschitz domains and 3-d star-shaped Lipschitz polyhedra given in \cite{ChSp:22}. The new BIOs $A_{k,\eta,\bZ}'$ and $A_{k,\eta,\bZ}$ have not been proved to be injective for all $k>0$. The formulations of  \cite{BuHi:05, BuSa:07, EnSt:07, EnSt:08} are invertible for all $k>0$ and the sum of a coercive operator and a compact operator, but only from ${H^{-1/2}(\Gamma)}\rightarrow {H^{1/2}(\Gamma)}$. We now give more detail on all these points.

The motivation for considering second-kind Laplace BIEs over first-kind Laplace BIEs is primarily based on the good conditioning of second-kind BIEs.
For Helmholtz BIEs, however, there is an additional complication compared to Laplace BIEs:~the BIEs in the exterior Helmholtz analogues of \eqref{eq:inteq1} and \eqref{eq:inteq2}, namely
\begin{equation}\label{eq:HBIE1}
S_k \partial_n^+ u = - \left(\frac{1}{2} I - D_k\right) g_D \quad\text{ and }\quad \left( \frac{1}{2} I + D'_k \right) \partial_n^+ u = H_k g_D,
\end{equation}
are not uniquely solvable for all values of $k$. Indeed,
$S_k$ is not injective when $k^2$ is an Dirichlet eigenvalue of $-\Delta$ in $\Omega_-$; this is because, from the first equation in \eqref{eq:HBIE1}, the Neumann trace of a Dirichlet eigenfunction of $-\Delta $ in $\Omega_-$ satisfies $S_k \partial_n^- u=0$.
Similarly, $\frac{1}{2} I + D'_k $ is not injective when $k^2$ is a Neumann eigenvalue of $-\Delta$ in $\Omega_-$; this is because
the operator $\frac{1}{2} I + D'_k $ appears in the indirect formulation of the interior Neumann problem (see Table \ref{tab:bies}) and thus cannot be injective when $k^2$ is a Neumann eigenvalue.

The standard remedy (going back to \cite{BrWe:65, Le:65, Pa:65} in the context of indirect BIEs and \cite{BuMi:71} in the context of direct BIEs) is to take a linear combination of the BIEs in \eqref{eq:HBIE1} and observe that, for $\eta\in \mathbb{C}\setminus\{0\}$, $\partial_n^+ u$ satisfies
\begin{equation}\label{eq:HBIE2}
A'_{k,\eta}\partial_n^+ u = B_{k,\eta} g_D
\end{equation}
with $A'_{k,\eta}$ and $B_{k,\eta}$ defined in \eqref{eq:opA}.
The mapping properties \eqref{eq:map2} imply that, for $|s|\leq 1/2$,
\begin{equation*}
A'_{k,\eta}: H^{s-1/2}(\Gamma)\rightarrow H^{s-1/2}(\Gamma) \quad\text{ and }\quad B_{k,\eta}: H^{s+1/2}(\Gamma)\rightarrow H^{s-1/2}(\Gamma),
\end{equation*}
 and one can then show that, for $k>0$ and $\Re\eta \neq 0$, $A'_{k,\eta}$ and $B_{k,\eta}$ are invertible; see, e.g, \cite[Theorem 2.27]{ChGrLaSp:12}.

Since $D'_k-D'$ and $S_k$ are both compact on ${L^2(\Gamma)}$ when $\Gamma$ is Lipschitz,
$A'_{k,\eta}$ is the sum of $\frac{1}{2} I + D'$ and a compact operator. Therefore the question of whether or not $A'_{k,\eta}$ is the sum of a coercive operator and a compact operator is equivalent to the analogous question for $\frac{1}{2} I +D'$. The results recalled in \S\ref{sec:openbook} imply that the answer to this question is yes if $\Gamma$ is $C^1$, $\Omega_-$ is a 2-d curvilinear polygon with $C^{1,\alpha}$ sides, or $\Gamma$ is a Lipschitz domain with sufficiently-small Lipschitz character; furthermore the answer to this question is no if $\Gamma$ is one of the 2- and 3-d Lipschitz domains or 3-d star-shaped Lipschitz polyhedra described in \cite{ChSp:22} (see \cite[\S1.2.1]{ChSp:22} for more detail).

The lack of a convergence theory for the Galerkin method applied to either $A'_{k,\eta}$ or $A_{k,\eta}$ on ${L^2(\Gamma)}$
motivated \cite{BuHi:05,BuSa:07, EnSt:07, EnSt:08} to introduce modifications of $A'_{k,\eta}$ and $A_{k,\eta}$.
For direct BIEs, the idea in these papers is to choose an operator $M: {H^{-1/2}(\Gamma)}\rightarrow {H^{1/2}(\Gamma)}$ so that
\begin{equation*}
\widetilde{A}_{k,\eta, M}:= M\left(\frac{1}{2} I + D'_k\right) - {\rm i} \eta S_k
\end{equation*}
maps ${H^{-1/2}(\Gamma)} \to {H^{1/2}(\Gamma)}$, is injective, and is the sum of a coercive operator and a compact operator (with this last property following from the facts that $S:{H^{-1/2}(\Gamma)}\rightarrow {H^{1/2}(\Gamma)}$ is coercive and $\widetilde{A}_{k,\eta, M}+{\rm i} \eta S_k$ is compact);
different choices of $M$ were then proposed and analysed in \cite{BuHi:05, BuSa:07, EnSt:07, EnSt:08}, and the corresponding interior boundary value problems proved to have a unique solution (so that the operator $\widetilde{A}_{k,\eta, M}$ is injective).
Convergence of the Galerkin method then follows from Part (b) of Theorem \ref{thm:Galerkin}. However, the BIE involving $\widetilde{A}_{k,\eta, M}$ is now a first-kind equation (in contrast to the second-kind equation \eqref{eq:HBIE2}), and so one expects the conditioning of the Galerkin linear systems to worsen as the meshes are refined (as discussed in \S\ref{sec:rationale}).

\subsection{Proofs of the results in \S\ref{sec:Helmholtz_state}} \label{sec:Helmholtz_prove}

\begin{proof}[Proof of Theorem \ref{thm:Hedp}]
The proof that the Neumann trace $\partial_n^+ u$ of the solution of the Helmholtz EDP satisfies the integral equation \eqref{eq:Hnew_ext} is very similar to the
analogous arguments for the Laplace integral equations; see the proofs of Theorems \ref{thm:Laplace_int} and \ref{thm:Laplace_ext} above.
As recalled in \S\ref{sec:Helm_diss}, the difference $A'_{k,\eta,\bZ}- A'_{E,\bZ,\alpha}$ is compact;
therefore, Part (iv) of Theorem \ref{thm:Laplace_ext} implies that, for any $\eta\in \mathbb{C}$, $A^\prime_{k,\eta,\bZ}$ is the sum of a coercive operator and a compact operator.
\end{proof}

\

\begin{proof}[Proof of Lemma \ref{lem:analyticFredholm}]
By the definitions of $A'_{E,\bZ,\alpha}$ \eqref{eq:opLEZ} and $A^\prime_{k,\eta,\bZ}$ \eqref{eq:AketaZpr_def},
\begin{align*}
A^\prime_{k,\eta,\bZ}&= A'_{E,\bZ,\alpha} + \big(K_{\bZ,k}'- K_\bZ\big) - {\rm i} \eta S_k - \alpha S,\\
&=A'_{E,\bZ,\alpha} \Big( I + (A'_{E,\bZ,\alpha})^{-1}\big( K_{\bZ,k}'- K_\bZ - {\rm i} \eta S_k - \alpha S\big)\Big).
\end{align*}
By Part (v) of Theorem \ref{thm:Laplace_ext}, $A'_{E,\bZ,\alpha}$ is invertible if $\alpha$ satisfies \eqref{eq:alpha_bound_boundary}. Furthermore, for all $k\in \mathbb{C}$ ($d\geq 3$ and odd) or $k\in \mathbb{C}\setminus \mathbb{R}^-$ ($d$ even),
$K_{\bZ,k}'- K_\bZ - {\rm i} \eta S_k - \alpha S$ is compact on ${L^2(\Gamma)}$; indeed, as discussed in the proof of Theorem \ref{thm:Hedp} above,
$S_k, S_0, D_k'-D',$ and $\nabla_\Gamma S_k  -\nabla_\Gamma S$ are all compact operators on ${L^2(\Gamma)}$. The result then follows from  analytic Fredholm theory; see, e.g., \cite[Theorem 8.26]{CoKr:98}.
\end{proof}

\

\begin{proof}[Proof of Theorem \ref{thm:star_in_oblique}]
The proof is similar to the proofs of Theorems \ref{thm:star_in_obliqueL} and \ref{thm:star_in_obliqueLE}.
If $u={\cal S}_k \phi$ with $\phi\in L^2(\Gamma)$ then, by the jump relations \eqref{eq:jump1}, the fact that $S_k:{L^2(\Gamma)}\rightarrow H^1(\Gamma)$ by \eqref{eq:map1} and the fact that ${\cal S}_k$ satisfies the same mapping properties as ${\cal S}$ in \eqref{eq:LPmap}, $u$ satisfies the oblique impedance BVP if and only if the boundary condition \eqref{eq:oblique} holds.
However, the jump relations \eqref{eq:jump1} imply that this equation is precisely \eqref{eq:BIE_ioip1}. On the other hand, if $u$ satisfies the oblique impedance problem then
$g_0 := \partial_n^- u - {\rm i}\eta\gamma^- u\in L^2(\Gamma)$. Recall the interior impedance BVP
\begin{equation*}
\Delta u +k^2 u =0 \quad\text{ in } \Omega_-, \qquad\partial_n^- u- {\rm i} \eta \gamma^- u = g \quad\text{ on }\Gamma,
\end{equation*}
 with $\eta\in\mathbb{R}\setminus\{0\}$ and $g=g_0$ has exactly one solution (see, e.g., \cite[Theorem 2.3]{ChGrLaSp:12}). Clearly this solution is $u$. Further, recalling first that $A^\prime_{k,k,\mathbf{n}}=A^\prime_{k,k}$, which is invertible by \cite[Theorem 2.27]{ChGrLaSp:12}, and second that the interior impedance problem is the special case of the interior oblique impedance problem when $\bZ=\mathbf{n}$, we see that $u = {\cal S}_k \phi$ with $\phi = (A^\prime_{k,k})^{-1}g_0\in L^2(\Gamma)$ by the first part of this theorem.
\end{proof}

\

\begin{proof}[Proof of Corollary \ref{cor:star_ioip_equiv}]
Surjectivity follows from Lemma~\ref{thm:star_in_oblique}.
Injectivity follows
provided that
$u={\cal S}_k\phi=0$ in $\Omega_-$ only if $\phi=0$.
Assume that $\phi\in L^2(\Gamma)$ and $u={\cal S}_k\phi$ is zero in $\Omega_-$. Then, by the jump relations \eqref{eq:jump1}, $u$ satisfies the Helmholtz exterior Dirichlet problem in $\Oe$ with zero Dirichlet data, and hence is zero by uniqueness of the solution of this BVP (see, e.g., \cite[Corollary 2.9]{ChGrLaSp:12}). The jump relations \eqref{eq:jump1} then imply that  $\phi = \partial_n^-u-\partial_n^+u=0$, as required.
\end{proof}

\

\begin{proof}[Proof of Theorem \ref{thm:Hoblique}]

(i) In this situation, $A^\prime_{k,\eta,\bZ}$ is coercive (not just coercive up to a compact perturbation) by \cite[Theorem 1.1]{SpChGrSm:11};
see Remark \ref{rem:star}.

(ii) This is proved in \cite[Theorem 5.6]{Gi:17} by showing that the corresponding $A^\prime_{k,\eta,\bZ}$ is injective;
for related results about this $A^\prime_{k,\eta,\bZ}$ see \cite[Appendix A]{GiChLaMo:21}.

(iii)
Let $\hbar:=k^{-1}$ be the semiclassical parameter.
In the notation of \cite[\S4]{GaLaSp:21}, the boundary condition \eqref{eq:oblique} then becomes
\begin{equation*}
\mathcal{N}\left(\frac{1}{{\rm i}} \hbar\partial_n u\right) - \mathcal{D} u = \frac{1}{{\rm i}} \hbar g,
\end{equation*}
where $\mathcal{N}:= \bZ \cdot \mathbf{n}$ and $\mathcal{D}= - \bZ \cdot\nabla_\Gamma/{\rm i} + c_1$; see \cite[Equation 4.1]{GaLaSp:21}.
We now check that the hypotheses of \cite[Theorem 4.6]{GaLaSp:21} hold, and then the result then follows from \cite[Theorem 4.6]{GaLaSp:21} (in fact, this theorem proves that, under these conditions, the interior oblique derivative problem is well-posed, and the bound on the solution in terms of the data has the same $k$ dependence as the bound for the standard interior impedance problem).
In the notation of \cite{GaLaSp:21}, $\sigma(\mathcal{N})=  \bZ \cdot \mathbf{n}$, $\sigma(\mathcal{D})= - \bZ\cdot \xi'+ c_1$, $m_1=0$, and $m_0=1$.
The first/second equation in \cite[Equation 4.2]{GaLaSp:21} holds because of the first/second equation in \eqref{eq:GLS1}, respectively.
Then, \cite[Equations 4.3, 4.10, 4.11, and 4.12]{GaLaSp:21} all hold because $m_0=m_1+1$.
\end{proof}

\

\begin{proof}[Proof of Theorem \ref{thm:inverse_formula_H}]
This is very similar to the proof of Theorem \ref{thm:inverse_formula_L} (the Laplace analogue), using
Lemma \ref{thm:star_in_oblique} in place of  Theorem \ref{thm:star_in_obliqueLE}.
\end{proof}

\begin{appendix}

\section{The vector field $\mathbf{Z}$}\label{app:Lipschitz}

This appendix contains (i) the construction of a vector field $\mathbf{Z}$ satisfying \eqref{eq:c} (\S\ref{sec:Z}) and (ii) the proof of the extension theorem Lemma \ref{lem:Lipschitz} (\S\ref{sec:Z2}).

\subsection{Constructing the vector field $\widetilde{\bZ}$}\label{sec:Z}

\begin{lemma}\label{lem:vf}
If $\Omega_-$ is a bounded Lipschitz open set then there exists a real-valued $\widetilde{\bZ}\in (C_{\rm{comp}}^\infty(\mathbb{R}^d))^d$ and $c>0$ such that $\bZ:= \widetilde{\bZ}|_\Gamma$ satisfies
\begin{equation}\label{eq:c2}
\bZ(\mathbf{x}) \cdot \mathbf{n}(\mathbf{x}) \geq c \quad \text{ for almost every } \mathbf{x}\in \Gamma.
\end{equation}
\end{lemma}
The construction of such a $\widetilde{\bZ}$ can be found in, e.g., \cite[Lemma 1.5.1.9]{Gr:85}, \cite[Proof of Lemma 1.3]{Ne:67}.

We now spell out the construction from \cite{Gr:85} (in slightly modified form), calculate bounds on the $\widetilde{\bZ}$ that we construct, and then show how both the constant $\alpha$ (required to satisfy \eqref{eq:alpha_bound_boundary}/\eqref{eq:alpha_bound}) and $\|A'_{I,\bZ,\alpha}\|_{{\LtG\rightarrow \LtG}}$ depend on constants in the construction of $\widetilde{\bZ}$ (see \eqref{eq_alphadef} and \eqref{eq:Abound} below).

Suppose that $\mathbf{x}_m\in \mathbb{R}^d$ and $a_m>0$, for $m=1, \ldots , M$, are chosen so that
\begin{equation} \label{eq:cover}
\Gamma \subset \bigcup_{m=1}^M B_{a_m}(\mathbf{x}_m),
\end{equation}
and so that, for each $m=1,...,M$, $\Gamma_m:= \Gamma\cap B_{a_m}(\mathbf{x}_m)$ is, in some rotated coordinate system $0\tilde x_1...\tilde x_d$, the graph of a Lipschitz continuous function with Lipschitz constant $L_m$ such that the unit vector $\widetilde{\bZ}_m$ in the $\tilde x_d$ direction points out of $\Omega$. It is clear from the definition of a Lipschitz open set (see, e.g., \cite[Definition 3.28]{Mc:00}) that it is always possible to find such $\mathbf{x}_m$ and $a_m$. Note that, for almost all $\mathbf{x}\in \Gamma_m$,
\begin{equation} \label{eq:Zn}
\mathbf{n}(\mathbf{x})\cdot \widetilde{\bZ}_m \geq c_m := (1+L_m^2)^{-1/2}.
\end{equation}

Suppose that $\theta_m\in C^{0,1}(\mathbb{R}^d)$  with $\mbox{supp}(\theta_m)\subset B_{a_m}(\mathbf{x}_m)$, for $m=1,...,M$, and that
\begin{equation} \label{eq:theta}
\sum_{m=1}^M \theta_m(\mathbf{x}) \geq 1, \quad \mathbf{x}\in \Gamma.
\end{equation}
Functions $\theta_1,...,\theta_M$ satisfying \eqref{eq:theta} exist by, e.g., \cite[Theorem 2.17]{Gr:09}, and indeed can be chosen so that equality rather than inequality holds in \eqref{eq:theta} and so that each $\theta_m\in C_{\rm comp}^\infty(\mathbb{R}^d)$, in which case we say that $(\theta_1,...,\theta_M)$ is a partition of unity for $\Gamma$ subordinate to the cover $(B_{a_1}(\mathbf{x}_1),...,B_{a_M}(\mathbf{x}_m))$. We construct explicitly functions satisfying \eqref{eq:theta} below; a choice of $\widetilde{\bZ}$ that satisfies \eqref{eq:c2} is then
\begin{equation} \label{eq:Zdef}
\widetilde{\bZ}(\mathbf{x}) := \sum_{m=1}^M \theta_m(\mathbf{x}) \widetilde{\bZ}_m,
\end{equation}
since, for almost all $\mathbf{x}\in \Gamma$,
\begin{equation}\label{eq:Zc}
\mathbf{n}(\mathbf{x}) \cdot \widetilde{\bZ}(\mathbf{x}) = \sum_{m=1}^M \theta_m(\mathbf{x}) \mathbf{n}(\mathbf{x}) \cdot \widetilde{\bZ}_m \geq \sum_{m=1}^M c_m\theta_m(\mathbf{x}) \geq c := \min_{m=1,...,M} c_m.
\end{equation}
Observe that, since $\widetilde{\bZ}_m$, $m=1, \ldots, M$, are unit vectors, for $\mathbf{x}\in\mathbb{R}^d$,
$$
|\widetilde{\bZ}(\mathbf{x})| \leq \Theta(\mathbf{x}) := \sum_{m=1}^M \theta_m(\mathbf{x}),
\quad\text{ and also } \quad
\nabla \cdot \widetilde{\bZ}(\mathbf{x}) = \sum_{m=1}^M \nabla \theta_m(\mathbf{x}) \cdot \widetilde{\bZ}_m,
$$
so that
\begin{equation}\label{eq:noise}
|\nabla \cdot \widetilde{\bZ}(\mathbf{x})| \leq \Theta^\prime(\mathbf{x}) := \sum_{m=1}^M |\nabla \theta_m(\mathbf{x})|.
\end{equation}
Furthermore, for unit vectors $\mathbf{a}$ and ${\bf b}$, writing $\widetilde{\bZ}_m$ in terms of its components as $\widetilde{\bZ}_m=(Z_{m,1},...,Z_{m,d})$,
$$
\partial_j Z_i a_ib_j = \sum_{m=1}^M \partial_j \theta_m Z_{m,i} a_ib_j = \sum_{m=1}^M (\nabla \theta_m \cdot {\bf b}) \; (\widetilde{\bZ}_m \cdot \mathbf{a}),
\quad\text{so that}\quad
|\partial_j Z_i a_ib_j |\leq \Theta^\prime.
$$
Therefore, with $D \widetilde{\bZ}$ the matrix with $(i,j)$th element $\partial_i Z_j$ and $\|\cdot\|_2$ denoting the matrix $2$-norm,
$
\|D \widetilde{\bZ}\|_2 \leq \Theta^\prime.
$
Combining this with \eqref{eq:noise}, we see that
the inequality for $\alpha$ \eqref{eq:alpha_bound}
holds if
\begin{equation} \label{eq_alpha_bound2}
2\alpha \geq 3 \|\Theta^\prime\|_{L^\infty(\mathbb{R}^d)}.
\end{equation}

We now choose a specific form for the partition-of-unity functions $\theta_m$, and hence obtain a bound on
$\|\Theta^\prime\|_{L^\infty(\mathbb{R}^d)}.$
The right hand side of \eqref{eq:cover} is still a cover for $\Gamma$ if we reduce the size of the balls slightly, i.e., \eqref{eq:cover} implies that, for some  $\mu \in (0,1)$,
\begin{equation} \label{eq:cover1}
\Gamma \subset \bigcup_{m=1}^M B_{\mu a_m}(\mathbf{x}_m).
\end{equation}
Assuming that \eqref{eq:cover1} holds for some $\mu\in (0,1)$, one way of constructing the functions $\theta_m$ to satisfy \eqref{eq:theta} is to set
\begin{equation} \label{eq:thetdef}
\theta_m(\mathbf{x}) := \chi\left(\frac{|\mathbf{x}-\mathbf{x}_m|}{a_m}\right)
\end{equation}
with $\chi\in C^{0,1}[0,\infty)$ chosen so that $0\leq \chi(t)\leq 1$, for $t\geq 0$, $\chi(t)=1$, for $0\leq t\leq \mu$, while $\chi(t)=0$ for $t\geq 1$.
(Observe that, with these choices of $\theta_m$, $\widetilde{\bZ}$ is no longer smooth, but only Lipschitz.)
Then, for $\mathbf{x}\in \mathbb{R}^d$,
\begin{equation} \label{eq:zbound}
|\widetilde{\bZ}(\mathbf{x})| \leq M^*,
\end{equation}
where $M^*\in \mathbb{N}$ is the smallest integer such that every $\mathbf{x}\in \Gamma$ is in at most $M^*$ balls $B_{a_m}(\mathbf{x}_m)$. Furthermore, for $\mathbf{x}\in \mathbb{R}^d$,
$$
\nabla \theta_m(\mathbf{x}) = \chi^\prime\left(\frac{|\mathbf{x}-\mathbf{x}_m|}{a_m}\right) \, \frac{\mathbf{x}-\mathbf{x}_m}{a_m|\mathbf{x}-\mathbf{x}_m|},
$$
so that, if
$$
a := \min_{m=1,...,M} a_m,
$$
then, by \eqref{eq:noise},
\begin{equation}\label{eq:Thetaprimebound}
\|\Theta^\prime\|_{L^\infty(\mathbb{R}^d)} \leq \frac{M^* \|\chi^\prime\|_\infty}{a}.
\end{equation}
 Thus \eqref{eq_alpha_bound2} holds if
\begin{equation} \label{eq_alpha_bound3}
2\alpha \geq  \frac{3M^*\|\chi^\prime\|_\infty}{a}, \quad \mbox{i.e., \,if} \quad 2\alpha \geq \frac{3M^*}{(1-\mu)a}
\end{equation}
if we make the simple choice that
\begin{equation} \label{eq:chidef}
\chi(t) = (1-t)/(1-\mu) \quad \mbox{for} \quad  \mu\leq t\leq 1.
\end{equation}
With $A'_{I,\bZ,\alpha}$ defined by \eqref{eq:opLIZ}, assume that
\begin{equation} \label{eq_alphadef}
\alpha :=  \frac{3M^*\|\chi^\prime\|_\infty}{2a}.
\end{equation}
Then, by Theorem \ref{thm:Laplace_int} combined with \eqref{eq_alpha_bound2} and \eqref{eq:Thetaprimebound}, $A'_{I,\bZ,\alpha}$ is coercive with coercivity constant $c$ given by \eqref{eq:Zc}. Furthermore, by \eqref{eq:zbound},
$A'_{I,\bZ,\alpha}$ is bounded with
\begin{equation} \label{eq:Abound}
\N{A'_{I,\bZ,\alpha}}_{{\LtG\rightarrow \LtG}} \leq M^* \left(\frac{1}{2} + \N{D'}_{{\LtG\rightarrow \LtG}} + \N{\nabla_\Gamma S}_{{\LtG\rightarrow \LtG}}\right) +  \frac{3\|\chi^\prime\|_\infty}{2a} \N{S}_{{\LtG\rightarrow \LtG}}.
\end{equation}

\subsection{Proof of Lemma \ref{lem:Lipschitz}}\label{sec:Z2}

\begin{proof}
By the Kirszbraun theorem \cite{Ki:34, Va:45},
$\bZ\in (C^{0,1}(\Gamma))^d$ can be extended to a function $\bZ_{\rm ext}\in (C^{0,1}(\mathbb{R}^d))^d$ with the same (non-zero) Lipschitz constant.

Let $R>0$ be such that $\overline{\Omega_-}\subset B_R$. For $a>1$ and $R^*>aR$, let
\begin{equation*}
\chi(r):=
\begin{cases}
1, & R<r <aR,\\
(R^*-r)/(R^*-aR), & aR<r<R^*,\\
0, & r> R^*,
\end{cases}
\end{equation*}
and let
\begin{equation*}
\widetilde{\bZ}_{\rm ext} (\mathbf{x}):=
\begin{cases}
\bZ_{\rm ext}(\mathbf{x}), & |\mathbf{x}|\leq R,\\
\bZ_{\rm ext}(R\widehat{\mathbf{x}})\,\chi(|\mathbf{x}|), & |\mathbf{x}|> R,
\end{cases}
\end{equation*}
where $\widehat \mathbf{x}:= \mathbf{x}/|\mathbf{x}|$, so that $\mathrm{supp}\, \widetilde{\bZ}_{\rm ext} = B_{R^*}$.

Let $L$ be the Lipschitz constant of $\bZ_{\rm ext}$.  We now prove that if $a$ and $R^*$ are both large enough then
\begin{equation}\label{eq:Lip}
\big|\widetilde{\bZ}_{\rm ext}(\mathbf{x}) - \widetilde{\bZ}_{\rm ext}(\mathbf{y})\big|\leq L |\mathbf{x}-\mathbf{y}| \quad\text{ for all } \mathbf{x},\mathbf{y} \in \mathbb{R}^d,
\end{equation}
i.e.~the Lipschitz constant of $\widetilde{\bZ}_{\rm ext}$ does not exceed that of $\bZ_{\rm ext}$.

First observe that there exists $a_0>1$ such that if $a\geq a_0$ and at least one of $|\mathbf{x}|$ and $|\mathbf{y}|$ are $\geq aR$ then
\begin{equation}\label{eq:largea}
|R\widehat{\mathbf{x}}-R\widehat{\mathbf{y}}| \leq \frac{1}{2} |\mathbf{x}-\mathbf{y}|.
\end{equation}
We now prove that if $a\geq a_0>1$ and
\begin{equation}\label{eq:R*}
R^*\geq aR + 2 \N{\bZ_{\rm ext}}_{L^\infty(B_R)}/L.
\end{equation}
then \eqref{eq:Lip} holds.

\paragraph{Case 1: $|\mathbf{x}|\leq R, |\mathbf{y}|\leq R$.} Since $\widetilde{\bZ}_{\rm ext}=\bZ_{\rm ext}$ in $B_R$, \eqref{eq:Lip} follows from the definition of $L$.

\paragraph{Case 2: $|\mathbf{x}|\leq R, R\leq |\mathbf{y}|\leq aR$.} The key point here is that $|\mathbf{x}-\mathbf{y}|\geq |\mathbf{x}- R\widehat{\mathbf{y}}|$ so that
\begin{equation*}
\big|\widetilde{\bZ}_{\rm ext}(\mathbf{x}) - \widetilde{\bZ}_{\rm ext}(\mathbf{y})\big|= \big|\bZ_{\rm ext}(\mathbf{x}) - \bZ_{\rm ext}(R\widehat{\mathbf{y}})\big|\leq L |\mathbf{x}- R\widehat{\mathbf{y}}|\leq L |\mathbf{x}-\mathbf{y}|.
\end{equation*}

\paragraph{Case 3: $|\mathbf{x}|\leq R, aR\leq |\mathbf{y}|\leq R^*$.}
Now
\begin{align*}
\big|\widetilde{\bZ}_{\rm ext}(\mathbf{x}) - \widetilde{\bZ}_{\rm ext}(\mathbf{y})\big|&= \big|\bZ_{\rm ext}(\mathbf{x}) - \bZ_{\rm ext}(R\widehat{\mathbf{y}})\chi(|\mathbf{y}|)\big|\\
&\leq \big|\bZ_{\rm ext}(\mathbf{x}) - \bZ_{\rm ext}(R\widehat{\mathbf{y}})\big| +\big|\bZ_{\rm ext}(R\widehat{\mathbf{y}})\big|(1- \chi(\mathbf{y}))\\
&\leq L\big|\mathbf{x}-R\widehat{\mathbf{y}}\big| +\N{\bZ_{\rm ext}}_{L^\infty(B_R)}\frac{ |\mathbf{y}|-aR}{R^*-aR}.
\end{align*}
If $a\geq a_0$, then
\begin{equation*}
\big|\mathbf{x}-R\widehat{\mathbf{y}}\big|\leq\big|R\widehat{\mathbf{x}}-R\widehat{\mathbf{y}}\big| \leq \frac{1}{2} |\mathbf{x}-\mathbf{y}\big|,
\end{equation*}
by \eqref{eq:largea}. Since $|\mathbf{y}|-aR \leq |\mathbf{x}-\mathbf{y}|$,  \eqref{eq:Lip} then follows if, with this choice of $a$,  $R^*$ is then large enough so that
\begin{equation}\label{eq:Rlarge}
\frac{\N{\bZ_{\rm ext}}_{L^\infty(B_R)}}{R^*-aR}\leq \frac{L}{2},
\end{equation}
which is equivalent to the inequality \eqref{eq:R*}.

\paragraph{Case 4: $|\mathbf{x}|\leq R, |\mathbf{y}|\geq R^*$.}
Now
\begin{equation*}
\big|\widetilde{\bZ}_{\rm ext}(\mathbf{x})- \widetilde{\bZ}_{\rm ext}(\mathbf{y})\big|= \big|\bZ_{\rm ext}(\mathbf{x})\big| \leq \N{\bZ_{\rm ext}}_{L^\infty(B_R)}.
\end{equation*}
Since $L(R^*-R)\leq L|\mathbf{x}-\mathbf{y}|$, \eqref{eq:Lip} follows if
\begin{equation}\label{eq:R*1}
\N{\bZ_{\rm ext}}_{L^\infty(B_R)}\leq L(R^*-R)\quad\text{ i.e.,} \quad
R^* \geq R + \N{\bZ_{\rm ext}}_{L^\infty(B_R)}/L,
\end{equation}
which is ensured by \eqref{eq:R*} since $a\geq a_0>1$.

\paragraph{Case 5: $R\leq |\mathbf{x}|\leq aR, R\leq |\mathbf{y}|\leq aR$.}

Similar to Case 2, the key point here is that $|\mathbf{x}-\mathbf{y}|\geq |R\widehat{\mathbf{x}}- R\widehat{\mathbf{y}}|$ so that
\begin{equation*}
\big|\widetilde{\bZ}_{\rm ext}(\mathbf{x}) - \widetilde{\bZ}_{\rm ext}(\mathbf{y})\big|= \big|\bZ_{\rm ext}(R\widehat{\mathbf{x}}) - \bZ_{\rm ext}(R\widehat{\mathbf{y}})\big|\leq L |R\widehat{\mathbf{x}}- R\widehat{\mathbf{y}}|\leq L |\mathbf{x}-\mathbf{y}|.
\end{equation*}

\paragraph{Case 6: $R\leq |\mathbf{x}|\leq aR, aR\leq |\mathbf{y}|\leq R^*$.}
Now
\begin{align*}
\big|\widetilde{\bZ}_{\rm ext}(\mathbf{x}) - \widetilde{\bZ}_{\rm ext}(\mathbf{y})\big|&= \big|\bZ_{\rm ext}(R\widehat{\mathbf{x}}) - \bZ_{\rm ext}(R\widehat{\mathbf{y}})\chi(|\mathbf{y}|)\big|\\
&\leq \big|\bZ_{\rm ext}(R\widehat{\mathbf{x}}) - \bZ_{\rm ext}(R\widehat{\mathbf{y}})\big| +\big|\bZ_{\rm ext}(R\widehat{\mathbf{y}})\big|(1- \chi(\mathbf{y}))\\
&\leq L\big|R\widehat{\mathbf{x}}-R\widehat{\mathbf{y}}\big| +\N{\bZ_{\rm ext}}_{L^\infty(B_R)}\frac{|\mathbf{y}|-aR}{R^*-aR}.
\end{align*}
Combining \eqref{eq:largea} and the inequality $|\mathbf{y}|-aR\leq |\mathbf{x}-\mathbf{y}|$, we see that \eqref{eq:Lip} holds if $R^*$ is chosen so that \eqref{eq:Rlarge} holds, which is equivalent to the inequality \eqref{eq:R*}.

\paragraph{Case 7: $R\leq |\mathbf{x}|\leq aR,|\mathbf{y}|\geq R^*$.}
Similarly  to Case 4,
\begin{equation*}
\big|\widetilde{\bZ}_{\rm ext}(\mathbf{x})- \widetilde{\bZ}_{\rm ext}(\mathbf{y})\big|= \big|\bZ_{\rm ext}(R\widehat{\mathbf{x}})\big| \leq \N{\bZ_{\rm ext}}_{L^\infty(B_R)}.
\end{equation*}
Since $L(R^*-aR)\leq L|\mathbf{x}-\mathbf{y}|$, \eqref{eq:Lip} follows if
\begin{equation*}
\N{\bZ_{\rm ext}}_{L^\infty(B_R)}\leq L(R^*-aR),
\end{equation*}
which is ensured by \eqref{eq:R*}.

\paragraph{Case 8: $aR\leq |\mathbf{x}|\leq R^*, aR\leq |\mathbf{y}|\leq R^*$.}
Now
\begin{align*}
\big|\widetilde{\bZ}_{\rm ext}(\mathbf{x}) - \widetilde{\bZ}_{\rm ext}(\mathbf{y})\big|&= \big|\bZ_{\rm ext}(R\widehat{\mathbf{x}})\chi(|\mathbf{x}|) - \bZ_{\rm ext}(R\widehat{\mathbf{y}})\chi(|\mathbf{y}|)\big|\\
&\leq \big|\bZ_{\rm ext}(R\widehat{\mathbf{x}}) - \bZ_{\rm ext}(R\widehat{\mathbf{y}})\big||\chi(|\mathbf{x}|) | +\big|\bZ_{\rm ext}(R\widehat{\mathbf{y}})\big|\big|\chi(|\mathbf{x}|)- \chi(|\mathbf{y}|)\big|\\
&\leq L\big|R\widehat{\mathbf{x}}-R\widehat{\mathbf{y}}\big| +\N{\bZ_{\rm ext}}_{L^\infty(B_R)}\frac{|\mathbf{x}-\mathbf{y}|}{R^*-aR}.
\end{align*}
Using \eqref{eq:largea}, we see that \eqref{eq:Lip} holds if $R^*$ is chosen so that \eqref{eq:Rlarge} holds, which is equivalent to the inequality \eqref{eq:R*}.

\paragraph{Case 9: $aR\leq |\mathbf{x}|\leq R^*, |\mathbf{y}|\geq R^*$.} Now
\begin{equation*}
\big|\widetilde{\bZ}_{\rm ext}(\mathbf{x})- \widetilde{\bZ}_{\rm ext}(\mathbf{y})\big|= \big|\widetilde{\bZ}_{\rm ext}(\mathbf{x})\big|=\big|\bZ_{\rm ext}(R\widehat{\mathbf{x}})\big|\frac{R^*-|\mathbf{x}|}{R^*-aR}\leq \N{\bZ_{\rm ext}}_{L^\infty(B_R)}\frac{R^*-|\mathbf{x}|}{R^*-aR}.
\end{equation*}
Since $(R^*-|\mathbf{x}|)\leq |\mathbf{x}-\mathbf{y}|$, it is sufficient to prove that the right-hand side of this last expression is $\leq L(R^*-|\mathbf{x}|)$,
and this is ensured by  the inequality \eqref{eq:R*}.

\paragraph{Case 10: $|\mathbf{x}|\geq R^*, |\mathbf{y}|\geq R^*$.} In this case, \eqref{eq:Lip} holds since $\widetilde{\bZ}_{\rm ext}(\mathbf{x})=\widetilde{\bZ}_{\rm ext}(\mathbf{y})=\mathbf{0}$.
\end{proof}

\section{Recap of mapping properties of layer potentials and boundary integral operators}\label{app:B}

Recall that the single-layer potential ${\cal S}\phi$ is defined by \eqref{eq:SLP}. For $\phi\in L^2(\Gamma)$, define the double-layer potential ${\cal D}\phi$ by
\begin{equation}\label{eq:DLP}
{\cal D}\phi(\mathbf{x}):= \int_\Gamma \pdiff{\Phi(\mathbf{x},\mathbf{y})}{n(\mathbf{y})}\phi(\mathbf{y})\, {\rm d} s(\mathbf{y}) \quad\text{ for } \mathbf{x} \in \mathbb{R}^d\setminus\Gamma.
\end{equation}
For $\chi\in C^\infty_{\rm comp}(\mathbb{R}^d$) and $|s|\leq 1/2$,
\begin{equation}\label{eq:LPmap}
\chi {\cal S} : H^{s-1/2}(\Gamma)\rightarrow H^{s+1}(\mathbb{R}^d) \quad\text{ and }\quad
\chi {\cal D} : H^{s+1/2}(\Gamma)\rightarrow H^{s+1}(\Omega^{\pm}).
\end{equation}
With $S,D,D',$ and $H$ defined by \eqref{eq:bio1} and \eqref{eq:bio2},
for all $|s|\leq 1/2$,
\begin{subequations}\label{eq:map}
\begin{align}\label{eq:map1}
S: H^{s-1/2}(\Gamma)\rightarrow H^{s+1/2}(\Gamma), \quad D: H^{s+1/2}(\Gamma)\rightarrow H^{s+1/2}(\Gamma),\\
D': H^{s-1/2}(\Gamma)\rightarrow H^{s-1/2}(\Gamma), \quad H: H^{s+1/2}(\Gamma)\rightarrow H^{s-1/2}(\Gamma).\label{eq:map2}
\end{align}
\end{subequations}
The results in \eqref{eq:LPmap} and \eqref{eq:map} for $|s|=1/2$ (which then imply the results for $|s|<1/2$ by interpolation) are consequences of the results in \cite{CoMcMe:82}, \cite{Ve:84}, and \cite{JeKe:95}; see, e.g., \cite[Theorems 2.15 and 2.16 and Corollary A.8]{ChGrLaSp:12}. (Note that the results in \eqref{eq:LPmap} for $|s|<1/2$ can also be obtained from mapping properties of the Newtonian potential and Green's integral representation, with the results in \eqref{eq:map} then following from results about the trace map; see \cite{Co:88}, \cite[Theorem 6.11]{Mc:00}.)

\begin{lemma}[Invertibility of $S:{L^2(\Gamma)}\rightarrow H^1(\Gamma)$ when $\Gamma$ is Lipschitz]
 \label{lem:SLP}
If $\Gamma$ is Lipschitz and either $d=3$, or $d=2$ and $a\neq\mathrm{Cap}_\Gamma$, then $S:L^2(\Gamma)\to H^{1}(\Gamma)$ is bounded and invertible.
\end{lemma}

\begin{proof}[References for the proof]
The boundedness is \eqref{eq:map1} above with $s=1/2$.
The invertibility is proved in  \cite[Theorem 5.1]{Ve:84} for $d=3$ and \cite[Theorem 4.11]{Ve:84} for $d=2$.
Note that \cite{Ve:84} assumes for simplicity that $\Omega_-$ and $\Gamma$ are connected, but it is clear that this implies that the result holds when $\Gamma$ is the boundary of any bounded Lipschitz open set
(with this result for $d=3$ contained in \cite[Theorem 4.1]{MitreaD:97}). Indeed, in this case, $\Gamma$ and $\mathbb{R}^d\setminus \Gamma$ each have finitely-many connected components, and the results of \cite{Ve:84} for the case when $\Gamma$ is connected imply that $S$ is Fredholm of index zero as an operator $L^2(\Gamma)\mapsto H^{1}(\Gamma)$. Further, $S:L^2(\Gamma)\to H^1(\Gamma)$ is injective since $S$ is invertible as an operator from $H^{-1/2}(\Gamma)$ to $H^{1/2}(\Gamma)$ \cite[Corollary 8.13, Theorem 8.16]{Mc:00}.
\end{proof}

\section{Recap of harmonic-analysis results}
\label{app:HA}

In this appendix we recap results on the behaviour of solutions to Laplace's or Poisson's equation near the boundary of the domain. For simplicity, these results are stated for a bounded Lipschitz domain $D$ with boundary $\partial D$. Analogues of the results then hold with $D=\Omega_-$ and $D=\Oe$, where in the latter case spaces such as $H^1(D)$ become $H^1_{\rm loc}(\Oe)$ (since these results do not assume any particular behaviour at infinity).

\begin{theorem}\mythmname{\cite[\S\S5.1.2, 5.2.1]{Ne:67}, \cite[Theorem 4.24]{Mc:00}}
\label{thm:Necas}
If $u\in H^1(D)$ and $\Delta u \in L^2(D)$, then $\partial_n u \in L^2(\Gamma)$ iff $\gamma u \in H^1(\Gamma)$.
\end{theorem}

Given $\mathbf{x}\in \Gamma$, let $\Theta(\mathbf{x})$ be the non-tangential approach set to $\mathbf{x}$ from $D$ defined, for some sufficiently large $C>1$, as in \eqref{eq:ntapproach2}. Given $u\in C^2(D)$ with $\Delta u=0$, let the non-tangential maximal function of $u$,  $u^*$, be defined by \eqref{eq:ntmax}, and let the non-tangential limit of $u$, $\widetilde{\gamma}u$, be defined by \eqref{eq:ntlim}.

\begin{theorem}\mythmname{\cite[Corollaries 5.5 and 5.7]{JeKe:95}}\label{thm:JeKe:95}
Let $u \in C^2(D)$ with $\Delta u=0$.

(i) $u^*\in L^2(\partial D)$ iff $u\in H^{1/2}(D)$.

(ii) $u^*\in L^2(\partial D)$ implies that $\widetilde{\gamma} u \in L^2(\partial D)$.

(iii) $(\nabla u)^*\in (L^2(\partial D))^d$ iff $u\in H^{3/2}(D)$.

(iv) $(\nabla u)^*\in (L^2(\partial D))^d$ implies that $\widetilde{\gamma} (\nabla u) \in (L^2(\partial D))^d$.
\end{theorem}

\begin{lemma}\mythmname{\cite[Lemma A.9]{ChGrLaSp:12}}\label{lem:A9}
If $u\in H^s(D)$ with $s>1/2$ and $\Delta u=0$, then $\widetilde{\gamma}u = \gamma u $.
\end{lemma}

\begin{lemma}
\label{lem:A10}
Let $u \in C^2(D)\cap H^1(D)$ with $\Delta u=0$.
Then $u\in H^{3/2}(D)$ iff $\partial_n u\in L^2(\partial D)$ and $\gamma u \in H^1(\partial D)$. Furthermore, if $u\in H^{3/2}(D)$ then, almost everywhere on $\partial D$,
\begin{equation}\label{eq:gradutrace}
\widetilde{\gamma}(\nabla u) = \mathbf{n} \,\partial_n u + \nabla_{\partial D}(\gamma u).
\end{equation}
\end{lemma}

\begin{proof}
The forward implication and the trace result \eqref{eq:gradutrace} are proved in \cite[Lemma A.10]{ChGrLaSp:12}.
For the reverse implication, assume that $\partial_n u \in L^2(\partial D)$ and $\gamma u \in H^1(\partial D)$. Let $v:= {\cal S}\phi$ where the single-layer potential ${\cal S}$ is defined by \eqref{eq:SLP} (with $\Gamma$ replaced by $\partial D$) with $a\neq {\rm Cap}_{\partial D}$ when $d=2$. Let $\phi:= S^{-1}\gamma u$, so that $\phi\in L^2(\partial D)$ by Lemma \ref{lem:SLP}. By the first jump relation in \eqref{eq:jump1} $\gamma v= \gamma u$, so that $v=u$ by uniqueness of the interior Dirichlet problem of Definition \ref{def:idp}.
Since $\phi\in L^2(\partial D)$, $(\nabla u)^*\in (L^2(\partial D))^d$ by \cite[Theorem 1.6]{Ve:84} (see also \cite[Chapter 15, Theorem 5]{MeCo:00}), and thus
$u\in H^{3/2}(D)$ by Part (iii) of Theorem \ref{thm:JeKe:95}.
\end{proof}

\begin{corollary}\label{cor:CoDa}
The space $V(D)$ defined by \eqref{eq:V} is equal to $\{ v : v \in H^{3/2}(D), \, \Delta v \in L^2(D)\}$.
\end{corollary}

\begin{proof}
With $\Phi$ defined by \eqref{eq:fund} and $f\in L^2(D)$, let
\begin{equation*}
{\cal N} f(\mathbf{x}) := \int_D \Phi(\mathbf{x},\mathbf{y}) \, f(\mathbf{y})\, {\rm d} \mathbf{y}, \quad\text{ for } \mathbf{x} \in D;
\end{equation*}
i.e., ${\cal N}$ is a Newtonian potential. Recall that ${\cal N}: L^2(D)\rightarrow H^2(D)$ by, e.g., \cite[Theorem 6.1]{Mc:00}, and $\Delta( {\cal N} f)=-f$.

Given $v \in V(D)$, observe that ${\cal N}(\Delta v) \in H^2(D)$, $\partial_n {\cal N}(\Delta v) = \mathbf{n} \cdot \gamma(\nabla {\cal N}(\Delta v))\in L^2(\partial D)$, and $\gamma {\cal N}(\Delta v) \in H^1(\partial D)$; to see this last point observe that Theorem \ref{thm:Necas} implies that
\begin{equation*}
\gamma: \big\{ u \, :\, u\in H^1(D),\Delta u \in L^2(D), \partial_n u \in L^2(\partial D) \big\} \rightarrow H^1(\partial D),
\end{equation*}
and thus, in particular, $\gamma : H^2(D) \rightarrow H^1(\partial D)$. Therefore ${\cal N}(\Delta v)\in V(D)$.
Now let $\widetilde{v}: = v + {\cal N}(\Delta v)$. Then $\widetilde{v}\in V(D)$ with $\Delta \widetilde{v}=0$ and hence $\widetilde{v}\in C^2(D)$ by elliptic regularity. Therefore $\widetilde{v}\in H^{3/2}(D)$ by Lemma \ref{lem:A10}. The result that $v\in H^{3/2}(D)$ then follows
since $v= \widetilde{v} - {\cal N}(\Delta v)$ and ${\cal N}(\Delta v)\in H^2(D)$.

The reverse inclusion is proved similarly:~given $v\in H^{3/2}(D)$ with $\Delta v\in L^2(D)$, define $\widetilde{v}$ as before.
Since ${\cal N}(\Delta v)\in H^2(D)$, $\widetilde{v}\in H^{3/2}(D)$. Since  $\Delta \widetilde{v}=0$, $\widetilde{v}\in C^2(D)$ by elliptic regularity, and then  $\partial_n \widetilde{v}\in L^2(\partial D)$ and $\gamma \widetilde{v}\in H^1(\Gamma)$ by Lemma \ref{lem:A10}; thus $\widetilde{v} \in V(D)$. The result that $v\in V(D)$ then follows from the definition of $\widetilde{v}$ and the fact that ${\cal N}(\Delta v)\in H^2(D)\subset V(D)$
\end{proof}

\

We also need the following results in $L^p(\partial D)$ for $p\neq 2$ (as opposed to the $L^2$-based results above).

\begin{theorem}\mythmname{\cite[Theorem 5.6.1]{Me:18}}\label{lem:Medkova}
Suppose $u \in C^2(D)$ with $\Delta u=0$ and $(\nabla u)^*\in (L^p(\partial D))^d$ for some $1<p<\infty$. Then $u^* \in L^p(\partial D)$, $\widetilde{\gamma} u \in L^p(\partial D)$, and
$\widetilde{\gamma} (\nabla u) \in (L^p(\partial D))^d$.
\end{theorem}

\begin{corollary}\label{cor:Lp}
Suppose $u \in C^2(D)$ with $\Delta u=0$ and $(\nabla u)^*\in (L^p(\partial D))^d$ for some $2<p<\infty$. Then $\gamma u \in L^p(\partial D)$, $\partial_n u\in L^p(\partial D)$, and $\nabla_{\partial D}(\gamma u) \in (L^p(\partial D))^d$.
\end{corollary}

\begin{proof}
Since $(\nabla u)^*\in (L^2(\partial D))^d$, $u\in H^{3/2}(D)$ by Part (iii) of Theorem \ref{thm:JeKe:95}. Further, $\widetilde{\gamma} u \in L^p(\partial D)$ and
$\widetilde{\gamma} (\nabla u) \in (L^p(\partial D))^d$ by Theorem \ref{lem:Medkova}. Also, $\gamma u =\widetilde{\gamma} u$ by Lemma \ref{lem:A9}, and thus $\gamma u \in L^p(\partial D)$. By Lemma \ref{lem:A10}, \eqref{eq:gradutrace} holds, and thus $\partial_n u = \mathbf{n} \cdot \widetilde{\gamma} (\nabla u)\in L^p(\partial D)$ and
$\nabla_{\partial D}(\gamma u) \in (L^p(\partial D))^d$.
\end{proof}

\section{Proofs of Theorems \ref{thm:BVPequiv} and \ref{thm:BVPequivH}}\label{app:BVPequiv}

\begin{proof}[Proof of Theorem \ref{thm:BVPequiv}]
We prove the result for the IDP when $d=3$; the proof for the EDP is very similar.
The proof for the IDP when $d=2$ is also similar, with use of \cite[Theorem 5.15.2]{Me:18} replaced by use of \cite[Theorem 5.15.3]{Me:18}.

If $u$ is the solution of the IDP in the sense of Definition \ref{def:idp} then, since $u\in H^1(\Omega_-)$, $\widetilde{\gamma}^- u = \gamma^- u$ by Lemma \ref{lem:A9} and thus $\widetilde{\gamma}^- u =g_D$. Furthermore, since $u\in H^{1/2}(\Omega_-)$, $u^* \in{L^2(\Gamma)}$ by Part (i) of Theorem \ref{thm:JeKe:95}. Finally, by elliptic regularity $u\in C^2(D)$. Therefore $u$ is a solution of the IDP in the sense of Definition \ref{def:idp2}

To prove the converse,
let  $v := (-{\cal D}+{\cal S})\phi$ for $\phi\in {L^2(\Gamma)}$, with ${\cal D}$ the double-layer potential defined by \eqref{eq:DLP} and ${\cal S}$ the single-layer potential defined by \eqref{eq:SLP}.
Now $\widetilde{\gamma}^+ v = (\frac{1}{2} I - D + S)\phi$, where we have used that (i) $\widetilde{\gamma}^- {\cal D}\phi = (-\frac{1}{2} I + D)\phi$
by \cite[Theorem 1.10]{Ve:84} (similarly to \eqref{eq:jumpCald}) and (ii) $\widetilde{\gamma}^+ {\cal S}\phi = \gamma^+ {\cal S}\phi=S\phi$ by Lemma \ref{lem:A9} and the first jump relation in \eqref{eq:jump1}.
Since $(\frac{1}{2} I - D+S):{L^2(\Gamma)}\rightarrow {L^2(\Gamma)}$ is invertible by  \cite[Theorem 5.15.2]{Me:18} \footnote{Note that the operator $K$ in \cite{Me:18} equals minus our $D$ (see \cite[\S5.3]{Me:18}).},
if $\phi := (\frac{1}{2} I - D+S)^{-1}g_D$, then $v$ is a solution to the IDP of Definition \ref{def:idp2}.
The solution of this BVP is unique by \cite[Page 41]{Ca:85} and \cite[Lemma 3.7]{Ve:84} \footnote{ \cite[Lemma 3.7]{Ve:84} justifies how the uniqueness argument of \cite[Theorem 2.3]{FaJoRi:78} for $C^1$ domains also holds for Lipschitz domains.},
 and thus $v=u$.
 Arguing as in the proof of \cite[Theorem 2.27]{ChGrLaSp:12}, one can show that $(-\frac{1}{2} I + D+S):{H^{1/2}(\Gamma)}\rightarrow {H^{1/2}(\Gamma)}$ is invertible -- this follows by proving that $(-\frac{1}{2} I + D'+S):{H^{-1/2}(\Gamma)}\rightarrow {H^{-1/2}(\Gamma)}$ is invertible, which in turn follows since $-\frac{1}{2} I + D':{H^{-1/2}(\Gamma)}\rightarrow {H^{-1/2}(\Gamma)}$ is Fredholm of index zero by \cite[Theorem 2.25]{ChGrLaSp:12}, $S:{H^{-1/2}(\Gamma)}\rightarrow {H^{-1/2}(\Gamma)}$ is compact by \eqref{eq:map1},
 and $-\frac{1}{2} I + D'+S$ is injective by uniqueness of the Laplace exterior Robin problem. Therefore, if $g_D\in {H^{1/2}(\Gamma)}$, then $\phi \in {H^{1/2}(\Gamma)}$, and $u\in H^1(\Omega_-)$ by the mapping properties of ${\cal D}$ and ${\cal S}$ in \eqref{eq:LPmap}. Finally, $\gamma^- u = \widetilde{\gamma}^-u$ by Lemma \ref{lem:A9}, and thus $\gamma u = g_D$ and $u$ is the solution of the IDP in the sense of Definition \ref{def:idp}.
\end{proof}

\

\begin{proof}[Proof of Theorem \ref{thm:BVPequivH}]
The proof is similar to the proof of Theorem \ref{thm:BVPequiv}. The main difference is that now we define $v:= ({\cal D}_k - {\rm i} k {\cal S}_k)\phi$, where the Helmholtz single- and double-layer potentials ${\cal S}_k$ and ${\cal D}_k$ are defined by \eqref{eq:SLP} and \eqref{eq:DLP} with $\Phi$ replaced by $\Phi_k$. Now $\widetilde{\gamma}^+ v = (\frac{1}{2} I + D_k - {\rm i} k S_k)\phi$, where we have used that (i) $\widetilde{\gamma}^+ {\cal D}\phi = (\frac{1}{2} I + D_k)\phi$ by \cite[\S4]{ToWe:93}, \cite[Page 111]{ChGrLaSp:12}, and (ii) $\widetilde{\gamma}^+ {\cal S}\phi = \gamma^+ {\cal S}\phi =S_k \phi$ by Lemma \ref{lem:A9} and the Helmholtz analogue of the first jump relation in \eqref{eq:jump1}.
The proof then follows the same steps as in the proof of Theorem \ref{thm:BVPequiv}, using uniqueness of the Helmholtz EDP in the sense of Definition \ref{def:Hedp2} (see \cite[Theorem 5.6, Part (ii)]{ToWe:93}) and the fact that $A_{k,k}:= \frac{1}{2} I + D_k - {\rm i} k S_k :H^{s+1/2}(\Gamma)\rightarrow H^{s+1/2}(\Gamma)$ is bounded and invertible for $|s|\leq 1/2$ by \cite[Theorem 2.27]{ChGrLaSp:12}.
\end{proof}

\end{appendix}

\section*{Acknowledgements}

EAS was supported by the UK Engineering and Physical Sciences Research Council grants EP/F06795X/1, EP/1025995/1, and EP/R005591/1, and SCW by grant EP/F067798/1.

\footnotesize{
\bibliographystyle{plain}