---
title: Potential Fredholm Neural Network (PFNN)
url: https://www.emergentmind.com/topics/potential-fredholm-neural-network-pfnn
type: topic
---

# Potential Fredholm Neural Network (PFNN)

to=arxiv_search.search  天天中彩票app  彩神争霸怎么json
{"query":"\"Potential Fredholm Neural Network\" OR \"Fredholm Neural Networks\" elliptic PDEs boundary integral", "max_results": 10, "sort_by": "relevance"}
to=arxiv_search.search  大发快三大小单双  լինգjson
{"query":"PFNN penalty-free neural network boundary value problems complex geometries", "max_results": 10, "sort_by": "relevance"}
to=arxiv_search.search  ส่งเงินบาทไทยjson
{"query":"Fredholm integral equations training neural networks tensor train Ritz-Galerkin", "max_results": 10, "sort_by": "relevance"}
Potential Fredholm Neural Network (PFNN) denotes a class of deep neural architectures for elliptic partial differential equations in which the network is derived directly from a boundary integral formulation and a fixed-point scheme from potential theory. In the formulation developed for two-dimensional linear and semi-linear elliptic problems, a PFNN combines a Fredholm Neural Network (FNN) that solves the boundary integral equation for an unknown boundary density with a further layer that evaluates the corresponding double-layer potential inside the domain. Its layers, widths, weights, biases, and hyperparameters are determined by the fundamental solution, quadrature rules, and the fixed-point iteration, rather than by generic end-to-end training of internal parameters [2507.06038].

## 1. Terminology and scope

Within recent arXiv literature, the acronym **PFNN** is not unique. It refers both to **Potential Fredholm Neural Network** in the boundary-integral literature and to **Penalty-Free Neural Network** in variational neural PDE solvers on complex geometries. The distinction is substantive rather than stylistic: the former is rooted in potential theory and Fredholm integral equations, while the latter is rooted in weak formulations and exact architectural enforcement of essential boundary conditions [2507.06038] [2004.06490].

| Acronym usage | Meaning | Representative source |
|---|---|---|
| PFNN | Potential Fredholm Neural Network | [2507.06038] |
| PFNN | Penalty-Free Neural Network | [2004.06490] |
| PFNN-2 | Domain decomposed penalty-free variant | [2205.00593] |

In the Potential Fredholm sense, the term “potential” refers to the boundary-integral method and layer-potential representations of elliptic PDEs. This should not be conflated with the “potential” or energy interpretation used in variational neural solvers, nor with function-space Fredholm formulations in which a parameter function acts as a “potential” over parameter space [2303.05262]. A recurring misconception is therefore to treat all PFNN papers as belonging to one method family; the arXiv record instead shows two distinct lines of development that share an acronym but differ in mathematical formulation and computational workflow.

## 2. Boundary-integral and potential-theoretic formulation

The PFNN developed for elliptic PDEs addresses two-dimensional Poisson, modified Helmholtz, and semi-linear elliptic equations on a bounded domain $\Omega \subset \mathbb{R}^2$, with detailed treatment of Dirichlet boundary conditions. For Poisson,
$$
\begin{cases}
\Delta u(x) = \psi(x), & x \in \Omega,\\
u(x) = f(x), & x \in \partial\Omega,
\end{cases}
$$
the solution is represented by a double-layer potential plus a volume potential,
$$
u(x)=\int_{\partial\Omega}\beta(y)\,\frac{\partial\Phi(x,y)}{\partial n_y}\,d\sigma_y
+\int_\Omega \Phi(x,y)\,\psi(y)\,dy,
$$
where $\Phi$ is the fundamental solution and $\beta$ is an unknown boundary density. Imposing the boundary condition through the jump relation produces a Fredholm integral equation of the second kind for $\beta$ on $\partial\Omega$ [2507.06038].

For Poisson, the fundamental solution is
$$
\Phi(x,y)=\frac{1}{2\pi}\ln|x-y|,
$$
while for $\mathcal L=\Delta-\lambda$ the modified Helmholtz fundamental solution is
$$
\Phi(x,y)=-\frac{1}{2\pi}K_0(\sqrt{\lambda}\,|x-y|),
$$
with $K_0$ the modified Bessel function of the second kind. In the Poisson case, the boundary density satisfies
$$
\beta(x^*) = 2\Big(f(x^*)-\int_\Omega \Phi(x^*,y)\psi(y)\,dy\Big)
-2\int_{\partial\Omega}\beta(y)\,\frac{\partial\Phi(x^*,y)}{\partial n_y}\,d\sigma_y,
\qquad x^*\in\partial\Omega.
$$
This Fredholm equation is the object solved by the FNN block inside the PFNN [2507.06038].

