Numerical analysis of the adiabatic variable method for the approximation of the nuclear hamiltonian
Maday, Yvon ; Turinici, Gabriel
ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique, Tome 35 (2001), p. 779-798 / Harvested from Numdam

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.

Publié le : 2001-01-01
Classification:  65N25,  35P15,  81V55
@article{M2AN_2001__35_4_779_0,
     author = {Maday, Yvon and Turinici, Gabriel},
     title = {Numerical analysis of the adiabatic variable method for the approximation of the nuclear hamiltonian},
     journal = {ESAIM: Mathematical Modelling and Numerical Analysis - Mod\'elisation Math\'ematique et Analyse Num\'erique},
     volume = {35},
     year = {2001},
     pages = {779-798},
     mrnumber = {1863280},
     zbl = {0995.65112},
     language = {en},
     url = {http://dml.mathdoc.fr/item/M2AN_2001__35_4_779_0}
}
Maday, Yvon; Turinici, Gabriel. Numerical analysis of the adiabatic variable method for the approximation of the nuclear hamiltonian. ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique, Tome 35 (2001) pp. 779-798. http://gdmltest.u-ga.fr/item/M2AN_2001__35_4_779_0/

[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 0846.65056

[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 1470226

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

[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 1470227

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

[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 0642.35001

[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 0513.65079

[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 0165.10801

[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 0799.65112

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

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