What Trace Powers Reveal About Log-Determinants: Closed-Form Estimators, Certificates, and Failure Modes
Abstract
Computing 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 , natural when matrix powers are available. Classical moment-based approximations Taylor-expand around the arithmetic mean. This requires and diverges when . We work instead with the moment-generating function for normalized eigenvalues . Since , the log-determinant becomes -- the problem reduces to estimating a derivative at . Trace powers give at positive integers, but interpolating directly is ill-conditioned due to exponential growth. The transform compresses this range. Normalization by ensures . With these anchors fixed, we interpolate through consecutive integers and differentiate to estimate . 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; is tail-sensitive. This motivates guaranteed bounds. From the same traces we derive upper bounds on . Given a spectral floor , we obtain moment-constrained lower bounds, yielding a provable interval for . A gap diagnostic indicates when to trust the point estimate and when to report bounds. All estimators and bounds cost , independent of . For , 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}
}