Anchor-Based Gauge Map

• Anchor-based gauge maps are explicit linear reconstruction schemes that approximate global diffusion geometry using local anchor-to-node shortest-path distances.
• They employ a trilateration approach by solving a linear system with transformed distance measurements to efficiently estimate spectral coordinates.
• In practice, these methods achieve near-spectral performance in tasks like drug–drug interaction prediction while significantly reducing computational costs.

An anchor-based gauge map, or trilateration map, is an explicit linear reconstruction scheme that recovers truncated diffusion (spectral) coordinates of nodes in a graph from anchor-to-node shortest-path distance encodings. This method bridges two foundational positional encoding paradigms in graph machine learning: spectral coordinates from Laplacian or diffusion operators and anchor-based shortest-path distance encodings. Anchor-based gauge maps systematically approximate global diffusion geometry using local distance measurements to a fixed anchor set, yielding accuracy guarantees in both pointwise and matrix norms. In recent work, this framework has provided a scalable and theoretically rigorous alternative to expensive spectral eigendecompositions, showing practical success in molecular graph learning tasks such as drug–drug interaction (DDI) prediction (Yan et al., 8 Jan 2026).

1. Diffusion Geometry and Spectral Coordinates

Let $G=(V,E)$ be a connected, undirected graph of $n=|V|$ nodes. The normalized Laplacian operator is $L=I-(1/r)A$ for $r$-regular graphs ($A$: adjacency matrix). The diffusion operator at time $t>0$ is $K_t = e^{-tL}$, whose entries $k_t(u,v)$ define the graph’s heat kernel. The diffusion distance is

$$d_t(u,v)^2 = \sum_{j=1}^n e^{-2t\lambda_j}\left(\varphi_j(u)-\varphi_j(v)\right)^2 = k_{2t}(u,u) + k_{2t}(v,v) - 2k_{2t}(u,v),$$

where $(\lambda_j,\varphi_j)$ are the Laplacian eigenpairs. Truncating to the first $m$ nontrivial modes ($j=2,...,m+1$; $0=\lambda_1 < \lambda_2\leq \dots \leq\lambda_n$) yields the $m$-dimensional diffusion map:

$$\psi^{(m)}(v) := \left( e^{-t\lambda_2} \varphi_2(v), ..., e^{-t\lambda_{m+1}} \varphi_{m+1}(v) \right) \in \mathbb{R}^m.$$

The truncated diffusion distance is $$d^{(m)}(u,v) := \|\psi^{(m)}(u) - \psi^{(m)}(v)\|_2$$.

2. Anchor-Based Distance Encodings

Anchor-based methods use a set $S = \{a_1,\ldots,a_{m+1}\} \subset V$ of $m+1$ anchors (theoretical analysis uses i.i.d. uniform anchors on random regular graphs; practical implementations on molecular graphs use $k \gg m$ anchors via farthest-point sampling in the shortest-path metric). For node $v\in V$, the anchor distance vector is

$$\delta(v|S) := \left( \text{SPD}(a_1,v), ..., \text{SPD}(a_{m+1},v) \right)^T \in \mathbb{R}^{m+1},$$

where $\text{SPD}(a,v)$ is the shortest-path distance. Optionally, a monotone scalar transform $\nu:\mathbb{R}_{\geq 0}\rightarrow \mathbb{R}_{>0}$—such as $\nu(d)=e^{-d}$ or $\nu(d)=\log(1+d)$—is applied elementwise to obtain $\nu(\delta(v|S)) \in \mathbb{R}^{m+1}$. This yields a low-dimensional anchor-based positional encoding.

3. Construction of the Trilateration (Gauge) Map

Given, for each anchor $a_i$, its spectral coordinates $p_i:=\psi^{(m)}(a_i)\in\mathbb{R}^m$, the true spectral coordinate $z:=\psi^{(m)}(v)$ for node $v$ satisfies $m+1$ sphere equations:

$$\|z-p_i\|_2^2 = d^{(m)}(a_i,v)^2, \quad i=1,\ldots,m+1.$$

Let $r_i := d^{(m)}(a_i,v)$. Subtracting the equation for $i=m+1$ from the first $m$ yields the linear system

$A z = b(r)$

with

$$A = 2 \begin{bmatrix} (p_1-p_{m+1})^T \ \vdots\ (p_m-p_{m+1})^T \end{bmatrix} \in \mathbb{R}^{m\times m},$$

and

$$b_i(r) = \|p_i\|_2^2 - \|p_{m+1}\|_2^2 + r_{m+1}^2 - r_i^2, \quad i=1,\ldots,m.$$

Under affine-independence of the anchors, $A$ is invertible and

$z = \psi^{(m)}(v) = A^{-1}b(r).$

In practice, $d^{(m)}(a_i,v)$ is unavailable, and one substitutes transformed raw distances $\rho_i := \nu(\text{SPD}(a_i,v))$ into $b(\rho)$, defining the trilateration (gauge) map:

$$T(v) := A^{-1} b(\nu(\delta(v|S))) \in \mathbb{R}^m,$$

which approximates $\psi^{(m)}(v)$. This "trilateration map" reconstructs the node's spectral embedding from its anchor distances.

