---
title: 'Local Potential Problems: A Survey'
url: https://www.emergentmind.com/topics/local-potential-problems
type: topic
---

# Local Potential Problems: A Survey

“Local potential problems” is not a single standardized term; on arXiv it denotes several formally distinct research programs that share a local variational or local equilibrium structure. In one line of work, local minimizers of nonlocal interaction energies induce classical or fractional obstacle problems for the interaction potential \(W * \mu\). In another, the term refers to locally checkable optimization or equilibrium conditions on graphs, local ground-state prediction tasks in disordered systems, or localized value functions in minimax optimization. It also appears in mixed local–nonlocal elliptic equations with Hardy-type singular potentials, in coarse-grained molecular modeling through local-density-dependent interaction laws, in nonequilibrium superconductivity through the local pair potential \(\Delta(\mathbf r,t)\), and in conceptual density functional theory through unsuccessful attempts to define a local hardness from \(\delta\mu/\delta n(\mathbf r)\) [1406.4040] [2505.02927] [2507.12038] [2407.06763] [1310.1288] [1107.4249]. This suggests that the expression is best understood as a family of locality-based formulations rather than a single canonical object.

## 1. Core meanings and recurring formal structure

Across these literatures, the phrase designates a problem in which a global object is reconstructed from local data, local constraints, or local response. Typical representatives include the interaction potential \(\psi=W*\mu\) near points of \(\operatorname{supp}\mu\), the local cost \(c_v\) and global potential \(\operatorname{Pot}(G,\ell)=\sum_{v\in V}\phi_\ell(G_r(v),v)\) in distributed graph problems, the local single-bond error \(\mathcal E_{ij}\) in spin glasses, and localized minimax value functions \(V_\tau(x)\) or \(V_{\epsilon_0}(x)\) under coupled constraints [2507.12038] [2505.02927] [2510.03861].

| Domain | Local object | Representative formulation |
|---|---|---|
| Interaction energies | \(\psi=W*\mu\) | Local obstacle problem |
| Distributed optimization | \(c_v,\operatorname{Pot}(G,\ell)\) | LOP / GLOP |
| Disordered systems | \(\mathcal E_{ij}\), \(J_{ij}^c\) | Local hardness |
| Mixed PDE | Hardy potential, Riesz/Wolff potentials | Mixed local–nonlocal elliptic problem |
| Mesoscopic modeling | \(\lambda(V)\), \(V_i\) | Local-density-dependent potential |
| DFT | \(\eta(\mathbf r)\), \(\tilde\mu(\mathbf r)\) | Local hardness / local chemical potential |

A second recurring feature is that “local” rarely means merely “small scale.” In the obstacle-problem setting, locality arises because \(d_\infty\)-minimality permits only small spatial perturbations. In distributed graph problems, locality is the constant-radius checkability of a labeling. In disordered systems, local solvers are tested against global critical thresholds and avalanches. In mixed elliptic problems, local singular potentials coexist with nonlocal diffusions, so local behavior is still influenced by tails. In conceptual DFT, the failure of a local hardness definition is traced precisely to the asymptotic behavior of the universal functional, showing that a purportedly local derivative can be controlled by behavior “at infinity” [1406.4040] [2507.12038] [2505.02927] [2407.06763] [1107.4249].

## 2. Interaction energies as local obstacle problems

In the interaction-energy setting, the basic functional is
\[
E[\mu]=\frac12\iint_{\mathbb R^N\times\mathbb R^N}W(x-y)\,d\mu(x)\,d\mu(y),
\]
for \(\mu\in\mathcal P(\mathbb R^N)\), with \(W\) nonnegative, lower semicontinuous, and locally integrable. Local minimality is metric, defined with the \(d_\infty\) optimal transport distance, and for an \(\varepsilon\)-local minimizer the associated interaction potential
\[
\psi(x)=(W*\mu)(x)
\]
satisfies the local Euler–Lagrange inequality
\[
\psi(x_0)\le \psi(x)\quad\text{for a.e. }x\in B_\varepsilon(x_0),
\]
for every \(x_0\in\operatorname{supp}\mu\). After continuity is established, this becomes pointwise and strengthens to local flatness on the support:
\[
\psi(x)\ge \psi(x_0)\ \forall x\in B_\varepsilon(x_0),\qquad
\psi(x)=\psi(x_0)\ \forall x\in \operatorname{supp}\mu\cap B_\varepsilon(x_0).
\]
The local potential problem is then to identify \(\psi\) as a solution of a unilateral PDE with coincidence set \(\{\psi=\psi(x_0)\}\) [1406.4040].

