On the electronic path integral normal modes of the Meyer-Miller-Stock-Thoss representation of nonadiabatic dynamics

Résumé technique : Sur les modes normaux de l'intégrale de chemin électronique de la représentation de Meyer–Miller–Stock–Thoss pour la dynamique non adiabatique

Énoncé du problème
La simulation précise et efficace de la dynamique non adiabatique est essentielle pour comprendre le transfert de charge et d'énergie dans les systèmes complexes. Un critère critique de validité de telles méthodes est la conservation de la distribution de Boltzmann quantique (QBD). Alors que la dynamique de Matsubara parvient à conserver la QBD pour la dynamique nucléaire sur une surface unique en tronquant les modes normaux supérieurs de l'intégrale de chemin en temps imaginaire, son extension aux systèmes non adiabatiques (NA-Mats) présente un défi. La méthode NA-Mats existante applique la troncature de Matsubara aux variables nucléaires mais conserve la représentation de Meyer–Miller–Stock–Thoss (MMST) pour les variables électroniques sans troncature des modes normaux. Par conséquent, cette méthode ne parvient pas à conserver la QBD pour les systèmes généraux. Cela soulève une question fondamentale : la troncature des modes normaux supérieurs de l'intégrale de chemin des variables électroniques au sein de la représentation MMST permet-elle la conservation de la QBD, ou ces modes sont-ils fondamentalement différents de leurs homologues nucléaires ?

Méthodologie
Les auteurs étudient les propriétés des modes normaux électroniques au sein du cadre de mapping MMST en utilisant un système purement électronique afin d'isoler le comportement des degrés de liberté électroniques.

• Cadre théorique : L'étude utilise l'Hamiltonien MMST, qui projette un système quantique multi-niveaux sur des oscillateurs harmoniques couplés avec des variables de position ($q$) et de quantité de mouvement ($p$) électroniques. Les auteurs définissent les modes normaux électroniques ($\check{q}_k, \check{p}_k$) en utilisant la même matrice de transformation ($T$) que celle employée pour les modes de Matsubara nucléaires, diagonalisant ainsi la matrice de constante de ressort du polymère en anneau libre.
• Analyse du Liouvillien : Le Liouvillien électronique exact est dérivé dans la base des modes normaux. Pour un système avec des degrés de liberté nucléaires et électroniques découplés, les équations d'évolution des modes normaux sont montrées comme étant séparables et de forme identique aux équations d'évolution des perles (beads).
• Fonctions de corrélation : Les auteurs formulent les fonctions de corrélation de la transformée de Kubo généralisée (KT) en termes de variables de perles et de variables de modes normaux. Ils examinent spécifiquement si les observables (telles que les populations d'états) peuvent être exprimées en fonction d'un nombre fini de modes normaux de basse fréquence.
• Protocole de troncature : Les tests computationnels impliquent de tronquer l'intégrale de chemin de deux manières :
  1. Troncature des perles : Propagation de seulement les $j$ premières perles tout en maintenant les autres à $t=0$.
  2. Troncature des modes normaux : Propagation de seulement les $k$ modes normaux les plus bas tout en maintenant les modes supérieurs à $t=0$.
• Modèle de système : Un système électronique à deux niveaux avec une matrice de potentiel ($V$) asymétrique est simulé en unités réduites ($\beta=1$) avec un polymère en anneau de 8 perles.

Contributions clés et résultats
L'étude présente la première analyse des modes normaux électroniques dans la représentation MMST, produisant plusieurs résultats contre-intuitifs par rapport à la dynamique de Matsubara nucléaire conventionnelle :

1. Observables non finies : Contrairement aux observables nucléaires, qui dépendent souvent d'un nombre fini de modes de Matsubara de basse fréquence, les observables électroniques dans la représentation MMST sont des fonctions de tous les modes normaux. La matrice de transformation ne peut être éliminée lors de la troncature, ce qui signifie que l'observable dépend de l'ensemble complet des modes.
2. Absence de rétrécissement de la distribution : Dans la dynamique de Matsubara nucléaire, la distribution des modes normaux supérieurs se rétrécit, permettant leur intégration. Dans la représentation électronique MMST, les modes normaux supérieurs ne sont pas contraints par la distribution ; leur distribution ne se rétrécit pas, empêchant l'intégration des modes de haute fréquence sans perte de précision.
3. Échec de la conservation de la QBD pour les trajectoires individuelles : L'analyse algébrique et les résultats computationnels démontrent que la troncature de la dynamique dans les modes normaux MMST ne conserve pas la QBD pour une trajectoire unique. La conservation de la QBD dans le cas multi-perles nécessite de compléter l'intégralité de la boucle de l'intégrale de chemin ; la troncature rompt cette condition.
4. Artéfact de moyenne d'ensemble : Bien que la QBD semble être conservée lors de la moyenne sur un ensemble de trajectoires (en raison d'un effet de moyenne), elle n'est pas conservée pour les trajectoires individuelles. Cela distingue la méthode d'une véritable dynamique conservant la QBD.
5. Fonctions de corrélation inexactes : La troncature de la dynamique dans les modes normaux MMST conduit à des fonctions de corrélation inexactes, car tous les modes sont nécessaires pour obtenir une dynamique raisonnable.

Signification et affirmations
L'article conclut que les modes normaux de l'intégrale de chemin MMST ne sont pas optimaux pour développer une méthode de dynamique non adiabatique précise et conservant la QBD. Les auteurs soutiennent que les propriétés spécifiques requises pour que la dynamique de Matsubara conserve la QBD — spécifiquement le fait que les observables soient des fonctions d'un nombre fini de modes et le rétrécissement des distributions des modes supérieurs — ne sont pas satisfaites par les variables électroniques MMST. Par conséquent, l'application simple de la troncature de Matsubara aux variables électroniques dans le cadre NA-Mats ne produit pas une méthode qui conserve la QBD pour une trajectoire unique. Ce travail suggère que des représentations de mapping ou des métriques alternatives pourraient être nécessaires pour atteindre simultanément une méthode de dynamique non adiabatique satisfaisant les critères de conservation de la QBD, de reproduction des oscillations de Rabi et de mise à l'échelle classique.