Numerical analysis of the adiabatic variable method for the approximation of the nuclear hamiltonian

Maday, Yvon; Turinici, Gabriel

ESAIM: Modélisation mathématique et analyse numérique, Tome 35 (2001) no. 4, pp. 779-798.

De nombreux problèmes en chimie quantique portent sur le calcul d’états fondamentaux ou excités de molécules et conduisent à la résolution de problèmes aux valeurs propres. Une des difficultés majeures dans ces calculs est la très grande dimension des systèmes qui sont en présence lors des simulations numériques. En effet les modes propres recherchés sont fonctions de $3n$ variables où $n$ est le nombre de particules (électrons ou noyaux) de la molécule. Afin de réduire la dimension des systèmes à résoudre les chimistes multiplient les idées intéressantes qui permettent d’approcher le système complet. La méthode des variables adiabatiques entre dans ce cadre et nous présentons ici une étude mathématique rigoureuse de cette approximation. En particulier nous proposons un estimateur a posteriori qui pourrait permettre de vérifier l’hypothèse d’adiabaticité faite sur certaines variables ; des simulations numériques qui implémentent cet estimateur sont aussi présentées.

Many problems in quantum chemistry deal with the computation of fundamental or excited states of molecules and lead to the resolution of eigenvalue problems. One of the major difficulties in these computations lies in the very large dimension of the systems to be solved. Indeed these eigenfunctions depend on $3n$ variables where $n$ stands for the number of particles (electrons and/or nucleari) in the molecule. In order to diminish the size of the systems to be solved, the chemists have proposed many interesting ideas. Among those stands the adiabatic variable method; we present in this paper a mathematical analysis of this approximation and propose, in particular, an a posteriori estimate that might allow for verifying the adiabaticity hypothesis that is done on some variables; numerical simulations that support the a posteriori estimators obtained theoretically are also presented.

[1] J. Antihainen, R. Friesner and C. Leforestier, Adiabatic pseudospectral calculation of the vibrational states of the four atom molecules: Application to hydrogen peroxide. J. Chem. Phys. 102 (1995) 1270.

[2] M. Azaiez, M. Dauge and Y. Maday, Méthodes spectrales et les éléments spectraux. Institut de Recherche Mathématique de Rennes, Prépublications 1994-17 (1994).

[3] I. Babuška and C. Schwab, A posteriori error estimation for hierarchic models of elliptic boundary value problems on thin domains. SIAM J. Numer. Anal. 33 (1996) 241-246. | Zbl

[4] C. Bernardi and Y. Maday, Spectral methods, in Handbook of numerical analysis, Vol. V, Part 2, Ph. G. Ciarlet and J.L. Lions Eds., North-Holland, Amsterdam (1997). | MR

[5] C. Bernardi and Y. Maday, Approximations spectrales de problèmes aux limites elliptiques. Springer, Paris, Berlin, New York (1992). | MR | Zbl

[6] G. Caloz and J. Rappaz, Numerical analysis for nonlinear and bifurcation problems, in Handbook of numerical analysis, Vol. V, Part 2, Ph.G. Ciarlet and J.L. Lions Eds., North-Holland, Amsterdam (1997). | MR

[7] C. Canuto, M.Y. Hussaini, A. Quarteroni and T.A. Zang, Spectral methods in fluid dynamics. Springer, Berlin (1987). | MR | Zbl

[8] R. Dutray and J.L. Lions, Analyse mathématique et calcul numérique pour les sciences et les techniques, Tome 5. Masson, CEA, Paris (1984). | Zbl

[9] R. Friesner, J. Bentley, M. Menou and C. Leforestier, Adiabatic pseudospectral methods for multidimensional vibrational potentials. J. Chem. Phys. 99 (1993) 324.

[10] R. Kosloff, Time-dependent quantum-mecanical methods for molecular dynamics. J. Chem. Phys. 92 (1988) 2087.

[11] D. Kosloff and R. Kosloff, Fourier method for the time dependent Schrödinger equation as a tool in molecular dynamics. J. Comp. Phys. 52 (1983) 35. | Zbl

[12] C. Leforestier, Grid representation of rotating triatomics. J. Chem. Phys. 94 (1991) 6388.

[13] J.L. Lions and E. Magenes, Problèmes aux limites non-homogènes et applications. Dunod, Paris (1968). | Zbl

[14] R. Verfürth, A posteriori error estimates for non-linear problems. Finite element discretisations of elliptic equations. Math. Comp. 62 (1994) 445-475 | Zbl

[15] R. Verfürth, A review of a posteriori error estimates and adaptative mesh-refinement techniques. Wiley-Teubner, Stuttgart (1997). | Zbl

[16] K. Yamashita, K. Mokoruma and C. Leforestier, Theoretical study of the highly vibrationally excited states of $FH{F}^{-}$: Ab initio potential energy surface and hyperspherical formulation. J. Chem. Phys. 99 (1993) 8848.