For Newtonian repulsion, the decomposition \(W=V+W_a\) with \(-\Delta V=\delta_0\) yields
\[
-\Delta\psi=\mu-(\Delta W_a)*\mu.
\]
Fixing \(x_0\in\operatorname{supp}\mu\) and writing \(F=\Delta W_a *\mu\), the potential solves in \(B_\varepsilon(x_0)\) the local obstacle problem
\[
\psi\ge C_0,\qquad -\Delta\psi\ge -F,\qquad -\Delta\psi=-F\ \text{in }\{\psi>C_0\},
\]
with obstacle \(C_0=\psi(x_0)\) and \(\mu=-\Delta\psi+F\). For Riesz-type repulsion \(W(x)\sim -|x|^{2s-N}\) with \(s\in(0,1)\), the same structure becomes fractional:
\[
(-\Delta)^s\psi=\mu-F,
\]
and \(\psi\) solves a local fractional obstacle problem in which the complement condition holds in \(B_\varepsilon(x_0)\cap\{\psi>C_0\}\) and the exterior datum is prescribed on \(\mathbb R^N\setminus B_\varepsilon(x_0)\) [1406.4040].

The regularity theory is dictated by the repulsion strength at the origin. Under Newtonian repulsion, compactly supported \(\varepsilon\)-minimizers satisfy \(\psi\in C^{1,1}(\mathbb R^N)\) and \(\mu=\rho(x)\,dx\) with \(\rho\in L^\infty(\mathbb R^N)\). If \(\Delta W_a\in W^{1,1}\), then \(\rho\in BV_{\mathrm{loc}}(\mathbb R^N)\), and if \(F=\Delta W_a*\mu>0\) near \(\partial\operatorname{supp}\mu\), then \(\operatorname{supp}\mu\) is a set of locally finite perimeter. For more singular-than-Newtonian repulsion, \(\psi\in C^{1,\gamma}\) for every \(\gamma<s\) and \(\rho\in C^\alpha(\mathbb R^N)\) for every \(\alpha<1-s\). The coincidence set of the obstacle problem identifies, up to negligible sets, with \(\operatorname{supp}\mu\), so free-boundary regularity transfers to the geometry of minimizers [1406.4040].

## 3. Landscape formulations: local hardness, overlap potentials, and calm local minimax

In disordered spin systems, “local potential problems” are studied through exact local optimization on subsystems. The paper on local optimization in EA and SK spin glasses introduces the local single-bond solver: one selects a target bond \(\langle i_0j_0\rangle\), cuts out a finite neighborhood, computes the exact subsystem ground state, and compares \(\sigma^{\rm sub}_{i_0}\sigma^{\rm sub}_{j_0}\) with the full-system ground-state bond product. The disorder-averaged error rate is
\[
\mathcal E_{ij}=\frac{1-[\sigma_i\sigma_j\sigma_i^{\rm sub}\sigma_j^{\rm sub}]}{2}.
\]
This error decays as a power law in subsystem size for the 2D and 3D EA models,
\[
\mathcal E_{ij}\sim \left(\frac{\ell_{\mathcal E}}{L_{\rm sub}}\right)^\kappa,
\]
with fitted parameters \(\ell_{\mathcal E}\approx0.20,\ \kappa\approx0.685\) in 2D and \(\ell_{\mathcal E}\approx 6\times 10^{-4},\ \kappa\approx0.18\) in 3D. The key organizing notion is local hardness: how large a local subsystem must be so that a local solver predicts a local observable with a prescribed small average error. That hardness is controlled not by \(P\) versus NP-hard global complexity, but by bond-specific critical thresholds \(J_{ij}^c\), gapless avalanche-like excitations, and zero-energy droplets. The local error obeys a stretched-exponential relation in distance to criticality,
\[
\mathcal E_{ij}\sim k_J e^{-\alpha |J_{ij}-J_{ij}^c|^\beta},
\]
and changes in critical thresholds decay algebraically with distance,
\[
[|\Delta J^c(r)|]\sim \left(\frac{\ell_c}{r}\right)^{a_c},
\]
with \(a_c\approx1.26\) in 2D and \(a_c\approx1.3\) in 3D [2505.02927].

