Dhruva Karkada, James B. Simon, Yasaman Bahri and Michael R. DeWeese (2025), Advances in Neural Information Processing Systems 38 (NeurIPS 2025), pp. 4507–4540 — ARXIV-2502.09863 (first posted 14 February 2025 as Solvable Dynamics of Self-Supervised Word Embeddings and the Emergence of Analogical Reasoning)

Key takeaways

  • The target is a deviation from independence. With tied weights and the quartic Maclaurin approximation of the word2vec loss, the objective is (1/4)Σ G_ij (W Wᵀ − M*)²_ij + const, where M_ij = (Ψ⁺P_ij − Ψ⁻P_iP_j) / (½(Ψ⁺P_ij + Ψ⁻P_iP_j)) and G_ij = Ψ⁺P_ij + Ψ⁻P_iP_j (Proposition 2). M agrees with PMI to third order in its argument (Appendix D.1) but is bounded, |M*_ij| ≤ 2 when the two reweightings are equal, so least squares no longer over-weights rarely co-occurring pairs.
  • Which low-rank solution is learned. When G is constant (Setting 3.1: symmetric Ψ⁺, Ψ⁻ with Ψ⁺P_ij + Ψ⁻P_iP_j = g), the problem is unweighted factorisation, and the global minima are the top-d eigenvectors of M* scaled by the square roots of their eigenvalues, up to a right orthogonal transformation (Theorem 1, by completing the square and Eckart–Young). For G of rank one the minimiser is also known (Proposition 2); for general G it is weighted low-rank approximation, which is NP-hard, and the paper does not treat it.
  • The order in which it is learned. From aligned initialisation each singular value follows Saxe et al.’s sigmoid, s_k²(t) = s_k²(0)λ_k e^{λ_k t} / (λ_k + s_k²(0)(e^{λ_k t} − 1)), realised at τ_k = (1/λ_k) ln(λ_k/s_k²(0)) (Lemma 3.1, proved). From vanishing random initialisation, the embeddings first align silently with the eigenbasis and then follow the same steps (Result 3). Result 3 is derived by discarding off-diagonal couplings and checking that they stay small; the authors call it a conjecture to be made rigorous. Early stopping and small d therefore both limit the number of eigen-features, not their character.
  • It is a good proxy for word2vec. On 2 billion Wikipedia tokens (V = 10,000, d = 200): Google analogies 68.0% for word2vec, 65.1% for the trained quadratic model, 66.3% for explicit factorisation of M*, 50.6% for truncated PPMI and 8.4% for truncated PMI. The leading principal directions of word2vec and the quadratic model are the same interpretable topics. An ablation shows the clean stepwise separation comes from the reweighting, not from the quartic approximation.
  • Analogy structure, empirically. Difference vectors within a semantic class (male − female, verb − past tense) concentrate on a few eigen-features, and the spectrum of their Gram matrix looks like a Marchenko–Pastur bulk plus one spike, the mean difference vector. The spike’s signal-to-noise ratio first rises and then falls as rank is added, and its maximum tracks analogy accuracy across classes. This part is observation with a manual fit, not theorem.

Standing in the record

Filed on 2026-10-09 at the owner’s request, with no stated context; it came in a batch with arXiv 2509.24914, 2504.12916, 2205.10343 and 2410.17770. It was read on its own merits (NOTE-660), and the reading produces THEORY-182.

It carries anthology-candidate: it is a theory of the training dynamics of an ML algorithm, published at NeurIPS, of the kind the anthology holds alongside word2vec (ANTH-LIT-604) and Levy and Goldberg (ANTH-LIT-612). It stays here for now beside the record’s spectral self-supervised line (LIT-242, THEORY-001, THEORY-009), which it extends to language data.