Hamiltonian System Approach to Distributed Spectral Decomposition in Networks
For researchers needing high-resolution spectral decomposition of large networks, this work offers a novel physics-inspired distributed approach, though it is incremental as it adapts existing physical analogies and numerical methods.
The paper develops distributed algorithms for computing eigenvalues and eigenvectors of graph matrices with higher resolution, modeling spectral computation as Hamiltonian/Lagrangian systems and using symplectic integrators for stability. Numerical simulations on real-world networks demonstrate effectiveness.
Because of the significant increase in size and complexity of the networks, the distributed computation of eigenvalues and eigenvectors of graph matrices has become very challenging and yet it remains as important as before. In this paper we develop efficient distributed algorithms to detect, with higher resolution, closely situated eigenvalues and corresponding eigenvectors of symmetric graph matrices. We model the system of graph spectral computation as physical systems with Lagrangian and Hamiltonian dynamics. The spectrum of Laplacian matrix, in particular, is framed as a classical spring-mass system with Lagrangian dynamics. The spectrum of any general symmetric graph matrix turns out to have a simple connection with quantum systems and it can be thus formulated as a solution to a Schrödinger-type differential equation. Taking into account the higher resolution requirement in the spectrum computation and the related stability issues in the numerical solution of the underlying differential equation, we propose the application of symplectic integrators to the calculation of eigenspectrum. The effectiveness of the proposed techniques is demonstrated with numerical simulations on real-world networks of different sizes and complexities.