---
title: Spectrum Projection Score (SPS)
url: https://www.emergentmind.com/topics/spectrum-projection-score-sps
type: topic
---

# Spectrum Projection Score (SPS)

Spectrum Projection Score (SPS) denotes two distinct projection-based quantities in the arXiv literature. In "Minimum Spectral Connectivity Projection Pursuit" [1509.01546], SPS is the second smallest eigenvalue of the graph Laplacian of projected data, used to find low-dimensional projections that maximize the separability of a binary partition of an unlabelled dataset. In "Spectrum Projection Score: Aligning Retrieved Summaries with Reader Models in Retrieval-Augmented Generation" [2508.05909], SPS is a supervision-free score defined from the residual of a summary representation outside a reader model’s principal subspace, used to rank retrieved summaries in retrieval-augmented generation. Both uses are projection-centric and spectrum-based, but they differ in mathematical object, optimization target, and application domain.

## 1. Dual usage of the term

The 2015 usage, developed by Hofmeyr et al., belongs to projection pursuit and spectral graph theory; the 2025 usage belongs to retrieval-augmented generation and inference-time control for large language models. The identical name therefore does not identify a single standardized metric across fields, but two separate constructions with different semantics and different computational pipelines [1509.01546] [2508.05909].

| Setting | SPS definition | Operational role |
|---|---|---|
| Projection-pursuit clustering | $\mathrm{SPS}(V)\equiv \lambda_2(L(V))$ | Minimize over projection subspaces to obtain a binary split into two weakly connected pieces |
| Retrieval-augmented generation | $\mathrm{SPS}(x)=\|(I-P)x\|_2$ | Rank candidate summaries; lower SPS indicates stronger alignment with the reader’s principal subspace |

This terminological overlap is a common source of confusion. In the first case, the score is attached to a projection matrix $V$ and a similarity graph on projected data. In the second, the score is attached to a candidate summary represented inside a reader LLM.

## 2. SPS in minimum spectral connectivity projection pursuit

In "Minimum Spectral Connectivity Projection Pursuit" [1509.01546], let $X=\{x_1,\dots,x_n\}\subset\mathbb{R}^d$ and let $V\in\mathbb{R}^{d\times r}$ be an orthonormal projection matrix with $r<d$. Writing $p_i=V^\top x_i$ and $P=[p_1\ \dots\ p_n]\in\mathbb{R}^{r\times n}$, the method forms a weighted similarity graph on the projected points with

$$
W_{ij}=k(\|p_i-p_j\|/\sigma),
$$

where $k:\mathbb{R}_+\to\mathbb{R}_+$ is a smooth, decreasing kernel and $\sigma>0$ is the scale or bandwidth parameter. The degree matrix is $D=\operatorname{diag}(d_1,\dots,d_n)$ with $d_i=\sum_j W_{ij}$, and the unnormalized graph Laplacian is

$$
L(V)=D-W.
$$

The Spectrum Projection Score of $V$ is the second smallest eigenvalue of $L(V)$:

$$
\mathrm{SPS}(V)\equiv \lambda_2(L(V)).
$$

The intended interpretation is explicit: $\lambda_2$ measures how well connected the graph is, and minimizing $\lambda_2$ seeks a projection that bisects $X$ into two weakly connected pieces. The corresponding optimization problem is

$$
\min_{V\in\mathbb{R}^{d\times r}} \lambda_2(L(V))
\quad\text{subject to}\quad
\|V^\top V-I\|=0.
$$

For $r=1$, this reduces to finding a unit vector $v\in\mathbb{R}^d$ with $\|v\|=1$ that minimizes $\lambda_2$ of its induced Laplacian. An equivalent formulation adds a soft orthogonality penalty,

$$
\min_V \lambda_2(L(V))+\omega\sum_{i\ne j}(v_i^\top v_j)^2,
$$

with $\omega$ large, to keep columns of $V$ nearly orthonormal.

The paper characterizes this objective as non-convex and non-smooth. Its significance lies in replacing variance-maximizing or reconstruction-based projections with a projection chosen directly for unsupervised separability.

## 3. Theoretical link to maximum-margin structure and computational method

A central theoretical result is the connection between minimum spectral connectivity and maximum-margin separation. In the univariate case, if $v^*(\sigma)=\arg\min_{\|v\|=1}\lambda_2(L(v;\sigma))$ and the maximum-margin hyperplane normal $v_0$ is unique, then

$$
\lim_{\sigma\to 0^+} v^*(\sigma)=v_0.
$$

The paper states that, as $\sigma\to 0$, the kernel becomes sensitive only to very small inter-point gaps, and proves that the optimal univariate projection based on spectral connectivity converges to the normal vector of the maximum-margin hyperplane through the data. Lemmas giving upper and lower spectral bounds in terms of the largest gap underpin this result [1509.01546].