A technical difficulty is the jump behavior of the double-layer potential near the boundary. The PFNN construction therefore uses a boundary-regular representation in which a boundary point $x^*\in\partial\Omega$ with the same angular coordinate as $x\in\Omega$ is introduced, and the potential is rewritten in terms of $\beta(y)-\beta(x^*)$. For Poisson,
$$
u(x)=\int_{\partial\Omega}(\beta(y)-\beta(x^*))\Big(\frac{\partial\Phi}{\partial n_y}(x,y)-\frac{1}{4\pi}\Big)\,d\sigma_y
+\frac{1}{2}\beta(x^*)+\int_{\partial\Omega}\frac{1}{4\pi}\beta(y)\,d\sigma_y
+\int_\Omega \Phi(x,y)\psi(y)\,dy.
$$
For Helmholtz, an analogous identity uses $\delta\Phi(x,y)=\Phi(x,y)-\Phi(x^*,y)$ and modified coefficients involving $\lambda$ [2507.06038]. This boundary-regular reformulation is central: it removes the singular jump term from the integral part and transfers it into explicit additive terms that are directly encoded in the network.

## 3. Network architecture and fixed-point interpretation

The PFNN consists of two coupled components. The first is a Fredholm NN that solves the boundary integral equation for $\beta$ by unrolling a fixed-point iteration; the second is a potential block that evaluates the PDE solution $u(x)$ from the approximate boundary density. In the linear boundary-integral setting, the fixed-point step is written abstractly as
$$
f(x)=(\mathcal T f)(x)=g(x)+\int_{\mathcal A}K(x,z)f(z)\,dz,
$$
and approximated by Krasnosel’skii–Mann iteration,
$$
\beta^{(k+1)}=(1-\kappa)\beta^{(k)}+\kappa\,\mathcal T\beta^{(k)},\qquad 0<\kappa\le 1.
$$
The PFNN depth therefore has a direct numerical interpretation: the number of hidden layers equals the number of fixed-point iterations, and the layer width equals the number of quadrature nodes on the boundary [2507.06038].

On a boundary grid $\{z_1,\dots,z_N\}$, the first hidden layer encodes the source term through
$$
W_1=(\kappa g(z_1),\dots,\kappa g(z_N))^\top,\qquad b_1=0.
$$
For inner layers $m=2,\dots,M-1$, the weight matrix has entries
$$
W_m(i,j)=
\begin{cases}
K(z_i,z_i)\Delta z+(1-\kappa), & i=j,\\
K(z_i,z_j)\Delta z, & i\neq j,
\end{cases}
$$
and the bias is
$$
b_m=(\kappa g(z_1),\dots,\kappa g(z_N))^\top.
$$
The activation in this block is linear. This is a literal network encoding of the discretized integral operator and its relaxation parameter, not a metaphorical analogy to iteration [2507.06038].

The additional “potential” part of the PFNN is what distinguishes it from a generic Fredholm NN. After the boundary-density approximation $\beta_M$ has been obtained, an extra hidden layer computes the vector $\beta(y_j)-\beta(x^*)$ via
$$
W_{M+1}=I_{N\times N},\qquad b_{M+1}=(-\beta(x^*),\dots,-\beta(x^*))^\top.
$$
The output layer then evaluates the regularized double-layer potential. In the Poisson case, the output weights are
$$
W_O=(\mathcal D\Phi(x,y_1)\Delta\sigma_y,\dots,\mathcal D\Phi(x,y_N)\Delta\sigma_y)^\top,
$$
where
$$
\mathcal D\Phi(x,y)=\frac{\partial\Phi}{\partial n_y}(x,y)-\frac{\partial\Phi}{\partial n_y}(x^*,y),
$$
and the output bias contains the explicit jump and volume terms,
$$
b_O=\frac12\beta(x^*)+\int_{\partial\Omega}\beta(y)\frac{\partial\Phi(x^*,y)}{\partial n_y}\,d\sigma_y
+\int_\Omega \Phi(x,y)\psi(y)\,dy.
$$
The modified Helmholtz case uses the corresponding formula with $\delta\Phi$ and $\lambda$-dependent coefficients [2507.06038].

This architecture makes PFNN explainable in a strong numerical-analysis sense. Each layer corresponds to a known operation: application of the boundary kernel, fixed-point relaxation, subtraction of the boundary trace, and potential evaluation. Hyperparameters are likewise numerical parameters: depth is iteration count, width is quadrature resolution, and the weights are discretized kernel evaluations. The closely related “Fredholm Neural Networks” framework formalized this iteration-as-network construction for linear and nonlinear Fredholm equations of the second kind and later served as the basis of the PFNN extension to elliptic PDEs [2408.09484].

## 4. Theoretical properties and error decomposition

