Dynamical perturbation theory for eigenvalue problems

by   Maseim Kenmoe, et al.

Many problems in physics, chemistry and other fields are perturbative in nature, i.e. differ only slightly from related problems with known solutions. Prominent among these is the eigenvalue perturbation problem, wherein one seeks the eigenvectors and eigenvalues of a matrix with small off-diagonal elements. Here we introduce a novel iterative algorithm to compute these eigenpairs based on a reformulation of the eigenvalue problem as an algebraic equation in complex projective space. We show from explicit and numerical examples that our algorithm outperforms the usual Rayleigh-Schrödinger expansion on three counts. First, since it is not defined as a power series, its domain of convergence is not a priori confined to a disk in the complex plane; we find that it indeed usually extends beyond the standard perturbative radius of convergence. Second, it converges at a faster slower rate than the Rayleigh-Schrödinger expansion, i.e. fewer iterations are required to reach a given precision. Third, the (time- and space-) algorithmic complexity of each iteration does not increase with the order of the approximation, allowing for higher precision computations. Because this complexity is merely that of matrix multiplication, our dynamical scheme also scales better with the size of the matrix than general-purpose eigenvalue routines such as the shifted QR or divide-and-conquer algorithms. Whether they are dense, sparse, symmetric or unsymmetric, we confirm that dynamical diagonalization quickly outpaces LAPACK drivers as the size of matrices grows; for the computation of just the dominant eigenvector, our method converges order of magnitudes faster than the Arnoldi algorithm implemented in ARPACK.


page 4

page 25

page 26


A fast iterative algorithm for near-diagonal eigenvalue problems

We introduce a novel iterative eigenvalue algorithm for near-diagonal ma...

Convergence of Eigenvector Continuation

Eigenvector continuation is a computational method that finds the extrem...

Iterative Refinement of Schur decompositions

The Schur decomposition of a square matrix A is an important intermediat...

Simple exponential acceleration of the power iteration algorithm

Many real-world problems rely on finding eigenvalues and eigenvectors of...

A multipoint perturbation formula for eigenvalue problems

Standard perturbation theory of eigenvalue problems consists of obtainin...

Solving the Parametric Eigenvalue Problem by Taylor Series and Chebyshev Expansion

We discuss two approaches to solving the parametric (or stochastic) eige...

Fast optimization of viscosities for frequency-weighted damping of second-order systems

We consider frequency-weighted damping optimization for vibrating system...

Please sign up or login with your details

Forgot password? Click here to reset