Papers
Topics
Authors
Recent
Search
2000 character limit reached

Spectrum Projection Score (SPS)

Updated 17 July 2026
  • Spectrum Projection Score (SPS) is a projection-based metric with dual formulations: one defining the second eigenvalue of a graph Laplacian for clustering and the other computing residual norms for summary alignment in RAG.
  • In clustering, SPS minimizes the second smallest eigenvalue to bisect data into weakly connected subsets and employs microcluster approximations to achieve significant speedups with minimal accuracy loss.
  • In retrieval-augmented generation, SPS measures how well a candidate summary aligns with a reader model’s principal subspace, leading to improved performance in QA tasks.

Spectrum Projection Score (SPS) denotes two distinct projection-based quantities in the arXiv literature. In "Minimum Spectral Connectivity Projection Pursuit" (Hofmeyr et al., 2015), 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" (Hu et al., 8 Aug 2025), 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 LLMs. The identical name therefore does not identify a single standardized metric across fields, but two separate constructions with different semantics and different computational pipelines (Hofmeyr et al., 2015, Hu et al., 8 Aug 2025).

Setting SPS definition Operational role
Projection-pursuit clustering SPS(V)λ2(L(V))\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 SPS(x)=(IP)x2\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 VV 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" (Hofmeyr et al., 2015), let X={x1,,xn}RdX=\{x_1,\dots,x_n\}\subset\mathbb{R}^d and let VRd×rV\in\mathbb{R}^{d\times r} be an orthonormal projection matrix with r<dr<d. Writing pi=Vxip_i=V^\top x_i and P=[p1  pn]Rr×nP=[p_1\ \dots\ p_n]\in\mathbb{R}^{r\times n}, the method forms a weighted similarity graph on the projected points with

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

where k:R+R+k:\mathbb{R}_+\to\mathbb{R}_+ is a smooth, decreasing kernel and SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_20 is the scale or bandwidth parameter. The degree matrix is SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_21 with SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_22, and the unnormalized graph Laplacian is

SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_23

The Spectrum Projection Score of SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_24 is the second smallest eigenvalue of SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_25:

SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_26

The intended interpretation is explicit: SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_27 measures how well connected the graph is, and minimizing SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_28 seeks a projection that bisects SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_29 into two weakly connected pieces. The corresponding optimization problem is

VV0

For VV1, this reduces to finding a unit vector VV2 with VV3 that minimizes VV4 of its induced Laplacian. An equivalent formulation adds a soft orthogonality penalty,

VV5

with VV6 large, to keep columns of VV7 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.

A central theoretical result is the connection between minimum spectral connectivity and maximum-margin separation. In the univariate case, if VV8 and the maximum-margin hyperplane normal VV9 is unique, then

X={x1,,xn}RdX=\{x_1,\dots,x_n\}\subset\mathbb{R}^d0

The paper states that, as X={x1,,xn}RdX=\{x_1,\dots,x_n\}\subset\mathbb{R}^d1, 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 (Hofmeyr et al., 2015).

The computational method is based on generalized gradient-based descent. To evaluate SPS at X={x1,,xn}RdX=\{x_1,\dots,x_n\}\subset\mathbb{R}^d2, the procedure projects points in X={x1,,xn}RdX=\{x_1,\dots,x_n\}\subset\mathbb{R}^d3, builds X={x1,,xn}RdX=\{x_1,\dots,x_n\}\subset\mathbb{R}^d4 in X={x1,,xn}RdX=\{x_1,\dots,x_n\}\subset\mathbb{R}^d5, forms X={x1,,xn}RdX=\{x_1,\dots,x_n\}\subset\mathbb{R}^d6 and X={x1,,xn}RdX=\{x_1,\dots,x_n\}\subset\mathbb{R}^d7, and computes X={x1,,xn}RdX=\{x_1,\dots,x_n\}\subset\mathbb{R}^d8 and the associated eigenvector X={x1,,xn}RdX=\{x_1,\dots,x_n\}\subset\mathbb{R}^d9 in VRd×rV\in\mathbb{R}^{d\times r}0. When VRd×rV\in\mathbb{R}^{d\times r}1 is a simple eigenvalue, standard matrix-perturbation gives

