English

What Trace Powers Reveal About Log-Determinants: Closed-Form Estimators, Certificates, and Failure Modes

Machine Learning 2026-01-21 v1 Machine Learning

Abstract

Computing logdet(A)\log\det(A) for large symmetric positive definite matrices arises in Gaussian process inference and Bayesian model comparison. Standard methods combine matrix-vector products with polynomial approximations. We study a different model: access to trace powers pk=\tr(Ak)p_k = \tr(A^k), natural when matrix powers are available. Classical moment-based approximations Taylor-expand log(λ)\log(\lambda) around the arithmetic mean. This requires λ\AM<\AM|\lambda - \AM| < \AM and diverges when κ>4\kappa > 4. We work instead with the moment-generating function M(t)=\E[Xt]M(t) = \E[X^t] for normalized eigenvalues X=λ/\AMX = \lambda/\AM. Since M(0)=\E[logX]M'(0) = \E[\log X], the log-determinant becomes logdet(A)=n(log\AM+M(0))\log\det(A) = n(\log \AM + M'(0)) -- the problem reduces to estimating a derivative at t=0t = 0. Trace powers give M(k)M(k) at positive integers, but interpolating M(t)M(t) directly is ill-conditioned due to exponential growth. The transform K(t)=logM(t)K(t) = \log M(t) compresses this range. Normalization by \AM\AM ensures K(0)=K(1)=0K(0) = K(1) = 0. With these anchors fixed, we interpolate KK through m+1m+1 consecutive integers and differentiate to estimate K(0)K'(0). However, this local interpolation cannot capture arbitrary spectral features. We prove a fundamental limit: no continuous estimator using finitely many positive moments can be uniformly accurate over unbounded conditioning. Positive moments downweight the spectral tail; K(0)=\E[logX]K'(0) = \E[\log X] is tail-sensitive. This motivates guaranteed bounds. From the same traces we derive upper bounds on (detA)1/n(\det A)^{1/n}. Given a spectral floor rλminr \leq \lambda_{\min}, we obtain moment-constrained lower bounds, yielding a provable interval for logdet(A)\log\det(A). A gap diagnostic indicates when to trust the point estimate and when to report bounds. All estimators and bounds cost O(m)O(m), independent of nn. For m{4,,8}m \in \{4, \ldots, 8\}, this is effectively constant time.

Keywords

Cite

@article{arxiv.2601.12612,
  title  = {What Trace Powers Reveal About Log-Determinants: Closed-Form Estimators, Certificates, and Failure Modes},
  author = {Piyush Sao},
  journal= {arXiv preprint arXiv:2601.12612},
  year   = {2026}
}