Daniel J. Korchinski, Dhruva Karkada, Yasaman Bahri and Matthieu Wyart (2025), Advances in Neural Information Processing Systems 38 (NeurIPS 2025), pp. 68995–69023 — ARXIV-2505.18651, DOI-10.52202/085713-2075

Key takeaways

  • The model. A word is α ∈ {−1, +1}^d, one entry per binary attribute (masculine/feminine, royal/not, singular/plural). Attribute k has incidence p_k and a signal strength s_k ∈ [0, 1], and P(i, j) = P(i)P(j) ∏_k P^(k)(α_i^k, α_j^k) with a 2 × 2 factor P^(k) fixed by normalisation up to s_k (Eqs. 3–4). In it the ratio postulate of earlier accounts, p(χ|king)/p(χ|queen) = p(χ|man)/p(χ|woman) for every context word χ, holds exactly. The noiseless model assumes it rather than deriving it.
  • The ratio matrix M is a Kronecker product (Theorem, proved by induction in A1): its eigenvectors are tensor products of the 2 × 2 factors’ eigenvectors and its eigenvalues the products of theirs. For weak signals the spectrum falls into bands: a constant mode, d attribute modes v_k(i) ∝ α_i^k with eigenvalue ∝ s_k, then about d²/2 modes ∝ s_k s_k′ that encode products of attributes, and so on. Embedding M gives exact analogies only if the s_k are small and narrowly spread and K ≤ d + 1, so that the kept eigenvectors are affine in the attributes. With broadly spread s_k the bands interleave and accuracy rises and then falls with K.
  • The PMI is bilinear in the attributes. log M = δ11ᵀ + ADAᵀ + Aη1ᵀ + 1ηᵀAᵀ, with A the words-by-attributes matrix (Eq. 12). Its rank is at most d + 1, its eigenvectors are affine in the attributes, and a spectral embedding of it satisfies W_A − W_B + W_C = W_D whenever α_D = α_A − α_B + α_C, whatever the s_k (Results 1–4). This is the paper’s account of why PMI targets beat M (observation iii) and why accuracy saturates once K reaches d. Appendix A5 keeps rank ≤ d + 1 and the affine eigenvectors with pairwise-correlated attributes; the eigenvectors then mix attributes.
  • Robustness, argued and simulated. Symmetric i.i.d. noise of scale σ in the PMI has spectral norm about 2σ2^(d/2), against semantic eigenvalues of order 2^d, so analogies survive until K ≈ 2^(d/2)/σ noise modes are admitted (§7, A4, Fig. 10). Pruning the vocabulary to m ≫ d random words keeps the top-d spectrum up to a factor m/2^d (Marchenko–Pastur, §8). Setting to zero the PMI of every pair that differs only in one attribute (all king–queen, man–woman pairs) is a sparse perturbation of norm about 2^(d/2), and the same relation survives in the embedding (§9, Fig. 3). Every step from small perturbation to stable eigenvectors goes through Weyl’s eigenvalue inequality alone.
  • Wikipedia. Spectral embeddings of M and of log(M + 10⁻²) over 10,000 words, scored on Mikolov et al.’s analogies, reproduce observations (ii)–(iv) qualitatively. The PMI’s spectrum is broad and close to log-normal, which the paper uses against Arora et al.’s isotropic latent space (ANTH-LIT-613). Zeroing the PMI of an analogy family’s own pairs costs little accuracy (Fig. 3b). No attribute is identified in the data, and no trained word2vec or GloVe model is examined.

Standing in the record

Filed on 2026-10-09 at the owner’s request: LIT-860 (Karkada, Korchinski, Nava, Wyart and Bahri 2026, read in NOTE-664) builds on it, and neither record held it. LIT-860 extends its binary-attribute model to continuous attributes, and its Appendix D (Theorem 5) combines the two. The paper was read on its own merits (NOTE-665), and the reading produces THEORY-185.

It sits between two readings the record already has. LIT-855 (THEORY-182) says what a word2vec-like model factorises, M*, which agrees with the PMI to third order; this paper says what geometry the PMI has when co-occurrence is generated by independent binary attributes. LIT-855’s analogy account is a fitted spiked random matrix; this one is derived, but for an idealised spectral embedding, not a trained one. LIT-860 (THEORY-183) does for continuous latent variables what this paper does for binary ones.