Theoretical Analysis of Thermodynamic Matrix Inversion: First-order Equivalence to Preconditioned Gradient Descent and Implications for Analog Computing

Abstract

Recent research has demonstrated the possibility of exploiting the thermodynamics of coupled electrical oscillators to implement computational tasks such as matrix inversion. While physical implementations rely on thermal noise to drive equilibration, we show that the underlying dynamics reduce to a deterministic iterative algorithm. Building on the framework of Aifer et al., we analyze the moment evolution of the Ornstein-Uhlenbeck process governing thermodynamic symmetric positive definite (SPD) matrix inversion. We prove that to a first-order approximation, the covariance dynamics are mathematically identical to preconditioned gradient descent on the Frobenius norm of the residual $\tilde{A}^{-1}A-I$. This equivalence demonstrates that thermal fluctuations, while essential for physical thermodynamic hardware, are algorithmically redundant for convex problems with a single global minimum. We validate the resulting algorithm against Thermox (Duffield et al.), a stochastic thermodynamic simulator, achieving speedups exceeding 100,000-fold while remaining competitive with the Newton-Schulz iteration. We also demonstrate acceleration through Schur complement techniques. These results establish a rigorous link between analog thermodynamic computing, statistical physics, and deterministic optimization methods.