Full Quantum Work Statistics for Non-Homogeneous Many-Body Systems

Cet article établit un cadre fondé sur les principes premiers utilisant la théorie de la fonctionnelle de la densité dépendante du temps thermique pour calculer les statistiques complètes du travail quantique et les moments du travail dissipé dans les systèmes à plusieurs corps en interaction, démontrant son pouvoir prédictif dans l'analyse du croisement entre l'isolant de Mott et l'isolant de bande au sein du modèle de Hubbard.

L'essentiel

La vue d'ensemble : Mesurer le « coût » de l'agitation d'un système quantique

Imaginez que vous avez une piste de danse très complexe et bondée (un système quantique) où tout le monde se tient la main et bouge en synchronisation. C'est un « système à corps multiples ». Maintenant, imaginez que vous poussez soudainement le mur de la pièce ou que vous changez le tempo de la musique (une force externe). Les danseurs vont trébucher, se cogner les uns contre les autres et finir par s'installer dans un nouveau rythme.

L'énergie perdue pendant ce trébuchement — la « friction » de la piste de danse — est appelée travail dissipé. Dans le monde quantique, il ne s'agit pas d'un glissement fluide ; c'est un événement chaotique et saccadé, rempli de fluctuations aléatoires.

Ce document présente une nouvelle carte de haute précision (un cadre mathématique) pour prédire exactement quelle quantité d'énergie est perdue et à quel point cette perte sera chaotique, sans avoir besoin de simuler chaque danseur individuellement.

Le problème : La « boîte noire » du chaos quantique

Pendant longtemps, les scientifiques avaient deux manières d'étudier ces systèmes :

1. La voie phénoménologique : Ils devinaient comment le système réagirait en se basant sur des règles générales, comme dire : « Il fait généralement plus chaud quand on le pousse. » C'est comme deviner la météo en regardant le ciel sans thermomètre. C'est utile, mais pas très précis.
2. La voie exacte : Ils essayaient de calculer le mouvement de chaque particule. Pour un système de milliards de particules, c'est comme essayer de compter chaque grain de sable sur une plage pendant un ouragan. C'est informatiquement impossible.

Les auteurs cherchaient une solution « entre les deux » (le juste milieu) : une méthode suffisamment précise pour voir les détails, mais assez simple pour être réellement exécutable sur un ordinateur.

La solution : L'astuce de l'ombre chinoise

Les auteurs ont utilisé une technique appelée Théorie de la fonctionnelle de la densité dépendante du temps thermique (thTDDFT).

Considérez le véritable système quantique complexe comme un spectacle de marionnettes géant et complexe avec des milliers de marionnettes qui interagissent. Il est trop difficile de suivre chaque fil et chaque articulation.

• L'astuce : Au lieu de suivre les vraies marionnettes, ils créent un spectacle d'ombres chinoises. Ce spectacle d'ombres est beaucoup plus simple (c'est un système de particules non interagissantes), mais il est mathématiquement conçu pour projeter la même ombre exacte (densité) sur le mur que le système complexe réel.
• Le bénéfice : En étudiant l'ombre simple, ils peuvent comprendre ce que fait le système complexe. Ils n'ont pas besoin de connaître les secrets de chaque interaction ; ils ont juste besoin de savoir comment l'« ombre » se déplace.

La découverte clé : Diviser la « friction »

Le document fait une distinction subtile entre deux types de « friction » ou de perte d'énergie :