VRd×rV\in\mathbb{R}^{d\times r}2

after which the chain rule is used to backpropagate to VRd×rV\in\mathbb{R}^{d\times r}3. The reported per-step complexity is VRd×rV\in\mathbb{R}^{d\times r}4.

To reduce the quadratic cost, the paper proposes a microcluster approximation. A fast VRd×rV\in\mathbb{R}^{d\times r}5-means clustering such as BIRCH or K-means produces microclusters VRd×rV\in\mathbb{R}^{d\times r}6 with centers VRd×rV\in\mathbb{R}^{d\times r}7 and cardinalities VRd×rV\in\mathbb{R}^{d\times r}8, where VRd×rV\in\mathbb{R}^{d\times r}9. The full graph Laplacian is then approximated by an r<dr<d0 weighted Laplacian r<dr<d1 with

r<dr<d2

and with degree r<dr<d3 on the diagonal. The approximation error is bounded by

r<dr<d4

where r<dr<d5 is the radius of microcluster r<dr<d6 and r<dr<d7. The paper states that for small r<dr<d8 relative to r<dr<d9 this bound is tiny. Computing pi=Vxip_i=V^\top x_i0 of the pi=Vxip_i=V^\top x_i1 Laplacian costs pi=Vxip_i=V^\top x_i2, approximate gradients cost pi=Vxip_i=V^\top x_i3, and in practice pi=Vxip_i=V^\top x_i4 yields pi=Vxip_i=V^\top x_i5–pi=Vxip_i=V^\top x_i6 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 pi=Vxip_i=V^\top x_i7 clusters, the binary projection-pursuit split is applied recursively in a top-down fashion. At each node corresponding to a subset pi=Vxip_i=V^\top x_i8, the method finds the best projection pi=Vxip_i=V^\top x_i9 and then bipartitions P=[p1  pn]Rr×nP=[p_1\ \dots\ p_n]\in\mathbb{R}^{r\times n}0 by the sign or 2-means on the second eigenvector of its Laplacian. Splitting continues on the largest, or most connected, leaf until P=[p1  pn]Rr×nP=[p_1\ \dots\ p_n]\in\mathbb{R}^{r\times n}1 leaves are obtained. The resulting divisive tree permits different nodes to use different projection subspaces P=[p1  pn]Rr×nP=[p_1\ \dots\ p_n]\in\mathbb{R}^{r\times n}2 (Hofmeyr et al., 2015).

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 P=[p1  pn]Rr×nP=[p_1\ \dots\ p_n]\in\mathbb{R}^{r\times n}3 samples, dimension up to approximately P=[p1  pn]Rr×nP=[p_1\ \dots\ p_n]\in\mathbb{R}^{r\times n}4, and clusters P=[p1  pn]Rr×nP=[p_1\ \dots\ p_n]\in\mathbb{R}^{r\times n}5 up to P=[p1  pn]Rr×nP=[p_1\ \dots\ p_n]\in\mathbb{R}^{r\times n}6. 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 P=[p1  pn]Rr×nP=[p_1\ \dots\ p_n]\in\mathbb{R}^{r\times n}7, such as P=[p1  pn]Rr×nP=[p_1\ \dots\ p_n]\in\mathbb{R}^{r\times n}8, microclustered SCPP runs P=[p1  pn]Rr×nP=[p_1\ \dots\ p_n]\in\mathbb{R}^{r\times n}9–Wij=k(pipj/σ),W_{ij}=k(\|p_i-p_j\|/\sigma),0 faster than the full Wij=k(pipj/σ),W_{ij}=k(\|p_i-p_j\|/\sigma),1 method while preserving more than Wij=k(pipj/σ),W_{ij}=k(\|p_i-p_j\|/\sigma),2 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 (Hu et al., 8 Aug 2025). 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 Wij=k(pipj/σ),W_{ij}=k(\|p_i-p_j\|/\sigma),3 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

