Single-Layer Acoustic Potential

• The single-layer acoustic potential is defined via the Helmholtz operator’s fundamental solution and serves as a boundary integral tool for solving Neumann problems.
• It employs Hölder regularity and negative exponent spaces to handle distributional normal derivatives, allowing solutions with infinite energy.
• Its mapping and jump relations provide a robust framework for nonvariational analysis, aiding applications in scattering and inverse acoustic problems.

The single layer acoustic potential is a boundary integral construction fundamental to the nonvariational theory of the Helmholtz operator $(\Delta + k^2)$ in a bounded open set $\Omega \subset \mathbb{R}^n$ with boundary of class $C^{1,\alpha}$ for some $\alpha \in (0,1)$. Unlike traditional variational approaches based on energy integrals, the single layer potential is formulated in spaces of H\"older regularity and is central to the analysis of the Neumann problem when boundary data are distributional normal derivatives, i.e., elements of the space $V^{-1,\alpha}(\partial\Omega)$. This framework accommodates solutions whose Dirichlet energy may be infinite, thus lying outside classical Sobolev settings (Cristoforis, 15 Apr 2025).

1. Fundamental Solution and Definition

The building block of the single layer acoustic potential is the fundamental solution $\Phi_k$ to the Helmholtz operator in $\mathbb{R}^n$ for $n \ge 2$. $\Phi_k$ satisfies

$$(\Delta + k^2)\,\Phi_k = \delta_0 \qquad \text{in }\mathcal{D}'(\mathbb{R}^n),$$

where $k \in \mathbb{C}$ and $\delta_0$ denotes the Dirac measure at the origin. The explicit forms are: $\Phi_k(x) = \begin{cases} \displaystyle \frac{e^{ik|x|}}{4\pi|x|}, & n=3, \[2ex] \displaystyle \frac{i}{4} H_0^{(1)}(k|x|), & n=2, \end{cases}$ or, in general,

$$\Phi_k(x) = \frac{i}{4} \left(\frac{k}{2\pi|x|}\right)^{\frac{n-2}{2}} H_{\frac{n-2}{2}}^{(1)}(k|x|),$$

with $H_\nu^{(1)}$ the Hankel function of the first kind. The acoustic single layer potential associated with a density $\mu$ on $\partial\Omega$ is

$$S_k[\mu](x) = \int_{\partial\Omega} \Phi_k(x-y)\, \mu(y)\, d\sigma(y), \qquad x \in \mathbb{R}^n \setminus \partial\Omega.$$

The interior and exterior traces,

$$v^+(x) = \lim_{\Omega \ni z \to x} S_k[\mu](z), \qquad v^-(x) = \lim_{(\mathbb{R}^n \setminus \overline{\Omega}) \ni z \to x} S_k[\mu](z),$$

yield a common extension on smooth enough boundaries.

2. Function Spaces and Their Properties

The theoretical foundation relies on function spaces tailored for boundary regularity and distributional derivatives:

• Schauder–H\"older spaces: For $m \in \mathbb{N}_0$,

$$C^{m,\alpha}(\overline{\Omega}) = \{ u \in C^m(\overline{\Omega}) : [D^m u]_{\alpha,\overline{\Omega}} < \infty \},$$

with norm

$$\|u\|_{C^{m,\alpha}(\overline{\Omega})} = \sum_{|\beta| \le m} \sup_{x \in \Omega} |D^\beta u(x)| + \sum_{|\beta|=m} \sup_{x \neq y} \frac{|D^\beta u(x) - D^\beta u(y)|}{|x-y|^\alpha}.$$

• Negative exponent Hölder spaces: The space $C^{-1,\alpha}(\Omega)$ consists of distributions $f$ on $\Omega$ that can be written as

$$f = f_0 + \sum_{j=1}^n \partial_{x_j} f_j, \qquad f_0, f_1, \dots, f_n \in C^{0,\alpha}(\overline{\Omega}),$$

with norm defined by the infimum over all such decompositions. $C^{-1,\alpha}(\Omega)$ is a Banach space continuously embedded in $\mathcal{D}'(\Omega)$.

• Solution space:

$$C^{1,\alpha}(\Omega) = \{ u \in C^1(\overline{\Omega}) : \Delta u \in C^{0,\alpha}(\overline{\Omega}) \},$$

with norm $$\|u\|_{C^{1,\alpha}(\Omega)} = \|u\|_{C^1(\overline{\Omega})} + \|\Delta u\|_{C^{0,\alpha}(\overline{\Omega})}$$.

• Boundary-distribution space: $V^{-1,\alpha}(\partial\Omega)$ is defined as

$$V^{-1,\alpha}(\partial\Omega) = \{ \mu_0 + S_{2,+}^T[\mu_1] : \mu_0, \mu_1 \in C^{0,\alpha}(\partial\Omega) \},$$

where $S_{2,+}$ is the Steklov–Poincaré map and $S_{2,+}^T$ its transpose in the duality $(V^{-1,\alpha},\,C^{1,\alpha})$. The norm on $V^{-1,\alpha}(\partial\Omega)$ is induced by

$$\|\mu\|_{V^{-1,\alpha}(\partial\Omega)} = \inf_{\mu = \mu_0 + S_{2,+}^T[\mu_1]} \left( \|\mu_0\|_{C^{0,\alpha}} + \|\mu_1\|_{C^{0,\alpha}} \right).$$

This is a Banach space with $$C^{0,\alpha}(\partial\Omega) \hookrightarrow V^{-1,\alpha}(\partial\Omega) \hookrightarrow (C^{1,\alpha}(\partial\Omega))'$$.

