Linear response and exact hydrodynamic projections in Lindblad equations with decoupled Bogoliubov hierarchies

Cet article étudie une classe d'équations de Lindblad pour des fermions sans spin présentant des hiérarchies BBGKY découplées, permettant de mapper leur dynamique sur des équations de Schrödier non hermitiennes afin d'obtenir des projections hydrodynamiques exactes et des fonctions de réponse linéaire.

L'essentiel

Imaginez que vous observez une foule de personnes dans une grande place. Normalement, si vous lancez une balle, elle rebondit, les gens bougent, et tout devient chaotique. Mais dans le monde quantique, les règles sont différentes. Cette recherche, menée par Patrik Penc et Fabian Essler, explore ce qui se passe quand cette foule de particules (des fermions sans spin) interagit non seulement entre elles, mais aussi avec un « bruit » extérieur, comme une pluie de bruit blanc qui perturbe leur mouvement.

Voici une explication simple de leur travail, en utilisant des images de la vie quotidienne.

1. Le Problème : Le Chaos Quantique et le Bruit

En physique, décrire comment un système complexe évolue quand il est perturbé par son environnement est souvent un cauchemar mathématique. C'est comme essayer de prédire la trajectoire de chaque goutte de pluie dans une tempête tout en tenant compte de chaque vent.

Les scientifiques utilisent une équation appelée équation de Lindblad pour décrire ce genre de systèmes « ouverts » (qui échangent de l'énergie avec l'extérieur). Le problème ? Ces équations sont généralement si complexes qu'on ne peut les résoudre qu'avec des approximations grossières ou des supercalculateurs.

2. La Découverte : Une Structure Cachée (La Hiérarchie Découplée)

L'idée géniale de cette équipe est de se concentrer sur une classe spéciale de modèles où les règles du jeu sont un peu magiques. Ils ont découvert que, dans ces modèles, la complexité se « découple ».

L'analogie de l'escalier :
Imaginez que la physique de ce système est comme un grand escalier.

• Normalement, pour comprendre le mouvement d'une personne au 10ème étage, vous devez connaître la position de tout le monde sur les étages 1 à 9, et tout est mélangé. C'est la « hiérarchie BBGKY » (un terme technique pour dire que tout dépend de tout).
• Dans les modèles étudiés ici, l'escalier est magique : les étages sont séparés par des murs étanches. Pour savoir ce qui se passe au 2ème étage (deux particules), vous n'avez pas besoin de regarder les étages supérieurs ou inférieurs. Vous pouvez résoudre le problème étage par étage, comme si chaque niveau était un petit puzzle indépendant.

3. Le Tour de Magie : Transformer le Temps

Leur deuxième astuce est de transformer le temps. Habituellement, on étudie comment les choses évoluent dans le temps réel. Ici, ils ont trouvé un moyen de transformer cette équation de mouvement en une équation de Schrödinger imaginaire.

L'analogie du film inversé :
Imaginez que vous filmez une tasse de café qui refroidit. Si vous jouez le film à l'envers, la chaleur semble remonter dans la tasse. C'est ce qu'ils font mathématiquement. Ils transforment l'évolution temporelle du système en un problème de mécanique quantique où le temps est « imaginaire ». Cela permet d'utiliser des outils mathématiques puissants (comme ceux utilisés pour les systèmes intégrables) pour obtenir des réponses exactes, sans approximation.

4. Les Résultats Clés : Ce qu'ils ont appris

A. La Diffusion (Le Mouvement de la Foule)

Dans ces systèmes, si vous avez une conservation de nombre de particules (comme une foule où personne ne rentre ni ne sort), les particules finissent par se disperser de manière diffusive.

• L'image : Imaginez une goutte d'encre dans un verre d'eau. Au début, elle est concentrée. Avec le temps, elle s'étale doucement.
• La découverte : Les auteurs ont pu calculer exactement comment cette « goutte » d'information se propage. Ils ont trouvé que, même si le système est complexe, il existe des « modes hydrodynamiques » (des façons spécifiques dont l'information voyage) qui dominent le comportement à long terme. Ils ont donné la formule exacte de cette propagation, qui ressemble à une courbe en cloche qui s'élargit avec le temps.

B. La Réponse Linéaire (Le Test de Résonance)

Ils ont aussi étudié ce qui se passe si on donne un petit « coup de pouce » au système (une perturbation) alors qu'il est déjà en train d'évoluer.

• L'image : C'est comme si vous tapotiez doucement une cloche qui est déjà en train de résonner dans une pièce bruyante.
• La découverte : Ils ont pu prédire exactement comment le système réagit à ce petit coup. Ils ont montré que le « bruit » (la dissipation) a pour effet d'effacer les détails fins et les pics nets de la réponse, lissant tout comme si on passait un filtre sur une photo. Cela aide à comprendre comment les matériaux réels se comportent dans des conditions de laboratoire réalistes.