4. Theoretical Guarantees on Reconstruction

On $r$-regular graphs, the connection between shortest-path and diffusion distances is formalized as follows. For $u,v$ with $\text{SPD}(u,v)\leq R$:

$$|d_t(u,v)-\chi(\text{SPD}(u,v))| \leq \epsilon_{\text{geom}}$$

(coupling argument), and

$$|d^{(m)}(u,v)-d_t(u,v)| \leq \epsilon_{\text{tail}}$$

(from truncation error). Thus, Theorem 2 states that for all $u,v$ with $\text{SPD}(u,v)\leq R$ and any strictly increasing $\nu$:

$$|d^{(m)}(u,v)-\nu(\text{SPD}(u,v))| \leq \epsilon_{\text{geom}}+\epsilon_{\text{tail}}.$$

Injecting $\nu(\text{SPD})$ into the linear system yields

$$\|\psi^{(m)}(v)-T(v)\|_2 \leq \|A^{-1}\|_{\text{op}}\cdot O(\sqrt{m})\cdot (4\max_{d\leq R}\nu(d) + 2) \cdot (\epsilon_{\text{geom}}+\epsilon_{\text{tail}})$$

with high probability (Theorem 5: Pointwise Reconstruction), and for the distance matrices $D_{\text{SPD}}, D^{(m)}$:

$$\|D^{(m)}-\nu(D_{\text{SPD}})\|_F \leq n(m+1)O(\epsilon_{\text{geom}}+\epsilon_{\text{tail}})$$

(Theorem 6: Frobenius-Gap). Thus, the trilateration map recovers leading $m$ diffusion modes up to controlled errors. The conditioning of $A$ governs stability, with heavy-tailed $\log\text{cond}(A)$ but only $\sim$30% exceeding $\text{cond}=10^3$ on molecular graphs.

5. Practical Implementation and Performance

On DrugBank molecular graphs ($|V|\in[15,200]$), $k=32$ anchors per molecule are selected by farthest-point sampling. The preferred transform is $\nu(d)=\exp(-d)$ (ablation supports its superiority). For exact theoretical trilateration, $m=8$ and $m+1=9$ anchors are sufficient, but diagnostics and Nyström approximations use the full $k$ vectors. The Nyström kernel approximation $K_t\approx \widetilde{K}=C(W+\lambda I)^{-1}C^T$ employs $C_{v,i} = \nu(\text{SPD}(v,a_i))$, $W_{i,j} = \nu(\text{SPD}(a_i,a_j))$, and Tikhonov regularization ($\lambda=10^{-6}$).

Diagnostic results on 80 graphs include:

• Relative Frobenius error of $K$ vs $\widetilde{K}$: $0.024 \pm 0.021$ (median $0.020$)
• Coordinate MSE after Procrustes ($\psi^{(m)}$ vs $T$): $3.9\times 10^{-4}$
• Pearson correlation ($d_t$ vs $\nu(\text{SPD})$): $0.988$ average

For DDI prediction (shared GNP backbone), various positional encoding schemes were benchmarked:

Method (DrugBank)	AUROC / F1
NoPE (no PE)	0.890 / 0.820
DE ($k$=16, $\nu(d)=e^{-d}$)	0.976 / 0.928
LapPE ($m$=8 Lap eigvecs)	0.980 / 0.934
RWSE (5-step walk)	0.892 / 0.812
HKS (5 heat times)	0.863 / 0.787

The anchor-based distance encoding with Nyström closely matches full spectral performance and consistently exceeds the no-encoding baseline. Similar results hold on ChCh-Miner. Ablation studies confirm $\nu(d)=\exp(-d)$ and increased $k$ improve or maintain performance.

6. Connection to Nyström Approximation and Scalability

The same anchor set used for gauge mapping naturally enables a Nyström-type approximation of the full diffusion kernel $K_t$, leveraging anchor-to-node affinities. Here, anchor-based encodings not only reconstruct coordinates via trilateration but also generate a scalable approximate kernel with strong empirical alignment to spectral embeddings. On a 155-node DrugBank molecule, DE-Nyström (with $k=32$) matches the exact Gaussian-kernel eigenmap with mean nodewise error $\approx 7.8\times 10^{-3}$ and maximum $\approx 1.15\times 10^{-1}$.

A plausible implication is that anchor-based gauge maps, coupled with the Nyström scheme, provide a principled and computationally tractable alternative to eigendecomposition in large-scale graph learning, especially where spectral coordinates are costly or impossible to obtain exactly.

7. Summary and Significance

The anchor-based gauge map, as formalized by Yan et al. (Yan et al., 8 Jan 2026), delivers a rigorous and practical solution for approximating spectral positional signals from local distance measurements to anchors. The explicit trilateration reconstruction guarantees recovery of leading diffusion modes with bounded error both pointwise and in matrix norms. On molecular graphs, this yields nearly equivalent downstream performance to full Laplacian spectral encodings in DDI prediction tasks, while offering orders-of-magnitude gains in scalability. These findings establish anchor-based distance encodings, interpreted through the lens of diffusion geometry, as a powerful and flexible tool for geometric deep learning on graphs.