A related landscape formulation appears in the overlap analysis of shortest paths and shortest path trees. There the Franz–Parisi potential is used as a constrained free-energy profile over an overlap parameter \(r\). For shortest paths, the overlap-gap property and the Franz–Parisi potential both indicate that local search fails in the path landscape, whereas for shortest path trees they indicate that local search should succeed. The tree Gibbs measure is
\[
\mu_{G_n,\beta}(T)=\frac{1}{Z_{G_n,\beta}}
\exp\!\Big(-\beta\log\log n\sum_{v\in\overline V_n} d_T(1,v)\Big),
\]
and the Franz–Parisi potential \(\mathcal F_\beta^{\mathrm{FP}}(r)\) is defined by constraining the overlap with a typical Gibbs sample. In the tree formulation, the overlap structure is continuous enough that there is no ensemble overlap-gap property, while in the path formulation local search encounters a barrier. Li and Schramm’s earlier lower bound against stable algorithms for shortest paths is thereby reinterpreted as a statement about one unfavorable landscape among two polynomial-time-equivalent formulations [2511.18666].

In minimax optimization with coupled constraints, the same local-potential language is used for a localized value function around a candidate equilibrium. The basic problem is
\[
\min_{x\in X}\max_{y\in Y(x)} f(x,y),
\]
with \(Y(x)\) depending on \(x\). A calm local minimax point \((\bar x,\bar y)\) is defined by the existence of a radius function \(\tau(\delta)\to0\) with \(\tau(\delta)\le \kappa\delta\) such that, for all small \(\delta\),
\[
f(\bar x,y)\le f(\bar x,\bar y)\le \max_{y'\in Y(x)\cap B_{\tau(\delta)}(\bar y)} f(x,y').
\]
This is equivalent to local optimality of a localized value function \(V_\tau(x)\) and inner calmness of the localized argmax map \(S_\tau(x)\). The resulting first- and second-order conditions are stated in terms of tangent cones, graphical derivatives \(DY(\bar x,\bar y)(u)\), and second subderivatives, and they specialize to KKT-type conditions under MSCQ, RS, MFCQ, RCRCQ, or affine/polyhedral assumptions [2510.03861].

## 4. Distributed graph-theoretic local potential problems

In distributed graph algorithms, a local potential problem is formalized as a locally checkable labeling endowed with a nonnegative local potential
\[
\phi : L(\Sigma_V^{\mathrm{in}}\times\Sigma_V^{\mathrm{out}},\,\Sigma_E^{\mathrm{in}}\times\Sigma_E^{\mathrm{out}},\,r,\Delta)\to \mathbb R_{\ge0},
\]
and a global potential
\[
\operatorname{Pot}(G,\ell)=\sum_{v\in V}\phi_\ell(G_r(v),v).
\]
A labeling is locally optimal when no node can change only its own output and incident half-edges in a way that strictly decreases its local cost and the potential in its radius-\(r\) neighborhood. This is the LOP framework; GLOPs extend it by allowing a node to relabel all nodes in its radius-\(r\) neighborhood, and every GLOP reduces in \(O(1)\) rounds to an LOP. The main algorithmic result is that on bounded-degree graphs every LOP admits a randomized LOCAL algorithm in \(O(\log^6 n)\) rounds w.h.p., and hence every local potential problem in this sense has deterministic and randomized complexity \(\log^{O(1)}n\). In particular, the deterministic complexity of locally optimal cut is settled to \(\log^{\Theta(1)}n\) [2507.12038].

A complementary classification is available on directed cycles for general local optimization problems of min-sum, max-sum, min-max, or max-min type. An opt LCL is specified by \((\Gamma,r,c,\mathrm{aggr},\mathrm{obj})\), where \(c:\Gamma^{r+1}\to \mathbb R_{\ge0}\cup\{\bot\}\) is a local cost or utility on radius-\(r\) neighborhoods and the global objective aggregates these local quantities by \(\sum\), \(\min\), or \(\max\). For every such problem and every constant approximation ratio \(\alpha\), the complexity on directed cycles falls into exactly one of four classes: \(O(1)\) deterministic and randomized; \(\Theta(\log^* n)\) deterministic and \(O(1)\) randomized; \(\Theta(\log^* n)\) in both models; or \(\Theta(n)\) in both models. Moreover, the complexity class and an asymptotically optimal algorithm can be found automatically by a centralized meta-algorithm [2602.13046].