The computational method is based on generalized gradient-based descent. To evaluate SPS at $V$, the procedure projects points in $O(ndr)$, builds $W_{ij}=k(\|p_i-p_j\|/\sigma)$ in $O(n^2r)$, forms $D$ and $L=D-W$, and computes $\lambda_2(L)$ and the associated eigenvector $u$ in $O(n^2)$. When $\lambda_2$ is a simple eigenvalue, standard matrix-perturbation gives

$$
d\lambda_2=u^\top\, dL\, u,
$$

after which the chain rule is used to backpropagate to $V$. The reported per-step complexity is $O(n^2r+ndr+n^2)\approx O(n^2r)$.

To reduce the quadratic cost, the paper proposes a microcluster approximation. A fast $m$-means clustering such as BIRCH or K-means produces microclusters $\{C_1,\dots,C_m\}$ with centers $c_j$ and cardinalities $n_j$, where $m\ll n$. The full graph Laplacian is then approximated by an $m\times m$ weighted Laplacian $L^c$ with

$$
B_{ij}=\sqrt{n_i n_j}\,k(\|V^\top c_i-V^\top c_j\|/\sigma),
$$

and with degree $N_{ii}=\sum_k n_k\,k(\cdots)$ on the diagonal. The approximation error is bounded by

$$
\frac{|\lambda_2(L)-\lambda_2(N-B)|}{\lambda_2(L)}
\le
\max_{i\ne j}
\left\{
1-\frac{k(D_{ij}-\rho_i-\rho_j)}{k(D_{ij})},
\frac{k(D_{ij}+\rho_i+\rho_j)}{k(D_{ij})}-1
\right\},
$$

where $\rho_j$ is the radius of microcluster $C_j$ and $D_{ij}=\|V^\top c_i-V^\top c_j\|$. The paper states that for small $\rho$ relative to $\sigma$ this bound is tiny. Computing $\lambda_2$ of the $m\times m$ Laplacian costs $O(m^2)$, approximate gradients cost $O(m^2r)$, and in practice $m\approx 200$ yields $50$–$100\times$ speedups with negligible loss in clustering quality.

## 4. Hierarchical divisive clustering and empirical behavior of SCPP

The SPS objective is embedded in a divisive hierarchical model to obtain more than two clusters. To obtain $K$ clusters, the binary projection-pursuit split is applied recursively in a top-down fashion. At each node corresponding to a subset $C\subset X$, the method finds the best projection $V_C$ and then bipartitions $C$ by the sign or 2-means on the second eigenvector of its Laplacian. Splitting continues on the largest, or most connected, leaf until $K$ leaves are obtained. The resulting divisive tree permits different nodes to use different projection subspaces $V$ [1509.01546].

The empirical evaluation covers Optical Digits, Pen Digits, Multiple-Feature Digits, Satellite, Image Segmentation, Breast-Cancer, Synthetic Chart, Isolet, Dermatology, Yeast, Smartphone Activity, Yale Faces, and Phoneme, with up to approximately $10\,000$ samples, dimension up to approximately $1\,200$, and clusters $K$ up to $26$. The competing methods are SC, SC\_PC / SC\_IC, DRSC, and iSVR\_G. The evaluation metrics are Purity and Normalized Mutual Information.

The reported result is that SCPP, the SPS-based method, attains the highest Purity and NMI on over half the data sets and is never far behind the best on the rest. Across 14 benchmarks the average relative Purity and NMI of SCPP exceed all competitors. On large $n$, such as $n\approx 10\,000$, microclustered SCPP runs $10$–$100\times$ faster than the full $n^2$ method while preserving more than $98\%$ of clustering accuracy. The paper summarizes this as state-of-the-art clustering performance across a wide range of data sets.

## 5. SPS in retrieval-augmented generation and xCompress

In the 2025 RAG formulation, SPS is introduced as a lightweight, supervision-free metric for measuring the semantic alignment of a retrieved summary with a reader model’s hidden representation [2508.05909]. The motivation is that token-level perplexity and single mean-pooled embeddings tend to favor “typical” or centrally clustered tokens rather than informative boundary tokens. SPS replaces this point-based view with an area-based view.

Let $D$ be the hidden dimension of the reader’s penultimate layer. From either the input-embedding matrix or a large bank of hidden representations, one forms

$$
W\in\mathbb{R}^{D\times M}.
$$

The singular value decomposition is

$$
W=U\,\Sigma\,V^\top.
$$

The method retains the top $r$ singular values whose cumulative energy is at least $95\%$ of $\operatorname{trace}(\Sigma^2)$, with corresponding left singular vectors $U_p\in\mathbb{R}^{D\times r}$. The projection matrix onto the principal subspace is

