skip to main content

A low-storage filter diagonalization method for quantum eigenenergy calculation or for spectral analysis of time signals

Mandelshtam, Vladimir A. ; Taylor, Howard S.

The Journal of chemical physics, 1997-03, Vol.106 (12), p.5085-5090 [Periódico revisado por pares]

Texto completo disponível

Citações Citado por
  • Título:
    A low-storage filter diagonalization method for quantum eigenenergy calculation or for spectral analysis of time signals
  • Autor: Mandelshtam, Vladimir A. ; Taylor, Howard S.
  • É parte de: The Journal of chemical physics, 1997-03, Vol.106 (12), p.5085-5090
  • Descrição: A new version of the filter diagonalization method of diagonalizing large real symmetric Hamiltonian matrices is presented. Our previous version would first produce a small set of adapted basis functions by applying the Chebyshev polynomial expansion of the Green’s function on a generic initial vector χ. The small Hamiltonian, H, and overlap, S, matrices would then be evaluated in this adapted basis and the corresponding generalized eigenvalue problem would be solved yielding the desired spectral information. Here in analogy to a recent work by Wall and Neuhauser [J. Chem. Phys. 102, 8011 (1995)] H and S are computed directly using only the Chebyshev coefficients cn=〈χ|Tn(Ĥ)|χ〉, calculation of which requires a minimal storage if the Ĥ matrix is sparse. The expressions for H and S are analytically simple, computationally very inexpensive and stable. The method can be used to obtain all the eigenvalues of Ĥ using the same sequence {cn}. We present an application of the method to a realistic quantum dynamics problem of calculating all bound state energies of H3+ molecule. Since the sequence {cn} is the only input required to obtain all the eigenenergies, the present method can be reformulated for the problem of spectral analysis of a real symmetric time signal defined on an equidistant time grid. The numerical example considers a model signal C(tn)=∑kdk cos(tnωk) generated by a set of N=100 000 frequencies and amplitudes, (ωk,dk). It is demonstrated that all the ωk’s and dk’s can be obtained to very high precision using the minimal information, i.e., 200 000 sampling points.
  • Idioma: Inglês

Buscando em bases de dados remotas. Favor aguardar.