3. Regularity and Continuity of the Single Layer Potential

The mapping properties of the single layer operator are formalized as follows:

For every $\mu \in V^{-1,\alpha}(\partial\Omega)$,

$v^+(x) = S_k[\mu](x), \ x \in \Omega$

satisfies

$$v^+ \in C^{0,\alpha}(\overline{\Omega}), \qquad (\Delta + k^2) v^+ = 0 \quad \text{in} \ \Omega,$$

thus $v^+ \in C^{1,\alpha}(\Omega)$, with

$$\|v^+\|_{C^{1,\alpha}(\Omega)} \leq C \|\mu\|_{V^{-1,\alpha}(\partial\Omega)}.$$

The constant $C$ depends only on $\Omega, n, \alpha, k$. For the exterior, $v^-(x) = S_k[\mu](x)$ for $x$ in the exterior domain belongs to $$C^{0,\alpha}_{\text{loc}}(\mathbb{R}^n \setminus \Omega)$$. These continuity and mapping results rely on decomposing $\mu$ into $C^{0,\alpha}(\partial \Omega)$ and $S_{2,+}^T$ terms, with the classical and distributional cases handled via Green identities and Schauder theory for each summand.

4. Jump Relations and Double Layer Operator

The jump relations determine the discontinuity in the normal derivative of $S_k[\mu]$ across the boundary $\partial \Omega$. The acoustic double-layer operator is defined as

$$(W_k\varphi)(x) = \int_{\partial\Omega} \partial_{\nu_y} \Phi_k(x-y) \, \varphi(y) \, d\sigma(y), \qquad x \notin \partial\Omega,$$

and its transpose $W_k^*$ acts on $V^{-1,\alpha}(\partial\Omega)$ in duality with $C^{1,\alpha}(\partial\Omega)$.

For every $\mu \in V^{-1,\alpha}(\partial\Omega)$,

$$\partial_\nu^+ S_k[\mu] = \left( -\frac{1}{2} I + W_k^* \right)[\mu], \qquad \partial_\nu^- S_k[\mu] = \left( \frac{1}{2} I + W_k^* \right)[\mu] \quad \text{on } \partial\Omega.$$

$W_k^*$ is a bounded operator on $V^{-1,\alpha}(\partial\Omega)$. Furthermore, for $0<\beta<\alpha$, the embedding $$V^{-1,\alpha}(\partial\Omega)\hookrightarrow V^{-1,\beta}(\partial\Omega)$$ is compact, and the mapping $W_k^* : V^{-1,\alpha} \to V^{-1,\alpha}$ is compact as a consequence.

Table: Main Operators and Spaces

Object	Definition/Formula	Properties
$\Phi_k$	Fundamental solution of $(\Delta + k^2)$	$C^2$ away from $0$, explicit formula above
$S_k[\mu]$	$$\int_{\partial\Omega} \Phi_k(x-y) \mu(y) d\sigma(y)$$	Maps $V^{-1,\alpha}$ to $C^{1,\alpha}(\Omega)$
$W_k$, $W_k^*$	Double-layer operator, its transpose in duality	$W_k^*$ bounded and compact on $V^{-1,\alpha}$
$V^{-1,\alpha}$	$$\{\mu_0 + S_{2,+}^T[\mu_1] : \mu_0,\mu_1 \in C^{0,\alpha}\}$$	Banach, intermediary for distributional boundary data

5. Nonvariational Solution of the Neumann Problem

The construction addresses the Neumann problem

$$(\Delta + k^2)u = 0 \quad \text{in } \Omega, \qquad \partial_\nu u = g \quad \text{on } \partial\Omega,$$

for $g \in V^{-1,\alpha}(\partial\Omega)$. The single-layer ansatz

$u(x) = S_k[\mu](x), \quad x\in\Omega$

leads to the boundary integral equation

$$\left( -\frac{1}{2} I + W_k^* \right)[\mu] = g \quad \text{on} \ \partial\Omega.$$

As $W_k^*$ is compact on $V^{-1,\alpha}$, $-\frac{1}{2}I + W_k^*$ is a Fredholm operator of index zero. The Fredholm alternative in the duality $(V^{-1,\alpha}, C^{1,\alpha})$ gives existence and uniqueness of $\mu \in V^{-1,\alpha}(\partial\Omega)$, and hence of $u = S_k[\mu]$.

A distinguishing feature is that the solution $u$ may not reside in $H^1(\Omega)$: its Dirichlet integral $\int_\Omega |\nabla u|^2$ can be infinite. Thus, the approach exists outside energy (variational) frameworks and traces historically to examples of H\"older-continuous harmonic functions with infinite energy by Prym and Hadamard.

6. Significance and Perspectives

The nonvariational theory of the single layer acoustic potential extends boundary integral methods for the Helmholtz equation to boundary data and solutions of lower regularity, beyond classical $H^1$ or $L^2$ theories. This enables rigorous treatment of Neumann problems where the data are distributional normal derivatives and the solutions may exhibit infinite energy. The flexible framework circumvents variational constraints and applies to $C^{1,\alpha}$ domains for arbitrary dimension $n \ge 2$, with all mapping and jump properties derived in compatible H\"older–Schauder scales (Cristoforis, 15 Apr 2025).

A plausible implication is the applicability of these results to scattering, inverse problems, and boundary element methods in nonstandard regularity settings, facilitating rigorous analysis where variational techniques are inadequate. The results foundationally support advances in mathematical acoustics and PDE theory where general boundary regularity and nonvariational data are encountered.