Méthode de Ritz
La méthode de Ritz est un procédé de calcul d'une solution approchée dans le cas de certains problèmes aux limites.
Elle porte le nom de Walther Ritz, bien qu'elle soit aussi communément appelée méthode de Rayleigh-Ritz et méthode de Galerkine.
En mécanique quantique, un système de particules peut être décrit en termes de « fonctionnelle énergétique » ou opérateur hamiltonien, dont les valeurs propres mesurent l'énergie du système de particules. On attache souvent une grande importance à la connaissance de l'état fondamental (état de plus basse énergie) d'un système, notamment en théorie de la fonctionnelle de la densité, parce que l'état fondamental permet de déterminer les valeurs moyennes d'autres observables du système. En pratique, il est souvent impossible de déterminer toutes les configurations (i.e. les fonctions d'onde propres de l'hamiltonien). Il est plus important de pouvoir trouver des méthodes numériques permettant d'approcher l'état de plus basse énergie.
La méthode Ritz peut être utilisée pour atteindre cet objectif. Dans le langage des mathématiques, il s'agit de la méthode des éléments finis utilisée pour le calcul des vecteurs propres et des valeurs propres d'un système hamiltonien.
Discussion
Comme avec d'autres méthodes variationnelles, on utilise une fonction d'onde approchée, , décomposée sur une base de fonctions, parfois appelée fonction d'onde de test ou ansatz. Cette forme de fonction d'onde est choisie pour respecter les conditions aux limites (et toute autre contrainte physique). La fonction d'onde ansatz, décrite par une combinaison linéaire de fonctions de la base choisie, contient un ou plusieurs paramètres réglables, qui sont calculés pour minimiser l'énergie et ainsi trouver la configuration de plus basse énergie.
On peut montrer que l'énergie de l'état fondamental, , satisfait une inégalité :
le membre de droite de cette inégalité est parfois appelé quotient de Rayleigh. La fonction d'onde ansatz donnera toujours un quotient de Rayleigh supérieur à l'énergie de l'état fondamental.
Si la fonction d'onde ansatz est connue pour être orthogonale à l'état fondamental, alors elle fournira une limite pour l'énergie d'un certain état excité.
La fonction ansatz (parfois dite de Ritz) est une combinaison linéaire de N fonctions de base connues , paramétré par des coefficients inconnus :
Avec un hamiltonien connu, sa valeur moyenne peut s'écrire sous la forme
Les fonctions de base ne sont généralement pas orthogonales, de sorte que la matrice de recouvrement S a des éléments non diagonaux non nuls. Soit ou (la conjugaison du premier) l'ensemble des coefficients à calculer pour minimiser la valeur attendue. Par exemple, en écrivant que les dérivées partielles de par rapport aux sont nulles, on obtient l'égalité suivante pour tout k = 1, 2, . . ., N :
ce qui conduit à un ensemble de N équations séculaires :
Dans les équations ci-dessus, l'énergie et les coefficients sont inconnus. Il s'agit d'un ensemble homogène d'équations linéaires par rapport aux , qui a une solution lorsque le déterminant des coefficients de ces inconnues est nul :
ce qui n'est vrai que pour N valeurs de . De plus, puisque l'hamiltonien est un opérateur hermitien, la matrice H est également hermitienne et les valeurs de sont réelles. La valeur la plus basse parmi (i=1,2,. .,N), , sera la meilleure approximation de l'état fondamental pour les fonctions de base utilisées. Les énergies restantes sont des estimations des énergies d'état excité. Une approximation de la fonction d'onde de l'état i peut être obtenue en trouvant les coefficients à partir de l'équation séculaire correspondante.
Relation avec la méthode des éléments finis
Dans le langage de la méthode des éléments finis, la matrice est précisément la matrice de rigidité de l'hamiltonien dans l'espace des éléments linéaires par morceaux, et la matrice est la matrice de masse . Dans le langage de l'algèbre linéaire, la valeur est une valeur propre de l'hamiltonien discrétisé, et le vecteur est un vecteur propre discrétisé.
Voir aussi
- Méthode de Rayleigh-Ritz
- Théorie de Sturm-Liouville
- Espace de Hilbert
- Méthode de Galerkine
- Théorie des perturbations
Sources
Articles scientifiques
- Walter Ritz (1909) "Über eine neue Methode zur Lösung gewisser Variationsprobleme der mathematischen Physik" Journal für die Reine und Angewandte Mathematik, vol. 135, pages 1 – 61. Disponible en ligne sur : http://gdz.sub.uni-goettingen.de/no_cache/dms/load/img/? IDDOC=261182 .
- JK MacDonald, "Approximations successives par la méthode de – Ritz", Phys. Rev. 43 (1933) 830 Disponible en ligne sur : http://journals.aps.org/pr/abstract/10.1103/PhysRev.43.830
Liens externes
- Méthode Ritz dans l' Encyclopédie des mathématiques
- Gander et Wanner, « From Euler, Ritz, and Galerkin to Modern Computing », SIAM Review, vol. 54, no 4, , p. 627–666 (DOI 10.1137/100804036)
Portail de la chimie
Portail de la physique