A new Lanczos-based algorithm for simulating high-frequency two-dimensional electron spin resonance spectra.
The Lanczos algorithm (LA) is a useful iterative method for the reduction of a large matrix to tridiagonal form. It is a storage efficient procedure requiring only the preceding two Lanczos vectors to compute the next. The quasi-minimal residual (QMR) method is a powerful method for the solution of linear equation systems, Ax = b. In this report we provide another application of the QMR method: we incorporate QMR into the LA to monitor the convergence of the Lanczos projections in the reduction of large sparse matrices. We demonstrate that the combined approach of the LA and QMR can be utilized efficiently for the orthogonal transformation of large, but sparse, complex, symmetric matrices, such as are encountered in the simulation of slow-motional 1D- and 2D-electron spin resonance (ESR) spectra. Especially in the 2D-ESR simulations, it is essential that we store all of the Lanczos vectors obtained in the course of the LA recursions and maintain their orthogonality. In the LA-QMR application, the QMR weight matrix mitigates the problem that the Lanczos vectors lose orthogonality after many LA projections. This enables substantially more Lanczos projections, as required to achieve convergence for the more challenging ESR simulations. It, therefore, provides better accuracy for the eigenvectors and the eigenvalues of the large sparse matrices originating in 2D-ESR simulations than does the previously employed method, which is a combined approach of the LA and the conjugate-gradient (CG) methods, as evidenced by the quality and convergence of the 2D-ESR simulations. Our results show that very slow-motional 2D-ESR spectra at W-band (95 GHz) can be reliably simulated using the LA-QMR method, whereas the LA-CG consistently fails. The improvements due to the LA-QMR are of critical importance in enabling the simulation of high-frequency 2D-ESR spectra, which are characterized by their very high resolution to molecular orientation.
Department of Chemistry, National Tsing Hua University, Hsinchu 30013, Taiwan. ywchiang@mx.nthu.edu.tw
Yun-Wei
YW
Chiang
Jack H
JH
Freed
D000465:Algorithms; D004578:Electron Spin Resonance Spectroscopy; D008958:Models, Molecular; D011789:Quantum Theory
Science. 2001 Jan 12;291(5502):266-9
11253218
Magn Reson Chem. 2005 Nov;43 Spec no.:S256-66
16235203
Biophys J. 2003 Apr;84(4):2619-33
12668470
J Chem Phys. 2008 Apr 28;128(16):165106
18447510
Adv Protein Chem. 2003;63:243-90
12629973
Biophys J. 2006 Jun 15;90(12):4452-65
16565045
Curr Opin Struct Biol. 2006 Oct;16(5):644-53
16949813
J Am Chem Soc. 2009 Feb 25;131(7):2597-605
19191603
J Phys Chem B. 2007 Sep 27;111(38):11260-70
17760438
J Phys Chem B. 2008 Sep 4;112(35):11014-27
18698714
J Phys Chem B. 2008 May 8;112(18):5755-67
18412413