$$
P=U_pU_p^\top.
$$

For a candidate summary of $n$ tokens, the reader produces penultimate hidden states $h_1,\dots,h_n\in\mathbb{R}^D$. The summary “area” vector is the element-wise max-pooled vector

$$
x=\max_{i=1,\dots,n} h_i \in \mathbb{R}^D.
$$

SPS is then defined by decomposing $x$ into its component inside the principal subspace and its residual outside it:

$$
\mathrm{SPS}(x)=\|(I-P)x\|_2.
$$

A lower value means that the summary’s salient envelope lies mostly within the reader’s core directions.

This score is integrated into xCompress, which sits between the retriever or compressor LLM and the final reader LLM. For each query, the system retrieves top-$N$ documents, compresses them into an initial summary $S_0$, and computes two simple norms on $S_0$’s representation in the reader:

- $L2\_mean=\|\operatorname{mean\_pool}(S_0\_\text{hidden\_states})\|_2$
- $L1\_max=\|\operatorname{max\_pool}(S_0\_\text{hidden\_states})\|_1$

with ratio $R=L2\_mean/L1\_max$. If $R>\tau$, the system skips sampling and feeds $S_0$ to the reader. Otherwise it samples $K$ diverse summaries. For text-to-text, the sampling uses stochastic decoding with temperature $=1.0$ and repetition\_penalty $=1.2$; for text-to-embedding, it injects small Gaussian probes, keeps $M$ probes with the smallest semantic-deviation scores $S\_\text{probe}$, plus the original. SPS is computed for each candidate, and the summary with the lowest SPS is selected for final answer generation.

The experimental validation covers HotpotQA, 2WikiMultiHop, Natural Questions, TriviaQA, and Musique, using LLaMA-3.1-8B-Instruct, Gemma3-12B-Instruct, Qwen3-8B, and Mistral-7B. The paper reports that PPL and LongPPL have near-zero or negative Pearson correlations with downstream Exact Match and F1, whereas SPS shows strong positive correlations and higher AUROC. For example, on HotpotQA, SPS achieves $\operatorname{PCC}(\mathrm{EM})=0.643$, $\operatorname{PCC}(\mathrm{F1})=0.753$, and AUROC of approximately $0.55$. In end-to-end QA on HotpotQA with LLaMA-3.1, naïve retrieval-to-reader gives EM/F1 of $19.6/28.9$, CompAct summary gives $34.0/43.2$, and xCompress + SPS gives $37.6/47.9$, corresponding to a $+3.6$ EM and $+4.7$ F1 improvement over CompAct. In text-to-embedding with Mistral-7B + xRAG, xCompress + SPS improves EM/F1 on HotpotQA from $5.2/16.6$ to $7.6/20.1$.

The ablations report that max-pooling outperforms mean-pooling and last-token pooling by 2–4 points EM/F1, penultimate-layer embeddings give the best downstream performance, $95\%$ variance retention is optimal, and performance saturates at $K=5$. Practical details include precomputing the SVD or PCA once per reader model, an $O(D\cdot r)$ projection cost per candidate at inference, adaptive norm-guided filtering that skips sampling for approximately $30\%$ of queries and cuts $70\%$ of sampling cost with minimal impact on accuracy, greedy decoding for the final answer, and an overhead of approximately $10$–$20\%$ over a single-pass RAG pipeline on standard GPU hardware such as NVIDIA A100.

## 6. Comparative interpretation and disambiguation

The two SPS definitions are mathematically unrelated beyond a broad reliance on projection and spectrum. In the clustering formulation, SPS is a second-smallest-eigenvalue objective on a graph Laplacian constructed from pairwise similarities in a projected space. In the RAG formulation, SPS is an $\ell_2$ residual norm after projecting a max-pooled summary representation onto a reader-specific PCA subspace [1509.01546] [2508.05909].

Their optimization logic also differs. The clustering method minimizes SPS over projection matrices $V$ to discover unsupervised binary partitions, and then composes these splits into a divisive hierarchy. The RAG method ordinarily holds the reader subspace fixed, computes SPS for each candidate summary, and selects the summary with the lowest score inside xCompress. In both cases, lower values are preferred, but “lower” means different things: weaker graph connectivity in the first case, and smaller residual outside the reader’s principal directions in the second.

A common misconception is that SPS denotes a single canonical score transferable across application areas. The published record summarized here does not support that interpretation. A more precise usage is to specify either the spectral-connectivity SPS of projection pursuit clustering or the reader-alignment SPS used in retrieval-augmented generation. This suggests that explicit citation is necessary whenever the acronym appears without context.

Source: https://www.emergentmind.com/topics/spectrum-projection-score-sps