Wij=k(pipj/σ),W_{ij}=k(\|p_i-p_j\|/\sigma),4

The singular value decomposition is

Wij=k(pipj/σ),W_{ij}=k(\|p_i-p_j\|/\sigma),5

The method retains the top Wij=k(pipj/σ),W_{ij}=k(\|p_i-p_j\|/\sigma),6 singular values whose cumulative energy is at least Wij=k(pipj/σ),W_{ij}=k(\|p_i-p_j\|/\sigma),7 of Wij=k(pipj/σ),W_{ij}=k(\|p_i-p_j\|/\sigma),8, with corresponding left singular vectors Wij=k(pipj/σ),W_{ij}=k(\|p_i-p_j\|/\sigma),9. The projection matrix onto the principal subspace is

k:R+R+k:\mathbb{R}_+\to\mathbb{R}_+0

For a candidate summary of k:R+R+k:\mathbb{R}_+\to\mathbb{R}_+1 tokens, the reader produces penultimate hidden states k:R+R+k:\mathbb{R}_+\to\mathbb{R}_+2. The summary “area” vector is the element-wise max-pooled vector

k:R+R+k:\mathbb{R}_+\to\mathbb{R}_+3

SPS is then defined by decomposing k:R+R+k:\mathbb{R}_+\to\mathbb{R}_+4 into its component inside the principal subspace and its residual outside it:

k:R+R+k:\mathbb{R}_+\to\mathbb{R}_+5

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-k:R+R+k:\mathbb{R}_+\to\mathbb{R}_+6 documents, compresses them into an initial summary k:R+R+k:\mathbb{R}_+\to\mathbb{R}_+7, and computes two simple norms on k:R+R+k:\mathbb{R}_+\to\mathbb{R}_+8’s representation in the reader:

  • k:R+R+k:\mathbb{R}_+\to\mathbb{R}_+9
  • SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_200

with ratio SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_201. If SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_202, the system skips sampling and feeds SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_203 to the reader. Otherwise it samples SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_204 diverse summaries. For text-to-text, the sampling uses stochastic decoding with temperature SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_205 and repetition_penalty SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_206; for text-to-embedding, it injects small Gaussian probes, keeps SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_207 probes with the smallest semantic-deviation scores SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_208, 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 SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_209, SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_210, and AUROC of approximately SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_211. In end-to-end QA on HotpotQA with LLaMA-3.1, naïve retrieval-to-reader gives EM/F1 of SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_212, CompAct summary gives SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_213, and xCompress + SPS gives SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_214, corresponding to a SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_215 EM and SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_216 F1 improvement over CompAct. In text-to-embedding with Mistral-7B + xRAG, xCompress + SPS improves EM/F1 on HotpotQA from SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_217 to SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_218.

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, SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_219 variance retention is optimal, and performance saturates at SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_220. Practical details include precomputing the SVD or PCA once per reader model, an SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_221 projection cost per candidate at inference, adaptive norm-guided filtering that skips sampling for approximately SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_222 of queries and cuts SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_223 of sampling cost with minimal impact on accuracy, greedy decoding for the final answer, and an overhead of approximately SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_224–SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_225 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 SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_226 residual norm after projecting a max-pooled summary representation onto a reader-specific PCA subspace (Hofmeyr et al., 2015, Hu et al., 8 Aug 2025).

Their optimization logic also differs. The clustering method minimizes SPS over projection matrices SPS(x)=(IP)x2\mathrm{SPS}(x)=\|(I-P)x\|_227 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.

Definition Search Book Streamline Icon: https://streamlinehq.com
References (2)

Topic to Video (Beta)

No one has generated a video about this topic yet.

Whiteboard

No one has generated a whiteboard explanation for this topic yet.

Follow Topic

Get notified by email when new papers are published related to Spectrum Projection Score (SPS).