The PFNN theory is built on two error layers: the error in approximating the boundary density by the FNN block, and the discretization error in evaluating the boundary and volume potentials. For a nonexpansive Fredholm operator, the underlying FNN approximation $f_M$ satisfies an explicit estimate of the form
$$
\|f-f_M\|_{\mathcal A}
\le
\frac{e^{1-q}}{1-q}
\left(\|\mathcal T g-g\|_{\mathcal A}+\frac{D(b-a)^2}{2n}\right)
e^{-(1-q)v_M},
$$
so the boundary-density error is controlled by the iteration count $M$, the quadrature resolution $n$, and the operator’s nonexpansive constant. In the constant-$\kappa$ case, $v_M=M\kappa$, making the dependence of accuracy on depth explicit [2507.06038].

The PFNN error bound for the PDE solution $u$ is then expressed in terms of the boundary-density approximation error $\|\beta-\beta_M\|_{\partial\Omega}$ and four discretization terms $D_1,D_2,D_3,D_4$ associated with the potential and boundary integrals. For interior points,
$$
\|u-\tilde u\|_\Omega
\le
\|\beta-\beta_M\|_{\partial\Omega}(\cdots)
+\|\beta_M\|_{\partial\Omega}D_1(\Omega)
+D_2(\Omega)+D_3(\Omega)+D_4(\Omega),
$$
while on the boundary,
$$
\|f-\tilde u\|_{\partial\Omega}
\le
\|\beta-\beta_M\|_{\partial\Omega}(\cdots)
+D_2(\partial\Omega)+D_4(\partial\Omega).
$$
The omitted coefficients are explicit integrals of $\mathcal D\Phi$, $\delta\Phi$, and $\partial\Phi/\partial n_y$. The structural point is that the total PFNN error is decomposed into iteration error from the boundary Fredholm solve and quadrature error from the potential evaluation; both contributions correspond directly to identifiable architectural components [2507.06038].

For semi-linear elliptic equations, the PFNN is embedded in an outer fixed-point loop, yielding a **Recurrent PFNN (RPFNN)**. If $T$ denotes the outer contraction map, then after $n$ outer iterations
$$
\|\tilde u_n-u\|_\Omega \le \epsilon_\Omega^{RPFNN}+q^n\|u_0-u\|_\Omega,
\qquad
\|\tilde u_n-f\|_{\partial\Omega} \le \epsilon_{\partial\Omega}^{RPFNN}+q^n\|u_0-u\|_{\partial\Omega},
$$
where the $\epsilon^{RPFNN}$ terms are precisely the inner PFNN discretization errors. This gives the semi-linear architecture a two-scale convergence structure: geometric decay from the outer contraction and explicit discretization control from the inner boundary-integral solver [2507.06038].

A distinctive claim of this framework is boundary fidelity. Because the PFNN is derived from a representation that already encodes the jump behavior of the double-layer potential, boundary conditions are respected by construction rather than through penalty terms. The earlier Fredholm NN paper showed a related Laplace-unit-disc example with maximum error about $2.27\times 10^{-7}$ for the potential-theoretic architecture, again emphasizing that the boundary treatment is inherited from the integral representation rather than imposed by an auxiliary loss term [2408.09484].

## 5. Forward and inverse formulations

For forward problems, PFNN is a deterministic solver. Given boundary data $f$ and source information $\psi$ or $F(x,u)$, the architecture is assembled from the PDE data, the boundary grid, quadrature rules, and the chosen number of iterations. Internal PFNN weights are not trained by generic backpropagation; they are precomputed from the Green’s function, quadrature, and fixed-point coefficients. This differentiates PFNN from PINN-style residual minimization, even though both are neural-network-based PDE solvers [2507.06038].

The reported forward experiments concern the unit disk. For Poisson,
$$
\Delta u = 2x_1,\qquad u|_{\partial\Omega}=0,
$$
with exact solution
$$
u(x)=\frac14 x_1(x_1^2+x_2^2-1),
$$
a boundary grid $N=1000$ and radial grid size $100$ were used, and both $L_\infty$ and MAE errors in the interior and on the boundary decreased monotonically with the number of layers $M$. For the modified Helmholtz example,
$$
\Delta u-u=-x_1^3+2x_2^2+6x_1-4,\qquad
u|_{\partial\Omega}=x_1^3+2x_1^2-2,
$$
with exact solution $u(x)=x_1^3-2x_2^2$, the PFNN again produced accurate solutions with very small boundary error. For the semi-linear equation
$$
\Delta u-e^u=-e^{1-x_1^2-x_2^2}-4,\qquad u|_{\partial\Omega}=0,
$$
with exact solution $u(x)=1-(x_1^2+x_2^2)$, an RPFNN with $12$ outer iterations, $100$ PFNN layers per iteration, a $120\times120$ polar grid, and $\lambda=1$ was used; the boundary error was again reported as near machine precision [2507.06038].