1. La partie « adiabatique » (l'étirement lent) : Imaginez que vous étirez lentement un élastique. Même si vous le faites très lentement, l'élastique résiste parce que sa forme change. C'est une perte d'énergie due au changement de forme du système, et non à cause du chaos.
2. La partie « non-adiabatique » (le claquement soudain) : Imaginez que vous cassez cet élastique d'un coup sec. La perte d'énergie ici provient des secousses et des transitions soudaines et chaotiques.

Les auteurs ont développé un moyen de séparer ces deux éléments. Ils ont montré que la partie « chaotique » (non-adiabatique) est directement liée à la manière dont le système répond à une poussée rapide (une « fonction de relaxation »). En utilisant leur méthode d'ombre chinoise, ils peuvent calculer cette fonction de réponse à partir de principes premiers (les lois fondamentales de la physique) plutôt que de deviner.

Le test : La piste de danse du « Modèle de Hubbard »

Pour prouver que leur carte fonctionne, ils l'ont testée sur un célèbre modèle théorique appelé le modèle de Hubbard.

• Le dispositif : Imaginez une ligne de danseurs (électrons) sur une grille. Ils peuvent sauter à l'emplacement suivant, mais si deux danseurs essaient de se tenir sur le même emplacement, ils subissent un « choc » (répulsion).
• L'expérience : Ils ont appliqué une poussée « alternée » (poussant les danseurs aux numéros impairs d'un côté et ceux aux numéros pairs de l'autre).
• Le résultat : À mesure qu'ils changeaient la force de la poussée et la température, le système passait entre différents « états de la matière » :
  • Isolant de Mott : Les danseurs sont bloqués sur place car ils ont peur de heurter leurs voisins.
  • Isolant de Bande : Les danseurs sont bloqués parce que le sol lui-même est incliné.
  • Isolant d'ordre de liaison (Bond-Order Insulator) : Un étrange terrain intermédiaire où les danseurs se regroupent par paires selon un motif spécifique.

Les auteurs ont découvert que leur méthode pouvait clairement voir les « signatures » de ces différentes phases dans la perte d'énergie. Par exemple, juste à la limite où le système passe d'une phase à une autre, la « friction » (perte d'énergie) augmentait de façon spectaculaire. Cela a confirmé que leur méthode peut détecter des changements subtils dans le monde quantique simplement en mesurant la quantité d'énergie gaspillée.

Pourquoi cela importe

Ce document n'invente pas une nouvelle batterie ou une nouvelle puce informatique. Il fournit plutôt un nouvel outil de mesure.

• Avant : Les scientifiques devaient deviner comment les systèmes quantiques se comporteraient lorsqu'ils étaient poussés, ou ils devaient attendre que les supercalculateurs plantent en essayant de faire les calculs.
• Maintenant : Ils disposent d'une recette fiable, basée sur les « principes premiers », pour calculer exactement l'énergie perdue et les fluctuations du système, même dans des systèmes quantiques complexes et encombrés.

Il comble le fossé entre la réalité désordonnée des particules en interaction et les mathématiques propres et solubles des systèmes d'« ombres », permettant aux scientifiques de prédire le coût thermodynamique des processus quantiques avec une haute précision.

Résumé technique

Résumé Technique : Statistiques Complètes du Travail Quantique pour les Systèmes à N-Corps Non Homogènes

Énoncé du Problème
L'article traite du défi que représente la caractérisation de la thermodynamique hors équilibre des systèmes quantiques à N-corps en interaction, en se concentrant spécifiquement sur les statistiques du travail dissipé lors de protocoles de pilotage à temps fini. Bien que la théorie de la réponse linéaire généralisée (GLRT) ait établi que la distribution complète du travail dissipé peut être exprimée via une fonction de relaxation, l'application pratique de ce cadre a été limitée. En l'absence de données microscopiques, la fonction de relaxation est souvent modélisée de manière phénoménologique, sacrifiant ainsi la précision prédictive. De plus, les approches existantes peinent à fournir une voie de premier principe pour reconstruire cette fonction pour les systèmes non homogènes et corrélés, où la cohérence quantique et les corrélations à N-corps jouent un rôle critique. Les auteurs visent à combler le fossé entre la théorie de la fonctionnelle de la densité thermique (thDFT) et les statistiques du travail hors équilibre afin de fournir un cadre prédictif et microscopique pour ces systèmes.

Méthodologie
Les auteurs développent un formalisme basé sur la Théorie de la Fonctionnelle de la Densité Temporelle Thermique à Réponse Linéaire (LR-thTDDFT). Le cœur de l'approche consiste à dériver la fonction de relaxation, $\psi(t)$, directement à partir de principes microscopiques plutôt que d'une modélisation phénoménologique.

1. Cadre Théorique :

   • Les statistiques du travail sont analysées dans le cadre de la GLRT, où tous les cumulants de la distribution du travail dissipé sont des fonctionnelles bilinéaires de la dérivée temporelle du protocole de pilotage, pondérées par la transformée de Fourier de la fonction de relaxation, $\tilde{\psi}(\omega)$.
   • La fonction de relaxation est décomposée en composantes adiabatiques ($\psi_{ad}$) et non adiabatiques ($\psi_{na}$). La partie adiabatique rend compte de l'irréversibilité due à la déformation temporelle du spectre d'énergie (indépendante de la vitesse de pilotage), tandis que la partie non adiabatique capture les véritables transitions dissipatives entre les niveaux d'énergie.
   • En utilisant les théorèmes de Runge-Gross et de van Leeuwen étendus aux températures finies, les auteurs établissent une correspondance biunivoque entre le potentiel externe dépendant du temps et la densité du système.
   • La composante non adiabatique de la fonction de relaxation est explicitement liée à la partie imaginaire de la fonction de réponse densité-densité, $\chi_{ij}^\beta(\omega)$, du système en interaction.
   • La fonction de réponse en interaction est obtenue à partir de la fonction de réponse de Kohn-Sham (KS) non interagissante via le formalisme de Gross-Kohn, qui incorpore le noyau Hartree-échange-corrélation (Hxc) à température finie, $f_{Hxc}^\beta(\omega)$.

2. Approximations :

   • Pour rendre le problème traitable, les auteurs emploient l'Approximation de la Densité Locale Adiabatique Thermique (thALDA). Dans cette approximation, le noyau Hxc est instantané dans le temps (indépendant de la fréquence) et local dans l'espace.
   • Ce choix néglige les effets de mémoire et les caractéristiques excitoniques au-delà des excitations particule-trou uniques, mais conserve la simplicité et l'efficacité computationnelle des fonctionnelles locales.

3. Système Modèle :

   • Le cadre est appliqué au modèle de Hubbard unidimensionnel soumis à un potentiel externe échelonné. Ce système est choisi pour sa riche diagramme de phases, incluant les phases d'isolant de Mott (MI), d'isolant de bande (BI) et d'isolant d'ordre de liaison (BOI).
   • L'étude analyse le système à travers le croisement Mott-vers-isolant de bande, en faisant varier l'intensité de l'interaction ($U$) et l'amplitude du potentiel échelonné ($v_0$).

Contributions Clés

• Dérivation de Premier Principe : L'article établit une voie directe pour calculer la fonction de relaxation à partir des densités microscopiques et des potentiels KS, éliminant ainsi le besoin de modélisation phénoménologique du noyau de dissipation.
• Décomposition Adiabatique/Non Adiabatique : Les auteurs fournissent une interprétation physique claire de la fonction de relaxation en séparant les contributions adiabatiques et non adiabatiques, liant la première à la susceptibilité statique et la seconde à la réponse dynamique.
• Statistiques du Travail Microscopiques : L'étude démontre comment calculer la distribution complète du travail quantique (spécifiquement ses premiers cumulants) pour des systèmes à N-corps en interaction en utilisant la thTDDFT, jetant un pont entre la DFT à l'équilibre et la thermodynamique hors équilibre.
• Étalonnage : L'approche est étalonnée par rapport à la diagonalisation exacte pour le dimère de Hubbard (modèle de Hubbard à deux sites), montrant un bon accord qualitatif et quantitatif, même dans les régimes fortement corrélés.

Résultats

• Signatures de Phase : Dans le modèle de Hubbard, le premier moment du travail dissipé (production d'entropie irréversible) présente un pic prononcé qui trace le croisement entre les régimes MI et BI, fournissant une signature thermodynamique de la région BOI. Cet accroissement est attribué à la sensibilité accrue des densités thermiques aux variations des paramètres de contrôle près des phases compétitives.
• Comportement des Cumulants : Les trois premiers cumulants de la distribution du travail dissipé ont été analysés en fonction de la durée de la trempe ($\tau$).
  • Le troisième cumulant s'annule dans la limite adiabatique ($\tau \to \infty$), ce qui est cohérent avec le fait que la distribution approche une forme gaussienne où seuls les deux premiers moments subsistent.
  • Le premier et le second cumulant conservent des signatures de la transition de phase (BOI vers BI) même dans le régime adiabatique, sous l'effet de la contribution adiabatique au travail irréversible.
  • Le système dans le régime d'isolant de bande atteint l'adiabaticité à des vitesses de trempe plus rapides par rapport aux phases de Mott ou de BOI.
• Relation d'Incertitude Thermodynamique (TUR) : L'étude confirme que la borne inférieure de la TUR ($F_W \geq 2/\beta$) n'est pas saturée en présence de fluctuations quantiques, même dans le régime de réponse linéaire, distinguant ainsi le comportement quantique des attentes classiques.
• Précision : Les résultats de la LR-thTDDFT pour le dimère de Hubbard montrent des erreurs relatives moyennes d'environ $1,36 \times 10^{-2}$, $2,76 \times 10^{-2}$ et $1,76 \times 10^{-2}$ pour le premier, le second et le troisième cumulant, respectivement, par rapport à la diagonalisation exacte.

Signification et Revendications
L'article affirme fournir un cadre microscopique et transférable pour la thermodynamique quantique dans les systèmes corrélés. En connectant directement les statistiques du travail hors équilibre aux densités thermiques et aux potentiels KS, l'approche offre un outil prédictif qui évite les hypothèses phénoménologiques. Les auteurs soulignent que, bien que l'implémentation actuelle utilise la thALDA (qui possède des limitations connues dans les régimes fortement corrélés dominés par des interactions non locales), le formalisme lui-même est exact dans la réponse linéaire et permet des améliorations systématiques grâce à des noyaux d'échange-corrélation plus raffinés, non locaux ou dépendants de la fréquence.

Ce travail met en évidence l'utilité des méthodes basées sur la DFT pour l'étude de la thermodynamique hors équilibre dans de grands systèmes (par exemple, les gaz atomiques ultra-froids, les matériaux de basse dimensionnalité) où les techniques à N-corps coûteuses en calcul (comme DMRG ou QMC) pourraient être prohibitives. Les auteurs notent avec modestie que si l'étude actuelle capture les signatures de phase et les tendances générales, elle ne tente pas d'extraire des lois d'échelle universelles (telles que l'échelle de Kibble-Zurek) près des points critiques, ce qui nécessiterait des noyaux plus avancés et une analyse dédiée. La contribution primaire est l'établissement du lien théorique et la démonstration de sa faisabilité pour l'extraction de l'irréversibilité thermodynamique à partir de modèles microscopiques.