Decay bounds and O(n) algorithms for approximating functions of sparse matrices

Electronic Transactions on Numerical Analysis

Volumn 28, Issue , 2007, Pages 16-39

  Benzi, Michele ^a   Razouk, Nader ^a  

^a EMORY UNIVERSITY   (United States)

Author keywords

Chebyshev polynomials; Decay rates; Density matrix; Determinant; Faber polynomials; Linear time algorithms; Matrix functions; Sparse and banded matrices; Trace

Indexed keywords

APPROXIMATION ALGORITHMS; CLUSTERING ALGORITHMS; DECAY (ORGANIC); POLYNOMIALS;

BANDED MATRICES; CHEBYSHEV POLYNOMIALS; DECAY RATE; DENSITY MATRIX; DETERMINANT; FABER POLYNOMIALS; LINEAR-TIME ALGORITHMS; MATRIX FUNCTIONS; TRACE;

MATRIX ALGEBRA;

EID: 54549085632     PISSN: None     EISSN: 10689613     Source Type: Journal    
DOI: None     Document Type: Article

References (79)

1. P. W. ANDERSON, Absence of diffusion in certain random lattices, Phys. Review, 109 (1958), pp. 1492-1505.
2. R. BAER AND M. HEAD-GORDON, Chebyshev expansion methods for electronic structure calculations on large molecular systems, J. Chem. Phys., 107 (1997), pp. 10003-10013.
3. _, Sparsity of the density matrix in Kohn-Sham density functional theory and an assessment of linear system-size scaling methods, Phys. Rev. Lett., 79 (1997), pp. 3962-3965.
4. Z. BAI, M. FAHEY AND G. H. GOLUB, Some large-scale matrix computation problems, J. Comput. Appl. Math., 74 (1996), pp. 71-89.
5. Z. BAI, M. FAHEY, G. H. GOLUB, M. MENON AND E. RICHTER, Computing partial eigenvalue sums in electronic structure calculations, Tech. Report SCCM-98-03, Stanford University, 1998.
6. Z. BAI AND G. H. GOLUB, Bounds for the trace of the inverse and the determinant of symmetric positive definite matrices, Ann. Numer. Math., 4 (1997), pp. 29-38.
7. C. BEKAS, E. KOKIOPOULOU, AND Y. SAAD, An estimator for the diagonal of a matrix, Appl. Numer. Math., to appear.
8. M. BENZI AND G. H. GOLUB, Bounds for the entries of matrix functions with applications to preconditioning, BIT, 39 (1999), pp. 417-438.
9. M. BENZI AND M. TŮMA, Orderings for factorized sparse approximate inverse preconditioners, SIAM J. Sci. Comput., 21 (2000), pp. 1851-1868.
10. L. BERGAMASCHI, M. CALIARI AND M. VIANELLO, Efficient approximation of the exponential operator for discrete 2D advection-diffusion problems, Numer. Linear Algebra Appl., 10 (2003), pp. 271-289.
11. L. BERGAMASCHI AND M. VIANELLO, Efficient computation of the exponential operator for large, sparse, symmetric matrices, Numer. Linear Algebra Appl., 7 (2000), pp. 27-45.
12. D. A. BINT, N. J. HIGHAM AND B. MEINI, Algorithms for the matrix pth root, Numer. Algorithms, 39 (2005), pp. 349-378.
13. A. BORIÇI, A Lanczos approach to the inverse square root of a large and sparse matrix, J. Comput. Phys., 162 (2000), pp. 123-131.
14. D. R. BOWLER AND M. J. GILLAN, Density matrices in O(N) electronic structure calculations: theory and applications, Comp. Phys. Comm., 120 (1999), pp. 95-108.
15. D. CALVETTI, S.-M. KIM AND L. REICHEL, Quadrature rules based on the Arnoldi process, SIAM J. Matrix Anal. Appl., 26 (2005), pp. 765-781.
16. M. CRAMER AND J. EISERT, Correlations, spectral gap and entanglement in harmonic quantum systems on generic lattices, New J. Phys., 8 (2006), Art. No. 71.
17. M. CRAMER, J. EISERT, M. B. PLENIO AND J. DRETSSIG, Entanglement-area law for general bosonic harmonic lattice systems, Phys. Rev. A, 73 (2006), Art. No. 012309.
18. J. H. CURTISS, Faber polynomials and the Faber series, Amer. Math. Monthly, 78 (1971), pp. 577-596.
19. P. I. DAVIES AND N. J. HIGHAM, Computing f(A)b for matrix functions f, in A. Boriçi, A. Frommer, B. Joo, A. Kennedy and B. Pendleton, eds., QCD and Numerical Analysis III, Lect. Notes Comput. Sci. Eng., Vol. 47, Springer-Verlag, Berlin, 2005, pp. 15-24.
20. N. DEL BUONO AND L. LOFEZ, A survey on methods for computing matrix exponentials in numerical schemes for ODEs, Lecture Notes in Comput. Sci., 2658 (2003), pp. 111-120.
21. N. DEL BUONO, L. LOPEZ AND R. PELUSO, Computation of the exponential of large sparse skew-symmetric matrices, SIAM J. Sci. Comput., 27 (2005), pp. 278-293.
22. S. DEMKO, W. F. MOSS AND P. W. SMITH, Decay rates for inverses of band matrices, Math. Comp., 43 (1984), pp. 491-499.
23. L. DIECI, B. MORINI AND A. PAPINI, Computational techniques for real logarithms of matrices, SIAM J. Matrix Anal. Appl., 17 (1996), pp. 570-593.
24. R. DIESTEL, Graph Theory, Springer-Verlag, Berlin, 2000.
25. V. L. DRUSKIN AND L. A. KNIZHNERMAN, Two polynomial methods for computing functions of symmetric matrices, Comput. Math. Math. Phys., 29 (1989), pp. 112-121.
26. _, Extended Krylov subspaces: approximation of the matrix square root and related functions, SIAM J. Matrix Anal. Appl., 19 (1998), pp. 755-771.
27. I. S. DUFF, A. M. ERISMAN, C. W. GEAR, AND J. K. REID, Sparsity structure and Gaussian elimination, SIGNUM Newsletter, Association for Computing Machinery, New York, 23 (1988), pp. 2-8.
28. M. EIERMANN AND O. ERNST, A restarted Krylov subspace method for the evaluation of matrix functions, SIAM J. Numer. Anal., 44 (2006), pp. 2481-2504.
29. S. W. ELLACOTT, Computation of Faber series with application to numerical polynomial approximation in the complex plane, Math. Comp., 40 (1983), pp. 575-587.
30. L. FOX AND I.B. PARKER, Chebyshev Polynomials in Numerical Analysis, Oxford University Press, London, 1968.
31. A. FROMMER AND V. SIMONCINI, Matrix functions, Tech. Report, Dipartimento di Matematica, Università di Bologna, 2006.
32. _, Stopping criteria for rational matrix functions of Hermitian and symmetric matrices, Tech. Report, Dipartimento di Maternatica, Università di Bologna, 2007.
33. D. GAIER, Lectures on Complex Approximation, Birkhäuser, Boston, 1987.
34. G. GALLI, Large-scale electronic structure calculations using linear scaling methods, Phys. Stat. Sol. B, 217 (2000), pp. 231-249.
35. E. GALLOPULOS AND Y. SAAD, Efficient solution of parabolic equations by Krylov suhspace methods, SIAM J. Sci. Stat. Comput., 13 (1992), pp. 1236-1264.
36. S. GOEDECKER, Decay properties of the finite-temperature density matrix in metals, Phys. Rev. B, 58 (1998), pp. 3501-3502.
37. _, Linear scaling electronic structure methods, Rev. Mod. Phys., 71 (1999), pp. 1085-1123.
38. G. H. GOLUB AND G. MEURANT, Matrices, moments, and quadratures, in Numerical Analysis 1993, D. F. Griffiths and G. A. Watson, eds., Pitman Research Notes in Mathematics, vol. 303, Essex, England, UK, 1994, pp. 105-156.
39. G. H. GOLUB AND C. F. VAN LOAN, Matrix Computations, 3rd edition, Johns Hopkins University Press, Baltimore, MD, 1996.
40. K. GRÖCHENIG AND M, LEINERT, Symmetry and inverse-closedness of matrix algebras and functional calculus for infinite matrices, Trans. Amer. Math. Soc., 358 (2006), pp. 2695-2711.
41. G. I. HAROREAVES AND N.J. HIGHAM, Efficient algorithms for the matrix cosine and sine, Numer. Algorithms, 40 (2005), pp. 383-400.
42. L. HE AND D. VANDERBILT, Exponential decay properties of Wannier functions and related quantities, Phys. Rev. Lett., 86 (2001), pp. 5341-5344.
43. N. J. HIOHAM, Functions of matrices, in L. Hogben, ed., Handbook of Linear Algebra, Chapman and Hall/CRC, Boca Raton, FL, 2006.
44. M. HOCHBRUCK AND C. LUBICH, On Krylov subspace approximations to the matrix exponential operator, SIAM J. Numer. Anal., 34 (1997), pp. 1911-1925.
45. _, Exponential integrators for large systems of differential equations, SIAM J. Sci. Comput., 19 (1998), pp. 1552-1574.
46. R. A. HORN AND C. R. JOHNSON, Topics in Matrix Analysis, Cambridge University Press, Cambridge, 1994.
47. I. C. F. IPSEN AND D. J. LEE, Determinant approximations, Tech. Report, Department of Mathematics, North Carolina State University, Raleigh, NC, 2006.
48. A. ISERLES, How large is the exponential of a banded matrix?, New Zealand J. Math., 29 (2000), pp. 177-192.
49. S. JAFFARD, Propriétés des matrices "bien localisées" près de leur diagonale et quelques applications, Ann. Inst. H. Poincaré Anal. Non Linéaire, 7 (1990), pp. 461-476.
50. C. S. KENNEY AND A. J. LAUB, Condition estimates for matrix functions, SIAM J. Matrix Anal. Appl., 10 (1989), pp. 191-209.
51. L. A. KNIZHNERMAN, Calculation of functions of unsymmetric matrices using Arnoldi's method, Comput. Math. Math. Phys., 31 (1991), pp. 5-9.
52. W. KOHN, Density functional and density matrix method scaling linearly with the number of atoms, Phys. Rev. Lett., 76 (1996), pp. 3168-3171.
53. 2 - I = 0, Math. Tables Other Aids Comp., 12 (1958), pp. 109-116.
54. D. J. LEE AND I. C. F. IPSEN, Zone determinant expansions for nuclear lattice simulations, Phys. Rev. C, 68 (2003), pp. 064003-1/8.
55. X.-P. LI, R. W. NUNES AND D. VANDERBILT, Density-matrix electronic-structure method with linear system-size scaling, Phys. Rev. B, 47 (1993), pp. 10891-10894.
56. L. LOPEZ AND A. PUGLIESE, Decay behaviour of functions of skew-symmetric matrices, in Proceedings of HERCMA 2005, 7th Hellenic-European Conference on Computer Mathematics and Applications, September 22-24, 2005, Athens, E. A. Lipitakis, ed., Electronic Editions LEA, Athens.
57. L. LOPEZ AND V. SIMONCINI, Analysis of projection methods for rational function approximations to the matrix exponential, SIAM J. Numer. Anal., 44 (2006), pp. 613-635.
58. A. I. MARKUSHEVICH, Theory of Functions of a Complex Variable, Vol. III, Prentice-Hall, Englewood Cliffs, NJ, 1967.
59. P. E. MASLEN, C. OCHSENFELD, C. A. WHITE, M. S. LEE AND M. HEAD-GORDON, Locality and sparsity of ab initio one-particle density matrices and localized orbitals, J. Phys. Chem. A, 102 (1998), pp. 2215-2222.
60. G. MEURANT, A review on the inverse of symmetric tridiagonal and block tridiagonal matrices, SIAM J, Matrix Anal. Appl., 13 (1992), pp. 707-728.
61. J. M. MLLAM AND G. E. SCUSERIA, Linear scaling conjugate gradient density matrix search as an alternative to diagonalization for first principles electronic structure calculations, J. Chem. Phys., 106 (1997), pp. 5569-5577.
62. C. MOLER AND C. VAN LOAN, Nineteen dubious ways to compute the exponential of a matrix, twenty-five years later, SIAM Rev., 45 (2003), pp. 3-49.
63. I. MORET AND P. NOVATI, The computation of functions of matrices by truncated Faber series, Numer. Funct. Anal. Optim., 22 (2001), pp. 697-719.
64. _, Interpolating functions of matrices on zeros of kernel polynomials, Numer. Linear Algebra Appl., 11 (2005), pp. 337-353.
65. R. NABBEN, Decay rates of the inverse of nonsymmetric tridiagonal and band matrices, SIAM J. Matrix Anal. Appl., 20 (1999), pp. 820-837.
66. M. K. NO AND J. PAN, Approximate inverse circulant-plus-eliagonal preconditioners for Toeplitz-plusdiagonal matrices, Tech. Report, Department of Mathematics, Hong Kong Baptist University, 2006.
67. A. REUSKEN, Approximation of the determinant of large sparse symmetric positive definite matrices, SIAM J. Matrix Anal. Appl., 23 (2002), pp. 799-818.
68. R. F. RINEHART, The equivalence of definitions of a matric function, Amer. Math. Monthly, 62 (1955), pp. 395-414.
69. Y. SAAD, Analysis of some Krylov subspace approximations to the matrix exponential operator, SIAM J. Numer. Anal., 29 (1992), pp. 209-228.
70. V. E. SACKSTEDER, Linear Algebra with Disordered Sparse Matrices that have Spatial Structure: Theory and Computation, PhD thesis, Department of Physics, Università degli Studi di Roma "La Sapienza," Rome, Italy, 2004.
71. _, O( N) algorithms for disordered systems, Numer. Linear Algebra Appl., 12 (2005), pp. 827-838.
72. N. SCHUCH, J. I. CIRAC AND M. M. WOLF, Quantum states on harmonic lattices, Comm. Math. Phys., 267 (2006), pp. 65-92.
73. V. I. SMIRNOV AND N. A. LEBEDEV, Functions of a Complex Variable: Constructive Theory, MIT Press, Cambridge, MA, 1968.
74. D. T. SMITH, Exponential decay of resolvents and discrete eigenfunctions of banded infinite matrices, J. Approx. Theory, 66 (1991), pp. 83-97.
75. P. K. SUETIN, Fundamental properties of Faber polynomials, Russian Math. Surveys, 19 (1964), p. 121-149.
76. H. TAL-EZER, Polynomial approximation of functions of matrices and applications, J. Sci. Comput., 4 (1989), pp. 25-60.
77. L. N. TREFETHEN, Numerical computation of the Schwarz-Christoffel transformation, SIAM J. Sci. Stat. Comput., 1 (1980), pp. 82-102.
78. L. N. TREFETHEN, M. CONTEDINI AND M. EMBREE, Spectra, pseudospectra, and localization for random bidiagonal matrices, Comm. Pure Appl. Math., 54 (2001), pp. 595-623.
79. L. N. TREFETHEN AND M. EMBREE, Spectra and Pseudospectra: The Behavior of Nonnormal Matrices and Operators, Princeton University Press, Princeton, NJ, 2005.