Inverse problems are handled differently. In the source-identification example for Poisson,
$$
\Delta u(x)=\psi(x),\qquad u|_{\partial\Omega}=f,
$$
the PFNN is kept fixed as the PDE solver, while a separate shallow neural model $\psi_\theta$ for the unknown source is trained from coarse interior measurements. The loss is
$$
L(\theta)=\frac1N\sum_{i=1}^N\big(\tilde u(x_i)-\hat u(x_i;\psi_\theta)\big)^2
+\lambda_{reg}\mathcal R(\theta),
$$
with Tikhonov regularization $\mathcal R(\theta)=\|\psi_\theta\|_2^2$. In the reported implementation, the source model had one hidden layer with $20$ nodes and $\tanh$ activation, the optimizer was Levenberg–Marquardt with $600$ iterations, and $\lambda_{reg}=10^{-12}$. On a test grid of size $50\times 50$, the mean training MSE was $8.45\times 10^{-7}$, the mean test MSE was $1.72\times 10^{-6}$, the mean test $L_\infty$ error was $7.73\times 10^{-3}$, and the boundary test MAE and $L_\infty$ were $1.95\times 10^{-15}$ and $5.26\times 10^{-15}$, respectively [2507.06038].

## 6. Relation to adjacent Fredholm and neural PDE frameworks

PFNN did not emerge in isolation. The immediate precursor is the Fredholm NN framework for linear and nonlinear Fredholm integral equations of the second kind, where a fully connected DNN with analytically specified weights reproduces fixed-point iterations and can be extended by a final integral layer to evaluate potential representations for elliptic PDEs [2408.09484]. A later generalization, **Fredholm Integral Neural Operators (FREDINOs)**, moved from fixed known kernels to learned contractive integral operators in arbitrary dimensions and explicitly incorporated PFNNs for a nonlinear elliptic PDE through a boundary integral equation formulation [2604.03034].

Two other nearby strands are conceptually related but mathematically distinct. First, the paper on “Fredholm integral equations for function approximation and the training of neural networks” recasts shallow-network training as the approximate least-squares solution of a Fredholm integral equation of the first kind, with a parameter function over parameter space learned by Ritz–Galerkin, Tikhonov regularization, and tensor-train methods [2303.05262]. Second, the LDNN work on nonlinear Volterra–Fredholm–Hammerstein integral equations uses a Legendre-based network and Gaussian quadrature collocation to minimize integral-equation residuals; it is not a PFNN in the potential-theoretic sense, but it shows another way of hard-wiring Fredholm structure into a neural architecture [2106.14320].

The most frequent confusion concerns the acronym collision with the **Penalty-Free Neural Network** literature. That PFNN uses two neural networks, a length factor function, and weak/energy formulations to enforce Dirichlet conditions without penalty terms for second-order boundary-value problems on complex geometries; PFNN-2 extends that line to non-self-adjoint and time-dependent PDEs with overlapping domain decomposition [2004.06490] [2205.00593]. The shared acronym should therefore not be taken as evidence of a shared architecture. Potential Fredholm PFNNs are boundary-integral fixed-point networks; Penalty-Free PFNNs are weak-form trial-space constructions.

## 7. Limitations and open directions

The current Potential Fredholm formulation is developed for two-dimensional elliptic PDEs with smooth boundaries and Dirichlet conditions, with examples on the unit disk. Extension to three dimensions is described as conceptually possible but computationally more expensive, and Neumann or Robin conditions would require modified integral formulations and jump relations. The semi-linear treatment is likewise restricted to monotone-iteration settings in which a contraction can be ensured [2507.06038].

The framework’s numerical strengths depend on boundary discretization and quadrature accuracy. For very complex geometries, both the construction of the boundary grid and the evaluation of the singular or weakly singular kernels become more delicate. This suggests a practical boundary between current PFNN implementations and more general geometric PDE solvers: PFNNs are strongest where boundary integral methods are already numerically advantageous, and less mature where geometry or coefficient structure makes boundary-element-style preprocessing expensive [2507.06038].

A plausible implication is that the later FREDINO program points toward a hybrid future for PFNNs. FREDINOs learn contractive linear and nonlinear integral operators while preserving fixed-point convergence, and in the PDE setting they learn a corrected potential kernel $\Phi_\theta$ within a PFNN-style boundary-integral architecture [2604.03034]. This suggests a broader interpretation of PFNN: not only as a deterministic encoding of a known boundary integral method, but also as an explainable scaffold into which learned kernels, operator corrections, or uncertainty-aware components can be inserted without abandoning the Fredholm and potential-theoretic structure.

Source: https://www.emergentmind.com/topics/potential-fredholm-neural-network-pfnn