C. Intégrabilité vs Chaos

Certains de leurs modèles sont « intégrables » (ils ont des règles cachées qui les rendent prévisibles, comme un jeu d'échecs parfait), d'autres ne le sont pas (plus chaotiques).

• Le résultat surprenant : Même si les règles sont différentes, le comportement à long terme (la diffusion) est très similaire. Cependant, la façon dont le système atteint cet état final est différente. Pour les modèles intégrables, ils ont pu utiliser des techniques de « Bethe Ansatz » (une méthode mathématique très sophistiquée) pour tout résoudre exactement. Pour les autres, ils ont dû utiliser des méthodes numériques, mais la structure « découplée » a rendu la tâche beaucoup plus facile que d'habitude.

En Résumé

Cette recherche est comme si on avait trouvé une clé universelle pour ouvrir des portes mathématiques qui étaient considérées comme verrouillées.

1. Ils ont identifié des systèmes où la complexité se simplifie naturellement (découplage).
2. Ils ont transformé un problème de temps réel difficile en un problème de temps imaginaire plus facile.
3. Ils ont obtenu des formules exactes pour décrire comment l'information se diffuse et comment le système réagit aux perturbations, même en présence de bruit.

C'est une avancée majeure car cela permet de prédire le comportement de systèmes quantiques ouverts (comme ceux utilisés dans les futurs ordinateurs quantiques ou les capteurs) avec une précision que l'on pensait impossible, sans avoir à faire de suppositions approximatives.

Résumé technique

1. Problématique et Contexte

L'étude de la dynamique des systèmes quantiques à plusieurs corps en présence de dissipation (systèmes ouverts) constitue un défi majeur en physique théorique. La dynamique est généralement décrite par une équation de Lindblad (LE). Cependant, l'analyse exacte de ces équations pour des systèmes à N corps est extrêmement difficile, la plupart des études se limitant à des approches perturbatives ou numériques.

L'article se concentre sur une classe spécifique d'équations de Lindblad pour des fermions sans spin, caractérisée par la découplage des hiérarchies BBGKY (Bogoliubov–Born–Green–Kirkwood–Yvon). Dans ces modèles, les équations du mouvement pour les opérateurs à $n$ particules ne dépendent pas des opérateurs à $n+1$ particules, ce qui permet de fermer les systèmes d'équations.

Les auteurs s'interrogent sur trois points principaux :

1. Peut-on obtenir des résultats analytiques exacts pour les projections hydrodynamiques (le comportement à long terme des opérateurs locaux) dans ces systèmes dissipatifs ?
2. Comment calculer les fonctions de réponse linéaire hors équilibre dans ce cadre ?
3. Quel est l'impact de l'intégrabilité (modèles Yang-Baxter intégrables vs non intégrables) sur la dynamique dissipative, sachant que les états stationnaires sont souvent des états à température infinie (matrice densité maximale) ?

2. Méthodologie

Les auteurs développent une approche basée sur la vectorisation des opérateurs dans l'image de Heisenberg.

• Transformation en équation de Schrödinger imaginaire : En vectorisant les opérateurs (en traitant les opérateurs de création et d'annihilation comme des opérateurs agissant sur un espace de Hilbert doublé), les équations du mouvement de Lindblad sont mappées sur une équation de Schrödinger dépendant du temps imaginaire :
  $$\frac{d}{dt}|\Psi(t)\rangle = \mathcal{H} |\Psi(t)\rangle$$
  où $\mathcal{H}$ est un Hamiltonien non hermitien agissant sur un espace de particules effectives (spin-up et spin-down dans la représentation vectorisée).
• Structure des Hamiltoniens : Pour les modèles étudiés, $\mathcal{H}$ prend la forme d'un modèle de fermions en 1D avec des termes de saut purement imaginaires et des interactions à courte portée.
• Résolution spectrale : La dynamique temporelle est obtenue en diagonalisant $\mathcal{H}$. Les auteurs résolvent l'équation aux valeurs propres pour les secteurs à 1, 2, 3 et 4 particules. Ils utilisent des ansatz d'ondes planes et, pour le modèle intégrable (Modèle I), des équations de Bethe ansatz.
• Analyse asymptotique : Le comportement à long terme est dominé par les modes propres dont la partie réelle de la valeur propre est la plus proche de zéro (modes hydrodynamiques diffusifs).

3. Contributions Clés et Résultats

A. Expressions exactes pour les opérateurs propres diffusifs

Les auteurs identifient et construisent explicitement les opérateurs propres du dual de Lindblad qui correspondent aux modes hydrodynamiques.

• Ces opérateurs sont des états liés "particule-trou" (particle-hole bound states) dans l'espace vectorisé.
• Leur amplitude décroît exponentiellement avec la distance entre les opérateurs de création et d'annihilation.
• La longueur de liaison tend vers zéro lorsque le moment total $p \to 0$.
• Les valeurs propres associées suivent une loi de diffusion : $\lambda(p) \approx -D p^2$ pour les petits moments, où $D$ est le coefficient de diffusion dépendant du modèle.

B. Projections hydrodynamiques exactes

L'article fournit des expressions fermées pour la projection des opérateurs locaux quadratiques en fermions (ex: $c^\dagger_{x_0} c_{y_0}$) sur les modes diffusifs à temps longs ($t \gg \gamma^{-1}$).

• Les amplitudes de ces projections sont exprimées en termes de fonctions hypergéométriques confluentes.
• À l'asymptote temporelle, ces amplitudes se réduisent à des lois de puissance décroissantes en $t$. L'exposant dépend de la parité et de la distance spatiale :
  $$\psi(t) \sim (Dt)^{-1 + X/2} \quad \text{ou} \quad (Dt)^{-2 + X/2}$$
  où $X$ et $Y$ sont des combinaisons linéaires des positions.
• Cela démontre que les opérateurs locaux acquièrent des "queues" (tails) hydrodynamiques, contrairement aux opérateurs portant un moment défini qui décroissent exponentiellement.

C. Fonctions de réponse linéaire hors équilibre

Les auteurs appliquent leur formalisme pour calculer la réponse linéaire après un "quench" quantique (changement soudain de paramètres) dans un état initial hors équilibre (ex: état fondamental d'une chaîne à sauts dimerisés).

• Ils montrent comment la dissipation affecte les formes de raies (lineshapes) des fonctions de réponse.
• Résultat majeur : La dissipation lisse les singularités de type "racine carrée" (typiques des systèmes unitaires intégrables) et fait émerger des signatures des modes diffusifs à basse fréquence.

D. Comparaison Intégrable vs Non-Intégrable

L'étude compare des modèles intégrables (Modèle I, lié au modèle de Hubbard avec interaction imaginaire) et non intégrables (Modèles II, III, IV).

• Intégrabilité : Dans le modèle intégrable, les états propres du Hamiltonien non hermitien peuvent être décrits par un ansatz de Bethe (strings $k-\Lambda$).
• Non-intégrabilité : Pour les modèles non intégrables, bien que les états propres ne soient pas de la forme Bethe, la structure de la hiérarchie BBGKY découplée permet toujours d'obtenir des résultats exacts pour les secteurs à deux particules. Les spectres numériques montrent que les modes hydrodynamiques existent également dans les modèles non intégrables, avec des coefficients de diffusion différents.

E. Comportement des opérateurs d'ordre supérieur

• Opérateurs cubiques : Ils décroissent exponentiellement en temps (pas de queues hydrodynamiques) en raison de leur parité impaire.
• Opérateurs quartiques : Ils peuvent présenter des queues hydrodynamiques. Pour le modèle intégrable, une analyse via Bethe ansatz est possible, tandis que pour les modèles non intégrables, le spectre est obtenu numériquement.

4. Signification et Implications

Ce travail est significatif pour plusieurs raisons :

1. Rigueur analytique : Il fournit l'une des rares solutions exactes pour la dynamique hors équilibre et les fonctions de réponse dans des systèmes quantiques ouverts à plusieurs corps, au-delà des approximations de champ moyen ou des simulations numériques limitées.
2. Compréhension de l'hydrodynamique dissipative : Il établit un lien explicite entre la structure spectrale des opérateurs de Lindblad (modes liés dans un Hamiltonien non hermitien) et les lois de puissance observées dans les fonctions de corrélation temporelle.
3. Rôle de l'intégrabilité : Il clarifie que l'intégrabilité dans les systèmes dissipatifs ne garantit pas l'existence de lois de conservation au sens unitaire (les états stationnaires sont à température infinie), mais impose une structure spécifique aux fonctions d'onde des modes hydrodynamiques.
4. Généralité : La méthode de découplage de la hiérarchie BBGKY est applicable en toute dimension et pour divers types de bruit, ouvrant la voie à l'étude de la transport et de la thermalisation dans des systèmes ouverts complexes.

En résumé, l'article démontre que la structure particulière de certaines équations de Lindblad permet de réduire un problème de dynamique quantique ouverte complexe à la résolution d'équations de Schrödinger à quelques corps non hermitiennes, offrant ainsi une fenêtre unique sur la physique hydrodynamique et la réponse linéaire des systèmes quantiques dissipatifs.