The broader LCL landscape shows that some local problems do separate stronger models, but trees are unusually rigid. One constructed LCL requires \(\Omega(\sqrt n)\) rounds with private randomness yet drops to \(O(\log n)\) with shared randomness, resolving the question of whether shared randomness ever helps with LCLs [2407.05445]. Another LCL, iterated GHZ, is solvable in one quantum-LOCAL round but requires \(\Omega(\Delta)\) rounds classically, giving the first super-constant distributed quantum advantage for a local problem [2411.03240]. By contrast, in trees many LCL classes have the same locality in deterministic LOCAL and randomized online-LOCAL; this yields corollaries that there is no distributed quantum advantage for those classes on trees, and that problems global in deterministic LOCAL remain global across all stronger models captured by randomized online-LOCAL [2409.13795].

## 5. Mixed local–nonlocal elliptic problems with singular potentials and measure data

In PDE, one major meaning of local potential problems concerns mixed local–nonlocal operators with Hardy-type singular terms. A model case is
\[
-\Delta u + (-\Delta)^s u - \gamma \frac{u}{|x|^2}=f \quad\text{in }\Omega,\qquad
u=0 \quad\text{in }\mathbb R^n\setminus\Omega,
\]
where \(\Omega\subset\mathbb R^n\) is bounded, \(0\in\Omega\), and \(0<s<1\). The natural energy space is
\[
X^{1,2}(\Omega)=\overline{C_c^\infty(\Omega)}^{\|\cdot\|_{X^{1,2}}},
\qquad
\|u\|_{X^{1,2}}^2=\int_{\mathbb R^n}|\nabla u|^2\,dx + C_{n,s}[u]_s^2.
\]
The decisive Hardy fact is that for the mixed norm the optimal Hardy constant remains the classical local one,
\[
\mathcal A_{s,n}=\Lambda_n=\frac{(n-2)^2}{4},
\]
and this value is not attained. Consequently, the critical potential is still \(|x|^{-2}\), not \(|x|^{-2s}\), and the solvability theory closely parallels the purely local Hardy problem. For \(f\in L^m(\Omega)\) with \(m\ge (2^*)'=\frac{2n}{n+2}\) and \(0<\gamma<\Lambda_n\), there is a unique weak solution in \(X^{1,2}(\Omega)\); for lower summability one obtains distributional or duality solutions, and for \(f\in L^1(\Omega)\) solvability is characterized by the condition \(\int_\Omega f\Phi_0<\infty\), where \(\Phi_0\) solves the mixed Hardy equation with right-hand side \(1\) [2407.06763].

A second PDE line studies mixed local–nonlocal \(p\)-Laplace type equations with measure data,
\[
-\operatorname{div}\mathcal A(x,Du)-\mathcal L_\Phi u=\mu,
\]
with \(\mu\) a finite signed Borel measure, \(\mathcal A\) a Carathéodory vector field satisfying \(p\)-growth and ellipticity, and \(K(x,y)\) merely measurable in the nonlocal term. The nonlinear potential theory is formulated through truncated Wolff and Riesz potentials,
\[
{\bf W}^{\mu}_{\beta,p}(x,R)=\int_0^R
\left[\frac{|\mu|(B_\rho(x))}{\rho^{n-\beta p}}\right]^{\!\frac1{p-1}}\frac{d\rho}{\rho},
\qquad
{\bf I}^{\mu}_{\beta}(x,R)=\int_0^R \frac{|\mu|(B_\rho(x))}{\rho^{n-\beta}}\frac{d\rho}{\rho},
\]
together with a nonlocal tail
\[
\operatorname{Tail}(f;x_0,r)=
\left(r^p \int_{\mathbb R^n\setminus B_r(x_0)}\frac{|f(x)|^{p-1}}{|x-x_0|^{n+sp}}\,dx\right)^{\!\frac1{p-1}}.
\]
The paper proves universal pointwise estimates for \(u\) and \(Du\) via Riesz and Wolff potentials, introduces a novel fractional maximum function that captures local and nonlocal features simultaneously, and shows that local oscillations of solutions are determined by local energy averages, the tail term, and nonlinear potentials of \(\mu\) [2510.13269].

A third line considers a quasilinear mixed operator with a mixed interpolated Hardy potential,
\[
\mathcal T(u)= -\Delta_p u + (-\Delta_p)^s u - \mu \frac{u^{p-1}}{|x|^{p\theta}},
\qquad \theta\in[s,1].
\]
The mixed Hardy inequality
\[
\int_\Omega \frac{|u|^p}{|x|^{p\theta}}\,dx
\le \frac{1}{\bar\mu(\theta)}
\left(\int_\Omega |\nabla u|^p\,dx +
\iint_{\mathbb R^{2N}}\frac{|u(x)-u(y)|^p}{|x-y|^{N+sp}}\,dx\,dy\right)
\]
controls the singular local potential by the mixed energy. The paper establishes a concentration–compactness principle in which concentration defects occur only at the singular point \(0\), combines it with Ricceri’s variational principle to obtain nontrivial weak solutions for subcritical nonlinearities, and applies the mountain pass theorem to the superlinear case [2606.29348].

## 6. Mesoscopic, superconducting, and conceptual-DFT formulations

In coarse-grained molecular modeling, a local potential problem arises when the pair interaction must depend on the local environment in order to encode mesoparticle compressibility. The basic ansatz is
\[
\mathcal U_{\mathrm{dd}}(r_{ij},V_{ij})=\mathcal U_{\mathrm{std}}(r_{ij}-\lambda_{ij}),
\qquad
\lambda_{ij}=\lambda(V_{ij}),\qquad
V_{ij}=\frac{V_i+V_j}{2},
\]
where \(V_i\) is a local volume associated with particle \(i\). Two choices are compared: spherical local density \(\rho_i^{\rm s}\), which depends on a kernel \(\chi(r)\) and a cutoff \(r_{\rm cut}\), and Voronoi-based density \(\rho_i^{\rm v}=m/V_i\), where \(V_i\) is the Voronoi cell volume. The Voronoi approach is parameter-free and satisfies \(\overline{\rho^{\rm v}}=\rho_{\rm th}\) exactly. In nitromethane, with one mesoparticle representing one thousand molecules, a quadratic \(\lambda(V)\) fitted against Hugoniot data reduces the RMS pressure error from \(F\approx 3.13\) GPa for \(p=0\) to \(F\approx0.034\) GPa for \(p=2\), and remains accurate on higher-density test states [1310.1288].

In nonequilibrium superconductivity, the local pair potential is
\[
\Delta(\mathbf r,t)=D_0\sum_n u_n(\mathbf r,t)v_n^*(\mathbf r,t)(1-2f_n),
\]
within the time-dependent Bogoliubov–de Gennes equations. The proposed mechanism is to use a control field \(V(\mathbf r,t)\) to localize Bogoliubov quasiparticle amplitudes near a target point \(x_0\), thereby enhancing \(\Delta(x_0,t_0)\) and the associated local effective transition temperature
\[
T_c(x_0,t_0)\approx \frac{\Delta(x_0,t_0)}{1.76k_B}.
\]
For a one-dimensional model with \(V(x,t)=x\,f(t)\) and a multi-frequency \(f(t)\), the paper reports about \(300\%\) enhancement of the local pair potential and local \(T_c\) compared with the initial state. The effect is transient and explicitly constrained by the weak-nonequilibrium regime of the TDBdG equations, pair breaking for \(\hbar\omega_{\max}>2|\Delta|\), decoherence, and phase slips in thin systems [1109.3815].

Conceptual DFT supplies an explicitly critical perspective. The traditional local hardness is defined by
\[
\eta(\mathbf r)=\left.\frac{\delta\mu}{\delta n(\mathbf r)}\right|_{v},
\]
and the analogous local chemical potential by \(\left.\delta E[n]/\delta n(\mathbf r)\right|_v\). The analysis shows that the fixed-\(v\) derivative is only a restricted derivative and is therefore ambiguous; if it is interpreted as a proper constrained derivative, it collapses to a constant value equal to the global hardness \(\eta\), and the local chemical potential likewise becomes the constant \(\mu\). If the external-potential constraint is dropped, the resulting derivative is tied to the asymptotic value of the second derivative of the universal functional,
\[
\frac{\delta\mu[n]}{\delta n(\mathbf r)}
=
\lim_{|\mathbf R|\to\infty}\frac{\delta^2 F[n]}{\delta n(\mathbf r)\,\delta n(\mathbf R)},
\]
which is then undermined by asymptotic divergence. The paper therefore concludes that the traditional approach is incapable of delivering a local hardness indicator and that the parallel local chemical-potential construction fails for the same reason [1107.4249].

These physical and chemical usages make clear that “local potential problems” can denote either a constructive modeling strategy or a diagnosis of conceptual failure. In the mesoparticle and superconducting settings, locality is engineered to make the potential responsive to environment or control. In conceptual DFT, the same localization program breaks down because the supposed local quantity is governed by restricted derivatives or asymptotic functional behavior rather than chemically meaningful sitewise response [1310.1288] [1109.3815] [1107.4249].

Source: https://www.emergentmind.com/topics/local-potential-problems